SciELO - Scientific Electronic Library Online

vol.53Curcumin prevents neuronal loss and structural changes in the superior cervical (sympathetic) ganglion induced by chronic sleep deprivation, in the rat modelHypoxia inducible factor-1 (HIF-1α) reduced inflammation in spinal cord injury via miR-380-3p/ NLRP3 by Circ 0001723 índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados




Links relacionados

  • En proceso de indezaciónCitado por Google
  • No hay articulos similaresSimilares en SciELO
  • En proceso de indezaciónSimilares en Google


Biological Research

versión impresa ISSN 0716-9760

Biol. Res. vol.53  Santiago  2020  Epub 05-Ago-2020 

Research Article

Microarray profiling and functional analysis of differentially expressed plasma exosomal circular RNAs in Graves' disease

Ying Sun1 

Wei Wang1 

Yuxiao Tang1 

Daping Wang1 

Liang Li2 

Min Na3 

Guantong Jiang4 

Qian Li5  * 

Shulin Chen1 

Jin Zhou1  *

1Department of Endocrinology, Affiliated Yantai Yuhuangding Hospital of Qingdao University Medical College, Yantai, Shandong, China.

2Department of Neurosurgery, Affiliated Zhongshan Hospital of Dalian University, Dalian, Liaoning, China.

3Department of Radiology, Dalian Sixth People's Hospital, Dalian, Liaoning, China.

4Department of Endocrinology and Metabolism, Binzhou Medical University Hospital, Binzhou, Shandong, China.

5Department of Scientific Research, Affiliated Yantai Yuhuangding Hospital of Qingdao University Medical College, Yantai, Shandong, China.



Circulating RNA (circRNA) regulates various bioactivities in cells. A better understanding of the exosomal circRNA can provide novel insights into the pathogenesis and treatment of Graves' disease (GD). We aimed to profile the differentially expressed circRNAs (DEcRs) in plasma exosomes of patients with GD and speculate and probe the functions of the DEcR by comprehensive bioinformatics analyses.


Serum exosomes were isolated from five primary GD patients and five healthy controls via ultracentrifugation. After verification with transmission electron microscopy, exosome samples were subjected to microarray profiling using human circRNA microarrays. Two up-regulated and two down-regulated DEcRs were selected for validation in plasma exosomes from 20 GD and 20 healthy control participants using reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR). The circRNA/microRNA/mRNA interaction network was then assembled and the analysis of the Gene Ontology and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways was utilized to predict the potential functions of the DEcR associated genes.


There were 15 DEcRs revealed in primary GD cases. The intronic circRNA hsa_circRNA_000102 was confirmed as an up-regulated component in plasma exosomes from patients with GD. The circRNA/microRNA/mRNA interaction network unveiled the most potential targeting microRNAs of hsa_circRNA_000102 and its associated genes. The functional analyses predicted involvement of hsa_circRNA_000102 associated genes in pathways of immune system activation, such as viral infection and interferon-beta signaling.


hsa_circRNA_000102 is a differentially up-regulated plasma exosomal circRNA in patients with GD. Our study highlights multiple pathways, particularly virus infection and interferon-beta signaling, for mediating immune activation in Graves' disease.

Keywords: Circular RNA; Graves' disease; Exosomes; hsa_circRNA_000102; Immune activation


Graves' disease (GD) is an autoimmune disorder characterized by typical clinical manifestations with unique association of diffuse goiter, ophthalmopathy, and pretibial myxedema [1]. GD is the most common cause of hyperthyroidism, identified in patients spanning the globe [2]. Timely diagnosis and appropriate treatments are critical for the prognosis of patients with GD. A recent observational study in Italy reported that the clinical phenotype of GD is becoming milder, likely due to increased awareness on thyroid disease and the availability of sensitive diagnostic tests, in conjunction with timely and effective treatments [3]. However, cases of extreme hyperthyroidism surface occasionally. Additionally, current treatments for GD, including antithyroid medications, radioactive iodine, and surgery (thyroidectomy), often result in adverse effects. Therefore, the importance of improving the understanding of the molecular and intermolecular interactions and the key signaling pathways affected by Graves' disease persists.

The pathogenesis of Graves' disease remains insufficiently elucidated. However, studies revealed GD as a B cell-mediated, T cell-dependent autoimmune disease with multiple types of immune cell involvement [4]. In patients with GD, overproduction of thyroid hormones and thyroid hyperplasia are often caused by stimulatory autoantibodies against the thyroid-stimulating hormone receptor (TSHR), also known as the thyroid-stimulating antibody (TSAb) [5]. Over the course of GD, communication between different immune cells is of great importance [6]. Aside from direct cell-to-cell contact, a novel and perhaps even more complex mode of immune regulation is provided by exosomes that are released by dendritic cells, B cells, T cells, and other cell types [7]. With the diameters ranging from 30 to 150 nm, exosomes are small membrane vesicles of endocytic origin, which are secreted by almost all cell types, displaying diverse functions dependent upon their origin [8]. Exosomes are rich in content and contain protein, lipid, and nucleic acid, which are involved in immune responses, antigen presentation, intercellular communication, protein and RNA transport and other physiological processes [9]. Current research into the genetic etiology of GD suggests that the “hidden” heredity of GD is likely associated with rare genomic variants including noncoding RNAs (ncRNAs), microRNAs (miRNAs), and/or epigenetic factors [10]. The exosomal circular RNAs (circRNAs), rather than those unprotected ncRNAs in blood, potentially represent novel bridges for communication of various immune cells in GD.

