MicroRNA Expression Profiling in Behçet’s Disease

Puccetti, Antonio; Pelosi, Andrea; Fiore, Piera Filomena; Patuzzo, Giuseppe; Lunardi, Claudio; Dolcino, Marzia

doi:https://doi.org/10.1155/2018/2405150

Journal of Immunology Research

On this page

Abstract Introduction Materials and Methods Results Discussion Data Availability Conflicts of Interest Authors’ Contributions Supplementary Materials References Copyright Related Articles

Special Issue

Immune System and Chronic Diseases 2018

View this Special Issue

Research Article | Open Access

Volume 2018 | Article ID 2405150 | https://doi.org/10.1155/2018/2405150

MicroRNA Expression Profiling in Behçet’s Disease

Antonio Puccetti,^1,2Andrea Pelosi,¹Piera Filomena Fiore,¹Giuseppe Patuzzo,³Claudio Lunardi,³and Marzia Dolcino³

Academic Editor: Andréia M. Cardoso

Received26 Jan 2018

Accepted07 Mar 2018

Published07 May 2018

Abstract

Background. Behçet’s disease (BD) is a chronic inflammatory multisystem disease characterized by oral and genital ulcers, uveitis, and skin lesions. MicroRNAs (miRNAs) are key regulators of immune responses. Differential expression of miRNAs has been reported in several inflammatory autoimmune diseases; however, their role in BD is not fully elucidated. We aimed to identify miRNA expression signatures associated with BD and to investigate their potential implication in the disease pathogenesis. Methods. miRNA microarray analysis was performed in blood cells of BD patients and healthy controls. miRNA expression profiles were analyzed using Affymetrix arrays with a comprehensive coverage of miRNA sequences. Pathway analyses were performed, and the global miRNA profiling was combined with transcriptoma data in BD. Deregulation of selected miRNAs was validated by real-time PCR. Results. We identified specific miRNA signatures associated with BD patients with active disease. These miRNAs target pathways relevant in BD, such as TNF, IFN gamma, and VEGF-VEGFR signaling cascades. Network analysis revealed several miRNAs regulating highly connected genes within the BD transcriptoma. Conclusions. The combined analysis of deregulated miRNAs and BD transcriptome sheds light on some epigenetic aspects of BD identifying specific miRNAs, which may represent promising candidates as biomarkers and/or for the design of novel therapeutic strategies in BD.

1. Introduction

Behçet’s disease (BD) is a rare and chronic multisystem disease characterized by a triple-symptom complex of recurrent oral aphthous ulcers, genital ulcers, and uveitis. Moreover, manifestations of vascular, articular, neurologic, urogenital, gastrointestinal, pulmonary, and cardiac involvement may occur. Hippocrates described BD in the fifth century BCE. In 1930, the Greek ophthalmologist Benediktos Adamantiades reported a patient with inflammatory arthritis, oral and genital ulcers, phlebitis, and iritis. The disease is named after the Turkish dermatologist Hulusi Behçet, who identified it in a patient in 1924 and published a description of the disease in 1937.

As virtually no unique histological or laboratory features have been identified to help in the diagnosis of the disease, clinical features are used to define and diagnose Behçet syndrome. An international study group on Behçet’s disease has recently revised the criteria for the classification/diagnosis of BD [1].

There are sporadic cases of BD all around the world, but it is most frequently seen along the ancient Silk Route because of its frequency in the Middle East and far-east Asia (prevalence of 14–20/100,000 inhabitants), and these regions have traditionally been considered the endemic areas for the condition [2].

BD is a sporadic disease, but a familial aggregation is well known. Carriers of HLA-B51/HLA-B5 have an increased risk of developing Behçet’s disease compared with noncarriers. HLA-B51 is the strongest associated genetic factor, and it has been shown to be more prevalent in Turkish, Middle Eastern, and Japanese populations, with a higher prevalence of Behçet’s disease in these populations [3]. Non-HLA genes also contribute to the development of BD [3]. Genome-wide association studies have shown that polymorphisms in genes encoding for cytokines, activator factors, and chemokines are associated with increased BD susceptibility. Among cytokines, IL-10 polymorphisms cause a reduction in the serum level of IL-10, an inhibitory cytokine that regulates innate and adaptive immune responses; on the other hand, IL-23 receptor polymorphism, which reduces the response to IL-23 stimulation, is associated with protection from BD [3–5]. Recent data reported also associations with CCR1, STAT4, and KLRC4 encoding for a chemokine receptor, a transcription factor implicated in IL-12 and IL-23 signaling and a natural killer receptor [6, 7]. Finally, susceptibility genes involved in the innate immune response to microbial exposure have recently been identified by Immunochip analysis [8].

Increased Th1, CD4⁺ and CD8⁺ T cell, γδ⁺ T cell, and neutrophil activities have been found both in the serum and in inflamed tissues of BD patients, suggesting the involvement of innate and adaptive immunity in the pathogenesis of BD [2, 9]. Studies on T lymphocytes have suggested a Th1-predominant response. Both CD4⁺ and CD8⁺ lymphocytes are higher in the peripheral blood, with increased levels of IL-2 and interferon- (INF-) γ cytokines [10]. The cytokine Th17 may also play an important role in the pathogenesis of the disease [2, 11]. This hypothesis is supported by the observation of high IL-21 and IL-17 levels in sera of patients affected by BD with neurologic involvement [12, 13]. Another study has reported a higher Th17/Th1 ratio in peripheral blood of patients with BD compared to healthy controls, and this ratio was higher in patients with uveitis or folliculitis compared with patients without these manifestations [14, 15].

MicroRNAs (miRNAs) are small noncoding RNAs that play an important role in the regulation of several biological processes through their interaction with cellular messenger RNAs [16]. Inflammatory responses have an impact on miRNA expression, regulating their biogenesis by altering the transcription and processing of precursor transcripts or influencing the stabilization of mature miRNAs [17, 18]. In recent years, the number of miRNAs implicated in immune system development and function has dramatically enhanced, and there has been widespread discussion of their potential use as therapeutics for immunological diseases [16]. Indeed, the aberrant expression of miRNAs frequently occurs in human diseases, including hematological disorders and autoimmunity [19, 20]. The concept that miRNAs participate in the pathogenesis of diseases, especially refractory diseases with unidentified mechanisms, might lead to a novel effective treatment. A number of studies have reported a differential expression of miRNAs in several inflammatory autoimmune diseases, such as in rheumatoid arthritis (RA), multiple sclerosis, systemic lupus erythematosus, psoriasis, and systemic sclerosis [21]. These studies highlighted a deep implication of miRNAs as regulatory molecules in autoimmunity and the intriguing possibility to use miRNAs as disease biomarkers in these immunological disorders.

