Front. Cell Dev. Biol.Frontiers in Cell and Developmental BiologyFront. Cell Dev. Biol.2296-634XFrontiers Media S.A.75656010.3389/fcell.2021.756560Cell and Developmental BiologyOriginal ResearchIdentification of lncRNA and mRNA Expression Profile in Relapsed Graves’ DiseaseYao et al.lncRNAs and GDYaoQiuming1SongZhenyu2WangBin1JiaXi1SongRonghua1*ZhangJinan1*
Department of Endocrinology, Shanghai University of Medicine and Health Sciences Affiliated Zhoupu Hospital, Shanghai, China
Ovarian Cancer Program, Department of Gynaecologic Oncology, Zhongshan Hospital, Fudan University, Shanghai, China
Edited by:Marina I. Arleevskaya, Kazan State Medical Academy, Russia
Reviewed by:Lin Liao, Shandong University, China
Jing Li, China Medical University, China
Shuang-Xia Zhao, Shanghai Ninth People’s Hospital, China
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Background: Graves’ disease (GD) is a common autoimmune disease, and its pathogenesis is unclear. Studies have found that the occurrence of GD is related to the immune disorder caused by the interaction of genetic susceptibility and environmental factors. The CD4+ T cell subset is closely related to the immune disorder of GD. LncRNAs are RNA molecules with a length of more than 200 nt and are involved in a variety of autoimmune diseases. However, the roles of lncRNAs in recurrent GD are still elusive. The purpose of this study is to identify lncRNA and mRNA expression profile in relapsed Graves’ disease.
Method: CD4+ T cells from 12 recurrent GD and 8 healthy controls were collected for high-throughput sequencing. The gene-weighted co-expression network analysis (WGCNA) was used to construct the co-expression module relevant to recurrent GD, and the key genes in the module were verified by RT-PCR.
Results: There are 602 upregulated lncRNAs and 734 downregulated lncRNAs in CD4+ T cells in recurrent GD patients compared with the healthy controls. The module most relevant to GD recurrence was constructed using WGCNA, and the key genes in the module were verified by RT-PCR. We found that the expression of RPL8, OAS2, NFAT5, DROSHA, NONHSAT093153.2, NONHSAT118924.2, and NONHSAT209004.1 was significantly decreased in GD group (p < 0.001, p < 0.001, p < 0.01, p < 0.05, p < 0.001, p < 0.05, and p < 0.01, respectively).
Conclusion: LncRNAs are closely related to the recurrence of GD. For the first time, we constructed the expression profile of lncRNAs and mRNAs in CD4+ T cells in recurrent GD patients.
lncRNAsrelapsed GDWGCNANONHSAT093153.2NONHSAT209004.1No. 81873636Yalong River Joint Fund10.13039/501100019490Introduction
Graves’ disease, also known as toxic diffuse goiter, is characterized by the production of antibodies against thyroid stimulating hormone receptors (TRAb), leading to the hypertrophy and hyperfunction of the thyroid follicular cells (Morshed et al., 2012). GD is the most common cause of hyperthyroidism, and its incidence is about 0.5% (Brent, 2008). Hyperthyroidism caused by GD also increases the risk of atrial fibrillation, congestive heart failure, and miscarriage in pregnant women (Iddah and Macharia, 2013). However, the pathogenesis of GD is unclear, and its therapeutic effect is not satisfactory. At present, there are three main treatment methods for GD, including anti-thyroid drugs, radioactive iodine (RAI), and surgical resection. Each of the three methods has its own advantages and disadvantages. Compared with RAI and surgery, the disadvantage of anti-thyroid drugs mainly includes lower remission rate of hyperthyroidism and higher recurrence rate of patients with high TRAb titer. In addition, anti-thyroid drugs can also cause side effects such as skin rash, joint pain, agranulocytosis, and liver toxicity (Kotwal et al., 2018). Although radioactive iodine therapy has a higher cure rate, it has more chances of causing permanent hypothyroidism. Moreover, RAI is contraindicated in pregnant and lactating women and patients with active thyroid eye disease (Kotwal et al., 2018). Although surgical removal of the thyroid can quickly improve the symptoms of hyperthyroidism, patients need to take thyroid hormones throughout their lives, and the operation itself can cause a variety of complications such as hypoparathyroidism, recurrent laryngeal nerve injury, and neck hematoma. Therefore, it is urgent to further study the pathogenesis of GD and, on this basis, develop etiological treatment methods to reduce the serious harm of GD to public health.
Long non-coding RNAs (lncRNAs) are new members with more than 200 nucleotides in length of the non-coding RNAs (Caley et al., 2010). Although lncRNAs do not encode any protein products, they can regulate the gene expression at the transcriptional, post-transcriptional, and epigenetic level (Caley et al., 2010; Guttman et al., 2011). LncRNAs are also involved in functionally distinct biological and physiological processes such as chromatin remodeling, RNA junction, and protein transport (Mercer et al., 2009). Several studies have shown that lncRNAs are associated with autoimmune diseases, such as Crohn’s disease (Qiao et al., 2013), systemic lupus erythematosus (SLE) (Zhang Y. et al., 2021), and rheumatoid arthritis (Zhang J. et al., 2021), but the association of lncRNAs with relapsed GD remains unclear. Here, this study explored the potential roles of lncRNAs in relapsed GD.
Materials and MethodsSubjects
Forty-six relapsed GD patients and 33 age- and sex-matched normal healthy controls (NC) were enrolled from Zhoupu Hospital. Among them, 12 GD patients and 8 healthy controls were collected for lncRNA and mRNA sequencing, and the rest of subjects were recruited for the subsequent validation. Relapsed GD was diagnosed based on recurrence of clinical symptoms, elevated free triiodothyronine(FT3) or free thyroxine(FT4), suppressed thyroid-stimulating hormone(TSH), and positive anti-thyrotropin receptor antibody (TRAb) after a 12- to 18-month course of ATD treatment. We also detected the antibody against thyroglobulin (TGAb) or thyroid peroxidase (TPOAb) level of relapsed GD patients. Individuals without any acute or chronic autoimmune or allergic or infectious diseases or any acute or chronic visceral diseases were recruited as healthy controls or normal controls(NC). The study was approved by the Ethics Committee of Zhoupu Hospital. All subjects signed an informed consent form.
CD4<sup>+</sup> T Cell Isolation
Firstly, peripheral blood mononuclear cells (PBMCs) of all subjects were isolated from freshly collected venous blood by lymphocyte separation medium (Sigma-Aldrich) according to the manufacturer’s instruction. Then, the human CD4 Micro Beads (Miltenyi Biotec, Germany) was used to purificate CD4+ T cells from fresh PBMCs. The CD4+ T cells with a purity greater than 95% were used for further research. We calculated the purity of CD4+ T cells by flow cytometer (BD Biosciences, USA).
Differentially Expressed Gene Screening
CD4+ T cells were isolated from PBMCs of 12 relapsed GD patients and 8 healthy controls. Then, we added 1 ml TRIzol to the CD4+ T cells. Total RNA was extracted from the CD4+ T cells using the TRIzol reagent (Takara) according to the manufacturer’s protocol. The samples were preserved at −80°C for further analysis. These samples were then subjected to mRNA and lncRNA-seq on the Illumina HiSeq platform following the standard procedures. The raw data were cleaned to obtain the reads with high quality. The clean reads with high quality were then aligned to the reference genome. Subsequently, the differentially expressed lncRNAs and mRNAs between GD patients and NC were screened in the expressing data using the “edeg” R package. The significantly changed genes with p < 0.05 and log2 fold change (FC)≥1 were considered as differentially expressed genes. Gene ontology (GO) enrichment and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis was conducted by R package.
Construction of Co-Expression Module
We constructed the co-expression module of the differentially expressed lncRNAs and mRNAs (p < 0.05 and average expression level >1) between relapsed GD group and healthy controls by the R package “WGCNA”. We chose 5 as the soft-thresholding power, and 30 was chosen as the minimum number of modules.
Hub Gene Identification
In the module-trait correlation analysis, the genes with gene significance greater than 0.4 and module group members (MM) greater than 0.9 were considered as hub genes, which are significantly associated with clinical phenotypes.
Statistical Analysis
Software R3.5.3. was used to perform WGCNA analysis. The comparison between relapsed GD group and healthy controls was analyzed using non-parametric test. A p value less than 0.05 was considered statistically significant.
ResultsThe Expression Profile of lncRNAs and mRNAs in Relapsed GD CD4<sup>+</sup> T Cells
To explore the crucial role of lncRNAs and mRNAs associated with the recurrence and development of GD, we performed RNA-Seq to detect the expression profile of lncRNAs and mRNAs in GD and NC group. Totally, we found that 1336 lncRNAs and 266 mRNAs were significantly differentially expressed between relapsed GD patients and healthy controls. Of the identified lncRNAs, 602 lncRNAs were significantly upregulated and 734 lncRNAs were significantly downregulated in the CD4+ T cells of GD patients.
Of those detected mRNAs, 128 mRNAs were upregulated and 138 mRNAs were downregulated in CD4+ T cells of GD patients. Hierarchical cluster analyses displayed lncRNA and mRNA expression profile in two groups (Figures 1A,C). Volcano plot analyses were also performed to visualize the differentially expressed lncRNAs and mRNAs (Figures 1B,D).
Expression profile of lncRNAs and mRNAs in CD4+ T cells of relapsed GD patients. (A) Hierarchical clustering of differentially expressed lncRNAs between GD group (n = 12) and normal controls (NC) (n = 8). (B) Volcano plots of lncRNA expression levels between two groups. (C) Hierarchical clustering of differentially expressed mRNAs between GD group (n = 12) and normal controls (NC) (n = 8). (D) Volcano plots of mRNA expression levels between two groups. Each column represents a sample, and each row indicates one gene. Red indicates those genes with relatively high expression level, and blue indicates those genes with relatively low expression level.
GO Analysis and Pathway Analysis
We conducted Gene ontology (GO) and KEGG pathway enrichment analyses to further explore the function of those differentially expressed genes. The GO analysis found that differentially expressed genes identified were mainly enriched in hemoglobin complex, oxygen transporter activity, excitatory postsynaptic potential, and positive regulation of mitotic nuclear division (Figure 2A). KEGG pathway analysis revealed that those genes were mainly enriched in glycine, serine, and threonine metabolism, complement and coagulation cascades, and cell adhesion molecules (CAMs) (Figure 2B).
GO and KEGG pathway analysis in relapsed GD CD4+ T cells. (A) GO analysis of differentially expressed genes between GD group and contols. According to biological process (circle), cellular component (triangle), and molecular function (square). (B) KEGG pathway analysis for differentially expressed mRNAs.
WGCNA Analysis
As shown in Figure 3, we constructed a total of 13 co-expression modules by WGCNA analysis (Figure 3A). Moreover, these constructed modules were independent of each other (Figure 3B). Figure 3C shows an eigengene dendrogram and adjacency heatmap.
WGCNA revealed gene co-expression modules in the CD4+ T cells of relapsed GD patients. (A) Clustering dendrograms of mRNAs and lncRNAs. Each color represents a co-expression module. (B) Network heatmap plot in the co-expression modules. (C) Eigengene dendrogram and eigengene adjacency heatmap.
As displayed in Figure 4, module-trait correlations showed that five modules were related to GD, including red, salmon, brown, yellow, and tan module. Interestingly, all these modules are also related to TSH and TRAb. While brown, yellow, red, salmon, and purple modules were associated with FT3, three modules including brown, yellow, and salmon were related to FT4. The brown, yellow, and red modules were associated with TPOAb. Figure 5 shows module significance values of co-expression modules associated with each phenotype, including GD (Figure 5A), FT3 (Figure 5B), FT4 (Figure 5C), TSH (Figure 5D), TPOAb (Figure 5E) and TRAb (Figure 5F). Figure 6 shows the scatterplots of gene significance for GD (Figure 6A) and TRAb (Figure 6B) vs. MM in different modules.
Heatmap of the correlation between each coexpression module and relapsed GD and different phenotypes. The phenotypes mainly include sex, FT3, FT4, TSH, TGAb, TPOAb, and TRAb. (The correlation coefficient and corresponding p value were shown in each cell.)
Module significance values of co-expression modules related to different phenotypes. (A) GD; (B) FT3; (C) FT4; (D), TSH; (E) TPOAb; (F) TRAb. Each color indicated one coexpression module.
The scatterplots of different modules highly related to GD and TRAb. (A) GD; (B) TRAb. The correlation coefficient and p value of module membership vs. gene significance were shown.
Functional Annotation of Key Co-Expression Modules
Figure 7 shows the results of GO analysis about genes in different module. Genes in brown module were mainly enriched in negative regulation of gene expression and epigenetic and chromatin silencing (Figure 7A); genes in yellow module were enriched in viral transcription, viral gene expression, and nuclear-transcribed mRNA catabolic process (Figure 7B); genes in red module were mainly enriched in inositol phosphate-mediated signaling, histone H3-K4 methylation, calcineurin-NFAT signaling cascade, and calcineurin-mediated signaling (Figure 7C); the salmon module was mainly enriched in type I interferon signaling pathway, response to virus, etc (Figure 7D).
GO analysis of genes in coexpression modules associated with relapsed GD. (A) Brown module; (B) Yellow module; (C) Red module; (D) Salmon module.
KEGG analysis (Figure 8) found that the genes in brown module were involved in alcoholism, systemic lupus erythematosus, and neutrophil extracellular trap formation (Figure 8A); genes in yellow module were involved in ribosome, coronavirus disease—COVID-19, RNA transport, etc. (Figure 8B). The genes in red module were involved in VEGF signaling pathway, viral carcinogenesis, and Yersinia infection (Figure 8C), and the genes in salmon module were involved in hepatitis C, influenza A, measles, human papillomavirus infection, and biosynthesis of cofactors (Figure 8D).
KEGG pathway analysis of genes in coexpression modules associated with relapsed GD. (A) Brown module; (B) Yellow module; (C) Red module; (D) Salmon module.
Validation of Hub Genes
The gene with gene significance value greater than 0.4 and module membership greater than 0.9 is considered as the hub gene. We verified the expression of eight mRNAs and six lncRNAs of interest by PCR, including RPL8 in brown module, PARP9, RSAD2, OAS2, USP18, and IFIH1 in salmon module, NFAT5 and DROSHA in red module, and NONHSAT093153.2, NONHSAT209004.1, NONHSAT101116.2, NONHSAT161865.1, NONHSAT118924.2, and NONHSAT077537.2 in tan module. As shown in Figure 9, our results showed that the expression of RPL8, OAS2, NFAT5, DROSHA, NONHSAT093153.2, NONHSAT118924.2, and NONHSAT209004.1 was significantly decreased in GD group (p < 0.001, p < 0.001, p < 0.01, p < 0.05, p < 0.001, p < 0.05, and p < 0.01, respectively). However, there was no significant difference in the expression level of PARP9, RSAD2, OAS2, USP18 A, NONHSAT101116.2, NONHSAT077537.2 (data not shown), and NONHSAT161865.1 (data not shown) between the two groups (all p > 0.05.)
Validation of key lncRNAs and mRNAs by qRT-PCR. The expression of RPL8, OAS2, NFAT5, DROSHA, NONHSAT093153.2, NONHSAT118924.2, and NONHSAT209004.1 was significantly decreased in GD group (p < 0.001, p < 0.001, p < 0.01, p < 0.05, p < 0.001, p < 0.05, and p < 0.01, respectively). No significant difference was found in the expression level of PARP9, RSAD2, OAS2, USP18 A, and NONHSAT101116.2 between GD and normal control (NC) groups. ∗p < 0.05, ∗∗p < 0.01, ∗∗∗ or ∗∗∗∗ p < 0.001, ns, not significant.
Discussion
Considering the high recurrence rate of GD after treatment, research on clarifying the pathogenesis of GD is an important but challenging task. The lymphocyte infiltration results in the destruction of thyroid tissues and amplifies the extent of autoimmune response. Among the lymphocytes, CD4+ T cells play an important role in the pathogenesis of GD, which mainly includes Th1, Th2, Th17, Th22, Tfh cells, and Treg. Emerging studies have shown that the imbalance between Th1 and Th2 cells leads to GD (Zemmour et al., 2017). Recent studies have revealed that the abnormal expression of Th17 (Su et al., 2020; Zake et al., 2021), Th22 (Peng et al., 2013; Vitales-Noyola et al., 2017), Tfh cells (Liu et al., 2018), and Treg (Chen et al., 2021a) is associated with GD pathogenesis. The above findings illustrated that the dysfunction of CD4+ T cells plays a vital role in the development of GD. Nevertheless, the underlying mechanisms of CD4+ T cell dysfunction need to be further clarified. In the present study, we generated a signature profile of numerous lncRNAs and mRNAs in CD4+ T cells of relapsed GD patients compared with healthy controls by high-throughput sequencing technologies.
We obtained a total of 13 co-expression modules by WGCNA analysis. Among them, five modules including brown, yellow, tan, red, and salmon module were the main modules involved in GD, containing 144, 105, 49, 90, and 43 genes, respectively.
Currently, widely used GO analysis is very powerful in classifying various biological entities into functional related groups (Rue-Albrecht et al., 2016). In the present study, we also used GO analysis to study the biological functions of genes in the five modules.
Our results showed that the genes in brown module were mainly enriched in negative regulation of gene expression and epigenetic and chromatin silencing; the genes in red module were mainly enriched in inositol phosphate-mediated signaling, histone H3-K4 methylation, calcineurin-NFAT signaling cascade, and calcineurin-mediated signaling; genes in salmon module were mainly enriched in type I interferon signaling pathway, response to virus, etc. The genes in yellow module were enriched in viral transcription, viral gene expression, etc. These findings suggest that multiple biological processes are involved in the pathogenesis of relapsed GD.
Among the identified hub genes, we found the expression of three lncRNAs (NONHSAT093153.2, NONHSAT118924.2, and NONHSAT209004.1) and four mRNAs (RPL8, OAS2, NFAT5, and DROSHA) were significantly downregulated in the relapsed GD patients, suggesting that these genes are involved in the occurrence of recurrent GD.
NONHSAT093153.2, NONHSAT118924.2, and NONHSAT209004.1 were firstly investigated in the relapsed GD patients. RPL8, a member of ribosomal proteins, is a component of the 60S ribosomal subunit in eucaryotic cells (Sun et al., 2015). It has been reported that RPL8 was related to multiple sclerosis (MS) and was a potential biomarker of MS(Chen et al., 2021b). RPL8 has not been reported in GD, and our study suggests that RPL8 was significantly decreased in GD and is worthy of further study. NFAT5 is a member of the Rel family of transcriptional factors (Lopez-Rodriguez et al., 1999). Recent emerging studies have reported the role of NFAT5 in the development and activation of macrophages and T cells (Lee et al., 2019). NFAT5 can induce the activation of pathogenic pro-inflammatory macrophages and pathogenic Th17 cells (Choi et al., 2016; Aramburu and López-Rodríguez, 2019). Numerous studies found that increased expression of NFAT5 was involved in inflammatory and autoimmune diseases (Choi et al., 2017; Choi et al., 2018). OAS2 is a potential new sensitive biomarker, which can predict the activity and severity of psoriasis, and can evaluate the clinical treatment efficacy (Zhou et al., 2020). OAS2 can also be considered as biomarker gene for systemic lupus erythematosus (SLE) diagnosis (Fang et al., 2021). OAS family genes including OAS2 were revealed to be closely related to lupus nephritis (Cao et al., 2020). Drosha is RNase III enzyme necessary for most miRNA biogenesis. Study has found that the Drosha polymorphism was associated with GD development (Saeki et al., 2016). Our results showed that Drosha expression was significantly decreased in relapsed GD patients.
Although IFIH1, RSAD2, and PARP9 were found to be associated with a variety of autoimmune or inflammatory disease development, such as SLE, RA, Sjögren’s syndrome (SS), type 1 diabetes (T1D), and AITD (Frommer and Kahaly, 2021; Zedan et al., 2021), we did not find that these genes were differentially expressed between the recurrent GD group and the normal group. Ubiquitin-specific peptidase 18 (USP18) plays a crucial role in the development of Th17 cells and can regulate the differentiation and function of Treg cells (Yang et al., 2021). In our study, no significant difference was found in the USP18 expression level between relapsed GD group and NC.
The present study also has some limitations. Firstly, we did not further explore the molecular mechanism of the hub genes in relapsed GD. Secondly, the number of samples we recruited to verify gene expression was too small, because blood samples from patients with recurrent GD are very difficult to collect. Thus, the hub gene expression and the potential role of them in GD still need to be further investigated and validated in more samples.
In summary, this study is the first to explore the coexpression gene networks including lncRNAs and mRNAs related to relapsed GD through WGCNA analysis with large sample size. Our study mainly finds involvement of the key gene co-expression modules, functional biological pathways, and hub genes in the development of relapsed GD. Although the potential mechanism of functional pathways and hub genes in relapsed GD still needs to be further investigated, these initial and innovative findings provide new insights into the pathogenesis of relapsed GD undoubtedly.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI [accession: PRJNA763124].
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Zhoupu Hospital. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
QY, ZS, BW, and JZ contributed for the experimental planning; QY and ZS performed the experiments; QY and ZS analyzed the data; QY wrote the main manuscript; RS and JZ collected samples; JZ revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (No. 81873636), Shanghai Medical Key Specialty (No. ZK2019C09) and Shanghai Health Medical College Clinical Research Center (20MC20200002) and Talent Youth Cultivation Plan of Pudong New District (No. PWRq2020-11).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
The authors would like to thank all participants in the study.
ReferencesAramburuJ.López-RodríguezC. (2019). Regulation of Inflammatory Functions of Macrophages and T Lymphocytes by NFAT5. Front. Immunol.10, 535. 10.3389/fimmu.2019.00535BrentG. A. (2008). Graves' Disease. N. Engl. J. Med.358 (24), 2594–2605. 10.1056/nejmcp0801880CaleyD. P.PinkR. C.TrujillanoD.CarterD. R. F. (2010). Long Noncoding RNAs, Chromatin, and Development. The Scientific World Journal.10, 90–102. 10.1100/tsw.2010.7CaoY.MiX.WangZ.ZhangD.TangW. (2020). Bioinformatic Analysis Reveals that the OAS Family May Play an Important Role in Lupus Nephritis. J. Natl. Med. Assoc.112 (6), 567–577. 10.1016/j.jnma.2020.05.006ChenX.HouH.QiaoH.FanH.ZhaoT.DongM. (2021). Identification of Blood-Derived Candidate Gene Markers and a New 7-gene Diagnostic Model for Multiple Sclerosis. Biol. Res.54 (1), 12. 10.1186/s40659-021-00334-6ChenZ.LiuY.HuS.ZhangM.ShiB.WangY. (2021). Decreased Treg Cell and TCR Expansion Are Involved in Long-Lasting Graves' Disease. Front. Endocrinol.12, 632492. 10.3389/fendo.2021.632492ChoiS. Y.LeeH. H.LeeJ. H.YeB. J.YooE. J.KangH. J. (2016). TonEBP Suppresses IL-10-mediated Immunomodulation. Sci. Rep.6, 25726. 10.1038/srep25726ChoiS. Y.LimS. W.SalimiS.YooE. J.Lee-KwonW.LeeH. H. (2018). Tonicity-Responsive Enhancer-Binding Protein Mediates Hyperglycemia-Induced Inflammation and Vascular and Renal Injury. Jasn29 (2), 492–504. 10.1681/asn.2017070718ChoiS.YouS.KimD.ChoiS. Y.KwonH. M.KimH.-S. (2017). Transcription Factor NFAT5 Promotes Macrophage Survival in Rheumatoid Arthritis. J. Clin. Invest.127 (3), 954–969. 10.1172/jci87880FangQ.LiT.ChenP.WuY.WangT.MoL. (2021). Comparative Analysis on Abnormal Methylome of Differentially Expressed Genes and Disease Pathways in the Immune Cells of RA and SLE. Front. Immunol.12, 668007. 10.3389/fimmu.2021.668007FrommerL.KahalyG. J. (2021). Type 1 Diabetes and Autoimmune Thyroid Disease-The Genetic Link. Front. Endocrinol.12, 618213. 10.3389/fendo.2021.618213GuttmanM.DonagheyJ.CareyB. W.GarberM.GrenierJ. K.MunsonG. (2011). lincRNAs Act in the Circuitry Controlling Pluripotency and Differentiation. Nature477 (7364), 295–300. 10.1038/nature10398IddahM. A.MachariaB. N. (2013). Autoimmune Thyroid Disorders. ISRN Endocrinol.2013, 509764. 10.1155/2013/509764KotwalA.StanCurrentM.TreatmentsFuture. (2018). Current and Future Treatments for Graves' Disease and Graves' Ophthalmopathy. Horm. Metab. Res.50 (12), 871–886. 10.1055/a-0739-8134LeeN.KimD.KimW.-U. (2019). Role of NFAT5 in the Immune System and Pathogenesis of Autoimmune Diseases. Front. Immunol.10, 270. 10.3389/fimmu.2019.00270LiuY.YuanX.LiX.CuiD.XieJ. (2018). Constitutive Changes in Circulating Follicular Helper T Cells and Their Subsets in Patients with Graves' Disease. J. Immunol. Res.2018, 8972572. 10.1155/2018/8972572Lopez-RodriguezC.AramburuJ.RakemanA. S.RaoA. (1999). NFAT5, a Constitutively Nuclear NFAT Protein that Does Not Cooperate with Fos and Jun. Proc. Natl. Acad. Sci.96 (13), 7214–7219. 10.1073/pnas.96.13.7214MercerT. R.DingerM. E.MattickJ. S. (2009). Long Non-coding RNAs: Insights into Functions. Nat. Rev. Genet.10 (3), 155–159. 10.1038/nrg2521MorshedS. A.LatifR.DaviesT. F. (2012). Delineating the Autoimmune Mechanisms in Graves' Disease. Immunol. Res.54 (1-3), 191–203. 10.1007/s12026-012-8312-8PengD.XuB.WangY.GuoH.JiangY. (2013). A High Frequency of Circulating Th22 and Th17 Cells in Patients with New Onset Graves' Disease. PLoS One8 (7), e68446. 10.1371/journal.pone.0068446QiaoY. Q.HuangM. L.XuA. T.ZhaoD.RanZ. H.ShenJ. (2013). LncRNA DQ786243 Affects Treg Related CREB and Foxp3 Expression in Crohn's Disease. J. Biomed. Sci.20, 87. 10.1186/1423-0127-20-87Rue-AlbrechtK.McGettiganP. A.HernándezB.NalpasN. C.MageeD. A.ParnellA. C. (2016). GOexpress: an R/Bioconductor Package for the Identification and Visualisation of Robust Gene Ontology Signatures through Supervised Learning of Gene Expression Data. BMC Bioinformatics17, 126. 10.1186/s12859-016-0971-3SaekiM.WatanabeM.InoueN.TokiyoshiE.TakuseY.ArakawaY. (2016). DICER and DROSHA Gene Expression and Polymorphisms in Autoimmune Thyroid Diseases. Autoimmunity49 (8), 514–522. 10.1080/08916934.2016.1230846SuX.YinX.LiuY.YanX.ZhangS.WangX. (2020). Gut Dysbiosis Contributes to the Imbalance of Treg and Th17 Cells in Graves' Disease Patients by Propionic Acid. J. Clin. Endocrinol. Metab.105 (11), dgaa511. 10.1210/clinem/dgaa511SunL.LiJ.YanB. (2015). Gene Expression Profiling Analysis of Osteosarcoma Cell Lines. Mol. Med. Rep.12 (3), 4266–4272. 10.3892/mmr.2015.3958Vitales-NoyolaM.Ramos-LeviA. M.Martínez-HernándezR.Serrano-SomavillaA.Sampedro-NuñezM.González-AmaroR. (2017). Pathogenic Th17 and Th22 Cells Are Increased in Patients with Autoimmune Thyroid Disorders. Endocrine57 (3), 409–417. 10.1007/s12020-017-1361-yYangL.JingY.KangD.JiangP.LiN.ZhouX. (2021). Ubiquitin-specific Peptidase 18 Regulates the Differentiation and Function of Treg Cells. Genes Dis.8 (3), 344–352. 10.1016/j.gendis.2020.03.004ZakeT.KalereI.Upmale‐EngelaS.SvirskisS.GersoneG.SkestersA. (2021). Plasma Levels of Th17‐associated Cytokines and Selenium Status in Autoimmune Thyroid Diseases. Immun. Inflamm. Dis.9 (3), 792–803. 10.1002/iid3.433ZedanM. M.AttiaZ. R.Abd El AzeemR. A.MutawiT. M.El ShehawyA. S.BakrA. (2021). Genetic Polymorphisms in Genes Involved in the Type I Interferon System (IFIH1/MDA-5, TNFAIP3/A20, and STAT4): Association with SLE Risk in Egyptian Children and Adolescents. JirVol. 14, 3349–3358. 10.2147/jir.s309008ZemmourD.PratamaA.LoughheadS. M.MathisD.BenoistC. (2017). Flicr, a Long Noncoding RNA, Modulates Foxp3 Expression and Autoimmunity. Proc. Natl. Acad. Sci. USA114 (17), E3472–E3480. 10.1073/pnas.1700946114ZhangJ.GaoF. F.XieJ. (2021b). LncRNA linc00152/NF-kappaB Feedback Loop Promotes Fibroblast-like Synovial Cells Inflammation in Rheumatoid Arthritis via Regulating miR-103a/TAK1 axis and YY1 Expression. Immun. Inflamm. Dis.9, 681-693. 10.1002/iid3.417ZhangY.XieL.LuW.LiJ.LiY.ShaoY. (2021a). LncRNA MIAT Enhances Systemic Lupus Erythematosus by Upregulating CFHR5 Expression via miR-222 Degradation. cejoi46 (1), 17–26. 10.5114/ceji.2021.105242ZhouY.WangP.YanB.-X.ChenX.-Y.LandeckL.WangZ.-Y. (2020). Quantitative Proteomic Profile of Psoriatic Epidermis Identifies OAS2 as a Novel Biomarker for Disease Activity. Front. Immunol.11, 1432. 10.3389/fimmu.2020.01432