The first circular RNA (circRNA) was identified as ‘scrambled exons’ more than 40 years ago [11], but it was only recently more and more evidence of their expression in various diseases have been emerging. Unlike traditional linear RNA, which contains 5' and 3' ends, circRNAs have unique structures resulting from a 3' to 5' end-joining event (backsplicing). Meanwhile, two different models of circRNA biogenesis have been described, the lariat or exon skipping model and the direct backsplicing model Meanwhile, two different models of circRNA biogenesis have been described, the lariat or exon skipping model and the direct backsplicing model. [12]. CircRNAs can be derived from exon, intron, untranslated or intergenic regions of the genome, and a large amount of them have been identified in mouse and human by high-throughput sequencing and bioinformatics analyses [13]. CircRNA is produced by special selective shear and is abundant in the cytoplasm of eukaryotic cells. It expresses with tissue, timing and disease specificity [14]. Though to the multiple functions of circRNAs recently inferred, the most important is the involvement of circRNAs in post-transcriptional regulation, similar to the roles of an endogenous RNA or miRNA sponge in competitively inhibiting RNA/miRNA transcriptional regulation [15, 16]. In addition, circRNA can also regulate variable shear or transcription and modify the expression of its related genes [17]. Therefore, circRNAs are expected to provide a new molecular basis of prediction, diagnosis and treatment of complex diseases. However, the means and extent to which exosomal circRNAs play a role in the pathogenesis of GD remains completely unknown.

In this study, the implications of circRNA in GD and exosomal circular RNAs were profiled in the serum samples from patients with and without GD through human circRNA microarray. Comprehensive bioinformatics analyses were conducted to explore and probe the functions of the differentially expressed circRNAs.


Validation of study cohorts and isolated exosomes

This study involved two cohorts of subjects: the discovery cohort and the validation cohort. The clinical characteristics of the participants in the discovery cohort are summarized in Additional file 1: Table S2. The mean age of patients with Graves' disease was 36.6 ± 11.5 years as compared to mean age of 37.8 ± 9.6 years among controls (P = 0.86). As expected, Graves' disease patients had significantly elevated FT3, elevated FT4, and suppressed TSH levels in comparison to controls (Additional file 1: Table S2). Notably, the concentration of plasma TRAb (TSH receptor antibodies, also called TSAb), the specific thyroid autoantibody produced in GD (29.49 ± 14.77 IU/L), was remarkably higher in GD patients than in healthy control subjects (0.55 ± 0.32 IU/L). In the validation cohort (Additional file 1: Table S3), there is no difference in the age of patients with Graves' disease (37.55 ± 12.18 years) and control groups (37.15 ± 10.89 years) (P = 0.91). TSH(GD group: 0.0145 ± 0.0294 IU/ ml vs. control group: 1.851 ± 0.8749 IU/ml), FT4(GD group: 90.66 ± 10.71 pmol/L vs. control group: 16.32 ± 2.06 pmol/L), FT3(GD group: 38.2 ± 10.83 pmol/L vs. control group: 4.85 ± 0.93 pmol/L), had a significant differences between two groups(P < 0.0001). There are significant differences in TRAb between GD (26.13 ± 12.9 IU/L) and healthy group (0.71 ± 0.46 IU/L) (P < 0.0001). The median TgAb level of GD and healthy control groups were 1887(127.6-4000) IU/ml and 20.82(10-114.1) IU/ml, respectively (P < 0.0001). Next, we verified isolated exosomes through observation of their morphology by transmission electron microscopy. As shown in Additional file 1: Figure S1, these exosomes displayed irregular spheres with a clearly defined and relatively intact membrane. They have diameters ranging from 30 nm to 100 nm, distinguishable from other much larger extracellular vesicles like microvesicles and apoptotic bodies.

Identification of differentially expressed plasma exosomal circRNAs from patients with GD and healthy control subjects

The plasma exosomes isolated from the discovery cohort were subjected to the microarray profiling of the differentially expressed circRNAs. A number of the 13,000 total circRNAs detected, developed an altered expression. As shown in the scatter plot, denoting the means of the expression values of the two groups (Fig. 1a), many circRNAs, from patients with GD, exhibit more than a twofold up-regulation or down-regulation of plasma exosomes compared to healthy control subjects. There are 15 significantly differentially expressed circRNAs that met the criteria of exhibiting fold changes > 1.2 and P-values < 0.05, as shown in the volcano plot representation of differentially expressed transcripts and their statistical significance (Fig. 1b). Among them, six circRNAs were up-regulated and nine circRNAs were down-regulated (Table 1).