As far as BD concerns, little is known about miRNA expression; moreover, no high-throughput miRNA expression studies have been conducted to identify miRNAs specifically associated with the disease and no study has been so far performed which combines the analysis of blood microRNAs with transcriptional profiles in patients with BD.

In the present study, we performed a miRNA microarray analysis on peripheral blood mononuclear cells (PBMCs) of BD patients. We identified specific miRNA signatures associated with patients with active BD. These deregulated miRNAs target signaling pathways typically implicated in BD pathogenesis, such as TNF, interferon gamma, and VEGF and VEGFR signaling cascades.

The modular analysis of differentially expressed genes in BD revealed pathogenetically relevant networks that are possibly targeted by the identified miRNAs. This study sheds light on some aspects of BD pathogenesis identifying deregulated miRNAs as promising candidates for the discovery of disease biomarkers and/or as molecular tools for designing novel therapeutic strategies in BD.

2. Materials and Methods

2.1. Patients

A group of 6 subjects with BD was used for the gene array study. All the patients attended the Unit of Autoimmune Diseases at the University Hospital in Verona, Italy.

All patients fulfilled the International Criteria for Behçet Disease (ICBD): oral aphthosis, genital ulcers, and ocular lesions were each given 2 points, whereas 1 point was assigned to each of skin lesions, vascular manifestations, and neurological manifestations. A patient scoring 4 points or above was classified as having BD [22, 23]. At enrollment, none of the patients had active infections or was affected by malignancies.

The clinical features of the patients are reported in Table 1 that also includes a description of the BD patients selected for the gene array study.

A written informed consent was obtained from all the participants to the study. The study was approved by the local Ethical Committee of the Azienda Ospedaliera Universitaria of Verona, Verona, Italy. All investigations have been conducted according to the principles expressed in the Helsinki declaration.

2.2. Microarray Analysis

Blood samples were collected in BD Vacutainer K₂EDTA tubes using a 21-gauge needle. Peripheral blood mononuclear cells (PBMC) were obtained upon stratification on Lympholyte® cell separation density gradient (Cedarlane, Burlington, Canada). Total RNA extraction from PBMC was performed with miRNeasy Mini Kit following the manufacturer’s protocol (Qiagen GmbH, Hilden, Germany). RNA concentration and purity were evaluated by spectrophotometric analysis (NanoDrop 2000; Thermo Fisher Scientific, Wilmington, DE, USA), and a further check for RNA integrity was performed with 2100 Bioanalyzer (Agilent Genomics, Santa Clara, CA, USA) before microarray hybridization. Sample hybridization and scanning were performed as recommended by the Affymetrix (Affymetrix; Thermo Fisher Scientific Inc.) supplied protocols, by the Cogentech Affymetrix Microarray Unit (Campus IFOM-IEO, Milan, Italy), and Affymetrix GeneChip® miRNA 4.0 (Affymetrix; Thermo Fisher Scientific Inc., Waltham, MA, USA) was used. The GeneChip miRNA 4.0 arrays contain more than 30,000 probes including those encoding for 2578 mature human miRNAs, according to Sanger miRBase v.20.

The arrays were analyzed employing the Transcriptome Analysis Console (TAC) 4.0 software (Applied Biosystems, Foster City, CA USA, by Thermo Fisher Scientific, Waltham, MA, USA). The Signal Space Transformation- (SST-) Robust Multiarray Average (RMA) algorithm was applied to background-adjust, normalize, and log-transform signal intensity.

Relative expression levels of each microRNA were validated applying a one-way analysis of variance (ANOVA) and false discovery rate (FDR) correction (). MicroRNAs with an expression level of at least 1.5-fold different in the test sample versus control sample were analyzed.

Targeted genes of significantly modulated miRNAs were identified using the integrative database for human microRNA target predictions mirDIP [24]. All the source filters and a very high confidence class were applied for our analyses.

