Ziyi Chen
Yufeng Liu2,3,4†- 1Department of Endocrinology, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’an, China
- 2MOE Key Laboratory for Intelligent Networks & Networks Security, School of Electronic and Information Engineering, Xi’an Jiaotong University, Xi’an, China
- 3Genome Institute, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’an, China
- 4BioBank, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’an, China
- 5Department of Endocrinology, Zhejiang Provincial People’s Hospital, Hangzhou, China
- 6Precision Medicine Center, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’an, China
Graves’ disease (GD) is a T cell-mediated organ-specific autoimmune disorder. GD patients who have taken anti-thyroid drugs (ATDs) for more than 5 years with positive anti-thyroid stimulating hormone receptor autoantibodies value were defined as persistent GD (pGD). To develop novel immunotherapies for pGD, we investigated the role of T cells in the long-lasting phase of GD. Clinical characteristics were compared between the pGD and newly diagnosed GD (nGD) (N = 20 respectively). Flow cytometric analysis was utilized to determine the proportions of Treg and Th17 cells (pGD, N = 12; nGD, N = 14). T cell receptor sequencing (TCR-seq) and RNA sequencing (RNA-seq) were also performed (pGD, N = 13; nGD, N = 20). Flow cytometric analysis identified lower proportions of Th17 and Treg cells in pGD than in nGD (P = 0.0306 and P = 0.0223). TCR-seq analysis revealed a lower diversity (P = 0.0025) in pGD. Specifically, marked clonal expansion, represented by an increased percentage of top V-J recombination, was observed in pGD patients. Interestingly, pGD patients showed more public T cell clonotypes than nGD patients (2,741 versus 966). Meanwhile, RNA-seq analysis revealed upregulation of the inflammation and chemotaxis pathways in pGD. Specifically, the expression of pro-inflammatory and chemotactic genes (IL1B, IL13, IL8, and CCL4) was increased in pGD, whereas Th17 and Treg cells associated genes (RORC, CARD9, STAT5A, and SATB1) decreased in pGD. Additionally, TCR diversity was negatively correlated with the expression of pro-inflammatory or chemotactic genes (FASLG, IL18R1, CCL24, and CCL14). These results indicated that Treg dysregulation and the expansion of pathogenic T cell clones might be involved in the long-lasting phase of GD via upregulating chemotaxis or inflammation response. To improve the treatment of pGD patients, ATDs combined therapies, especially those aimed at improving Treg cell frequencies or targeting specific expanded pathogenic TCR clones, are worth exploring in the future.
Introduction
Graves’ disease (GD) is an organ-specific autoimmune disease that arises due to the breakthrough of tolerance to thyroid stimulating hormone receptor (TSHR). GD has an annual incidence of 20 to 50 cases per 100,000 persons, and females are more susceptible than males (1). This disease is manifested as weight loss, fatigue, heat intolerance, tremor, and palpitations (1). After anti-thyroid drugs (ATDs) withdrawal (12–18 months), 40–50% of GD patients achieve durable remission with euthyroidism, while others will be trapped in long-lasting disease phase (1, 2). Those patients who have been taking ATDs treatment for more than 5 years with positive anti-thyroid stimulating hormone receptor autoantibodies (TRAb) value are classified as persistent GD (pGD) cases, namely refractory GD (rGD) patients in previous studies (3). Although radioactive iodine or thyroidectomy are alternative treatments for pGD patients, either poses a risk of permanent hypothyroidism. Therefore, many pGD patients prefer to continue on ATDs for decades to maintain euthyroidism (4), and this continuous treatment imposes enormous financial and psychological burdens. In addition, based on that GD involved a progressive accumulation of changes in the thyroid glands from short and long clinical course cases (5), there are potential perpetuating factors that maintain the autoimmune response. Considering the unsatisfactory treatment effects of current pGD patients, there is an urgent need to develop novel therapies targeting the underlying accumulated immune factors.
T cells are crucial to the pathogenesis of GD in that they not only infiltrate the thyroid with auto-reactive capacity but also stimulate B cell differentiation toward antibody-secreting plasma cells (6, 7). Once exposed to a specific antigen, naive CD4+ T cells differentiate into different subtypes, including pro-inflammatory T helper 17 (Th17) cells and immunosuppressive regulatory T cells (Tregs) (8). In our previous studies, the presence of decreased Treg proportions and increased Th17/Treg ratios in both animal models and GD patients were confirmed (9, 10), which is consistent with prior achievements (11, 12).
Besides GD, T cell subtypes may also participate in pGD. The -31C/T polymorphism in the IL1B gene producing high quantities of IL-1β (a Th17 cell inducer) together with an elevated peripheral percentage of Th17 cells are found (13, 14) in pGD. Additionally, the -3279AA genotype of the Forkhead box protein P3 (Foxp3) gene, which leads to defective transcription of FOXP3 via the loss of binding to certain transcription factors (E47 and C-Myb), is present in 11.3% of pGD patients but absent in remittent patients (15). It has been reported that FOXP3 controls Treg function through co-operation with nuclear factor of activated T cells (NFAT), and deficiency of the FOXP3 gene impairs the suppressor function of Treg cells, which often exacerbates the progression of autoimmune disease (16). Taken together, the results of the above studies indicate that Th17 and Treg cells potentially contribute to the long-lasting phase of GD.
Not only T cell subtypes, but clonal expansion of T cells is also involved in autoimmune diseases (17). T cell receptor (TCR), the diversity of which is a prerequisite for the adaptive immune system to respond to a wide variety of antigens, and the structure of the TCR β gene is a useful tool for identifying clonal expansions of T cells (17, 18). With respect to GD, restricted usage of TCR Vα and TCR Vβ was observed in T cells infiltrating the thyroid together with extrathyroidal sites (19), which emphasizes the crucial role of the TCR repertoire in the pathogenesis of GD. Apart from GD, TCR repertoire is also associated with persistence in other autoimmune diseases. In a study of autologous hematopoietic stem cell transplantation (HSCT) for poor-prognosis multiple sclerosis, non-responders had less TCR diversity than responders early in the course of HSCT (20). In addition, analysis of the CD8+ TCR repertoire in alemtuzumab-treated patients with relapsing-remitting multiple sclerosis (RRMS) who develop a secondary autoimmune disease (vitiligo) revealed greatly increased clonality and reduced repertoire diversity in those patients compared to treatment-naive patients with RRMS (21). Thus, we speculate that TCR repertoire may be involved in the long-lasting phase of GD.
In this study, we investigated the T cell immune mechanisms of the long-lasting phase in GD in terms of the T cell subset frequency (Th17 and Treg cells) and TCR repertoire. The results suggested that dysregulation of Tregs and pathogenic T cell clonal expansion may be involved in the persistence of GD by upregulating chemotaxis or inflammation response. For clinical treatment of pGD patients, precise therapies aimed at improving Treg cell frequencies or targeting specific expanded pathogenic TCR clones are promising approaches for future treatments.
Materials and Methods
Study Subjects
We enrolled a total of 20 pGD patients who had taken ATDs for more than 5 years with positive anti-thyroid stimulating hormone receptor autoantibodies (TRAb) value and 34 newly diagnosed GD (nGD) patients who had taken ATDs less than 3 months. These patients were diagnosed according to the guidelines of the American Thyroid Association (ATA). All participants were recruited from the Department of Endocrinology, First Affiliated Hospital of Xi’an Jiaotong University. This study was approved by the Ethic Committee of the First Affiliated Hospital of Xi’an Jiaotong University. We obtained consent from each patient after explaining the purpose of our study. Among these patients, specific clinical characteristics of 20 pGD and 20 nGD patients are presented in Table 1. Flow cytometry was performed in 12 pGD and 14 nGD patients. Samples from 13 pGD and 20 nGD patients were subjected to RNA sequencing (RNA-seq) and T cell receptor sequencing (TCR-seq) analysis. The clinical trial registration number is ChiCTR-IPR-16009305.
Table 1 Clinical characteristic of enrolled pGD and nGD patients.
Sample Preparation
Nine milliliters of peripheral blood was drawn from each patient and collected in an EDTA anticoagulant tube (BD Biosciences, USA). Next, peripheral blood mononuclear cells (PBMCs) were acquired utilizing Ficoll Paque Plus (GE Healthcare, USA) gradient cell separation according to the manufacturer’s instructions. For RNA-seq and TCR-seq analysis, PBMCs were placed in TRIzol (Ambion, USA), and mRNA was extracted according to the manufacturer’s instructions. Then, the quality and integrity of the total RNA were assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, USA) and 1.2% agarose gel electrophoresis; the RNA concentration was measured using a NanoDrop 2000 spectrophotometer (NanoDrop Technologies, USA).
Flow Cytometry
The concentration of isolated PBMCs was adjusted to 107/ml, and 100 μl was used to conduct flow cytometry. Before Th17 cell staining, PBMCs were stimulated with Cell Stimulation Cocktail (Invitrogen, USA) for 6 h at 37°C and 5% CO2. Surface staining including CD3, CD8, CD25, and CD4 (all from BD Biosciences, USA) was performed at 4°C for 30 min. After permeabilization with the Transcription Factor Buffer Set (BD Biosciences, USA) according to the manufacturer’s guidelines, anti-Foxp3 and IL17A antibodies (eBiosciences, USA) were incubated at 4°C for 45 min. Four-color flow cytometric analysis was performed using a FACSCalibur flow cytometer (BD Biosciences, USA).
TCR Sequencing and Analysis
In this study, HTBI primers and Arm-PCR from iRepertoire were used to construct the libraries as previously reported (22). After gel purification, the PCR product was subjected to high-throughput sequencing using the Illumina HiSeq 2000 platform. We amplified all CDR3 sequences present in 33 samples obtained from pGD and nGD patients and then filtered the sequence reads. After the high-quality paired reads were merged using COPE and FqMerger (BGI, China), the results designated as contigs were subsequently aligned to reference TCR Vβ/Dβ/Jβ gene sequences (http://www.imgt.org/download/GENE-DB/) using BLAST. The TCRβ CDR3 regions were identified within the sequencing reads according to the definition established by the International ImMunoGeneTics (IMGT) collaboration. In the end, we identified 62 Vβ and 14 Jβ segments. These results were used to analyze the Vβ and Jβ usage of the CDR3 amino acid clonotypes in each sample. Clones with a frequency of more than 0.1% were considered highly expanded clones (HECs). All TCR-seq raw sequencing data are available through the NCBI SRA accession PRJNA634746.
RNA Sequencing and Analysis
mRNA was enriched and cleaved into 300-bp fragments using fragmentation buffer. Then, the fragments were used as templates for double-stranded cDNA (dscDNA) synthesis using the SuperScript dscDNA synthesis kit (Invitrogen, USA) according to the manufacturer’s instructions. The dscDNA with indexed adapters was purified and enriched by PCR for 15 cycles and then assessed using a TBS-380 fluorometer, a NanoDrop 2000 spectrophotometer, and an Agilent 2100 Bioanalyzer. The prepared libraries were sequenced using an Illumina HiSeq 3000 sequencer.
For all raw reads, the low-quality bases and adapter sequences were removed, and then the method based on fragments per kilobase of exon per million fragments mapped (FPKM) was used to calculate gene expression. To identify differentially expressed genes (DEGs) between the libraries, we used the DESeq2 package (23). Genes with a fold change ≥2 (|log2FC|>1) and a false discovery rate (FDR) <0.05 were considered significant. Two features were extracted from all genes of each group with an unsupervised principal component analysis (PCA) method. Heatmaps and volcano plots were plotted using the pheatmap and ggplot2 packages. We conducted functional enrichment analyses including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways with the “clusterProfiler” R package (24) and set P value <0.05 as the threshold. Biological processes (BP) and molecular function (MF) were included in the GO enrichment analysis. The R package “GOplot” (25) was used to graphically visualize the results of the GO analyses. The protein-protein interaction (PPI) network was constructed based on data obtained from the SPRING online database (https://string-db.org/) (a confidence score >0.4 was set as the cut-off value) and was displayed using Cytoscape version 3.8.0 (Degree ≥13 regarded as key genes). Pathways with a P-value ≤0.05 were considered significantly enriched. The R package “clusterProfiler” was used to perform GSEA. The gene sets in the Molecular Signatures Database (MSigDB) (26) were selected as the reference gene sets, and P value <0.01 was used as the threshold. All RNA-seq raw sequencing data are available through the NCBI SRA accession PRJNA634746.
Statistical Analysis
All of the data are presented as the mean, median, or median (range), depending on whether they followed a normal distribution. Statistical analysis was performed using R 3.6.2, SPSS (version 23, IBM, USA) and Prism (version 7, GraphPad, USA). Categorical variables were analyzed by the χ2 test. For continuous variables, Student’s t-test (normal distribution determined by Kolmogorov-Smirnov Test) or the Mann-Whitney U-test (non-normal distribution determined by Kolmogorov-Smirnov Test) was adopted. Two-sided P values <0.05 were considered statistically significant. The flow cytometry results were analyzed using FlowJo (version 10, Tree Star, USA). Power analysis was conducted using R 3.6.2 with pwr package 1.3-0 (27) and the results were exhibited in Table S1.
Results
Clinical Characteristics of pGD and nGD Patients
The clinical characteristics of the patients are shown in Table 1. There were 20 pGD patients (age: 38.30 ± 9.510; gender: 2 males and 18 females) and 20 nGD patients (age: 40.40 ± 10.758; gender: 4 males and 16 females). The average time for taking ATD among the pGD patients was 10.29 years. No significant differences were observed in smoking (including passive), thyroid nodules, and thyroid disease family history between the pGD and nGD patients. After more than 5 years of ATDs use, the FT4 and total thyroxine (TT4) levels significantly declined in the pGD patients compared with those observed in the nGD patients (FT4: 21.00 pmol/L vs 34.60 pmol/L, P = 0.036; TT4: 12.35 pmol/L vs 20.64 pmol/L, P = 0.003), while there were no differences in other thyroid-associated parameters, including TSH, anti-thyroglobulin autoantibodies (TGAb), anti-thyroid microsomal autoantibodies (TMAb), and anti-thyroid stimulating hormone receptor autoantibodies (TRAb). Regarding immune-related cells, while monocyte (MONO) and PBMC counts decreased in the pGD group (MONO: 0.43 × 109/L vs 0.31 × 109/L, P = 0.024; PBMC: 2.45 × 109/L vs 1.86 × 109/L, P = 0.011), no obvious differences were observed in the quantities of platelets, white blood cells, lymphocytes, basophils, or eosinophils between the groups. Additionally, there was no significant difference in cholesterol concentration between the two groups, which is consistent with our previous finding that high pretreatment cholesterol levels were not correlated with the efficiency of intravenous methylprednisolone in thyroid-associated ophthalmopathy (28).
Flow Cytometric Analysis Identified Decreased Frequencies of Th17 and Treg Cells in pGD Patients
Considering that Th17 and Treg potentially contribute to the persistence of GD, we conducted flow cytometry of Th17 cells (CD3+CD8−IL17A+, Figure S1A) and Treg cells (CD4+CD25+FOXP3+, Figure S1B) in 12 pGD and 14 nGD patients. Compared with the nGD group, the proportion of Th17 cells among CD4+ T cells was significantly decreased in the pGD group (0.4595 vs 0.9195%, P = 0.0307, Figures 1A, B). Moreover, the percentage of Tregs was significantly diminished in the pGD group (3.163 vs 4.466%, P = 0.0306, Figures 1C, D).
Figure 1 Flow cytometric analysis confirmed declined frequencies of Th17 and Treg cells in pGD patients. (A, C) Representative flow cytometry plots of (A) Th17 cells (IL17A+) gated in CD3+CD8− T cells and (C) Treg cells (CD25+Foxp3+) gated in CD4+ T cells from pGD (left) and nGD (right) patients respectively. The numbers denote the percentage of cells in each rectangle gate. (B, D) Whisker plots showed percentage of (B) Th17 cells and (D) Treg cells in CD4+ T cells in pGD (triangle, N = 12) and nGD (circle, N = 14) patients. Data were exhibited as mean with SD with individual values. pGD, persistent Graves’ disease; nGD, newly diagnosed Graves’ disease; SD, standard deviation.
Restricted TCR Clonotypes Were Revealed in pGD Compared With nGD Patients
Given the potential role of TCR repertoire in the long-lasting phase of GD, we decided to perform TCR-seq analysis of PBMCs from pGD and nGD patients (N = 13 and 20, respectively). Compared with the nGD patients, the pGD patients showed decreased Shannon entropy (an index to evaluate TCR diversity), which suggested the persistence of expanded clones in the TCR repertoire of pGD (P = 0.0025, Figure 2A). There were no differences in the V-J evenness (Gini coefficient) and frequency of HECs between the two groups. However, the diminished diversity was confirmed by the reduced quantity of HECs in pGD compared to nGD patients (P = 0.0006, Figure 2A). Lower diversity of the TCR repertoire may reflect more recognition of certain specific antigen and less responses to diversified antigens. Such restricted recognition implied T cell clones targeting these antigens undergo expansion, which is also called clonal expansion (29). Subsequently, we assessed the gene usage in these two distinct groups. Analysis of the frequency of TRBV and TRBJ segment suggested that TRBV28, TRBJ1-2, and TRBJ2-5 (P = 0.0021, 0.00098, and 0.022, respectively) were preferentially used in pGD patients, while TRBV20-1 and TRBJ2-1 (P = 0.00026 and <0.0001, respectively) were preferred in nGD patients (Figures 2B, C; Figure S2). In the pGD samples, the top V-J recombination was often observed as a large percentage, such as in pGD patients R11 (43.4%), R7 (21.3%), and R20 (15.6%) (Figure 2D). The clonal expansion in pGD patients potentially resulted from persistent exposure to the pathogenic auto-antigen. TCR clones shared by more than one individual are generally considered public. Compared with the nGD patients (N = 966), the sharing of public clones was increased among the pGD patients (N = 2,741), suggesting a selection advantage of public clones (Figure 2E). Therefore, compared to the nGD patients, a more skewed TCR repertoire, represented by less TCR diversity, a larger percentage of top V-J combinations, and more public clones, was observed in pGD patients.
Figure 2 Restricted TCR clonetype were revealed in pGD compared with nGD patients. (A) Whisker plots showed Shannon entropy, Gini, high expanded clone number, as well as high expanded clone frequency in pGD patients (red, N = 13) compared with nGD patients (blue, N = 20). Data were exhibited as mean with SD with individual values. (B, C) TCR (B) V segment and (C) J segment frequencies of reads in pGD patients (N = 13) compared to nGD patients (N = 20). Columns represent individual samples. The colors in the bars within each column indicate different V and J segments, and were denoted on the right. (D) TRBV and TRBJ gene segment usage and V-J recombination were illustrated by circos plots in sample pGD patients (upper) represented by P7, P20, P11 and nGD patients (bottom) represented by N5, N51, N31. The TRBV and TRBJ genes were clockwise arranged in the order of their frequency from low to high. A VJ recombination was illustrated by colored curved paths whose thickness represents their frequencies in TCR repertoires. (E) Flower plots presented the amount of the overlap (core) and individual (petal) in TCR clonotypes in pGD patients (upper, N = 13) and nGD patients (bottom, N = 20), respectively. Freq, frequency; pGD, persistent Graves’ disease; nGD, newly diagnosed Graves’ disease; SD, standard deviation; TCR, T cell receptor.
Transcriptional Profiling of PBMCs in pGD and nGD Patients
To further explore the potential function of decreased Treg cell frequency and marked TCR clonal expansion in pGD, we performed transcriptional profiling of PBMCs from 13 pGD and 20 nGD patients. The principal component analysis (PCA) score trajectory plot of the persistent group did not overlap with that of the newly diagnosed group, indicating that there were significant differences in the transcriptional profiles between the two groups (Figure 3A). Combined with the hierarchical cluster analysis of the global gene expression profiles, these results revealed the existence of discrimination between pGD and nGD patients (Figure S3).
Figure 3 Transcriptional profiling of PBMCs in pGD and nGD patients. (A) PCA plot showed obvious differences with all expressed genes in pGD and nGD group. Red triangle represented pGD and blue circle for nGD. (B) Volcano plot showed DEGs between pGD and nGD group. Blue dots represented downregulated genes (N = 375) and red dots were upregulated genes (N = 696). The texts denoted in the plot were names for corresponding genes. (C) The top 10 most significant KEGG pathways for all DEGs. Size of bubble represented counts of the DEGs in corresponding pathway and color of bubble meant adjust p value. (D, E) GO enrichment analysis of DEGs among (D) biological processes and (E) molecular functions. Circos plots showed the relationship between genes and GO terms. The top 10 most significant GO processes were shown in the plot retrospectively. (F) PPI network of potential key DEGs (degree ≧13). Size and color of each gene represented the its degree score. pGD, persistent Graves’ disease; nGD, newly diagnosed Graves’ disease; PCA, principal component analysis; DEGs, differentially expressed genes; KEGG, Kyoto Encyclopedia of Genes and Genomes; GO, Gene ontology; PPI, Protein-protein interaction.
Compared to the nGD patients, 1,071 differentially expressed genes (DEGs) were identified in the persistent group, including 375 downregulated genes and 696 upregulated genes (Figure 3B). For the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, cytokine-cytokine receptor interaction was most often enriched by DEGs (Figure 3C). Gene Ontology (GO) analysis of the DEGs revealed that regulation of adaptive immune response, immune activities associated with chemokines and pro-inflammatory IL-1β were enriched among biological processes (BP) (Figure 3D). Likewise, cytokine and chemokine related activities were the most significant among molecular functions (MF) in the GO analysis (Figure 3E). A protein-protein interaction (PPI) network of DEGs was constructed, and 103 genes were identified as potential hub genes (Degree ≧13), where those genes with highest degree were inflammation and chemotaxis related including IL6, TNF, CXCL8, and IL-1β (Figure 3F). Taken together, these results suggested that inflammation and chemotaxis are correlated with the long-lasting phase of GD.
Specific Gene Expression in pGD Patients
Subsequently, we performed a detailed analysis of specific gene expression in the persistent and newly diagnosed group. Compared with the nGD patients, pro-inflammatory genes, including IL1B and IL13, were upregulated in the persistent group (IL1B: P < 0.0001; IL13: P < 0.0001, Figure 4A), which was consistent with previous reports (12, 30). In addition, upregulation of chemokines (IL8 and CCL4) was observed in pGD patients (IL8: P < 0.0001; CCL4: P < 0.0001, Figure 4B). These results were in agreement with a previous report that IL8 and CCL4 are associated with relapse in autoimmune polychondritis (31).
Figure 4 Specific genes expression in pGD compared with nGD patients. (A–D) Whisker plots indicated normalized gene expression associated with (A) inflammation (B) chemotaxis (C) Th17 cells (D) Treg cells between pGD (red, N = 13) and nGD (blue, N = 20) samples. Data were exhibited as mean with SD with individual values. (E) The hierarchical clustering of 13 pGD and 20 nGD patients by differentially expressed immune related genes. The heatmap showed, for each gene, the scaled expression in each patient. pGD, persistent Graves’ disease; nGD, newly diagnosed Graves’ disease; SD, standard deviation; FPKM, fragments per kilobase of exon per million fragments mapped.
Apart from the increased expression of pro-inflammatory and chemokine genes, consistent with our flow cytometric results, the expression of Th17- and Treg-related genes was decreased in pGD compared with that observed nGD patients. Genes associated with Th17 cells, such as RORC (P = 0.0189), CARD9 (P < 0.0001), IL17RA (P = 0.0037), and IL17RC (P < 0.0001), were significantly downregulated in the persistent group (Figure 4C). In addition, although the master regulator gene FOXP3 in pGD patients presented no difference compared to the newly diagnosed group, the expression levels of genes that are crucial in Treg differentiation and maintenance, including SATB1 (P < 0.0001), STAT5A (P < 0.0001), and TNFRSF18 (P = 0.0263), were obviously decreased (Figure 4D). In addition to the specific genes observed above, samples from pGD patients expressed higher levels of chemotaxis and inflammation molecules, including TNF, IL-6, CXCL2, and CXCL3 (Figure 4E).
Robust Inflammatory and Cytokine Responses Pathways Are Prominent in pGD Patients With Lower FOXP3 and IKZF2 Expression
To further validate the crucial role of Treg cells in the pathogenesis of pGD and elucidate potentially associated functional pathways, we performed GSEA analyses to map the biological processes (Figure 5). FOXP3 and IKZF2 were previously reported as master genes for the differential and suppressive function of Treg cells and its subset naturally occurring Treg cells (nTregs) (32). Therefore, we divided the pGD patients into two groups based on the expression of these two genes to perform subsequent analyses. In pGD patients with lower FOXP3 or IKZF2 expression, inflammatory response, complement, TNFA signaling via NFKB in the hallmark gene set together with the cytokine-cytokine receptor interaction pathway in KEGG pathways were significantly enriched (Figures 5A–D). However, there was no obvious enrichment in pGD patients with lower RORC or CARD9 gene expression (key genes for Th17 cells) compared with those with higher expression. Taken together, the GSEA results suggested that Treg cells contribute to the long-lasting phase of GD by aggravating the inflammatory response.
Figure 5 Robust inflammatory and cytokine responses pathways were prominent in pGD patients with lower FOXP3 and IKZF2 expression. (A, C) Gene set enrichment analysis (GSEA) showed complement, inflammatory response, KRAS signaling up pathways in hallmark gene set were most enriched in pGD patients with (A) lower FOXP3 and (C) lower IKZF2 expression. (B, D) GSEA showed cytokine-cytokine receptor interaction pathway in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways was obviously enriched in pGD patients with (B) lower FOXP3 and (D) lower IKZF2 expression. pGD, persistent Graves’ disease.
Lower TCR Diversity Is Correlated With Inflammation and Chemotaxis in pGD Patients
To further investigate the potential role of the marked expansion of T cell clones in the long-lasting phase of GD, we conducted linear regression analysis between TCR diversity (as reflected by entropy) and specific gene expression among the pGD patients. The results revealed that TCR diversity exhibited a negative correlation with inflammation-associated genes, including FASLG (P = 0.0006), IL18R1 (P = 0.0155), IL5RA (P = 0.0203), and FGF2 (P = 0.0440) (Figure 6A). In addition, chemotaxis-related genes, including CCL24, CCL14, CCL2, and CCR8, were significantly negatively correlated with TCR diversity (P = 0.0003, 0.0061, 0.0060, and 0.0447, respectively; Figure 6B). In summary, lower TCR diversity was obviously associated with more robust inflammation and chemotaxis among pGD patients.
Figure 6 Lower TCR diversity was correlated with inflammation and chemotaxis in pGD. (A, B) The correlations of TCR Entropy with the expression of (A) inflammation (B) chemotaxis associated genes in pGD patients (N = 13) by linear regression analysis using Pearson method. The gray zone represented 95% confidence intervals. Texts denoted in the left bottom corner were the correlation coefficient (cor) and P value. pGD, persistent Graves’ disease; TCR, T cell receptor.
Discussion
In the present study we investigated the underlying perpetuating factors that maintain the autoimmune response in GD through flow cytometric, TCR-seq, and RNA-seq analyses, and our findings emphasized the essential role of T cells. To our knowledge, this is the first research that exploring the mechanism of pGD by TCR-seq and RNA-seq. We identified that decreased Treg frequency and a marked expansion of TCR clonotypes were correlated with the long-lasting phase of GD via potential aggravation of inflammation and chemotaxis, which shed new light on the treatment for pGD patients.
RNA-seq analysis of PBMCs revealed robust inflammation and chemotaxis in the pGD group compared with that observed in the nGD group. In agreement with our results, the ratio of intracellular TLR7 (iTLR7) and iTLR9 intensities in B cells, which can indicate an excessive release of pro-inflammatory cytokines, were enhanced in pGD patients (33). In addition, serum-derived exosomes from pGD patients stimulated the mRNA expression of pro-inflammatory IL-1β and tumor necrosis factor α (TNF-α) in PBMCs (34). In another study, the MIG rs2276886 AG genotype, which regulates the migration of Th17 cells by altering the affinity for C-X-C motif chemokine receptor 3 (CXCR3), was reduced in pGD patients (35). Mechanistically, robust inflammation in the thyroid can be boosted by immune cell infiltration resulting from overexpressed chemotactic molecules. Subsequently, such inflammation will aggravate the autoimmune reaction by amplifying apoptotic signals, cell activation, and tissue remodeling, and may result in the long-lasting phase of GD (36). Taken together, these findings suggest that intensive inflammation and chemotaxis are involved in the progressive evolution of GD.
Given that the primary effector molecular of Th17 cells, IL-17, is associated with inflammation and chemoattraction, we speculate that the robust inflammation and chemotaxis in pGD was mediated by Th17 cells (37). However, the decreases in Th17 cell frequency and related gene expression together with no obvious enrichment in the GSEA results of pGD patients with lower RORC or CARD9 refute our hypothesis. Consistent with our findings, the frequency of the suppressor of cytokine signaling 3 (SOCS3) rs4969170 AA genotype, which negatively regulates Th17 generation, was elevated in patients with pGD (38). Previous studies have shown an increase in Th17 cell-related factors in pGD patients compared with remittent patients (12, 13, 39). A possible explanation for this observation is that nGD patients, possessing higher Th17 cells frequencies and levels of the related cytokines CXCL10 and IL23 (40, 41), rather than those remittent ones, were used as controls in the present study, which further resulted in a decreased Th17 ratio after comparison. In that way, the divergence between our results and previous findings can be partially explained.
Furthermore, we identified that descending Treg-related genes and frequencies dominated in pGD patients through RNA-seq and flow cytometric analysis, which further deteriorate the inflammatory and chemotactic state as suggested by both the GSEA results and previous studies (42). Although achievements in Treg alteration for pGD patients are insufficient, studies concerning immune suppression regulation in pGD agree with our findings. The MIR125A rs12976445 CC genotype, which is correlated with reduced expression of miR-125a (a negative regulator of activated T cells), together with hsa-miR-98-5p (a negative regulator of immunosuppressive IL-10), were increased in pGD patients compared to that observed in remitted patients (43, 44). Interestingly, under inflammatory environments, loss of Foxp3 expression in Treg cells has been observed, which suggests that impaired Treg cells in pGD can be intrinsic and/or extrinsic (45). Combined with the potential role of inflammation and chemotaxis in the pathogenesis of pGD, the dysregulation of Treg cells is closely associated with the long-lasting phase of GD.
Another assumption of the origin of the excessive inflammatory and chemotactic response is the expansion of unique pathogenic T cell clones in pGD. Indeed, TCR-seq analysis revealed marked expansion of T cell clones and increased public clones in pGD relative to nGD patients. These public T cell clones detected in pGD patients are potentially TSHR-specific and pathogenic similar to experimental autoimmune encephalomyelitis (EAE)-specific T cells with a shared Vβ8.2 segment (46). Mechanistically, antigen-specific T cells can initiate a cascade of immune cell immigration into the antigen-loaded site and further induce inflammation mediated by IFN-γ and TNF-α (47). It is reasonable to infer that the persistent expansion of antigen specific T cell clones in pGD patients will result in accumulated inflammation, chemotaxis, and a refractory autoimmune condition. Such inference can be partially demonstrated by our findings in that the diversity index was positively correlated with the expression of specific pro-inflammatory and chemotactic genes. Collectively, marked clonal expansion of pathogenic T cell clones leads to the long-lasting phase of GD by eliciting more intensive inflammatory and chemotactic responses.
The current study has some limitations. First, the quantity of samples was insufficient, which may have caused our results to be inaccurate to some extent. It is expected that a larger number of samples could be enrolled to perform a more detailed analysis. Second, considering that thyroid is the primary lesion while inaccessible, PBMCs utilized in this research can only partially disclose the root cause of pGD. Third, given the inherent defects of RNA-seq analysis, some alterations that are prominent in certain cell subtypes and frequently repeated in other types may be ignored. Methods such as single-cell RNA-seq could eliminate such weaknesses, and the application of such methods to explore the long-lasting phase of GD is a promising approach.
In summary, marked expansion of pathogenic T cell clones in concert with dysregulated Treg cells contribute to the long-lasting phase of GD by presumably eliciting excessive inflammation and chemotaxis responses. For patients with pGD, precise therapies aimed at improving Treg cell frequencies or targeting specific expanded pathogenic TCR clones are promising approaches for future treatments.
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 and PRJNA634746 https://www.ncbi.nlm.nih.gov/bioproject/PRJNA634746.
Ethics Statement
The studies involving human participants were reviewed and approved by Ethic Committee of the First Affiliated Hospital of Xi’an Jiaotong University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
ZC, YW, and BS conceived and designed the study. BS, SH, and ZC screened participants for study entry. ZC, YL, and YW processed the data. SH, MZ, and ZC collected samples and data. ZC and SH performed FACS analysis. ZC, YW, YL, and BS wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Key R&D Program of China [grant NO. 2018YFC1311500 (BS)], National Science Foundation of China (NSFC) [grant NO. 81970679 (BS), 81500690 (YW)], Natural Science Foundation of Shaanxi Province [2018JM70990 (YW)], the Fundamental Research Funds for the Central Universities [1191329875 (YW)], China Postdoctoral Science Foundation [224646 (YW)], and Clinical research award of the first affiliated hospital of Xi’an Jiaotong University (XJTU1AF-CRF-2019-020).
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2021.632492/full#supplementary-material
Supplementary Figure 1 | Flow cytometry analysis of Th17 and Treg cells. (A, B) Gating of (A) Th17 cells and (B) Treg cells.
Supplementary Figure 2 | Global gene expression profiles from pGD and nGD samples. Pairwise Pearson correlations of expression of genes with top 500 mad (median absolute deviation) value between individual patient samples. Samples were organized by hierarchical clustering based on their Pearson coefficients across samples. The disease is indicated by red (pGD, N = 13) and blue (nGD, N = 20) squares on the top and left of the matrix. The patient number for each sample is indicated to the bottom of the matrix. pGD, persistent Graves’ disease; nGD, newly diagnosed Graves’ disease.
Supplementary Figure 3 | Specific TRBV and TRBJ segment usage in pGD and nGD samples. (A, B) Boxplots showed the frequency of (A) TRBV28, TRBJ1-2, TRBJ2-5, (B) TRBV20-1, TRBJ2-1 in all reads between pGD patients (red, N = 13) and nGD patients (blue, N = 20). Boxplot with plot center and box corresponding to median and extremum respectively and include individual data points. pGD, persistent Graves’ disease; nGD, newly diagnosed Graves’ disease.
References
1. Smith TJ, Hegedüs L. Graves’ Disease. N Engl J Med (2016) 375:1552–65. doi: 10.1056/NEJMra1510030
2. Kim YA, Cho SW, Choi HS, Moon S, Moon JH, Kim KW, et al. The Second Antithyroid Drug Treatment Is Effective in Relapsed Graves’ Disease Patients: A Median 11-Year Follow-Up Study. Thyroid (2017) 27:491–96. doi: 10.1089/thy.2016.0056
3. Yamamoto K, Itoh M, Okamura T, Kimura M, Yokoyama A, Yoshino Y, et al. Relative levels of the inflammatory cytokine TNFa and the soluble CD40 ligand profile in serum correlate with the thyrotoxic activity of Graves’ disease. Thyroid (2012) 22:516–21. doi: 10.1089/thy.2011.0222
4. Lang BH, Woo YC, Chiu KW. Two-year outcomes of single-session high-intensity focused ultrasound (HIFU) treatment in persistent or relapsed Graves’ disease. Eur Radiol (2019) 29:6690–98. doi: 10.1007/s00330-019-06303-8
5. Ruiz-Riol M, Barnils MP, Colobran OR, Pla AS, Borràs SF, Lucas-Martin A, et al. Analysis of the cumulative changes in Graves’ disease thyroid glands points to IFN signature, plasmacytoid DCs and alternatively activated macrophages as chronicity determining factors. J Autoimmun (2011) 36:189–200. doi: 10.1016/j.jaut.2011.01.002
6. Brown RS, Lombardi A, Hasham A, Greenberg DA, Gordon J, Concepcion E, et al. Genetic analysis in young-age-of-onset Graves’ disease reveals new susceptibility loci. J Clin Endocrinol Metab (2014) 99:E1387–91. doi: 10.1210/jc.2013-4358
7. Hasham A, Tomer Y. Genetic and epigenetic mechanisms in thyroid autoimmunity. Immunol Res (2012) 54:204–13. doi: 10.1007/s12026-012-8302-x
8. Borst J, Ahrends T, Bąbała N, Melief C, Kastenmüller W. CD4(+) T cell help in cancer immunology and immunotherapy. Nat Rev Immunol (2018) 18:635–47. doi: 10.1038/s41577-018-0044-0
9. Wang Y, Wu LP, Fu J, Lv HJ, Guan XY, Xu L, et al. Hyperthyroid monkeys: a nonhuman primate model of experimental Graves’ disease. J Endocrinol (2013) 219:183–93. doi: 10.1530/JOE-13-0279
10. Chen Z, Wang Y, Ding X, Zhang M, He M, Zhao Y, et al. The proportion of peripheral blood Tregs among the CD4+ T cells of autoimmune thyroid disease patients: a meta-analysis. Endocr J (2020) 67:317–26. doi: 10.1507/endocrj.EJ19-0307
11. McLachlan SM, Rapoport B. Breaking tolerance to thyroid antigens: changing concepts in thyroid autoimmunity. Endocr Rev (2014) 35:59–105. doi: 10.1210/er.2013-1055
12. Yuan Q, Zhao Y, Zhu X, Liu X. Low regulatory T cell and high IL-17 mRNA expression in a mouse Graves’ disease model. J Endocrinol Invest (2017) 40:397–407. doi: 10.1007/s40618-016-0575-9
13. Hayashi F, Watanabe M, Nanba T, Inoue N, Akamizu T, Iwatani Y. Association of the -31C/T functional polymorphism in the interleukin-1beta gene with the intractability of Graves’ disease and the proportion of T helper type 17 cells. Clin Exp Immunol (2009) 158:281–86. doi: 10.1111/j.1365-2249.2009.04034.x
14. Nanba T, Watanabe M, Inoue N, Iwatani Y. Increases of the Th1/Th2 cell ratio in severe Hashimoto’s disease and in the proportion of Th17 cells in intractable Graves’ disease. Thyroid (2009) 19:495–501. doi: 10.1089/thy.2008.0423
15. Inoue N, Watanabe M, Morita M, Tomizawa R, Akamizu T, Tatsumi K, et al. Association of functional polymorphisms related to the transcriptional level of FOXP3 with prognosis of autoimmune thyroid diseases. Clin Exp Immunol (2010) 162:402–06. doi: 10.1111/j.1365-2249.2010.04229.x
16. Lu L, Barbi J, Pan F. The regulation of immune tolerance by FOXP3. Nat Rev Immunol (2017) 17:703–17. doi: 10.1038/nri.2017.75
17. Turner SJ, Doherty PC, McCluskey J, Rossjohn J. Structural determinants of T-cell receptor bias in immunity. Nat Rev Immunol (2006) 6:883–94. doi: 10.1038/nri1977
18. Qi Q, Liu Y, Cheng Y, Glanville J, Zhang D, Lee JY, et al. Diversity and clonal selection in the human T-cell repertoire. Proc Natl Acad Sci U S A (2014) 111:13139–44. doi: 10.1073/pnas.1409155111
19. Heufelder AE, Wenzel BE, Scriba PC. Antigen receptor variable region repertoires expressed by T cells infiltrating thyroid, retroorbital, and pretibial tissue in Graves’ disease. J Clin Endocrinol Metab (1996) 81:3733–39. doi: 10.1210/jc.81.10.3733
20. Muraro PA, Robins H, Malhotra S, Howell M, Phippard D, Desmarais C, et al. T cell repertoire following autologous stem cell transplantation for multiple sclerosis. J Clin Invest (2014) 124:1168–72. doi: 10.1172/JCI71691
21. Ruck T, Pfeuffer S, Schulte-Mecklenbeck A, Gross CC, Lindner M, Metze D, et al. Vitiligo after alemtuzumab treatment: Secondary autoimmunity is not all about B cells. Neurology (2018) 91:e2233–7. doi: 10.1212/WNL.0000000000006648
22. Han Y, Liu X, Wang Y, Wu X, Guan Y, Li H, et al. Identification of characteristic TRB V usage in HBV-associated HCC by using differential expression profiling analysis. OncoImmunology (2015) 4:e1021537. doi: 10.1080/2162402X.2015.1021537
23. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol (2014) 15:550. doi: 10.1186/s13059-014-0550-8
24. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS (2012) 16:284–87. doi: 10.1089/omi.2011.0118
25. Walter W, Sanchez-Cabo F, Ricote M. GOplot: an R package for visually combining expression data with functional analysis. Bioinformatics (2015) 31:2912–14. doi: 10.1093/bioinformatics/btv300
26. Liberzon A, Birger C, Thorvaldsdóttir H, Ghandi M, Mesirov JP, Tamayo P. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst (2015) 1:417–25. doi: 10.1016/j.cels.2015.12.004
27. Parker RI, Hagan-Burke S. Useful effect size interpretations for single case research. Behav Ther (2007) 38:95–105. doi: 10.1016/j.beth.2006.05.002
28. Hu S, Wang Y, He M, Zhang M, Ding X, Shi B. Factors associated with the efficacy of intravenous methylprednisolone in moderate-to-severe and active thyroid-associated ophthalmopathy: A single-centre retrospective study. Clin Endocrinol (Oxf) (2019) 90:175–83. doi: 10.1111/cen.13855
29. Feng Z, Fang Q, Kuang X, Liu X, Chen Y, Ma D, et al. Clonal expansion of bone marrow CD8(+) T cells in acute myeloid leukemia patients at new diagnosis and after chemotherapy. Am J Cancer Res (2020) 10:3973–89. doi: 10.21203/rs.3.rs-50231/v1
30. Inoue N, Watanabe M, Morita M, Tatusmi K, Hidaka Y, Akamizu T, et al. Association of functional polymorphisms in promoter regions of IL5, IL6 and IL13 genes with development and prognosis of autoimmune thyroid diseases. Clin Exp Immunol (2011) 163:318–23. doi: 10.1111/j.1365-2249.2010.04306.x
31. Arnaud L, Mathian A, Haroche J, Gorochov G, Amoura Z. Pathogenesis of relapsing polychondritis: a 2013 update. Autoimmun Rev (2014) 13:90–5. doi: 10.1016/j.autrev.2013.07.005
32. Thornton AM, Korty PE, Tran DQ, Wohlfert EA, Murray PE, Belkaid Y, et al. Expression of Helios, an Ikaros transcription factor family member, differentiates thymic-derived from peripherally induced Foxp3+ T regulatory cells. J Immunol (2010) 184:3433–41. doi: 10.4049/jimmunol.0904028
33. Inoue N, Katsumata Y, Watanabe M, Ishido N, Manabe Y, Watanabe A, et al. Polymorphisms and expression of toll-like receptors in autoimmune thyroid diseases. Autoimmunity (2017) 50:182–91. doi: 10.1080/08916934.2016.1261835
34. Hiratsuka I, Yamada H, Munetsuna E, Hashimoto S, Itoh M. Circulating MicroRNAs in Graves’ Disease in Relation to Clinical Activity. Thyroid (2016) 26:1431–40. doi: 10.1089/thy.2016.0062
35. Akahane M, Watanabe M, Inoue N, Miyahara Y, Arakawa Y, Inoue Y, et al. Association of the polymorphisms of chemokine genes (IL8, RANTES, MIG, IP10, MCP1 and IL16) with the pathogenesis of autoimmune thyroid diseases. Autoimmunity (2016) 49:312–19. doi: 10.3109/08916934.2015.1134507
36. Gianoukakis AG, Khadavi N, Smith TJ. Cytokines, Graves’ disease, and thyroid-associated ophthalmopathy. Thyroid (2008) 18:953–58. doi: 10.1089/thy.2007.0405
37. Miossec P, Kolls JK. Targeting IL-17 and TH17 cells in chronic inflammation. Nat Rev Drug Discovery (2012) 11:763–76. doi: 10.1038/nrd3794
38. Kunisato T, Watanabe M, Inoue N, Okada A, Nanba T, Kobayashi W, et al. Polymorphisms in Th17-related genes and the pathogenesis of autoimmune thyroid disease. Autoimmunity (2018) 51:360–69. doi: 10.1080/08916934.2018.1534963
39. Li JR, Hong FY, Zeng JY, Huang GL. Functional interleukin-17 receptor A are present in the thyroid gland in intractable Graves disease. Cell Immunol (2013) 281:85–90. doi: 10.1016/j.cellimm.2013.02.002
40. Peng D, Xu B, Wang Y, Guo H, Jiang Y. A high frequency of circulating th22 and th17 cells in patients with new onset graves’ disease. PloS One (2013) 8:e68446. doi: 10.1371/journal.pone.0068446
41. He M, Wang Y, Wang J, Sui J, Ding X, Chen Z, et al. The potential markers involved in newly diagnosed graves’ disease and the development of active graves’ orbitopathy. Cytokine (2020) 127:154998. doi: 10.1016/j.cyto.2020.154998
42. Sharabi A, Tsokos MG, Ding Y, Malek TR, Klatzmann D, Tsokos GC. Regulatory T cells in the treatment of disease. Nat Rev Drug Discov (2018) 17:823–44. doi: 10.1038/nrd.2018.148
43. Inoue Y, Watanabe M, Inoue N, Kagawa T, Shibutani S, Otsu H, et al. Associations of single nucleotide polymorphisms in precursor-microRNA (miR)-125a and the expression of mature miR-125a with the development and prognosis of autoimmune thyroid diseases. Clin Exp Immunol (2014) 178:229–35. doi: 10.1111/cei.12410
44. Takuse Y, Watanabe M, Inoue N, Ozaki R, Ohtsu H, Saeki M, et al. Association of IL-10-Regulating MicroRNAs in Peripheral Blood Mononuclear Cells with the Pathogenesis of Autoimmune Thyroid Disease. Immunol Invest (2017) 46:590–602. doi: 10.1080/08820139.2017.1322975
45. Scheinecker C, Göschl L, Bonelli M. Treg cells in health and autoimmune diseases: New insights from single cell analysis. J Autoimmun (2020) 110:102376. doi: 10.1016/j.jaut.2019.102376
46. Fazilleau N, Delarasse C, Sweenie CH, Anderton SM, Fillatreau S, Lemonnier FA, et al. Persistence of autoreactive myelin oligodendrocyte glycoprotein (MOG)-specific T cell repertoires in MOG-expressing mice. Euro J Immunol (2006) 36:533–43. doi: 10.1002/eji.200535021
Keywords: Graves’ disease, persistent Graves’ disease, refractory Graves’ disease, regulatory T cell, T cell receptor sequencing, RNA sequencing
Citation: Chen Z, Liu Y, Hu S, Zhang M, Shi B and Wang Y (2021) Decreased Treg Cell and TCR Expansion Are Involved in Long-Lasting Graves’ Disease. Front. Endocrinol. 12:632492. doi: 10.3389/fendo.2021.632492
Received: 23 November 2020; Accepted: 15 February 2021;
Published: 12 April 2021.
Edited by:
Marian Elizabeth Ludgate, Cardiff University, United KingdomReviewed by:
Angela Lombardi, Albert Einstein College of Medicine, United StatesSyed A. Morshed, Icahn School of Medicine at Mount Sinai, United States
Copyright © 2021 Chen, Liu, Hu, Zhang, Shi and Wang. 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.
*Correspondence: Bingyin Shi, shibingy@126.com; Yue Wang, shelly1021@126.com
†These authors have contributed equally to this work