Fig. 1 Differential expression of circRNAs in plasma exosomes from patients with GD and healthy control subjects. a Scatter plots with the mean expression values of samples in each group. Raw junction reads for all samples were normalized by total mapped read numbers and were transformed by log2. Each dot represents a single circRNA. The exhibiting fold change is 2.0 of Line 1 and Line 3, whereas the dots on Line 2 indicate circRNAs with equal expression between GD group and healthy control group. b Volcano plot representation of differentially expressed transcripts and their statistical significance. The fold changes are represented in log2 scale as depicted on the x-axis, whereas the - log10 P-value is depicted on the y-axis. Transcripts with greater statistical significance are higher in the plot. The red circles represent circRNA genes that show differential expression with fold changes > 1.2 and P-values < 0.05 between the GD and normal control exosomes samples 

Table 1 The list of differential expressed circRNAs in plasma exosome samples from patients with GD and healthy control subjects 

circRNA ID P-value FC (abs) Regulation circRNA_type GeneSymbol
hsa_circRNA_102059 0.038876 1.38207 Up Exonic MED1
hsa_circRNA_405965 0.01488 1.884644 Up Exonic CCNT2-AS1
hsa_circRNA_008389 0.010443 1.941364 Up Exonic DNAJC11
hsa_circRNA_000102 0.015919 1.543085 Up Intronic AKNAD1
hsa_circRNA_402094 0.007404 1.466367 Up Exonic HKR1
hsa_circRNA_005133 0.04595 1.301035 Up Exonic DDR1
hsa_circRNA_405709 0.00715 1.203396 Down Exonic MIER2
hsa_circRNA_072697 0.026879 1.267682 Down Exonic PPWD1
hsa_circRNA_025500 0.03979 1.288536 Down Exonic DUSP16
hsa_circRNA_001067 0.038894 1.214824 Down Intronic RP11-351M8.1
hsa_circRNA_043898 0.019398 1.219508 Down Exonic EZH1
hsa_circRNA_004939 0.049734 1.231611 Down Sense overlapping IL4R
hsa_circRNA_104619 0.007381 1.20487 Down Exonic PRKDC
hsa_circRNA_407279 0.031666 1.277486 Down Exonic SYTL5
hsa_circRNA_001454 0.027548 1.267223 Down Exonic KIAA0922

FC fold change (absolute value)

hsa_circRNA_000102 is the validated differentially expressed plasma exosomal circRNA in GD

Among 15 significantly differentially expressed circRNAs, 12 are exonic, two are intronic, and one is a combination of both (Table 1). We examined the genes associated with the back splicing-mediated generation of these circRNAs, and found MED1, AKNAD1, PPWD1 and IL4R were related to nuclear co-receptor, transcription factor, mRNA splicing and cytokine receptor pathways, respectively. These molecules and pathways possess the ability for an immune activation, as observed in patients with GD. Therefore, the corresponding four circRNAs, including two up-regulated circRNAs (hsa_circRNA_000102, hsa_circRNA_102059) and two down-regulated circRNAs (hsa_circRNA_004939, hsa_circRNA_072697), were initially selected for further verification by RT-qPCR. As shown in Fig. 2, the consistency between RT-qPCR results and microarray analysis data had been demonstrated only for hsa_circRNA_000102. However, the expression levels of hsa_circRNA_102059, hsa_circRNA_004939, and hsa_circRNA_072697 did not show significant differences between plasma exosomal samples from patients with GD and healthy control subjects (Fig. 2a-c). Furthermore, it was observed that expression levels of hsa_circRNA_000102/β-actin(r = 0.5983, P = 0.0053), hsa_circRNA_000102/18srRNA(r = 0.5462, P = 0.0127), hsa_circRNA_000102/ GAPDH(r = 0.4526, P = 0.0451) correlate significantly with TRAb levels (Fig. 2d-f). Therefore, our next stage of bioinformatics analyses focused on delineating the intermolecular interactions and the key signaling pathways associated with hsa_circRNA_000102.

Fig. 2 Validation of selected candidate differentially expressed circRNAs by Real-time qPCR. Data represents the mean ± SD of three independent experiments.n = 10 for each group. *P < 0.05 

Prediction of the circRNA/microRNA/mRNA interaction network for hsa_circRNA_000102

Recent studies on circRNA/microRNA interactions have indicated that circRNAs play a key role in the regulation of gene expression by interacting with their target microRNAs (miRNAs) [18]. To find the potential target miRNAs and mRNA associated with GD, a circRNA/ microRNA/mRNA interaction network was constructed to predict the intermolecular interaction for the experimentally confirmed circRNA (hsa_circRNA_000102). An Arraystars homemade miRNA target prediction software based on TargetScan and miRanda was utilized (Fig. 3). There were nine miRNAs (hsa-miR-3151-5p, hsa-miR-1227-5p, hsa-miR-194-3p, hsa-miR-1296-3p, hsa-miR-3688-3p, hsa-miR-7112-3p, hsa-miR-8063, hsa-miR-4512 and hsa-miR-7848-3p) and numerous mRNAs predicted to be interacting with hsa_circRNA_000102 directly or indirectly.