Pathway enrichment analysis of miRNA gene targets and differentially expressed genes (DEGs) in Behçet’s disease was carried out using FunRich (http://www.funrich.org/) [25], and only Bonferroni-corrected enriched values ≤0.01 calculated by the hypergeometric test were considered.

Pathways enrichment analysis of DEGs in BD was also performed employing the Panther expression analysis tools (http://pantherdb.org/) [26].

2.3. Protein-Protein Interaction (PPI) Network Construction and Network Clustering

Differentially expressed genes (DEGs) in BD samples from our previous analyses (Puccetti et al. unpublished observations) were mapped to the STRING database (version 1.0; http://string-db.org/) [27] to detect protein-protein interactions (PPI) pairs validated by experimental studies. A score of ≥0.7 for each PPI pair was selected to construct the PPI network. For the topological analysis of the built network, Cytoscape software was used and network clustering analysis was performed with the MCODE plugin of Cytoscape, based on the thresholds of module [28].

2.4. Real-Time PCR

Mature miRNA expression was assayed by TaqMan® Advanced miRNA assay chemistry (Applied Biosystems, Foster City, CA, USA). Briefly, 10 ng of total RNA was reverse transcribed and preamplified with TaqMan Advanced miRNA cDNA synthesis kit following the manufacturer’s instructions (Applied Biosystems, Foster City, CA, USA). Preamplified cDNA was diluted 1/10 in nuclease-free water, and 5 μL of diluted cDNA for each replicate was loaded in PCR. 20 μL PCR reactions were composed by 2x Fast Advanced Master Mix and TaqMan Advanced miRNA assays for hsa-miR-143-3p (477912_mir), hsa-miR-199a-5p (478231_mir), and hsa-miR-4505 (477842_mir). The mean of Ct for hsa-miR-16-5p (477860_mir) and hsa-miR-26a-5p (477995_mir) expression was used to normalize miRNA expression. Real-time PCR was carried out in triplicate on a QuantStudio 6 Flex instrument (Applied Biosystems, Foster City, CA, USA). Expression values were reported as fold change with respect to healthy controls by the ΔΔCt method using QuantStudio Real-Time PCR system software versus 1.3.

3. Results

3.1. High-Throughput MicroRNA Expression Profiling in Peripheral Blood Mononuclear Cells of Behçet’s Disease

Since a global miRNA expression analysis in BD with an updated coverage of miRNA sequences has not been performed yet, we wanted to provide a careful description of miRNAs associated with BD by interrogating the transcription of a large amount of different microRNA sequences in BD PBMCs by microarray strategy.

Therefore, PBMCs derived from 6 patients with BD and 6 healthy age- and sex-matched donors were analyzed using a dedicated and high-density array with a coverage for more than 2500 human microRNA transcripts and all mature miRNA sequences in miRBase Release 20. The clinical features of patients included in the microarray study are reported in Table 1.

Microarray analysis revealed a high number (269) of modulated miRNAs that satisfied the FDR-corrected value criterion () and the fold change criterion (), showing a robust and statistically significant variation between BD and healthy control samples (Supplementary Table 1). Such a large number of modulated transcripts clearly reflect the high performance of the array in the detection of a wide range of microRNA sequences. We thereafter sharpened our analysis by selecting only modulated miRNAs annotated as “high confidence” in miRBase 21 (http://www.mirbase.org), and to make our results more informative, we further narrowed down the analysis to miRNAs for which gene targets were annotated in FunRich.

By these criteria, we selected 47 modulated miRNAs that are shown in Table 2. Interestingly, almost all (45/47) miRNAs were down-modulated with only two up-regulated microRNAs.

Selected miRNAs significantly deregulated in the microarray analysis were validated by real-time PCR in the entire series of patients analyzed (see Figure 1).

(a)

(b)

(c)

3.2. Pathway Enrichment Analysis of miRNAs Deregulated in Behçet’s Disease

In a second part of our analysis, we wanted to identify all the molecular pathways that were targeted by the selected miRNA performing a pathway enrichment analysis based on annotated gene targets in FunRich. The software allows to evaluate the miRNA regulatory effect and to identify controlled pathways based on predicted and/or validated miRNA-target interactions. We also applied a pathway enrichment analysis on the dataset of differentially expressed genes (DEGs) obtained in our previous study of gene expression profiling in BD (Puccetti et al., unpublished observations). In this study, we were able to select modulated genes that may play an important role in BD pathogenesis since they are involved in biological processes strongly connected to the typical features of the disease.

Despite the strong statistical stringency applied to the two datasets (), we obtained a high number of significantly enriched pathways both from the miRNA target datasets (277) and from the BD DEGs (164) (Supplementary Tables 2 and 3). Notably, we found that a large proportion of pathways from the BD DEGs was also enriched in the list of miRNA-validated targets (64%, 106/164; Supplementary Table 4), indicating that the selected miRNAs exert a strong impact on the molecular pathways altered in the disease. Moreover, the large number of enriched pathways clearly reflected the multisystemic involvement typical of Behçet’s disease. Figure 2 shows a graphical representation of selected commonly enriched pathways. Interestingly, the enriched categories were involved in vascular biology (i.e., glypican pathway, vascular endothelial growth factor, VEGF and VEGFR network, endothelin pathway, PAR1-mediated thrombin signaling events, thrombin/protease-activated receptor (PAR) pathway, EGFR-dependent endothelin signaling events, platelet-derived growth factor (PDGF) receptor signaling network, and urokinase-type plasminogen activator (uPA), and uPAR-mediated signaling) and in apoptosis (TRAIL, p53, and FAS signaling pathways). In addition, other relevant pathways enriched in the two datasets were related to the immune response (i.e., IL6-mediated signaling events, TCR and BCR signaling, calcineurin-regulated NFAT-dependent transcription in lymphocytes, and toll-like receptor cascades) and to the inflammatory response (i.e., tumour necrosis factor, TNF receptor, IFN-gamma, p38 MAP kinase, pathway, and IL1- and CXCR4-mediated signaling events.

3.3. Comparative Analysis of Selected miRNA Gene Targets and Differentially Expressed Genes in BD

To better define the role played by miRNAs in BD pathogenesis, we wanted to select miRNAs that were able to target genes modulated in BD. Therefore, we used a more sophisticated integrative database for human microRNA target predictions (mirDIP) [24] to obtain the lists of genes that were targeted, with a very high score class, by each of the selected miRNAs. Then, we compared the resulting target lists to genes differentially expressed in BD from our previous study (Puccetti et al., unpublished observations) and we observed that 65% of DEGs were targeted by the selected miRNAs. Interestingly, the vast majority of these DEGs showed an opposite modulation with respect to the relative targeting miRNA (Figure 3(a)), consistently with the typical role of miRNAs as negative regulator of gene expression (i.e., typically, up-regulated genes are targeted by down-modulated miRNAs and underexpressed genes are targeted by overexpressed miRNAs). All the above-mentioned miRNAs are presented in Figure 3(b) along with the compiled lists of their targets DEGs.

(a)

(b)

All these selected miRNAs targeted genes involved in biological processes implicated in the disease pathogenesis including apoptosis, immune response, inflammation, and vascular damage (Figure 4). Thus, we could identify microRNAs that may control the gene modulation involved in the disease pathogenesis. Table 3 shows all the targeted genes and their corresponding targeting miRNAs.

Selected miRNAs targeted many DEGs that sustained the inflammatory response typically associated with BD including TNF, IL1A, IL10, IL6R, CXCL2, CXCR4, TNFAIP3, OLR1, S100A8, HSP90B1, and CCL3 (see Figure 3(b) and Table 3). Interestingly, TNF (targeted by hsa-miR-181d-5p and -181a-5p), CCL3 (targeted by hsa-let-7e-5p, -7d-5p, and -7f-5p), IL10 (targeted by hsa-let-7f-5p, -7e-5p, -7d-5p, -miR-146a-5p, and -27b-3p), and IL1A (targeted by hsa-miR-505-3p, -181d-5p, and -181a-5p) have been detected at increased concentrations in sera or plasma of BD patients when compared to normal subjects [29–31]. Moreover, down-modulated miRNAs targeted DEGs involved in the adaptive immune response including genes that played a role in T and in B cell immune response. Among these genes, several miRNAs targeted CD28 (hsa-miR-143-3p, -27b-3p, -15b-5p, and -195-5p), CD4 (hsa-miR-181a-5p), ICOS (hsa-miR 27b-3p and -let-7f-5p, -7e-5p, and -7d-5p), CTLA4 (hsa-miR-143-3p), EGR1 (hsa-miR-143-3p, -146a-5p, -181a-5p, -199b-3p, and -199a-3p), and AKIRIN2 (hsa-miR-139-5p, -27b-3p, and -181a-5p). Interestingly, 18 miRNAs targeted genes of the TH-17 gene signature including CD28, CD4, ICOS, SOCS3, IL6ST, YY1, IL6R, TNF, CXCL2, and SOCS1 (Figure 3(b)).

In addition, selected miRNAs could also control DEGs involved in the innate immune response including, for example, NKTR (hsa-miR-181a-5p, -486-5p, -151-3p, and -27b-3p), KIR2DL4 (hsa-miR-146a-5p), STAT1 (hsa-miR-584-5p, -27b-3p, -128-3p, and -146a-5p), STAT2 (hsa-miR-143-3p), and DEGs belonging to the Toll-like receptor (TLR) pathway such as TLR2, IKBKB, JUN, REL, NAMPT, HSP90B1, HSPA4, and ARF3. In particular, we have to mention that TLR2 (targeted by hsa-miR-143-3p and -146a-5p) is thought to be up-regulated in BD patients [32, 33]. As many as 22 miRNAs targeted genes of the type I interferon signature like VEGFA, DDX3X, and the above-mentioned STAT1, STAT2, SOCS1, and IL10. The coexpression of miRNAs targeting DEGs involved in the type I IFN signaling and Th-17 related genes may reflect the presence of a synergy between IFN and Th17 pathways that is typical of autoimmune diseases [34–43], thus suggesting that an autoimmune mechanism can be involved in the BD pathogenesis.

Noteworthily, several miRNAs targeted all the eight genes involved in the JAK/STAT signaling pathway that we found up-regulated in our previous work (Puccetti et al., unpublished observations) including PIK3R1 and the above-mentioned STAT1, STAT2, IL6ST, IL6R, IL10, SOCS1, and SOCS3. The downmodulation of miRNAs targeting this pathway further supports the hypothesis of an autoimmune origin of BD since its activation is very frequently associated with autoimmune diseases [44]. Moreover, the JAK/STAT signaling pathway is active in CD4⁺ T cells of patients with BD [45].

Interestingly, modulated miRNAs targeted DEGs involved in the vascular damage associated with BD that is characterized by myointimal proliferation, fibrosis, and thrombus formation [46]. Indeed, several miRNAs targeted DEGs associated with blood coagulation (i.e., THBS1, F5, and LMAN1; see Table 3), a process whose alteration is typically associated with BD vasculitis. Noteworthily, 10 down-modulated miRNAs (hsa-let-7d/7e/7f-5p, hsa-miR-151-3p, -628-3p, -139-5p, 143-3p, -194-5p, -199a-3p, and -199b-3p) targeted THBS1, suggesting a crucial role of THBS1 in the pathogenesis of vascular damage. Other modulated miRNAs targeted DEGs that played a role in angiogenesis including MMP8, VEGFA, NR4A3, and NAPA. Moreover, NR4A3 is a transcription factor involved in vascular development [47] and NAPA promotes vascular endothelial- (VE-) cadherin localization at endothelial junctions [48]. Interestingly, increased serum levels of soluble VE-cadherins have been detected in BD patients [49]. Among the above-mentioned genes, we have to notice that VEGFA was predictively targeted by a high number (9) of down-modulated miRNAs, including hsa-miR-199a-3p, -15b-5p, -195-5p, -361-5p, -126-3p, -486-5p, -199a-5p, 199b-3p, and hsa-miR-139-5p. Other DEGs involved in vascular damage and targeted by modulated miRNAs were FOSL2, THBD, and PTX3. THBD, targeted by hsa-miR-139-5p, is increased in sera of BD patients compared with healthy subjects, and PTX3, targeted by hsa-miR-628-3p, is considered as a marker of small vessel vasculitis [50].

Moreover, DEGs involved in the apoptotic process were targeted, and among these, we can cite MCL1, DNAJB1, BCL2L11, ARHGDIA, ZNF331, and IER3 (Table 3). Finally, several modulated miRNAs targeted DEGs that played a role in cell proliferation, a process that can be induced in response to both apoptosis and to skin ulcer formation. Among these miRNAs, we can mention those that targeted BTG2, EREG, PRRC2, and APLP2. In particular, EREG (targeted by hsa-miR-199a-3p, -199b-3p, and -181a-5p and by hsa-miR-192-5p) and APLP2 (targeted by hsa-miR-199a-3p, -199b-3p, -139-5p, and -181a-5p) were involved in the proliferation of corneal epithelial cells and in corneal epithelial wound healing, respectively [51, 52]. Thus, this gene regulation may have a role in ocular manifestations of BD like keratitis.

3.4. Network Analysis of Genes Differentially Expressed in BD

We performed a network analysis in which the functional interactions between the protein products of modulated genes in BD were evaluated. By this approach, a protein-protein interaction (PPI) network comprising 171 genes (nodes) and 3272 pairs of interactions (edges) was constructed (Figure 5(a)). We then performed a clustering analysis to identify areas of densely interconnected nodes (clusters/modules; CL) that are predicted to be involved in common biological processes and to have a crucial role in the disease pathogenesis. We could detect six clusters that are graphically represented in Figures 5(b)–5(g). Comparing the list of miRNA targets to DEGs included in the six clusters (Supplementary Table 5), we found that, in each cluster, a large number of DEGs were targeted (Supplementary Table 6 and Figures 5(b)–5(g)). In particular, we observed that several of these genes were involved in immune and inflammatory responses including TNF (CL2), IL10 (CL2), and IL1A (CL1). Moreover, many of such genes played a role in B cell response, like EGR1 (CL1), or in T cell response like, for example, CTLA4, CD4, and MLL (CL1); DUSP2 (CL3); and CD28 (CL4) and NR4A2 (CL5). Other DEGs targeted in clusters were involved in TLR signaling including HSPA4, IKBKB, JUN, REL, S100A8, TLR2 (CL1), and ARF3 and NAMPT (CL2). Interestingly, in the six clusters, we also observed that various genes involved in angiogenesis and/or in vascular damage were targeted by selected miRNAs including THBS1, MMP8, VEGFA (CL1), FOSL2, THBD (CL2), HIPK1 and DDX6 (CL3), PGK1 (CL4), and ACTR2 (CL6).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Given the well-known biological significance of highly connected gene clusters, to gain insight on the most relevant signaling networks that were controlled by the deregulated miRNAs, we performed a pathway enrichment analysis on targeted genes that were present in the six clusters. All the above-mentioned enriched pathways are listed in the Supplementary Table 7, and Table 4 shows a selection of the most relevant enriched signaling networks.

The TLR and the JAK/STAT signaling networks, two pathways notoriously involved in immune response and frequently associated with autoimmunity, were enriched in CL1-targeted genes (CL1-TGs). Other pathways implicated in the immune response were enriched in several genes targeted in clusters including T cell activation (CL1-TGs, CL2-TGs, and CL4-TGs), B cell activation (CL1-TGs), and PI3K signaling (CL2-TGs and CL3-TGs). This pathway is a crucial element in the regulation of adaptive and innate immune response [53, 54] and is a key player in inflammatory response. Therefore, in recent years, highly selective inhibitors of PI3K have been developed for anti-inflammatory treatments [55]. An enrichment in signalings involved in the inflammatory response was also found in CL1-TGs (inflammation mediated by chemokine and cytokine and interleukin pathways) and in CL5-TGs (integrin pathway).

In CL1-TGs and in CL3-TGs, the oxidative stress response pathway was enriched. Interestingly, an increased oxidative stress has been described in BD and it has been correlated to the severe inflammatory and degenerative clinical manifestations of the disease [56]. Several networks clearly connected to vascular damage were enriched in cluster TGs including PDGF, angiogenesis (CL1-TGs and CL5-TGs), Wnt, blood coagulation, endothelin, VEGF (CL2-TGs), and cadherin (CL6-TGs) signaling pathway. VE-cadherins have a crucial role in endothelial barrier integrity, and noteworthily, in CL6 the ACTR2 gene was targeted, known to stabilize adherens junctions between endothelial cells of the vascular wall interacting with cadherins [57]. Moreover, the Wnt signaling pathway is critically involved in vascular biology [58]. Finally, in CL1-TGs and in CL2-TGs the apoptosis and the p53 pathways were enriched, respectively.

4. Discussion

A systematic analysis of miRNA expression profiles in BD has not been performed yet. The aim of our work is therefore to provide a compiled description of miRNAs associated with BD using an array able to query a very large number of transcripts in the attempt to dissect the possible regulatory effects exerted by these molecules on the molecular pathways relevant for the disease pathogenesis.

In a previous work, we investigated BD-associated transcriptional profiles by a gene expression analysis of PBMC derived from BD patients and identified a gene modulation strictly connected to the disease pathogenesis (Puccetti et al., unpublished observations). Moreover, we showed the presence of a type I interferon and Th-17 gene signature, which suggests an autoimmune origin of BD. In this study, we aimed to complement this gene expression analysis detecting modulated miRNAs that may target differentially expressed genes (DEGs) identified in our previous analysis.

To this purpose, we used a sophisticated target prediction system (mirDIP) [24] to obtain the list of miRNA targets, which was then compared to the list of DEGs identified in blood samples of BD patients. The comparison of the two analyses (gene expression and miRNAs) allows selecting within deregulated miRNAs only those related to the modulation of genes that would be effectively altered in the PBMCs of BD patients and might have therefore a pivotal role in the pathogenesis of the disease.

We observed that there was a good overlap (60%) between the selected miRNA targets and the BD-modulated genes, indicating that the majority of the identified miRNAs regulate genes differentially expressed in BD. Such target genes belonged to functional classes strictly connected to typical features of BD. Indeed, several miRNAs targeted proinflammatory genes such as TNF and IL1A. Moreover, transcripts involved in both adaptive and innate immune response were also targeted. In this regard, we have to mention that selected down-modulated miRNAs controlled several Th17 cell-associated genes and transcripts involved in type I interferon response that we found up-regulated in our previous analysis. This may suggest a loss of control in BD on two synergistic mechanisms typically associated with an autoimmune response. We moreover observed a downregulation of miRNAs that control members of TLRs and JAK/STAT pathways, two molecular signalings involved in autoimmune diseases [44, 59] that are also active in BD [32, 60, 61]. In addition, miRNAs target genes involved in angiogenesis, and in blood coagulation, two processes commonly associated with BD vasculitis were also down-modulated in BD samples.

Consistently with the tight correspondence between deregulated miRNA and DEGS in BD, we observed a good overlap (64%) between the pathways enriched in the deregulated miRNA targets and in the genes differently expressed in BD. This indicated that, globally, the pathways identified by the selected miRNAs reflected fairly well the gene regulation described in our gene expression study, suggesting that the criteria applied for the selection of miRNAs had effectively identified miRNAs able to explain the gene modulation which we had previously described.

Interestingly, we found that both in BD-miRNA and in BD-DEGs dataset, meaningful signaling networks including, for example, VEGF and VEGFR, endothelins, PDGF receptor, TNF receptor, and IL1- and IL6-mediated TCR and BCR pathways were enriched.

A pivotal task in molecular biology is to understand gene modulation in the context of biological networks. Indeed, proteins achieve their functions in the protein-protein-interacting network, and interestingly, the dynamics of such networks can be influenced by miRNAs [62]. We therefore wanted to identify significant relationships among modulated miRNAs and the PPI network in which the protein product of modulated genes in BD can be involved. Moreover, since it is known that the repressive effect of a miRNA may lead to more severe biological effects when it is exerted on proteins with more interacting partners [62], to highlight more efficient interrelations, we focused our attention on the most connected proteins of the network, inspecting the presence of miRNA targets inside the six clusters extracted from the PPI network. We found that the six clusters identified were extensively targeted by several modulated miRNAs, thus indicating that the deregulation of selected miRNAs may have a meaningful effect on the dynamics of a protein network that control the disease pathogenesis. Indeed, genes targeted in the six clusters played an important role in vascular biology (i.e., VEGFA, THBS, THBS1, etc.), in inflammation (i.e., TNF, IL1A, IL6R, CXCR4, etc.), and in immune response (i.e., CD28, CTLA4, EGR1, TLR2, etc.). Moreover, the analysis of pathways that were enriched in the target genes present in the clusters confirmed the essential role of these transcripts and their relative targeting miRNA in BD pathogenesis.

In conclusion, this work represents the first analysis performed on such a large number of miRNAs and integrated with the study of the profiles of gene expression in BD. The study allowed correlating the expression of miRNAs and the modulation of genes important for the pathogenesis of the disease. Using this approach, we have been able to identify the specific molecular pathways on which the regulation of these miRNAs may occur.

This study sheds light on some epigenetic aspects of BD identifying specific miRNAs, which may represent promising candidates for the identification of disease biomarkers and/or the design of novel therapeutic strategies in BD.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Antonio Puccetti, Claudio Lunardi, and Marzia Dolcino conceived and designed the experiments. Piera Filomena Fiore and Andrea Pelosi performed the experiments. Marzia Dolcino and Andrea Pelosi analyzed the data. Giuseppe Patuzzo selected the patients and contributed reagents. Marzia Dolcino wrote the paper with inputs from Claudio Lunardi and Antonio Puccetti. Antonio Puccetti, Andrea Pelosi, Claudio Lunardi, and Marzia Dolcino equally contributed to this paper.

Supplementary Materials

Supplementary 1. Supplementary Table 1: miRNAs significantly modulated in BD.

Supplementary 2. Supplementary Table 2: pathways enriched in genes targeted by modulated miRNAs in BD.

Supplementary 3. Supplementary Table 3: pathways enriched in differentially expressed genes in BD.

Supplementary 4. Supplementary Table 4: pathways enriched both in genes targeted by BD deregulated miRNAs and in BD differentially expressed genes.

Supplementary 5. Supplementary Table 5: differentially expressed genes in BD that were included in the six clusters.

Supplementary 6. Supplementary Table 6: differentially expressed genes in BD that were included in clusters that are targeted by modulated miRNAs in BD.

Supplementary 7. Supplementary Table 7: enriched pathways in target genes that are included in clusters.

References

International Team for the Revision of the International Criteria for Behçet’s Disease (ITR-ICBD), F. Davatchi, S. Assaad-Khalil et al., “The International Criteria for Behçet’s Disease (ICBD): a collaborative study of 27 countries on the sensitivity and specificity of the new criteria,” Journal of the European Academy of Dermatology and Venereology, vol. 28, no. 3, pp. 338–347, 2014.
View at: Publisher Site | Google Scholar
E. Alpsoy, “Behçet’s disease: a comprehensive review with a focus on epidemiology, etiology and clinical features, and management of mucocutaneous lesions,” The Journal of Dermatology, vol. 43, no. 6, pp. 620–632, 2016.
View at: Publisher Site | Google Scholar
L. T. Morton, D. Situnayake, and G. R. Wallace, “Genetics of Behçet’s disease,” Current Opinion in Rheumatology, vol. 28, no. 1, pp. 39–44, 2016.
View at: Publisher Site | Google Scholar
N. Mizuki, A. Meguro, M. Ota et al., “Genome-wide association studies identify IL23R-IL12RB2 and IL10 as Behçet’s disease susceptibility loci,” Nature Genetics, vol. 42, no. 8, pp. 703–706, 2010.
View at: Publisher Site | Google Scholar
E. F. Remmers, F. Cosan, Y. Kirino et al., “Genome-wide association study identifies variants in the MHC class I, IL10, and IL23R-IL12RB2 regions associated with Behçet’s disease,” Nature Genetics, vol. 42, no. 8, pp. 698–702, 2010.
View at: Publisher Site | Google Scholar
I. Fresko, “15th international congress on Behçet’s disease,” Clinical and Experimental Rheumatology, vol. 30, Supplement 72, no. 3, pp. S118–S128, 2012.
View at: Google Scholar
M. Takeuchi, D. L. Kastner, and E. F. Remmers, “The immunogenetics of Behçet’s disease: a comprehensive review,” Journal of Autoimmunity, vol. 64, pp. 137–148, 2015.
View at: Publisher Site | Google Scholar
M. Takeuchi, N. Mizuki, A. Meguro et al., “Dense genotyping of immune-related loci implicates host responses to microbial exposure in Behçet’s disease susceptibility,” Nature Genetics, vol. 49, no. 3, pp. 438–443, 2017.
View at: Publisher Site | Google Scholar
H. Direskeneli, H. Fujita, and C. A. Akdis, “Regulation of T_H17 and regulatory T cells in patients with Behçet disease,” Journal of Allergy and Clinical Immunology, vol. 128, no. 3, pp. 665-666, 2011.
View at: Publisher Site | Google Scholar
N. Sugi-Ikai, M. Nakazawa, S. Nakamura, S. Ohno, and M. Minami, “Increased frequencies of interleukin-2- and interferon-gamma-producing T cells in patients with active Behçet’s disease,” Investigative Ophthalmology & Visual Science, vol. 39, no. 6, pp. 996–1004, 1998.
View at: Google Scholar
K. Hamzaoui, “Th17 cells in Behçet’s disease: a new immunoregulatory axis,” Clinical and Experimental Rheumatology, vol. 29, Supplement 67, no. 4, pp. S71–S76, 2011.
View at: Google Scholar
G. Geri, B. Terrier, M. Rosenzwajg et al., “Critical role of IL-21 in modulating T_H17 and regulatory T cells in Behçet disease,” Journal of Allergy and Clinical Immunology, vol. 128, no. 3, pp. 655–664, 2011.
View at: Publisher Site | Google Scholar
S. Y. Na, M. J. Park, S. Park, and E. S. Lee, “Up-regulation of Th17 and related cytokines in Behçet’s disease corresponding to disease activity,” Clinical and Experimental Rheumatology, vol. 31, Supplement 77, no. 3, pp. 32–40, 2013.
View at: Google Scholar
J. Kim, J. A. Park, E. Y. Lee, Y. J. Lee, Y. W. Song, and E. B. Lee, “Imbalance of Th17 to Th1 cells in Behçet’s disease,” Clinical and Experimental Rheumatology, vol. 28, Supplement 60, no. 4, pp. S16–S19, 2010.
View at: Google Scholar
M. C. Mat, A. Sevim, I. Fresko, and Y. Tuzun, “Behçet’s disease as a systemic disease,” Clinics in Dermatology, vol. 32, no. 3, pp. 435–442, 2014.
View at: Publisher Site | Google Scholar
A. Mehta and D. Baltimore, “MicroRNAs as regulatory elements in immune system logic,” Nature Reviews Immunology, vol. 16, no. 5, pp. 279–294, 2016.
View at: Publisher Site | Google Scholar
Y. Bronevetsky and K. M. Ansel, “Regulation of miRNA biogenesis and turnover in the immune system,” Immunological Reviews, vol. 253, no. 1, pp. 304–316, 2013.
View at: Publisher Site | Google Scholar
J. T. Mendell and E. N. Olson, “MicroRNAs in stress signaling and human disease,” Cell, vol. 148, no. 6, pp. 1172–1187, 2012.
View at: Publisher Site | Google Scholar
A. Esquela-Kerscher and F. J. Slack, “Oncomirs - microRNAs with a role in cancer,” Nature Reviews Cancer, vol. 6, no. 4, pp. 259–269, 2006.
View at: Publisher Site | Google Scholar
K. M. Pauley, S. Cha, and E. K. L. Chan, “MicroRNA in autoimmunity and autoimmune diseases,” Journal of Autoimmunity, vol. 32, no. 3-4, pp. 189–194, 2009.
View at: Publisher Site | Google Scholar
Z. Qu, W. Li, and B. Fu, “MicroRNAs in autoimmune diseases,” BioMed Research International, vol. 2014, Article ID 527895, 8 pages, 2014.
View at: Publisher Site | Google Scholar
H. Yazici and Y. Yazici, “Criteria for Behçet’s disease with reflections on all disease criteria,” Journal of Autoimmunity, vol. 48-49, pp. 104–107, 2014.
View at: Publisher Site | Google Scholar
F. Davatchi, C. Chams-Davatchi, H. Shams et al., “Behcet’s disease: epidemiology, clinical manifestations, and diagnosis,” Expert Review of Clinical Immunology, vol. 13, no. 1, pp. 57–65, 2017.
View at: Publisher Site | Google Scholar
T. Tokar, C. Pastrello, A. E. M. Rossos et al., “MirDIP 4.1—integrative database of human microRNA target predictions,” Nucleic Acids Research, vol. 46, no. D1, pp. D360–D370, 2018.
View at: Publisher Site | Google Scholar
M. Pathan, S. Keerthikumar, C. S. Ang et al., “FunRich: an open access standalone functional enrichment and interaction network analysis tool,” Proteomics, vol. 15, no. 15, pp. 2597–2601, 2015.
View at: Publisher Site | Google Scholar
H. Mi and P. Thomas, “Panther pathway: an ontology-based pathway database coupled with data analysis tools,” Methods in Molecular Biology, vol. 563, pp. 123–140, 2009.
View at: Publisher Site | Google Scholar
A. Franceschini, D. Szklarczyk, S. Frankild et al., “String v9.1: protein-protein interaction networks, with increased coverage and integration,” Nucleic Acids Research, vol. 41, Database issue 1, pp. D808–D815, 2013.
View at: Publisher Site | Google Scholar
G. D. Bader and C. W. Hogue, “An automated method for finding molecular complexes in large protein interaction networks,” BMC Bioinformatics, vol. 4, no. 1, p. 2, 2003.
View at: Publisher Site | Google Scholar
B. Turan, H. Gallati, H. Erdi, A. Gürler, B. A. Michel, and P. M. Villiger, “Systemic levels of the T cell regulatory cytokines IL-10 and IL-12 in Bechçet’s disease; soluble TNFR-75 as a biological marker of disease activity,” The Journal of Rheumatology, vol. 24, no. 1, pp. 128–132, 1997.
View at: Google Scholar
H. T. E. Ozer, E. Erken, R. Gunesacar, and O. Kara, “Serum RANTES, MIP-1α, and MCP-1 levels in Behçet’s disease,” Rheumatology International, vol. 25, no. 6, pp. 487-488, 2005.
View at: Publisher Site | Google Scholar
N. Duzgun, E. Ayaslioglu, H. Tutkak, and O. T. Aydintug, “Cytokine inhibitors: soluble tumor necrosis factor receptor 1 and interleukin-1 receptor antagonist in Behçet’s disease,” Rheumatology International, vol. 25, no. 1, pp. 1–5, 2005.
View at: Publisher Site | Google Scholar
N. Seoudi, L. A. Bergmeier, E. Hagi-Pavli, D. Bibby, M. A. Curtis, and F. Fortune, “The role of TLR2 and 4 in Behçet’s disease pathogenesis,” Innate Immunity, vol. 20, no. 4, pp. 412–422, 2014.
View at: Publisher Site | Google Scholar
Y. Liu, H. Yin, M. Zhao, and Q. Lu, “TLR2 and TLR4 in autoimmune diseases: a comprehensive review,” Clinical Reviews in Allergy and Immunology, vol. 47, no. 2, pp. 136–147, 2014.
View at: Publisher Site | Google Scholar
R. A. Gordon, G. Grigoriev, A. Lee, G. D. Kalliolias, and L. B. Ivashkiv, “The interferon signature and STAT1 expression in rheumatoid arthritis synovial fluid macrophages are induced by tumor necrosis factor α and counter-regulated by the synovial fluid microenvironment,” Arthritis & Rheumatism, vol. 64, no. 10, pp. 3119–3128, 2012.
View at: Publisher Site | Google Scholar
A. Nzeusseu Toukap, C. Galant, I. Theate et al., “Identification of distinct gene expression profiles in the synovium of patients with systemic lupus erythematosus,” Arthritis & Rheumatism, vol. 56, no. 5, pp. 1579–1588, 2007.
View at: Publisher Site | Google Scholar
R. M. Thurlings, M. Boumans, J. Tekstra et al., “Relationship between the type I interferon signature and the response to rituximab in rheumatoid arthritis patients,” Arthritis & Rheumatism, vol. 62, no. 12, pp. 3607–3614, 2010.
View at: Publisher Site | Google Scholar
H. G. Raterman, S. Vosslamber, S. de Ridder et al., “The interferon type I signature towards prediction of non-response to rituximab in rheumatoid arthritis patients,” Arthritis Research & Therapy, vol. 14, no. 2, article R95, 2012.
View at: Publisher Site | Google Scholar
F. Moschella, G. F. Torelli, M. Valentini et al., “Cyclophosphamide induces a type I interferon–associated sterile inflammatory response signature in cancer patients’ blood cells: implications for cancer chemoimmunotherapy,” Clinical Cancer Research, vol. 19, no. 15, pp. 4249–4261, 2013.
View at: Publisher Site | Google Scholar
N. I. Maria, Z. Brkic, M. Waris et al., “MxA as a clinically applicable biomarker for identifying systemic interferon type I in primary Sjögren’s syndrome,” Annals of the Rheumatic Diseases, vol. 73, no. 6, pp. 1052–1059, 2014.
View at: Publisher Site | Google Scholar
G. Caignard, M. Lucas-Hourani, K. P. Dhondt et al., “The V protein of Tioman virus is incapable of blocking type I interferon signaling in human cells,” PLoS One, vol. 8, no. 1, article e53881, 2013.
View at: Publisher Site | Google Scholar
A. Ambrosi, A. Espinosa, and M. Wahren-Herlenius, “Il-17: a new actor in IFN-driven systemic autoimmune diseases,” European Journal of Immunology, vol. 42, no. 9, pp. 2274–2284, 2012.
View at: Publisher Site | Google Scholar
R. C. Axtell, B. A. de Jong, K. Boniface et al., “T helper type 1 and 17 cells determine efficacy of interferon-β in multiple sclerosis and experimental encephalomyelitis,” Nature Medicine, vol. 16, no. 4, pp. 406–412, 2010.
View at: Publisher Site | Google Scholar
Z. Brkic, O. B. Corneth, C. G. van Helden-Meeuwsen et al., “T-helper 17 cell cytokines and interferon type I: partners in crime in systemic lupus erythematosus?” Arthritis Research & Therapy, vol. 16, no. 2, article R62, 2014.
View at: Publisher Site | Google Scholar
K. Hirahara, D. Schwartz, M. Gadina, Y. Kanno, and J. J. O'Shea, “Targeting cytokine signaling in autoimmunity: back to the future and beyond,” Current Opinion in Immunology, vol. 43, pp. 89–97, 2016.
View at: Publisher Site | Google Scholar
A. Tulunay, M. G. Dozmorov, F. Ture-Ozdemir et al., “Activation of the jak/stat pathway in Behcet’s disease,” Genes & Immunity, vol. 16, no. 2, p. 176, 2015.
View at: Publisher Site | Google Scholar
N. Maruotti, F. P. Cantatore, B. Nico, A. Vacca, and D. Ribatti, “Angiogenesis in vasculitides,” Clinical and Experimental Rheumatology, vol. 26, no. 3, pp. 476–483, 2008.
View at: Google Scholar
Y. Zhao and D. Bruemmer, “NR4A orphan nuclear receptors: transcriptional regulators of gene expression in metabolism and vascular biology,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 30, no. 8, pp. 1535–1541, 2010.
View at: Publisher Site | Google Scholar
A. V. Andreeva, M. A. Kutuzov, R. Vaiskunaite et al., “Gα₁₂ interaction with αSNAP induces VE-cadherin localization at endothelial junctions and regulates barrier function,” The Journal of Biological Chemistry, vol. 280, no. 34, pp. 30376–30383, 2005.
View at: Publisher Site | Google Scholar
T. Chen, Z. P. Guo, N. Cao, S. Qin, M. M. Li, and R. Z. Jia, “Increased serum levels of soluble vascular endothelial-cadherin in patients with systemic vasculitis,” Rheumatology International, vol. 34, no. 8, pp. 1139–1143, 2014.
View at: Publisher Site | Google Scholar
F. Fazzini, G. Peri, A. Doni et al., “PTX3 in small-vessel vasculitides: an independent indicator of disease activity produced at sites of inflammation,” Arthritis & Rheumatism, vol. 44, no. 12, pp. 2841–2850, 2001.
View at: Publisher Site | Google Scholar
S. Morita, Y. Shirakata, A. Shiraishi et al., “Human corneal epithelial cell proliferation by epiregulin and its cross-induction by other EGF family members,” Molecular Vision, vol. 13, pp. 2119–2128, 2007.
View at: Google Scholar
J. Guo, G. Thinakaran, Y. Guo, S. S. Sisodia, and F. X. Yu, “A role for amyloid precursor-like protein 2 in corneal epithelial wound healing,” Investigative Ophthalmology & Visual Science, vol. 39, no. 2, pp. 292–300, 1998.
View at: Google Scholar
D. A. Fruman, H. Chiu, B. D. Hopkins, S. Bagrodia, L. C. Cantley, and R. T. Abraham, “The PI3K pathway in human disease,” Cell, vol. 170, no. 4, pp. 605–635, 2017.
View at: Publisher Site | Google Scholar
J. Tsalikis, D. O. Croitoru, D. J. Philpott, and S. E. Girardin, “Nutrient sensing and metabolic stress pathways in innate immunity,” Cellular Microbiology, vol. 15, no. 10, pp. 1632–1641, 2013.
View at: Publisher Site | Google Scholar
A. K. Stark, S. Sriskantharajah, E. M. Hessel, and K. Okkenhaug, “PI3K inhibitors in inflammation, autoimmunity and cancer,” Current Opinion in Pharmacology, vol. 23, pp. 82–91, 2015.
View at: Publisher Site | Google Scholar
E. D. Sezer, K. Aksu, O. Caglayan, G. Keser, G. Karabulut, and G. Ercan, “DNA damage and its relationship with other oxidative stress parameters in Behcet’s disease,” Rheumatology International, vol. 32, no. 1, pp. 217–222, 2012.
View at: Publisher Site | Google Scholar
C. Rajput, V. Kini, M. Smith et al., “Neural Wiskott-Aldrich syndrome protein (N-WASP)-mediated p120-catenin interaction with Arp2-actin complex stabilizes endothelial adherens junctions,” Journal of Biological Chemistry, vol. 288, no. 6, pp. 4241–4250, 2013.
View at: Publisher Site | Google Scholar
P. Meyer, “Calcium and cellular responses in hypertension,” Journal of Cardiovascular Pharmacology, vol. 12, Supplement 6, pp. S36–S38, 1988.
View at: Publisher Site | Google Scholar
L. Duffy and S. O'Reilly, “Toll-like receptors in the pathogenesis of autoimmune diseases: recent and emerging translational developments,” ImmunoTargets and Therapy, vol. 5, pp. 69–80, 2016.
View at: Publisher Site | Google Scholar
Y. Kirino, M. Takeno, R. Watanabe et al., “Association of reduced heme oxygenase-1 with excessive toll-like receptor 4 expression in peripheral blood mononuclear cells in Behçet’s disease,” Arthritis Research & Therapy, vol. 10, no. 1, article R16, 2008.
View at: Publisher Site | Google Scholar
A. Tulunay, M. G. Dozmorov, F. Ture-Ozdemir et al., “Activation of the JAK/STAT pathway in Behcet’s disease,” Genes and Immunity, vol. 16, no. 2, pp. 170–175, 2015.
View at: Publisher Site | Google Scholar
H. Liang and W. H. Li, “MicroRNA regulation of human protein–protein interaction network,” RNA, vol. 13, no. 9, pp. 1402–1408, 2007.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2018 Antonio Puccetti et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

2698

Downloads

1342

Citations