Systems biology of myasthenia gravis, integration of aberrant lncRNA and mRNA expression changes
BMC Medical Genomics volume 8, Article number: 13 (2015)
A novel class of transcripts, long non-coding RNAs (lncRNAs), has recently emerged as a key player in several biological processes, and important roles for these molecules have been reported in a number of complex human diseases, such as autoimmune diseases, neurological disorders, and various cancers. However, the aberrant lncRNAs implicated in myasthenia gravis (MG) remain unknown. The aim of the present study was to explore the abnormal expression of lncRNAs in peripheral blood mononuclear cells (PBMCs) and examine mRNA regulatory relationship networks among MG patients with or without thymoma.
Microarray assays were performed, and the outstanding differences between lncRNAs or mRNA expression were verified through RT-PCR. The lncRNAs functions were annotated for the target genes using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) biological pathway. The potential regulatory relationships between the lncRNAs and target genes were analyzed using the ‘cis’ and ‘trans’ model. Outstanding lncRNAs were organized to generate a TF-lncRNA-gene network using Cytoscape software.
The lncRNA and mRNA expression profile analysis revealed subsets of differentially expressed genes in MG patients with or without thymoma. A total of 12 outstanding dysregulated expression lncRNAs, such as lncRNA oebiotech_11933, were verified through real-time PCR. Several GO terms including the cellular response to interferon-γ, platelet degranulation, chemokine receptor binding and cytokine interactions were very important in MG pathogenesis. The chromosome locations of some lncRNAs and associated co-expression genes were demonstrated using ‘cis’ analysis. The results of the ‘trans’ analysis revealed that some TFs (i.e., CTCF, TAF1and MYC) regulate lncRNA and gene expression. The outstanding lncRNAs in each group were implicated in the regulation of the TF-lncRNA-target gene network.
The results of the present study provide a perspective on lncRNA expression in MG. We identify a subset of aberrant lncRNAs and mRNAs as potential biomarkers for the diagnosis of MG. The GO and KEGG pathway analysis provides an annotation to determine the functions of these lncRNAs. The results of the ‘cis’ and ‘trans’ analyses provide information concerning the modular regulation of lncRNAs.
Myasthenia gravis (MG) is a T and B cell-mediated autoimmune disease of the neuromuscular junction; cytokines and chemokines may play a crucial role in the pathogenesis and perpetuation of this disease [1-3]. Muscle weakness and fatigue, as hallmarks of MG, indicate improper signaling between T and B cells, resulting in the elicitation of an antibody-mediated autoimmune response against the acetylcholine receptor (AChR) located at neuromuscular junctions [4,5]. In most cases, autoantibodies against the AChR are detected. Recently, other targets, such as muscle-specific kinase (MuSK) and lipoprotein-related protein 4 (LRP4), have been described . However, the mechanisms that contribute to the development or pathogenesis of MG are highly complex. A puzzling yet interesting characteristic of MG is that many patients show thymic abnormalities, such as thymic germinal center hyperplasia, whereas other patients develop thymic tumors [6-9]. The thymus has long been considered to be closely associated with the pathogenesis of MG, and some studies have suggested that MG autoimmune reactions originate in the thymus [7,10]. The relationship between MG and thymoma has been repeatedly suggested; however, many questions and controversies still remain. The biological and clinical behavior of thymomatous MG following thymectomy has not been adequately investigated, and evidence is lacking concerning both the clinical course of MG and the prognosis of thymomas associated with MG following surgical management [7,11]. Therefore, additional studies concerning the role of the thymus gland in the pathogenesis of MG are required.
Long non-coding RNAs (lncRNAs) are a group of RNA transcripts that are more than 200 nucleotides in length and lack significant open reading frames (ORFs) [12-15]. lncRNAs have only been recently identified, and most of the functions of these molecules remain unclear. The long nucleotide chain of lncRNAs can either form a complex spatial structure that interacts with protein factors or provide a large region for the concurrent binding of many molecules that collectively participate in genomic imprinting, X-chromosome silencing, chromosome modification, intranuclear transport, and transcriptional activation and interference [13,14,16-18],thereby regulating cell growth, differentiation, development, senescence and death [19-21]. Although several lncRNAs molecules have been implicated in diverse processes and diseases [22,23], only a few examples of the regulation of the immune system through lncRNAs have been described [24-26].
In the present study, we performed an array of lncRNA chip assays on PBMCs of MG patients. The outstanding lncRNAs functions were annotated based on co-expression genes and the GO biological process. The relationships among the lncRNAs were revealed through ‘cis’ and ‘trans’ analyses. These results provide information for further studies on MG.
lncRNA and mRNA expression profile in MG
To investigate the expression levels of lncRNAs associated with MG with or without thymoma, lncRNA and mRNA microarray analyses were performed on the PBMCs of MG patients. After separating the signal from noise and performing a t-test, significant differences in lncRNAs and mRNA expression of up to two-fold (P < 0.05 and FDR < 0.05) were observed. These results are summarized in Table 1, and detailed information is provided in Additional file 1. The lncRNA and mRNA expression data were clustered using Cluster3.0, as shown in Figure 1. Using dendrogram-based methods for clustering, the samples were further separated into two subgroups through hierarchical clustering based on similar expression patterns, and the results indicate that the expression of these lncRNAs and mRNAs was significantly different in MG patients with or without thymoma. This difference could potentially distinguish the disease group from MG patients without thymoma or healthy individuals. We also calculated the common upregulation or downregulation of lncRNA or mRNA expression in MG patients with or without thymoma versus healthy controls. A total of 42 lncRNAs exhibited upregulated expression, and 93 lncRNAs exhibited downregulated expression (Figure 1D), whereas 80 mRNAs exhibited upregulated expression and 43 mRNAs exhibited downregulated expression in the two groups (Figure 1E). The detailed data are provided in Additional file 2. These results suggest that these lncRNAs and mRNAs may have common functions in the development of MG.
The lncRNA function annotation and statistical characteristics network
To further explore the function of lncRNAs in MG, the results of the lncRNA and mRNA chip analyses were subjected to Pearson’s correlation coefficient analysis, in which co-expression was considered at P > 0.8. To determine the level of lncRNA and mRNA co-expression, we divided the different lncRNAs into two subsets: the upregulated lncRNAs and downregulated lncRNAs. When the subset contained more than 100 different lncRNAs, we selected the top 100 most distinguished lncRNAs, and when the subset contained less than 100 lncRNAs, all lncRNAs in the subset were used. Because the file is extremely large, the information pertaining to lncRNA oebiotech_11933 is shown as a representative result. The results for the other lncRNAs are provided in Additional file 3.
lncRNA oebiotech_11933 exhibited the highest upregulation of lncRNA expression among MG patients with thymoma versus healthy controls. A total of 552 genes(e.g., G0S2) relative to lncRNA oebiotech_11933 exhibited P > 90% (Table 2). These genes were further analyzed using DAVID functional annotation software (http://david.abcc.ncifcrf.gov/ gene2gene.jsp). The lncRNA oebiotech_11933 function was annotated using GO and KEGG pathway analyses. Selecting the reliability prediction terms(according to the P-value and enrichment), a total of 20 enrichment GO terms were obtained. Table 3 indicates that lncRNA oebiotech_11933 was associated with transcriptional regulation, cell surface receptor-linked signal transduction, and immune responses, among other terms.
Moreover, lncRNA oebiotech_11933 KEGG pathways are listed in Table 4, including ‘cytokine-cytokine receptor interaction’, ‘MAPK signaling pathway’, ‘chemokine signaling pathway’, ‘NOD-like receptor signaling pathway’, and ‘Toll-like receptor signaling pathway’. Previous studies have reported that these pathways are associated with lymphocyte immune cell proliferation and cancer development pathways [27-30]. The results of the GO and KEGG pathway analyses consistently showed that lncRNA oebiotech_11933 is associated with lymphocyte immune cell proliferation and cancer development pathways.
Furthermore, the lncRNA co-expression genes in each group were also analyzed using DAVID functional annotation software. The results of the GO and KEGG pathway analyses for lncRNA oebiotech_11933 are presented. In addition, the aberrant lncRNA genes were subjected to GO and KEGG pathway analyses. We selected the top 100 and 200 reliability prediction terms(according to the P-value and enrichment) for co-expressed and aberrant lncRNA genes, respectively (Figure 2). In MG patients with thymoma versus healthy controls, the top 100 terms in the GO and KEGG pathway analyses were associated with ‘cellular response to interferon-γ’, ‘positive regulation of cytokine production’, ‘smooth muscle cell proliferation’, and ‘cytokine receptors’, among other terms (Figure 2A). The top 200 terms in the GO and KEGG pathway analyses were associated with cellular responses to interferon-γ (Figure 2B). The 7 most common GO terms are shown in Figures 2A and 2B. The terms ‘cellular response to interferon-γ’, ‘chemokine receptor binding’ and ‘positive regulation of cytokine production’ were the most enriched GO terms. In MG patients without thymoma versus healthy controls, GOCC platelet alpha granule and GOMF chemokine receptor binding were the most enriched among the top 200 GO terms. A total of 14 common GO terms are shown in the Figures 2C and 2D. ‘chemokine receptor binding’, ‘cytokine-cytokine receptor interaction’, and ‘platelet alpha granule’ were the most enriched GO terms. In the MG with thymoma group versus the MG without thymoma group, ‘chemokine receptor binding’ and ‘cytokine-cytokine receptor interaction’ were the most enriched among the top 200 GO terms. A total of 13 common GO terms are shown in Figures 2E and 2F. ‘chemokine receptor binding’ and ‘cytokine-cytokine receptor interaction’ were the most enriched GO terms. The results of the GO and KEGG pathway analyses confirmed that lncRNAs play important roles in the lymphocyte immune system, such as inflammation, cell differentiation and proliferation.
Validation of disturbed lncRNA expression
To verify the disruption of lncRNA expression in MG patients, real-time PCR was performed to examine the up-or downregulation of lncRNAs in each group. As shown in Figures 3A-3D, differences in the expression of 4 lncRNAs were detected in MG patients with thymoma compared with healthy controls: lncRNA oebiotech_11933 was the most elevated (19.96-fold higher expression), followed by lncRNA A_24_P927716 (14.75-fold higher expression), whereas lncRNA A_21_P0010030 and lncRNA A_21_P0002844 exhibited 3.94- and 6.43-fold lower expression, respectively. Figures 3E-3H show the differences in the expression of 4 other lncRNAs in MG patients without thymoma compared with healthy controls. Figures 3I-3L showthe differences in the expression of 4 other lncRNAs in the MG with thymoma group compared with the MG without thymoma group. These results were consistent with the results obtained from the microarray chip analysis.
‘Cis’ analysis of the expression of lncRNAs and adjacent co-expression genes
Evidence suggests that several lncRNAs regulate their own transcription in ‘cis’, as well as that of nearby genes, by recruiting remodeling factors to local chromatin . Therefore, we identified the chromosomal co-expression genes 300 kbp upstream and downstream of the differentially expressed lncRNAs to determine potential lncRNA ‘cis’ genes. We searched the genes in the chromosome 300 kbp upstream and downstream of the differentially expressed lncRNAs and identified the common lncRNA co-expression genes within these regions of the chromosome; these genes were considered potential ‘cis’ genes of the lncRNAs. A data reduction strategy was used to overlay the genes located upstream or downstream lncRNAs and co-expression genes. The results of the ‘cis’ analyses were shown in Additional file 4. In MG patients with thymoma versus healthy controls, 157 lncRNAs had 452 ‘cis’ genes. Among these, 4 lncRNAs (oebiotech_11658, oebiotech_12721, oebiotech_21725 and oebiotech_21831) had 8 cis genes, and the lncRNA with the highest aberrant expression, lncRNA oebiotech_11933, had 3 cis genes (C1orf74, G0S2 and TRAF3IP3). The aberrantly expressed lncRNA A_24_P927716 had 2 cis genes (ACSL1 and LOC731424). In MG patients without thymoma versus healthy controls, 36 lncRNAs had 127 ‘cis’ genes upstream or downstream of their chromosomal position: lncRNA oebiotech_22652 had the highest number of ‘cis’ genes and the aberrantly expressed lncRNA oebiotech_11933 had 3 ‘cis’ genes (C1orf74, G0S2, and TRAF3IP3). The lncRNA oebiotech_03926 had 2 ‘cis’ genes (BAALC and AZIN1). In MG patients with versus without thymoma, 58 lncRNAs had 271 ‘cis’ genes upstream or downstream of their chromosomal position: lncRNA oebiotech_22652 had the highest number of ‘cis’ genes, whereas lncRNA oebiotech_12244 had 6 ‘cis’ genes(PPBP, CXCL1, CXCL5, PF4V1, PF4, and C14orf45). The cis relationships of 5 outstanding dysregulated lncRNAs (lncRNA oebiotech_11658, oebiotech_12721, oebiotech_21831, oebiotech_11933, and oebiotech_22652) are shown in Figure 4. According to the ‘cis’ transcriptional regulation mechanism, the genes that regulate transcription in ‘cis’ are affected by other genes whose chromosome positions are not far apart. Thus, the data reduction strategy was used to overlay the genes located upstream or downstream of the lncRNAs and the co-expression genes.
The ‘trans’ mechanism of the aberrant lncRNAs and the construction of the TF-lncRNA-target gene network
Currently, the known ‘trans’ regulation mechanism involves the factors mediating chromatin regulation and transcription. Using previously described methods [32,33], we calculated the lncRNA co-expression genes of chromatin regulators and transcription factors(TF) in the Encyclopedia of DNA Elements(ENCODE)  to identify common genes involved in lncRNA regulation. The “TF-lncRNA” two-element network was generated using Cytoscape software. The “TF-lncRNA” network is large and complex because numerous aberrant lncRNAs are involved. Therefore, we selected the top 100 largest relationships with the “TF-lncRNA” network to generate a core network map (Additional file 5). Figure 5A shows the “TF-lncRNA” core network map for MG patients with thymoma versus healthy controls. The transcription factor CTCF modulated the expression of 72 lncRNAs, whereas the TF TAF1 modulated the expression of 24 lncRNAs, and the TF MYC modulated the expression of 4 lncRNAs. Figure 5B shows the “TF” core network map for MG patients without thymoma versus healthy controls. The transcription factor CTCF modulated the expression of 60 lncRNAs, whereas the TF MYC modulated the expression of 9 lncRNAs. The “TF-lncRNA” core network map for MG patients without thymoma versus MG patients with thymoma indicated that TF TAF1 modulated the expression of 8 lncRNAs (Figure 5C). The transcription factor CTCF modulated the expression of 63 lncRNAs, the TF TAF1 modulated the expression of 26 lncRNAs, and the TF MYC modulated the expression of 11 lncRNAs. These three maps provide a vivid picture of the relationship between the lncRNAs and transcription factors and generate additional information for future studies.
Based on the results of the lncRNA co-expression analysis, we added the target genes into the “TF-lncRNA” network to determine the “TF-lncRNA-target genes” three-element network relationship. The networks were large and complex; therefore, we selected the top 300 relationships among lncRNAs, TFs and target genes to produce the core network map (Additional file 6). Figure 6A shows the core TF-lncRNA-target gene relationship for MG patients with thymoma versus healthy controls, containing 8 lncRNAs with disrupted expression (lncRNAs oebiotech_24272, oebiotech_23755, oebiotech_18319, oebiotech_13727, oebiotech_08281, A_21_P0008564, A_21_P0006010 and A_21_P0001906), 27 target genes and 1 core TF CTCF in this core map. Figure 6B shows the core TF-lncRNA-target gene relationship in MG patients without thymoma versus healthy controls, containing 10 lncRNAs with disrupted expression (lncRNAsoebiotech_25571, oebiotech_25019, oebiotech_23755, oebiotech_20751, oebiotech_20171, oebiotech_16583, oebiotech_07061, oebiotech_02627, oebiotech_00144 and A_21_P0007083), 41 target genes and 1 core TF CTCF in this core map. Figure 6C shows the core TF-lncRNA-target gene relationship in MG patients with thymoma versus MG patients without thymoma, containing 8 lncRNA with disrupted expression (lncRNAs oebiotech_22642, oebiotech_18319, oebiotech_18319, oebiotech_09353, oebiotech_06898, oebiotech_00715, A_21_P0010245, and A_21_P0009360), 27 target genes and1 core TF CTCF in this core map. The relationship among these three elements could be visualized using these three maps. Similar to the results shown in Figure 6A, the only core TF CTCF association regulated the expression of 8 lncRNAs and 27 target genes. The 8 lncRNAs might also regulate the expression of the 27 target genes. As observed for “CTCF-lncRNA oebiotech_24272-SOST” in this map, target genes, such as SOST, were co-expression genes for lncRNA oebiotech_24272. The transcription factor CTCF may regulate the expression of lncRNA oebiotech_24272and target genes, such as SOST. Thus, these maps provided valuable information concerning transcription factors, lncRNAs and target genes.
The thymus plays a distinct role in the pathogenesis of different MG subtypes. Marx et al.  divided MG into subgroups (early onset MG, late onset MG, and thymoma-associated MG) based on age of disease onset, gender and genetic biases, antibody specificities and associated thymus pathologies. However, in the present study, we divided MG patients into two groups: MG with or without thymoma. The first aim was to identify lncRNA subgroups as biomarkers to distinguish MG patients with or without thymoma from healthy individuals. As shown in Table 1 and Figure 2, numerous lncRNA subgroups were identified, providing the first evidence of the involvement of lncRNAs in MG. We selected 12 lncRNAs (lncRNA oebiotech_11933, lncRNA A_24_P927716, lncRNA A_21_P0010030, lncRNA A_21_P0002844, lncRNA oebiotech_11933, lncRNA oebiotech_02627, lncRNA oebiotech_22482, lncRNA A_19_P00315959, lncRNA oebiotech_13222, lncRNA oebiotech_22652, and lncRNA oebiotech_16223) for further validation among the 34 MG patients with or without thymoma. In addition, we identified 42 upregulated and 93 downregulated lncRNAs and 80 upregulated and 43 downregulated mRNAs (Figure 1D and Additional file 2) exhibiting similar expression patterns, suggesting that these lncRNAs and mRNAs are potential biomarkers for the diagnosis of MG. Thus, the identification of a subset of aberrant lncRNAs and mRNAs is the first highlight of the present study.
Precisely annotating the functions of lncRNAs remains complex. Here, we interpret the lncRNA functions based on co-expression gene GO and pathway analyses. As shown in Table 4), lncRNA oebiotech_11933is associated with the MAPK signaling pathway, chemokine signaling pathway, cytokine-cytokine receptor interaction, NOD-like receptor signaling pathway, Toll-like receptor signaling pathway, and adipocytokine-signaling pathway. Colombara et al.  reported that the constitutive activation of p38 and ERK1/2 MAPKs in the epithelial cells of a myasthenic thymus leads to IL-6 and RANTES overexpression. Cytokines and chemokines regulate immune responses in inflammatory diseases. Naïve CD4+cells differentiate into T helper type(Th1, Th2, Th9, Th17, and Th22) or regulatory T cells (Treg) through the expression of cytokines and chemokines. Because of the responses to chemokines and the interactions of these molecules with other cells, cytokines and chemokines are likely important in the pathogenesis of MG. Cufi et al  conducted RT-PCR analyses demonstrating large increases in the expression of interferon (IFN)-I subtypes (IFN-α2, IFN-α8, IFN-ϖ and IFN-β) in thymoma-associated MG, but not in MG patients without thymoma in control thymuses, and the abnormal regulation of dsRNA-sensing molecules, with the increased expression of Toll-like receptor 3 (TLR3) and decreased expression of protein kinase R (PKR) and dsRNA helicases (RIG-I and MDA5) in MG patients with thymoma. These authors suggested that INF-I overexpression and the activation of innate immunity pathways in thymoma-associated MG may develop after a pathogen infection. Cufi et al [27,29] also observed the contribution of IFN-I and IFN-III subtypes in thymic changes associated with MG. The IFN-I and IFN-III subtypes, particularly IFN-β, specifically induce α-AChR expression in thymic epithelial cells (TECs), and IFN-β increases TEC death and the uptake of TEC proteins in dendritic cells. In parallel, these results showed that IFN-β increased the expression of the chemokines CXCL13 and CCL21 in TECs and lymphatic endothelial cells, and similar results were obtained in C57BL/6 mice. All of these studies were performed using thymus tissues. In the present study, our aim was to identify potential biomarkers for the diagnosis of MG; therefore, we performed our studies on PBMCs. We identified several cytokine and chemokine GO terms and pathways that are consistent with the aforementioned studies. It is reported that the cytokines/chemokines in PBMCs may originate from the thymus [7,10]. Uzawa et al.  measured the serum levels of 27 cytokines/chemokines in 47 anti-acetylcholine receptor antibody-positive patients with MG and 20 healthy controls to investigate the contribution of cytokines/chemokines in MG pathogenesis. These authors identified eight cytokines that were significantly changed among MG patients with thymoma (i.e., IL-4, IL-8, IL-15, eotaxin, macrophage inflammatory protein-1α, macrophage inflammatory protein-1β, VEGF and IL-1b). Cordiglieri et al.  reported that TLR4 stimulation in MG TECs increased CCL17 and CCL22 expression and induced the production of Th17-related cytokines. TLR4 signaling in the MG thymic milieu might affect cell-to-cell interactions, favoring autoreactive T-cell activation. These results and the results obtained in the present study suggest that lncRNA oebiotech_11933 could play an important role in the pathogenesis of MG.
The most enriched GO terms in the predicted target genes of the lncRNAs were ‘cell response to interferon-γ’, ‘platelet degranulation’, ‘chemokine receptor binding’ and ‘cytokine-cytokine interaction’ in the pathogenesis between MG patients with thymoma and MG patients without thymoma (Figure 3). Many studies have shown that the number of immune cells secreting IL-2, IL-4, IFN-γ, and IL-10 are significantly higher in patients with MG than in healthy individuals [38-40]. The results of the present study demonstrated that the most important GO term was ‘cell response to interferon-γ’, suggesting that aberrant lncRNAs may play an important role in regulating IFN-γ expression. Interferon γ is a soluble cytokine and a member of the type II interferon family with antiviral, immunoregulatory and anti-tumor properties as a potent activator of macrophages. Tuzun et al.  reported that MG patients with low plasma IL-6 and IFN-γ levels attained better clinical improvement following etanercept treatment. An increasing number of studies have shown that IFN-γ regulates T and B cell differentiation [42-44]. The GO and KEGG pathway analysis of lncRNA co-expression genes to determine the functions of lncRNAs is the second highlight in the present study.
The molecular regulation through lncRNAs remains unknown because the functions of lncRNAs vary . Indeed, lncRNAs have been identified using a variety of methods, and the number of specific lncRNAs shown to influence genomic functions, including roles in imprinting , enhancer functions [46,47], X chromosome inactivation , chromatin structure  and genomic rearrangements during the generation of antibody diversity , is increasing. Here, we utilized the ‘cis’ versus ‘trans’ regulatory mechanisms to obtain additional information on dysregulated lncRNAs. In the present study, we explored ‘cis’ regulatory relationships (Additional file 4). The outstanding lncRNAs are shown in Figure 4. The increased expression of cytokines and chemokines, such as IL-8, CXCL1, CXCL3, CXCL3, and CXCL5, was observed among the oebiotech_11658, oebiotech_12721, oebiotech_21831, oebiotech_11933, and oebiotech_22652 ‘cis’ genes, suggesting that these lncRNAs may regulate the expression of these cytokines and chemokines. The results of this ‘cis’ analysis provide additional information concerning the regulation mediated through lncRNAs and the biological processes in the pathogenesis of MG. We constructed the “TF-lncRNA” and “TF-lncRNA-target gene” network based on the results of the ‘trans’ analysis (Additional file 6). The core TF-lncRNA-target gene network (Figures 5 and 6) showed that TFs, including CTCF, TAF1and MYC, regulate lncRNA expression in MG. The CTCF gene is a member of the BORIS + CTCF gene family, encoding a transcriptional regulator protein with 11 highly conserved zinc finger (ZF) domains. This nuclear protein uses different combinations of the ZF domains to bind different DNA target sequences and proteins. Depending on the context of the site, this protein binds to the histone acetyltransferase (HAT)-containing complex and functions as a transcriptional activator or binds to the histone deacetylase (HDAC)-containing complex and functions as a transcriptional repressor. When the protein is bound to a transcriptional insulator element, the communication between enhancers and upstream promoters is blocked, thereby regulating imprinted expression [51,52]. Thus, the ‘trans’ analysis provides another method to interpret lncRNA function and the biological processes in the pathogenesis of MG. The results of the ‘cis’ and ‘trans’ analyses provide information concerning the modular regulation of lncRNAs, representing the third highlight of the present study.
In the present study, we identified a subset of aberrant lncRNAs to distinguish MG patients with or without thymoma compared with healthy individuals. The function and biological processes of lncRNAs in the pathogenesis of MG were determined according to co-expression gene GO and pathway annotations. The results of the ‘cis’ and ‘trans’ analyses provide information for future studies of aberrantly expressed lncRNAs in MG. These results provide support for future investigations of the pathogenesis of MG.
Patients and sample collection
MG Patient Demographics
A total of 34 MG patients examined at the Neurology Department of Xiangya Hospital and 13 healthy donors were recruited from May 2010 to March 2012. MG was diagnosed based on a combination of fluctuating muscle weakness with a positive neostigmine test or abnormal single-fiber EMG test. The MG patients were divided into two groups: patients with or without thymoma. The Myasthenia Gravis Foundation of America (MGFA) clinical classification was also used to identify MG subgroups. All MG patients did not receive immunomodulatory or immunosuppressive treatment and were not treated with thymectomy. Information pertaining to sex, age at onset, disease duration, MGFA clinical classification upon first visit to the hospital, additional autoimmune diseases, including thyroid disorders, systemic lupus erythematosus (SLE), rheumatoid arthritis, and Sjogren syndrome, and previous treatments, including acetylcholinesterase inhibitors, immunosuppressive drugs, prednisone, plasma exchange, IVIg and thymectomy, is summarized in Additional file 7. The samples MG1-8 and Normal1-4 were selected for lncRNA chip analysis, and all of the samples in Additional file 7 were used for the validation of lncRNA expression. All patients and healthy volunteers provided informed consent for sample collection and have signed informed consent forms. Collections and use of tissue samples were approved by the Ethical Review Committee of Xiangya Hospital. The PBMCs were isolated from approximately 10 ml of EDTA whole blood using a Ficoll-Paque™ PLUS density gradient, and then, 5x107 cells were resuspended in Trizol® (Invitrogen) and stored at -80°C.This study was approved by the ethical review committees of Xiangya Hospital.
RNA extract and the lncRNA chip array
Total RNA was extracted from PBMCs using Trizol® reagent (Invitrogen). Approximately 200 ng of total RNA from each sample was used for the lncRNA microarray analysis. lncRNA expression was analyzed from May 2012 to August 2013 using OE_Biotech Human lncRNA chip software V2.0(4*180K),containing 46,506 lncRNAs and 30,656 mRNAs collected from eight authoritative databases, including Agilent_ncRNA, lncRNAdb, Gencode V13, H-invDB, NONCODE v3.0, RefSeq, UCR and UCSC_lncRNAs Transcripts. The lncRNA chip experiments were conducted at the OEbiotech Corporation in Shanghai, P.R. China.
Difference lncRNAs and mRNA screen
Raw data from each array were first normalized using GeneSpring software (version 12.5) and subsequently analyzed using an unpaired t-test, with a P-value cut-off of 0.05 and a fold-change cut-off of 2.0.
Difference integration analysis (Venn analysis)
MG patients with or without thymoma were compared with each other and healthy controls. The common characteristic elements between the three groups were determined using Venn analysis.
Difference lncRNA and mRNA clustering analysis
Different lncRNAs and mRNAs were analyzed using Cluster 3.0 software, and the data were used to examine a series of parameters, such as log transform data, normalized genes and arrays, and hierarchical parameters of genes and arrays. The results were further analyzed using Tree View software. Green-yellow indicates low expression, and red indicates high expression.
lncRNA co-expression analysis and gene function annotation
The expression of different lncRNAs and mRNAs was analyzed using Pearson’s correlation coefficient. The absolute value of 0.8 was considered relevant, a value less than 0.8 represented a negative correlation, and a value greater than 0.8 represented a positive correlation. A P-value of less than 0.05 was considered statistically significant. The expression of genes encoding each differentially expressed lncRNA, the ontology classification of the co-expression genes based on gene annotation and summary information are available through DAVID (Database for Annotation, Visualization and Integrated Discovery). The predicted target genes were assigned to functional groups based on molecular function, biological processes and specific pathways. The lncRNA gene function was predicted based on the co-expression gene GO functional annotation, selecting (according to the P-value) the top 100 GO and 200 function reliability prediction terms. Statistical function annotation generated additional GO terms, and the most enriched terms might reflect potential lncRNA functions.
‘Cis’ analysis of lncRNAs and co-expression genes
The gene location for different lncRNAs on the chromosome was determined. Subsequently, the common lncRNA co-expression genes were intersected to identify the genes 300 kbp upstream or downstream of the lncRNAs as potential ‘cis’ genes. The schematic shows the chromosome location of the lncRNAs and cis genes.
lncRNA sequences were mapped to the genome in the Sanger database. Jemboss software was used to examine the alignment of lncRNA and putative transcription factor binding sequences. The genome browser database was used to build the network describing the relationships between transcription factors and lncRNAs. An adjacency matrix was implemented in Java according to the binding of lncRNA and transcription factors. The core transcription factor is the most important center in the network, with the highest degree of expression [53,54]. Pearson’s correlation analysis  was used to measure the regulatory ability of transcription factors by calculating the correlation between transcription factors and lncRNAs. The TF and lncRNA relationship was generated using Cytoscape software. The red circle indicates the enriched lncRNAs or TFs, and the larger size indicates increasing enrichment.
The TF-lncRNA-gene network was constructed based on the interactions of lncRNAs and target co-expression genes as previously described . The lncRNA co-expression genes and the compounds for chromatin regulation and transcription complexes were intersection elements. Intersection elements might participate in lncRNA-mediated gene regulation. The three groups were generated using Cytoscape software based on the “TF-lncRNA” two-element network. Because many lncRNAs are involved, the “TF-lncRNA” network is difficult and complex. Therefore, we selected the top100 closest associations with the “TF-lncRNA” network to produce the core network map.
Quantitative reverse transcription–polymerase chain reaction analysis
Total RNA was extracted using Trizol® reagent (Invitrogen). The primers for RT-PCR were designed based on the lncRNA sequences from the UCSC. The primers were synthesized and purified at Invitrogen (Shanghai, China). The RT reactions were performed using a cDNA synthesis kit (Bio-Rad, Hercules, CA). Real-time PCR was performed using the ABI StepOne Plus™ Multicolor. The qPCR cycle was 98°C for 2 min, followed by 40 cycles of 95°C for 15 sec and 60°C for 30 sec. A final melt-curve analysis (60–95°C) was included. The standard curve was produced with slopes at approximately -3.32 (~100% efficiency). The lncRNA PCR results were quantified using the 2ΔΔct method against GAPDH for normalization. The data represent the means of three experiments.
Real-time PCR primers:
Long non-coding RNAs
Kyoto encyclopedia of genes and genomes
Lipoprotein-related protein 4
Systemic lupus erythematosus
Open reading frames
Toll-like receptor 3
Protein kinase R
Thymic epithelial cells
Myasthenia Gravis Foundation of America
Encyclopedia of DNA Elements
Peripheral blood mononuclear cell
Berrih-Aknin S, Frenkian-Cuvelier M, Eymard B. Diagnostic and clinical classification of autoimmune myasthenia gravis. J Autoimmun. 2014;48–49:143–8.
Le Panse R, Berrih-Aknin S. Autoimmune myasthenia gravis: autoantibody mechanisms and new developments on immune regulation. Curr Opin Neurol. 2013;26(5):569–76.
Berrih-Aknin S, Le Panse R. Myasthenia gravis: a comprehensive review of immune dysregulation and etiological mechanisms. J Autoimmun. 2014;52:90–100.
Thiruppathi M, Rowin J, Li Jiang Q, Sheng JR, Prabhakar BS, Meriggioli MN. Functional defect in regulatory T cells in myasthenia gravis. Ann N Y Acad Sci. 2012;1274:68–76.
Matsui N, Nakane S, Saito F, Ohigashi I, Nakagawa Y, Kurobe H, et al. Undiminished regulatory T cells in the thymus of patients with myasthenia gravis. Neurology. 2010;74(10):816–20.
Fujii Y. Thymus, thymoma and myasthenia gravis. Surg Today. 2013;43(5):461–6.
Bouchikh M, El Malki HO, Ouchen F, Achir A, Benosman A. [Thymoma-associated myasthenia gravis: Clinical features and surgical results]. Rev Neurol. 2013;169(11):879–83.
Li J. Myasthenia gravis and recurrent thymoma. Ann Thorac Surg. 2011;92(1):406. author reply 406-407.
Cizeron-Clairac G, Le Panse R, Frenkian-Cuvelier M, Meraouna A, Truffault F, Bismuth J, et al. Thymus and Myasthenia Gravis: what can we learn from DNA microarrays? J Neuroimmunol. 2008;201-202:57–63.
Marx A, Pfister F, Schalke B, Saruhan-Direskeneli G, Melms A, Strobel P. The different roles of the thymus in the pathogenesis of the various myasthenia gravis subtypes. Autoimmun Rev. 2013;12(9):875–84.
Vachlas K, Zisis C, Rontogianni D, Tavernarakis A, Psevdi A, Bellenis I. Thymoma and myasthenia gravis: clinical aspects and prognosis. Asian Cardiovasc Thorac Ann. 2012;20(1):48–52.
Song X, Cao G, Jing L, Lin S, Wang X, Zhang J, et al. Analysing the relationship between lncRNA and protein-coding gene and the role of lncRNA as ceRNA in pulmonary fibrosis. J Cell Mol Med. 2014;18(6):991–1003.
Necsulea A, Soumillon M, Warnefors M, Liechti A, Daish T, Zeller U, et al. The evolution of lncRNA repertoires and expression patterns in tetrapods. Nature. 2014;505(7485):635–40.
Dong R, Jia D, Xue P, Cui X, Li K, Zheng S, et al. Genome-wide analysis of long noncoding RNA (lncRNA) expression in hepatoblastoma tissues. PLoS One. 2014;9(1):e85599.
Wan G, Hu X, Liu Y, Han C, Sood AK, Calin GA, et al. A novel non-coding RNA lncRNA-JADE connects DNA damage signalling to histone H4 acetylation. EMBO J. 2013;32(21):2833–47.
Froberg JE, Yang L, Lee JT. Guided by RNAs: X-inactivation as a model for lncRNA function. J Mol Biol. 2013;425(19):3698–706.
Engreitz JM, Pandya-Jones A, McDonel P, Shishkin A, Sirokman K, Surka C, et al. The Xist lncRNA exploits three-dimensional genome architecture to spread across the X chromosome. Science. 2013;341(6147):1237973.
Maury P, Moreau A, Hidden-Lucet F, Leenhardt A, Fressart V, Berthet M, et al. Novel SCN5A mutations in two families with "Brugada-like" ST elevation in the inferior leads and conduction disturbances. J Interv Card Electrophysiol. 2013;37(2):131–40.
Li J, Chen Z, Tian L, Zhou C, He MY, Gao Y, et al. LncRNA profile study reveals a three-lncRNA signature associated with the survival of patients with oesophageal squamous cell carcinoma. Gut. 2014;63(11):1700–10.
Tang H, Wu Z, Zhang J, Su B. Salivary lncRNA as a potential marker for oral squamous cell carcinoma diagnosis. Mol Med Rep. 2013;7(3):761–6.
Liu Q, Huang J, Zhou N, Zhang Z, Zhang A, Lu Z, et al. LncRNA loc285194 is a p53-regulated tumor suppressor. Nucleic Acids Res. 2013;41(9):4976–87.
Wapinski O, Chang HY. Long noncoding RNAs and human disease. Trends Cell Biol. 2011;21(6):354–61.
Sauvageau M, Goff LA, Lodato S, Bonev B, Groff AF, Gerhardinger C, et al. Multiple knockout mouse models reveal lincRNAs are required for life and brain development. eLife. 2013;2:e01749.
Li Z, Chao TC, Chang KY, Lin N, Patil VS, Shimizu C, et al. The long noncoding RNA THRIL regulates TNFalpha expression through its interaction with hnRNPL. Proc Natl Acad Sci U S A. 2014;111(3):1002–7.
Carpenter S, Aiello D, Atianand MK, Ricci EP, Gandhi P, Hall LL, et al. A long noncoding RNA mediates both activation and repression of immune response genes. Science. 2013;341(6147):789–92.
Hu G, Tang Q, Sharma S, Yu F, Escobar TM, Muljo SA, et al. Expression and regulation of intergenic long noncoding RNAs during T cell development and differentiation. Nat Immunol. 2013;14(11):1190–8.
Cufi P, Dragin N, Ruhlmann N, Weiss JM, Fadel E, Serraf A, et al. Central role of interferon-beta in thymic events leading to myasthenia gravis. J Autoimmun. 2014;52:44–52.
Cavalcante P, Serafini B, Rosicarelli B, Maggi L, Barberis M, Antozzi C, et al. Epstein-Barr virus persistence and reactivation in myasthenia gravis thymus. Ann Neurol. 2010;67(6):726–38.
Cufi P, Dragin N, Weiss JM, Martinez-Martinez P, De Baets MH, Roussin R, et al. Implication of double-stranded RNA signaling in the etiology of autoimmune myasthenia gravis. Ann Neurol. 2013;73(2):281–93.
Cufi P, Soussan P, Truffault F, Fetouchi R, Robinet M, Fadel E, et al. Thymoma-associated myasthenia gravis: On the search for a pathogen signature. J Autoimmun. 2014;52:29–35.
Guenzl PM, Barlow DP. Macro lncRNAs: a new layer of cis-regulatory information in the mammalian genome. RNA Biol. 2012;9(6):731–41.
Guttman M, Rinn JL. Modular regulatory principles of large non-coding RNAs. Nature. 2012;482(7385):339–46.
Khalil AM, Guttman M, Huarte M, Garber M, Raj A, Rivea Morales D, et al. Many human large intergenic noncoding RNAs associate with chromatin-modifying complexes and affect gene expression. Proc Natl Acad Sci U S A. 2009;106(28):11667–72.
Gerstein MB, Kundaje A, Hariharan M, Landt SG, Yan KK, Cheng C, et al. Architecture of the human regulatory network derived from ENCODE data. Nature. 2012;489(7414):91–100.
Colombara M, Antonini V, Riviera AP, Mainiero F, Strippoli R, Merola M, et al. Constitutive activation of p38 and ERK1/2 MAPKs in epithelial cells of myasthenic thymus leads to IL-6 and RANTES overexpression: effects on survival and migration of peripheral T and B cells. J Immunol. 2005;175(10):7021–8.
Uzawa A, Kawaguchi N, Himuro K, Kanai T, Kuwabara S. Serum cytokine and chemokine profiles in patients with myasthenia gravis. Clin Exp Immunol. 2014;176(2):232–7.
Cordiglieri C, Marolda R, Franzi S, Cappelletti C, Giardina C, Motta T, et al. Innate immunity in myasthenia gravis thymus: pathogenic effects of Toll-like receptor 4 signaling on autoimmunity. J Autoimmun. 2014;52:74–89.
Zhang B, Shen C, Bealmear B, Ragheb S, Xiong WC, Lewis RA, et al. Autoantibodies to agrin in myasthenia gravis patients. PLoS One. 2014;9(3):e91816.
Zheng K, Xu G, Lu X, Zhang J, Zhang P. Expression and polymorphisms of T cell immunoglobulin domain and mucin domain protein-1 in thymoma with or without myasthenia gravis. Oncol Lett. 2014;8(1):317–22.
Levin MH, Gopal PP, Galetta SL. Systemic inflammatory pseudotumor and myasthenia gravis. JAMA Ophthalmol. 2014;132(3):359–61.
Tuzun E, Meriggioli MN, Rowin J, Yang H, Christadoss P. Myasthenia gravis patients with low plasma IL-6 and IFN-gamma benefit from etanercept treatment. J Autoimmun. 2005;24(3):261–8.
Zhang L, Liu J, Wang H, Zhao C, Lu J, Xue J, et al. Double filtration plasmapheresis benefits myasthenia gravis patients through an immunomodulatory action. J Clin Neurosci. 2014;21(9):1570–4.
Yilmaz V, Tutuncu Y, Baris Hasbal N, Parman Y, Serdaroglu P, Deymeer F, et al. Polymorphisms of interferon-gamma, interleukin-10, and interleukin-12 genes in myasthenia gravis. Hum Immunol. 2007;68(6):544–9.
Yi JS, Guidon A, Sparks S, Osborne R, Juel VC, Massey JM, et al. Characterization of CD4 and CD8 T cell responses in MuSK myasthenia gravis. J Autoimmun. 2014;52:130–8.
Nagano T, Fraser P. Emerging similarities in epigenetic gene silencing by long noncoding RNAs. Mamm Genome. 2009;20(9–10):557–62.
Kim A, Zhao H, Ifrim I, Dean A. Beta-globin intergenic transcription and histone acetylation dependent on an enhancer. Mol Cell Biol. 2007;27(8):2980–6.
Ho Y, Elefant F, Liebhaber SA, Cooke NE. Locus control region transcription plays an active role in long-range gene activation. Mol Cell. 2006;23(3):365–75.
Lee JT. Lessons from X-chromosome inactivation: long ncRNA as guides and tethers to the epigenome. Genes Dev. 2009;23(16):1831–42.
Rinn JL, Kertesz M, Wang JK, Squazzo SL, Xu X, Brugmann SA, et al. Functional demarcation of active and silent chromatin domains in human HOX loci by noncoding RNAs. Cell. 2007;129(7):1311–23.
Krangel MS. T cell development: better living through chromatin. Nat Immunol. 2007;8(7):687–94.
Ong CT, Corces VG. CTCF: an architectural protein bridging genome topology and function. Nat Rev Genet. 2014;15(4):234–46.
Zuin J, Dixon JR, van der Reijden MI, Ye Z, Kolovos P, Brouwer RW, et al. Cohesin and CTCF differentially affect chromatin architecture and gene expression in human cells. Proc Natl Acad Sci U S A. 2014;111(3):996–1001.
Prieto C, Risueno A, Fontanillo C, De las Rivas J. Human gene coexpression landscape: confident network derived from tissue transcriptomic profiles. PLoS One. 2008;3(12):e3911.
Vermeirssen V, Barrasa MI, Hidalgo CA, Babon JA, Sequerra R, Doucette-Stamm L, et al. Transcription factor modularity in a gene-centered C. elegans core neuronal protein-DNA interaction network. Genome Res. 2007;17(7):1061–71.
Joung JG, Hwang KB, Nam JW, Kim SJ, Zhang BT. Discovery of microRNA-mRNA modules via population-based probabilistic learning. Bioinformatics. 2007;23(9):1141–7.
This work was supported through funding from the National Natural Science Foundation, China (grant numbers: 81301024), the Hunan Natural Science Foundation (grant numbers: 14JJ3033 and 20155JJ3136) and the China Postdoctoral Science Foundation (grant number:2014M550427).
The authors declare that they have no competing interests.
ZHL and YL conducted the majority of the experiments. XFL and MCL collected the biopsies and data. ZHL, BX, LQX and HY contributed to the study design and analyzed and interpreted the data. ZHL, YBL, BX and HY drafted the manuscript. ZHL and HY provided financial support. LQX and YL collected and assembled the data. All of the authors have approved the publication of the manuscript.
The difference lncRNAs and mRNAs.
The same tendency expression lncRNAs and mRNAs.
The dyexpression lncRNAs annotation.
The dyexpression lncRNAs cis analysis.
The lncRNAs and their core TF relationship.
The lncRNAs-TF-target gen coreNET relationship.
The samples information.
About this article
Cite this article
Luo, Z., Li, Y., Liu, X. et al. Systems biology of myasthenia gravis, integration of aberrant lncRNA and mRNA expression changes. BMC Med Genomics 8, 13 (2015). https://doi.org/10.1186/s12920-015-0087-z
- Long non-coding RNAs (lncRNAs)
- Myasthenia gravis(MG)
- Transcription factors