Fig. 3 The circRNA-miRNA-mRNA regulatory network of hsa_circRNA_000102. CircRNA, miRNA, and mRNA are indicated as spheres in brown, red, blue color, respectively. CircRNA circular RNA, miRNA microRNA, mRNA messenger RNA 

GO analysis and KEGG pathway analysis of hsa_circRNA_000102 associated genes

Predicted genes in the circRNA/microRNA/mRNA interaction network of hsa_circRNA_000102 were used for GO function analysis to annotate and speculate the potential functions of this circRNA. GO analysis was divided into three parts: biological process (BP), cell component (CC) and molecular function (MF) (Fig. 4). GO analysis of BP showed that hsa_circRNA_000102 was significantly linked to Poly(A)+ mRNA export from nucleus, cellular response to interferon-beta, mRNA-containing ribonucleoprotein complex export from nucleus, mRNA export from nucleus, RNA transport, nucleic acid transport, regulation of production of molecular mediator of immune response, establishment of RNA localization, and positive regulation of cytokine biosynthetic process. Significant GO cell component terms of hsa_circRNA_000102 showed that it was associated with endocytic vesicle lumen, endolysosome membrane, Cul2-RING ubiquitin ligase complex, endolysosome, high-density lipoprotein particle, lipoprotein particle, plasma lipoprotein particle, cytoplasmic vesicle lumen, vesicle lumen, and ubiquitin ligase complex. For molecular function, hsa_circRNA_000102 was associated with transition metal ion binding, receptor ligand activity, phosphatidic acid binding, receptor regulator activity, protein transporter activity, metallocarboxypeptidase activity, tumor necrosis factor receptor binding, chemoattractant activity, molecular function regulator, and Zinc ion binding. The top 10 KEGG pathways showed that hsa_circRNA_000102 associated genes might be involved in Herpes simplex infection, Influenza A, RNA transport, mRNA surveillance pathway, Ribosome biogenesis in eukaryotes, Hepatitis B, steroid hormone biosynthesis, necroptosis, NOD-like receptor signaling pathway and Huntington disease (Fig. 5).

Fig. 4 The top 10 significantly enriched terms of biological process by GO analysis of hsa_circRNA_000102 associated genes. GO analysis was divided into three parts: biological process, cell component and molecular function 

Fig. 5 The top 10 significant enriched pathways by KEGG pathway analysis of hsa_circRNA_000102 associated genes 


In this study, we profiled the plasma exosomal circRNAs from patients with primary GD and healthy control subjects by human circRNA microarray. Before circRNA labeling and hybridization, RNAse R is used to efficiently digesting nearly all linear RNA species. Although RNAse R treatment may not exclude all entities, it is one of few ways at current available for enriching circRNAs. Our initial screening resulted in 15 differentially expressed circRNAs. One of the four selected candidate circRNAs, hsa_circRNA_000102, was verified to be significantly enriched in the plasma exosomes of patients with GD. Our comprehensive bioinformatics analyses, including prediction of a circRNA/microRNA/mRNA interaction network for hsa_circRNA_000102 and GO and KEGG pathway annotations of hsa_circRNA_000102 associated genes, revealed potential pathways correlated to the pathogenesis of GD.

hsa_circRNA_000102 is a circular intronic RNA (ciRNA), a new type of circRNA in human cells that is derived from introns discovered by Zhang et al. and others [19]. These RNAs escape debranching and depend on consensus RNA elements near the 5' splice site and the branchpoint for proper processing [19]. Interestingly, ciRNAs partially accumulate to their sites of synthesis and colocalize with the elongation RNA polymerase II (Pol II). Moreover, the existence of ciRNAs might influence the transcription of their parent genes. For example, they are specifically associated with phosphorylated Pol II whereby a depletion of ciRNAs may lead to a significant reduction in the transcription of parent genes.

Therefore, the higher level of hsa_circRNA_000102 in plasma exosomes of patients with GD, indicates the potentially elevated expression of its parent gene, AKNAD1, which encodes a protein with a domain found in an AT-hook-containing transcription factor. Though being not yet well-characterized, AKNAD1 is considered to involve in biological processes of “microtubule-based movement” and “microtubule-based process'! Till now, the relation between gene AKNAD1 and circRNA has not been discussed. Only one research has identified a significant copy number variation in gene AKNAD1 which was validated to be associated with type-2 diabetes through stress of endoplasmic reticulum [20], suggesting a potential link of AKNAD1 and autoimmune diseases. However, more work on AKNAD1 in immune regulation and the underlying mechanisms on potential immune activation in GD is required to delineate the relationship between AKNAD1 expression and the pathogenesis of GD.

The roles of circRNAs have been identified in tumor signaling pathways, invasion and migration, as well as tumor immunity [20]. Previous research proposed that mRNAs, transcribed pseudogenes, and circRNAs communicate with and co-regulate each other through competition for binding to the shared miRNAs [21]. Studies have also revealed that the activity of circRNA is decided by a series of factors, including the abundance of miRNAs, RNA editing, changes in the RNA 3' untranslated region (UTR) and RNA binding proteins (RBPs) [22]. Therefore, interpreting the cross talks among circRNAs, miRNAs and mRNAs would provide novel insights into the mechanisms of gene regulation and immune cell interaction in Graves' disease. We predicted the circRNAs/miRNAs/mRNAs interaction network of hsa_circRNA_000102, and then performed GO enrichment analysis and KEGG pathway analysis of the particular circRNA associated genes. Evidence demonstrated that hsa_circRNA_000102 is related to Herpes simplex infection, Influenza A signaling and interferon-beta signaling.

Both genetic and environmental factors have been suggested to contribute to the epidemiology of Graves' disease [23]. It is not known whether the infection of a Herpes virus or Influenza A and Graves' disease could co-exist simply by chance or the former could trigger the latter. A previous study has detected Herpes virus DNA in post-operative thyroid tissue specimens of patients with autoimmune thyroid disease, including Graves' disease and Hashimoto thyroiditis [24]. The possible role of human Herpes viruses in the pathogenesis of autoimmune thyroid diseases is largely unknown and data from one study did not support their association [25]. However, another study indicated that Graves' disease is linked to infectious mononucleosis due to primary Epstein-Barr virus infection in patients [26]. This suggests that both Herpes virus and Influenza A virus can employ similar mechanisms as Epstein-Barr virus to trigger immune activation.

We speculate that the association between viral infection and GD is due to the efficient T cell activation. Viral infection potentially increases the efficiency of T cell activation through fast attraction of dendritic cells from peripheral to secondary lymphoid tissues. Furthermore, continuous activation of innate immunity against viral infection might also profoundly strengthen autoimmune disorders by promoting autoantigen release without the need for specific activation of autoreactive T cells. As previously reported [27, 28], thyroid autoimmunity and dysfunctions have been observed over the course of IFN-beta therapy. A patient with multiple sclerosis received IFN-beta-1 therapy for 22 months and ultimately developed a severe Graves' disease with persistently positive TRAb [27]. Since thyroidal side effects can be late to diagnose, patients undergoing long-term IFN-beta therapy should be closely monitored for thyroid hormones and antibodies throughout the treatment to prevent the possible occurrence of GD [27]. Taken together, although further studies are needed to clarify the relationship between infection of Herpes virus/Influenza A or IFN-beta signaling and GD, our study indicates that the impact of virus infection and IFN signaling should be taken into consideration during the prevention and/or treatment of GD.

The current investigation has several strengths and limitations. Its major strength is its novelty; to the best of our knowledge, the current study is the first one on profiling of plasma exosomal circRNAs in patients with GD. Our work also provides a feasible framework for searching and validating circular RNAs to better facilitate target identification and development of target-oriented therapeutics. The major limitation of the current study is its sample size of patients used for discovery and validation. Moreover, the patients recruited in this study were from a single institution and thus are subject to referral bias. Thirdly, there are some other limitations in this study. For examples, ribosomal RNAs were not removed, and the expression of AKNAD1 was not measured. In addition, no mock-treated samples were used during samples prep stage for generation of circRNA profile. Also, circRNA arrays we used is based on backsplice junctions predesigned for those already existing in database. Therefore, this profile is very likely to lack novel circRNAs that could have been discovered using alternative methods to generate circRNA profiles.


In this study, we performed human circRNA microarray in the plasma exosomes from patients with GD and healthy control subjects and verified hsa_circRNA_000102 as a differentially up-regulated plasma exosomal circRNA in GD. Our comprehensive bioinformatics analyses including construction of circRNA/ microRNA/mRNA interaction network, GO annotation and KEGG pathway analysis, highlighted multiple pathways, particularly the virus infection and interferon-beta signaling, in mediating immune activation in Graves' disease.

Materials and methods


Between August 2017 and October 2017, five newly diagnosed patients with GD (two males and three females) and five healthy control subjects were recruited by the Affiliated Yantai Yu-huangding Hospital of Qingdao University Medical, Shandong Province, China. The newly diagnosed GD patients had typical characteristics of hyperthyroidism: a palpable, soft, and diffuse goitre confirmed by ultrasound examination, increased thyroid radioiodine uptake, elevated serum free thyroxine (FT4), very low thyroid-stimulating hormone (TSH), and high levels of thyroid-stimulating antibody (TSAb). Age- and gender-matched healthy control subjects with normal thyroid function (TSH: 0.27-4.2 mIU/L), free of goitre and any family history of GD who meeting the exclusion criteria selected from healthy hospital staff were admitted into the study. Exclusion criteria included patients with liver disorder, renal disorder, clinical evidence of tumor, infection, any autoimmune diseases, history of intake of immunosuppressive drugs, such as mannan peptide, lentinan, interferon, corticosteroids, and cyclophosphamide, within 1 year of the study, pregnant, and lactating women. In the validation cohort, 20 patients with GD and 20 healthy control subjects were recruited adhering to the same inclusion and exclusion criteria. The research protocol was approved by Medical Ethics Committee of Affiliated Yantai Yu-huangding Hospital of Qingdao University Medical and the written informed consent was obtained from all participants.

Laboratory test

The levels of serum TSH (11731459, Roche), FT4 (07976836, Roche), FT3 (free triiodothyronine) (06437206, Roche), TSAb(04388780190, Roche), thyroid peroxidase antibody (TPOAb) (06368590190, Roche), and thyroglobulin antibody (TgAb) (06368697, Roche) were measured via chemiluminescent enzyme immunoassay specific kits (Cobas 6000, Roche).

Ultracentrifugation exosomes isolation and RNA extraction

Venous blood samples from subjects were collected in evacuated tubes with the anticoagulant Ethylenediamine tetraacetic acid dipotassium (K2EDTA) and stored at — 80 °C for exosome isolation. All blood samples were processed by centrifugation at 1000×g for 10 min. The plasma was immediately separated and transferred to a new RNase free tube. The liquid samples were centrifuged at 500×g for 5 min to eliminate cells. The supernatant was then transferred to a new polycarbonate tube and centrifuged at 2000×g for 10 min. Supernatants from the previously centrifuged samples were collected and centrifuged again at 10,000×g for 30 min to eliminate shed microvesicles (sMV, 200-1000 nm). The supernatants were then collected, filtered with 0.22 μm membrane filter (Merck Millipore), and centrifuged at 100,000×g for 2 h to obtain exosomes. For RNA isolation, an exosome pellet was washed once with 1×PBS and centrifuged at 100,000×g for 2 h. Exosomes were observed under a transmission electron microscope (TEM, TecnaiG2 Spirit 120 kV). Total RNA was extracted from exosomes using the TRIzol LS Reagent (10296028, Invitrogen Life Technologies, Carlsbad, CA, USA) according to the manufacturer's protocol.

Linear RNA digestion, circRNA labeling and hybridization

The exosomal RNAs extracted from the discovery cohort (five GD patients and five healthy control subjects) were used for microarray analysis, as previously described [29]. RNase R (RNR07250, Epicentre, Inc.) is used to digest nearly all linear RNA species and enrich circRNAs. The circRNA expression microarray slide (AS-S-CR-H-V2.0, Arraystar Human circRNA Arrays V2, 8 × 15 K; Arraystar Inc.) was used. Sample labeling and array hybridization were performed according to the manufacturer's protocol. The hybridized arrays were scanned using the Agilent Scanner G2505C (Agilent Technologies, USA).

CircRNA microarray data analysis

Raw data extraction of the scanned images was conducted using the Agilent Feature Extraction software (version A series of data processing, including quantile normalization using a log2 ratio, was performed using the R software limma package. The circRNAs containing flags in “P” or “M” in at least five out of 10 samples (defined by GeneSpring software) were retained for further differential analyses. Differentially expressed circRNAs displaying statistically significant differences between the GD and healthy control groups were identified using fold change cutoff and through Volcano Plot filtering. The “fold change” (i.e. the ratio of the group averages) between the groups for each circRNA was computed. A t-test was employed to determine the statistical significance of the differences measured. CircRNAs denoted as significantly differentially expressed, were those having fold changes > 1.2 and P-values < 0.05 between the GD and healthy control exosome samples. The package gplots and function heatmap2 in R software were used for mapping. The distance metric was set as Euclidean distance.

Real-time quantitative PCR (RT-qPCR) validation of microarray data

Total RNA from each exosome sample was quantified using the NanoDrop ND-1000 spectrophotometer. The reverse-transcription (RT) primers for circRNAs and β-actin were designed using Primer 5.0 and listed in Additional file 1: Table S1. Complementary DNA (cDNA) was synthesized using SuperScript™ III Reverse Transcriptase (18080-044, Invitrogen). CircRNAs transcripts were quantified using a SYBR Green Real-time PCR Master Mix (AS-MR-006-25, Arraystar, US). PCR was performed with the ViiA 7 Real-Time PCR System (Applied Biosystems). The thermal cycler conditions consisted of one cycle of denaturation at 95 °C for 10 min followed by 40 cycles of 95 °C for 10 s and 60 °C for 60 s. Melting curves were generated to ensure the specificity of the RT-qPCR. The expression of circRNAs was calculated by standard curve method using β-actin, 18srRNA and GAPDH as housekeeping genes. Then we compare the relationship among circRNAs and TRAb. All reactions were performed in triplicate.

Annotation for circRNA/miRNA Interaction

RNA transcripts crosstalk by competing for common microRNAs and have been termed as competing endogenous RNAs (ceRNAs) [30]. Studies identified micro-RNA response elements (MREs) as the foundation of this interaction [31]. CeRNAs include pseudogene transcripts, long non-coding RNAs (lncRNAs), circRNAs and mRNAs, all of whose transcripts can compete for the same MRE to regulate mutually. Potential microRNA targets were identified with home-made miRNA target prediction software based on TargetScan and miRanda [32, 33]. The ceRNA network was constructed through merging the common targeted miRNAs. Three conditions must exist for the ceRNA network to occur, as previously described, which include: well defined relative concentrations of ceRNAs and their microRNAs; an appropriate number of microRNAs that ceRNA can “sponge”; and knowledge that not all of the MREs on ceRNAs are equal [30]. Therefore, only the ceRNA-pair relations passing these measures were accepted as filtered. The co-expression network was created according to the correlation analysis between normalized intensity of the differential circRNAs, miRNAs and mRNAs. Pearson correlation coefficients over 0.95 were retained for the construction using CYTOSCAPE software (The Cytoscape Consortium, San Diego, CA).

Bioinformatics analysis

Gene ontology (GO) analysis was performed to explore the functional roles of circRNA-targeting genes in terms of biological processes, cellular components and molecular functions, as previously described [34]. Biological systems defined by the Kyoto Encyclopedia of Genes and Genomes (KEGG) (, was used to explore those pathways related to circRNA-targeting genes, as previously described [35].


SPSS software, version 17.0 (SPSS Inc.) was used for statistical analyses. Data are expressed as the mean ± SD (standard deviation) or median values with ranges, as indicated. An unpaired two-tailed Student's t-test, or a nonparametric Mann-Whitney test was used for the comparisons between groups. Welch's correction was added when the variances were significantly different between the two groups. Pearson's correlation coefficient was used to analyze correlation between hsa_circRNA_000102 and TRAb. A two-sided P value < 0.05 was considered to be statistically significant.


This project was supported by the National Natural Science Foundation of China (No. 81700695) and China Scholarship Council (No. 201909370036).

Availability of data and materials

All data generated or analysed during this study are included in this published article.

Ethics approval and consent to participate

The present study was approved by Medical Ethics Committee of Affiliated Yantai Yu-huangding Hospital of Qingdao University Medical and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.

Consent for publication

Not applicable.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


I would like to show my deepest gratitude to Nana Huang and Fangfang Wei for all their work. I thank all my friends for their encouragement and support.



Circulating RNA


Graves' disease


Differentially expressed circRNAs


Kyoto Encyclopedia of Genes and Genomes


Thyroid-stimulating hormone receptor


Thyroid-stimulating antibody


Noncoding RNAs




Free thyroxine


Free triiodothyronine


TSH receptor antibodies


Circular intronic RNA


RNA binding proteins


Thyroid peroxidase antibody


Thyroglobulin antibody


Endogenous RNAs


microRNA response elements


Long non-coding RNAs


1. McLachlan SM, Nagayama Y, Rapoport B. Insight into Graves' hyperthyroidism from animal models. Endocr Rev. 2005;26:800-32. [ Links ]

2. Vanderpump MP. The epidemiology of thyroid disease. Br Med Bull. 2011;99:39-51. [ Links ]

3. Bartalena L, Masiello E, Magri F, Veronesi G, Bianconi E, Zerbini F, et al. The phenotype of newly diagnosed Graves' disease in Italy in recent years is milder than in the past: results of a large observational longitudinal study. J Endocrinol Invest. 2016;39:1445-51. [ Links ]

4. Smith TJ, Hegedus L. Graves' disease. N Engl J Med. 2016;375:1552-65. [ Links ]

5. Rapoport B, Chazenbalk GD, Jaume JC, McLachlan SM. The thyrotropin (TSH) receptor: interaction with TSH and autoantibodies. Endocr Rev. 1998;19:673-716. [ Links ]

6. Ramos-Levi AM, Marazuela M. Pathogenesis of thyroid autoimmune disease: the role of cellular mechanisms. Endocrinol Nutr. 2016;63:421-9. [ Links ]

7. Lindenbergh MFS, Stoorvogel W. Antigen presentation by extracellular vesicles from professional antigen-presenting cells. Annu Rev Immunol. 2018;36:435-59. [ Links ]

8. Raposo G, Stoorvogel W. Extracellular vesicles: exosomes, microvesicles, and friends. J Cell Biol. 2013;200:373-83. [ Links ]

9. Tang MK, Wong AS. Exosomes: emerging biomarkers and targets for ovarian cancer. Cancer Lett. 2015;367:26-33. [ Links ]

10. RobertRich TF, William S, Harry S, Anthony F, Cornelia W. Clinical immunology: principles and practice. 5th ed. New York: Elsevier; 2018. [ Links ]

11. Nigro JM, Cho KR, Fearon ER, Kern SE, Ruppert JM, Oliner JD, et al. Scrambled exons. Cell. 1991;64:607-13. [ Links ]

12. Jeck WR, Sorrentino JA, Wang K, Slevin MK, Burd CE, Liu J, et al. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA. 2013;19:141-57. [ Links ]

13. Rybak-Wolf A, Stottmeister C, Glazar P, Jens M, Pino N, Giusti S, et al. Circular RNAs in the mammalian brain are highly abundant, conserved, and dynamically expressed. Mol Cell. 2015;58:870-85. [ Links ]

14. Li J, Yang J, Zhou P, Le Y, Zhou C, Wang S, et al. Circular RNAs in cancer: novel insights into origins, properties, functions and implications. Am J Cancer Res. 2015;5:472-80. [ Links ]

15. Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495:333-8. [ Links ]

16. Taulli R, Loretelli C, Pandolfi PP. From pseudo-ceRNAs to circ-ceRNAs: a tale of cross-talk and competition. Nat Struct Mol Biol. 2013;20:541-3. [ Links ]

17. Qu S, Yang X, Li X, Wang J, Gao Y, Shang R, et al. Circular RNA: a new star of noncoding RNAs. Cancer Lett. 2015;365:141-8. [ Links ]

18. Xiong DD, Dang YW, Lin P, Wen DY, He RQ, Luo DZ, et al. A circRNA-miRNA-mRNA network identification for exploring underlying pathogenesis and therapy strategy of hepatocellular carcinoma. J Transl Med. 2018;16:220. [ Links ]

19. Zhang Y, Zhang XO, Chen T, Xiang JF, Yin QF, Xing YH, et al. Circular intronic long noncoding RNAs. Mol Cell. 2013;51:792-806. [ Links ]

20. Zhong Y, Du Y, Yang X, Mo Y, Fan C, Xiong F, et al. Circular RNAs function as ceRNAs to regulate and control human cancer progression. Mol Cancer. 2018;17:79. [ Links ]

21. Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, et al. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495:384-8. [ Links ]

22. Cheng DL, Xiang YY, Ji LJ, Lu XJ. Competing endogenous RNA interplay in cancer: mechanism, methodology, and perspectives. Tumour Biol. 2015;36:479-88. [ Links ]

23. Hemminki K, Li X, Sundquist J, Sundquist K. The epidemiology of Graves' disease: evidence of a genetic and an environmental contribution. J Autoimmun. 2010;34:J307-13. [ Links ]

24. Thomas D, Liakos V, Michou V, Kapranos N, Kaltsas G, Tsilivakos V, et al. Detection of herpes virus DNA in post-operative thyroid tissue specimens of patients with autoimmune thyroid disease. Exp Clin Endocrinol Diabetes. 2008;116:35-9. [ Links ]

25. Noura AZ, Monem F. Are human herpes viruses associated with autoimmune thyroid disease? J Infect Dev Ctries. 2011;5:890-2. [ Links ]

26. Akahori H, Takeshita Y, Saito R, Kaneko S, Takamura T. Graves' disease associated with infectious mononucleosis due to primary Epstein-Barr virus infection: report of 3 cases. Intern Med. 2010;49:2599-603. [ Links ]

27. Rotondi M, Mazziotti G, Biondi B, Manganella G, Del Buono AD, Montella P et al. Long-term treatment with interferon-beta therapy for multiple sclerosis and occurrence of Graves' disease. J Endocrinol Invest. 2000;23:321-4. [ Links ]

28. Schwid SR, Goodman AD, Mattson DH. Autoimmune hyperthyroidism in patients with multiple sclerosis treated with interferon beta-1b. Arch Neurol. 1997;54:1169-90. [ Links ]

29. Li H, Li K, Lai W, Li X, Wang H, Yang J, et al. Comprehensive circular RNA profiles in plasma reveals that circular RNAs can be used as novel biomarkers for systemic lupus erythematosus. Clin Chim Acta. 2018;480:17-25. [ Links ]

30. Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PR A ceRNA hypothesis: the Rosetta stone of a hidden RNA language? Cell. 2011;146:353-8. [ Links ]

31. Phelps M, Coss C, Wang H, Cook M, Reproducibility Project: Cancer B, Reproducibility Project Cancer B. Registered report: coding-independent regulation of the tumor suppressor PTEN by competing endogenous mRNAs. Elife. 2016;5:e12470. [ Links ]

32. Huang M, Zhong Z, Lv M, Shu J, Tian Q, Chen J. Comprehensive analysis of differentially expressed profiles of lncRNAs and circRNAs with associated co-expression and ceRNA networks in bladder carcinoma. Oncotarget. 2016;7:47186-200. [ Links ]

33. Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS. MicroRNA targets in Drosophila. Genome Biol. 2003;5:R1. [ Links ]

34. Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, et al. The gene ontology (GO) database and informatics resource. Nucleic Acids Res. 2004;32:D258-61. [ Links ]

35. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27-30. [ Links ]

Supplementary information

Supplementary information accompanies this paper at

Additional file 1.

Additional figure and tables.

Received: September 18, 2019; Accepted: July 06, 2020; pub: July 29, 2020

* Correspondence:;

Authors' contributions

JZ and SC designed the research. YS, WW, YT and DW conducted the experiments, analyzed data, and wrote the manuscript. QL contributed to the study design and manuscript writing. MN edited and revised manuscript. LL analyzed data. GJ performed experiments. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Ying Sun, Wei Wang, Yuxiao Tang, Daping Wang and Jin Zhou contributed equally to this paper

Creative Commons License This is an Open Access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.