Abstract

During Graves’ disease (GD) treatment, Graves’ ophthalmopathy (GO) is often ignored because only mild ocular symptoms are present in early GD. Therefore, we performed isobaric tags for relative and absolute quantification (iTRAQ) analysis and measured relevant endocrine hormones to identify predisposing factors of GO. Serum samples from 3 patients with mild GD and GO and 3 patients with GD but without GO were analyzed by iTRAQ. Based on their clinical data, 60 patients with GD were divided into the GO-free and GO groups. All patients were followed up for 7 months. Their eye conditions and changes in related biochemical indexes were recorded. The iTRAQ results showed that RhoA expression was upregulated and correlated significantly with the tight junction pathway and immunity. The changes in FT3 and RhoA from baseline to 7 months, the FT3 and RhoA baseline levels, and the TRAb titer levels in patients with GD significantly differed between the groups. ELISA and western blotting for RhoA, TRAb, and FT3 in the serum samples from GO patients showed significant upregulation, as well as elevated serum RhoA and TRAb levels in the mild stage of GO. At 7 months, the serum RhoA and FT3 levels were elevated. RhoA is a potential biomarker for mild GO. In GD patients, if an elevated serum RhoA level is accompanied by an elevated TRAb or FT3 level, GO is highly likely to occur, even when obvious ocular symptoms are absent.

1. Introduction

Graves’ disease (GD) is an autoimmune thyroid disease and the main cause of hyperthyroidism in adults. Incidence rates are approximately 3% for females and 0.5% for males [1, 2]. Graves’ ophthalmopathy (GO) is the most common extrathyroidal manifestation. Approximately 25%-50% of GD patients have ocular involvement [3, 4]. Most GO patients have hyperthyroidism, though a small number of cases may be hypo- or euthyroid [57]. Approximately 90% of GO patients display symptoms of eyelid retraction, which is an initial symptom and sign of early disease [8, 9]. Approximately 3% to 5% of GO patients will develop severe disease, which involves severe inflammation, eye pain, and vision-threatening corneal ulcers or compressive optic neuropathy (CON) [3, 4, 10, 11]. The appearance of GO is often ignored due to the presence of only mild ocular symptoms in the early stage.

Proteomics is a large-scale, high-throughput, and systematic study of the composition and function of all proteins in a certain type of cell, tissue, or body fluid. High-quality protein separation and identification technology combined with full utilization of bioinformatics reflects changes in the composition and expression level of proteins as a whole, which helps to identify those related to the formation and development of diseases and to explore the process of disease occurrence and development. Biological mechanisms and clinically related research as well as research and development of new drug targets, among others, have been widely employed [1214]. In recent years, high-throughput mass spectrometry has been implemented and used in the study of screening serum proteins, revealing many molecules with biomarker values.

In the past few decades, the incidence of GO has shown an increasing trend. GO severely affects vision, leading to blindness and depression in most patients. Therefore, treatment early in the mild stage of the disease is likely to result in a better prognosis. The purpose of this study was to determine risk factors leading to the onset of GO, which is very important for its early diagnosis. We observed changes in levels of endocrine hormones in GD patients and discovered possible markers for early diagnosis of GO through proteomics, with a view toward early diagnosis of CON and prevention of visual impairment caused by it. We sought to provide new ideas for improving the diagnosis and treatment of GO and for evaluating prognosis.

2. Materials and Methods

2.1. Subjects

We collected clinical data and serum samples for GD patients admitted to the Tianjin Medical University General Hospital from June 2018 to October 2019. We randomly collected venous blood from 3 patients in the mild stage who were first diagnosed with GO complicated with GD (experimental group) and 3 diagnosed with GD without GO (control group) treated at the Tianjin Medical University Eye Hospital from January 2020 to October 2020. The present study was approved by the Tianjin Medical University Eye Hospital Foundation Institutional Review Board (REC No. 2020KY(L)-01) and followed the principles of the Declaration of Helsinki. Written informed consent was obtained from all study subjects. The study participants were not financially compensated.

Diagnostic criteria of GD were as follows [1, 2]: (1) clinical diagnosis of hyperthyroidism, (2) diffuse thyroid enlargement (confirmed by palpation and B-ultrasound), with some cases lacking goiter, (3) exophthalmos and other infiltrating eye signs, (4) anterior tibial mucin edema, and (5) TRAb, TSAb, and TPOAb positivity. Of the above criteria, items (1) and (2) were essential for diagnosis, and items (3), (4), and (5) were auxiliary factors for diagnosis.

GO diagnostic criteria were as follows [5]: (1) Mild GO usually involves one or more of the following symptoms: mild eyelid retraction (<2 mm), mild soft tissue damage,  mm of the same population and sex, no diplopia or intermittent diplopia, and corneal exposure which can be improved by a lubricant, and (2) Moderate to severe GO usually has the following two or more symptoms: eyelid  mm, moderate or severe soft tissue damage, exophthalmos ≥ the same population and sex normal value 3 mm, and persistent or intermittent diplopia.

The inclusion criteria for the experimental group were as follows: (1) initial diagnosis of GD with ophthalmopathy, (2) ophthalmopathy diagnosed as mild GO, (3) oral methimazole prescribed for GD treatment, and (4) no treatment for ophthalmopathy.

The inclusion criteria for the control group were as follows: (1) initial diagnosis of GD without ophthalmopathy and (2) oral methimazole prescribed for GD treatment.

Exclusion criteria were as follows: (1) hyperthyroidism caused by nodular goiter, thyroid cancer, high-functioning adenoma, and thyroiditis, (2) thyroidectomy, (3) primary hypothyroidism, (4) orbital inflammatory pseudotumors, orbital cellulitis, and orbital granuloma, among others, and (5) pregnancy.

The patients were divided into two groups based on changes to the eye condition: a GO-free (36 cases) group and a GO (24 cases) group. All patients were given methimazole for GD. According to EUGOGO [5], all patients in the GO group were evaluated according to the CAS. No patient scored higher than 3. Therefore, none of the patients were considered in the active disease phase. Because these patients were first diagnosed in the Endocrine Department, ocular manifestations of the disease may not have been obvious.

2.2. Isobaric Tags for Relative and Absolute Quantification (iTRAQ) Protein Spectrum Analysis

Peripheral venous blood was collected from patients with GO and GO-free disease and centrifuged as soon as possible to separate the serum, which was immediately stored at -80°C. High-abundance proteins were removed according to the instructions of the ProteoMiner™ kit (Bio-Rad Laboratories, USA). Protein quantification was carried out according to the BCA method, and SDS-PAGE was performed. The samples were reduced, alkylated, and enzymatically hydrolyzed. The isotope mass label iTRAQ reagent was used for labeling and high pH reverse-phase fractionation. Finally, the protein samples were analyzed using a nanoupgraded reverse-phase chromatography-TripleTOF™ 5600.

2.3. ELISA

Using a serum separator tube, the blood samples were allowed to clot for two hours at room temperature or overnight at 4°C before centrifugation for 20 minutes at . Freshly prepared serum was assayed immediately or stored in aliquots at -20°C or -80°C for later use. Experiments were performed according to the manufacturer’s protocol of the FT3 ELISA kit (LEVELOP, Wuxi, China) and TRAb ELISA kit (GTX, Shanghai, China).

2.4. Western Blotting

Whole blood samples were centrifuged at 4°C for 20 minutes. The serum samples were added to serum protein extracts and centrifuged at at 4°C for 10 minutes. Proteins (30 μg) were separated by SDS-PAGE and transferred to polyvinylidene fluoride (PVDF) membranes (Bio-Rad Laboratories, USA). The membranes were washed for 1-2 minutes in 0.1% Tween-20 (TBST) (Solarbio, Beijing, China), blocked (5% skimmed milk powder) at 37°C for 1 hour, and incubated with a primary antibody against RhoA (1 : 5000, Abcam, UK), overnight at 4°C. After rinsing with TBST three times, the membrane was incubated with a rabbit/mouse secondary antibody at room temperature for 1 hour and exposed.

2.5. Biochemical Index Detection

Blood was drawn from fasting patients in the morning. Chemiluminescence immunoassays using Abbott reagents and instruments were used to measure free thyroid function. Chemiluminescence immunoassays using German Siemens reagents and instruments were employed for TPOAb and TgAb detection. Electrochemiluminescence immunoassays using Roche reagents and instruments were used to detect TRAb. The reference value range was 2.63~5.70 pmol/L for FT3, 9.01~19.05 pmol/L for FT4, 0.350~4.940 μU/L for TSH, 0.00~1.75 IU/L for TRAb, 0.00~35.00 IU/mL for TPOAb, and 0.00~40.00 IU/mL for TgAb. Because the measurement range for TPOAb was 10~1000 IU/mL, a  IU/mL was considered to be 1000 IU/mL, similarly, as the range for TgAb was 20~3000 IU/mL, a  IU/mL was considered to be 3000 IU/mL.

2.6. Data Collection

We conducted a 7-month follow-up to record patient serum samples and clinical data from June 2018 to October 2019. Age, sex, treatment of thyroid disease, history of hyperthyroidism, and monthly levels of thyroid-related hormones (FT3, FT4, TSH, TRAb, TPOAb, and TgAb) were collected, related biochemical indexes and the time of occurrence of various eye signs were measured.

2.7. Statistical Analysis

Analyses were performed using SPSS version 23.0 and GraphPad Prism version 8.0 for Windows. All numerical data conforming to a normal distribution are expressed as the . Data not conforming to a normal distribution are presented as medians (interquartile ranges), and enumeration data are given as percentages. The Mann–Whitney test was used for comparison of age. Two-tailed Fisher’s exact tests were applied to compare enumeration data. Mauchly’s test and the Greenhouse-Geisser correction method were utilized to determine significant differences in thyroid hormone levels of continuous data between the groups. Spearman’s rank correlation coefficient was used to determine whether different hormone or antibody titer levels of the two groups correlated with the occurrence of GO. ROC curves were drawn to assess whether hormone or antibody levels predict GO pathogenesis at certain cutoff levels. The area under the ROC curve (AUC) was evaluated for diagnostic value. was considered statistically significant.

3. Results

3.1. Demographic and Ophthalmological Data

The GO and GO-free groups were similar in age at diagnosis and sex (Supplementary Table 1).

3.2. Results of iTRAQ

We randomly collected venous blood from 3 GO with GD patients (experimental group) and 3 GD patients without GO (control group) with matching sex and age and performed iTRAQ protein profile analysis. As Pearson coefficients between the control group and the experimental group samples were greater than 0.86, the data were significant and repeatable.

In the principal component analysis of the protein quantitative results, the control group and the experimental group aggregated obviously into different clusters, indicating that protein expression in the experimental group was significantly different from that in the control group (Figure 1(a)).

We detected 1552 proteins. Compared with the control group, the experimental group showed 229 differential proteins with expression . Among these, 129 proteins were upregulated and 100 downregulated in the experimental group. Among them, the most significantly upregulated proteins were YBX3, RAB3GAP2, APOC3, GCA, and CSRP2, the most significantly downregulated proteins were POLA1, SMAP, ETFB, RNMT, and MAGED4 (Figures 1(b) and 1(c)).

To better investigate the biochemical process of the occurrence and development of GO, we conducted Gene Ontology analysis of differential proteins, mainly from three aspects: cell component, molecular function, and biological process. GO analysis was performed on upregulated proteins. The functions of these targeted proteins in cell components mainly involved focal adhesions and cell substrate joints. Molecular function terms largely were cell adhesion molecule binding and actin binding. In terms of biological process, enrichment in such processes as actin polymerization and depolymerization, cell adhesion, and neutrophil-mediated immunity was observed. For downregulated differential proteins (Figure 1(e)), cell composition terms were mostly the cytoplasm and extracellular matrix. Molecular function mainly involved antigen binding and serine endopeptidase activity, and biological process terms were the protein activation cascade, complement activation, and acute inflammation.

Through KEGG analysis (Figures 1(f) and 1(g)), upregulated differential proteins participate in signaling pathways such as actin cytoskeleton regulation, platelet activation, pathogenic Escherichia coli infection, tight junctions, and transendothelial migration of leukocytes. Downregulated proteins are mainly involved in the complement and coagulation cascade, Staphylococcus aureus infection, and prion diseases. We also performed gene set enrichment analysis () on the quantitative data for differential proteins (Figures 2(a) and 2(b)). The results showed tight junction signaling pathways to be enriched and upregulated in GO, consistent with the KEGG results. The main enriched proteins were also similar to the KEGG results, including YBX3, TJP2, RhoA, RAC1, MYL12A, MYH9, F11R, and CTTN. These proteins are all upregulated in GO.

Pearson correlation analysis revealed a high correlation of expression of proteins involved in the GSEA tight junction pathway (Figure 2(c)). High correlations, such as for RhoA, RAC1, F11R, and TJP2, suggest that these proteins may be similar in function. Furthermore, protein-protein interaction analysis of differential proteins related to tight junctions showed that RhoA is at the center of the interaction (Figure 2(d)). Next, we used western blotting to detect RhoA in serum samples from the two groups (Figure 2(e)) and calculated the grayscale value (Figure 2(g)). ELISA was applied to detect FT3 and TRAb (Figure 2(f)). The results showed that RhoA, FT3, and TRAb levels were significantly higher in the GO group than in the GO-free group.

3.3. Changes in RhoA, FT3, FT4, and TSH Levels

In the GO-free group, 24 patients (66.7%) developed hypothyroidism, with reversion to normal thyroid function levels occurring in 12 (33.3%). In the GO group, 20 patients (83.3%) had hypothyroidism, and 4 experienced reversions to normal thyroid function levels (16.7%). Overall, changes in the thyroid hormone levels of the two groups of patients were recorded for 7 consecutive months (Figures 3(b)3(d)), and serum RhoA levels were assessed by ELISA (Figure 3(a)). We performed repeated-measurement data analysis of variance for FT3, FT4, TSH, and RhoA at different time points. Mauchly’s test of sphericity was performed on factors in each group, however, the data did not satisfy the sphere hypothesis (, Supplementary Table 2). Hence, the values were adjusted using the Greenhouse-Geisser correction method (Supplementary Tables 2 and 3). A correlation between repeated measurements of FT3, FT4, TSH, and RhoA data was detected. After eliminating the interaction of hormones and time, FT3, FT4, and RhoA levels in the GO-free and GO groups were statistically significant, though there were no significant differences in TSH. Nevertheless, monthly FT3, FT4, and RhoA levels were higher in the GO group than in the GO-free group.

3.4. Changes in Basal Levels, Endpoint Values, and Differences in Thyroid-Related Hormone and RhoA Levels

For the two groups, we calculated differences in hormone and RhoA levels between baseline and the follow-up endpoint to evaluate fluctuations (Table 1). Numerical data in accordance with a normal distribution are expressed as the , and the parameter -test was used. Data that did not conform to a normal distribution are presented as the median (range of quartiles), and the Mann–Whitney test was used. The difference between basal levels and endpoint values was statistically significant for RhoA and FT3 levels but not for FT4 or TSH.

3.5. Changes in Thyroid Antibody Titer at the First Diagnosis of GD

Because TRAb levels change slightly during short-term treatment and TPOAb and TgAb are not affected by antithyroid drugs, we compared antibody titers of the two groups of patients at first diagnosis of GD using the Mann–Whitney test. The difference in TRAb titer () was statistically significant (Figure 4(a)), whereas that for TPOAb and TgAb titers () was not (Figures 4(b) and 4(c)).

3.6. Correlation Analysis of Thyroid Hormone and Related Antibodies with GO

Thyroid hormone levels, antibody titers, and RhoA levels were analyzed using Spearman rank correlation analysis, revealing that basal RhoA levels correlated positively with the incidence of GO. In addition, the higher the basal RhoA level was, the more likely the patient was to develop GO. Moreover, the difference between the RhoA endpoint value and basal level correlated negatively with the incidence of GO. As the absolute value was greater when there was a smaller difference between the RhoA endpoint value and the basal level, such patients were likely to develop GO. In addition, basal FT3 levels correlated positively with GO incidence, where the greater the basal FT3 level was, the more prone the patient was to developing GO. The difference between the FT3 endpoint value and basal level correlated negatively, with a smaller difference between the FT3 endpoint value and basal level indicating a greater absolute value and likely GO development. TRAb titers at first diagnosis of GD also correlated positively with the incidence of GO, and patients with a large TRAb value were prone to develop GO. Conversely, FT4, TSH, TPOAb, and TgAb levels did not correlate with GO onset (Table 2).

3.7. Diagnostic Significance of RhoA, FT3, FT4, TSH, and Related Antibodies in Onset of GO

We used specificity as the horizontal axis and sensitivity as the vertical axis to generate a ROC curve to evaluate the diagnostic significance of thyroid hormones, RhoA, and related antibodies. The results demonstrated that basal FT3 levels, basal RhoA levels, and the difference between basal FT3 levels and FT3 endpoint values and the difference between basal RhoA and RhoA endpoint values had diagnostic values (the difference indicates the absolute value) (Figure 5(a)). When patients were diagnosed with GD, TRAb titers also had diagnostic values and were superior to FT3. However, FT4, TSH, TPOAb, and TgAb levels exhibited no diagnostic value (Figure 5(b) and Table 2).

4. Discussion

GO is an autoimmune disease and a common orbital disease that affects vision eye movement and appearance. GO is seriously detrimental to quality of life. The incidence rate of GO is approximately 16 females and 3 males per 100,000 people [11]. Most patients with GD are overlooked due to the mild ocular symptoms, and treatment tends to focus on GD, delaying the diagnosis and treatment of GO. Hence, we aimed to identify biomarkers for early GO diagnosis through iTRAQ analysis. We investigated GD patients with new-onset GO to determine predisposing factors influencing the occurrence of GO. We used iTRAQ and WB and identified that expression of RhoA was significantly increased in GO, as verified in 60 patients and found that expression of RhoA was significantly higher in the GD with GO group than in the GD group at the initial stage (basal , ).

Ras homolog family member A (RhoA) is a low-molecular-weight G protein, and serine/threonine kinase Rho-associated kinase (ROCK) is one of the main downstream effectors of RhoA [15]. Studies have shown that RhoA plays important roles in cell adhesion, migration, and transformation, which is consistent with the results of our KEGG analysis. For example, Wei et al. [16] studied an in vitro model of GO and found that RhoA is involved in these processes. The RhoA/ROCK signaling pathway is also involved in the differentiation of myofibroblasts induced by TGF-β. And in the pathogenesis of GO, orbital fibroblasts are the central target of the immune response and are involved in the process of inflammation and remodeling of orbital connective tissue and extraocular muscles. GO is caused by an autoimmune reaction of shared antigens between the thyroid and the orbit, mainly involving the connective tissue behind the bulb and extraocular muscles. Based on immunohistochemical experiments of orbital tissues of patients with GO [17, 18], Weetman et al. [19] reported infiltration of a large number of lymphocytes, neutrophils, serum cells, and mucopolysaccharide deposits, mainly in fat cells and intermuscular spaces, and that the lymphocytes were T cells. It has also been found that the Th1/Th2 balance participates in the progression of GO, with a shift toward Th1 dominance in the later stage [20]. RhoA is a ubiquitously expressed cytoplasmic protein belonging to the small GTPase family and plays a key regulatory role in innate and adaptive immunity. For example, many scientists have shown that Th0 cells differentiate into one of four main subgroups (Th1, Th2, Th17, and Treg) to perform functions when the initial T lymphocytes entering the periphery are activated by the dual signals of antigen-presenting cells [2123]. Tregs play the most important role in regulating self-reactive T cells and in preventing an excessive immune reaction, which is harmful to the host [21]. RhoA is one of the core proteins for T cell migration [23], and studies on T cell-specific knockout RhoA mice show that RhoA is vital for T cell proliferation, activation, and migration. Therefore, we speculate that RhoA affects the balance of Th1/Th2 cells by interrupting Th0 cell differentiation, which promotes GO pathogenesis.

With the gradual development of GO, severe exophthalmos can lead to incomplete eyelid closure, dry exfoliation of the corneal epithelium, exposure keratitis, corneal ulcers, and even corneal perforation. Hypertrophic extraocular muscles oppress the optic nerve, which can lead to optic neuropathy and even vision loss. ROCK inhibitors have been shown to promote corneal cell proliferation and repair [24, 25]. Okumura et al. [26] found that apoptosis can activate phosphorylation within the Rho/ROCK signaling pathway in corneal endothelial cells, and the loss of cell adhesion can induce apoptosis. ROCK inhibitors can counteract the loss of cell adhesion by activating the focal adhesion complex. Inhibition of the Rho/ROCK pathway can also promote axonal regeneration and functional repair of the central nervous system. Yu et al. [27] established a rabbit model of optic nerve injury and administered the Rho/ROCK inhibitor fasudil. The state of retinal ganglion cells was better in the model group than in the saline group or the dexamethasone group, suggesting that inhibition of the Rho/ROCK signaling pathway can promote optic nerve repair to some extent. Furthermore, we observed increases in TRAb in 60 patients with GD and GO compared with GD patients. GD is characterized by the presence of antibodies against TSHR [2830]. TSH-receptor antibodies (TRAbs) are divided into thyroid-stimulating (TSAbs), thyroid-blocking (TBAbs), and neutral (neutral Abs) antibodies. Studies have demonstrated a close clinical and temporal relationship between GD and GO, which suggests a common pathogenic antigen in the thyroid and orbit [4, 5, 31]. As TRAb is an autoantibody that causes GD [32], TRAb levels may have a certain clinical value in determining the relationship between GD and GO [33]. In our study, the average TRAb titer in the GO-free group at first diagnosis of GD was 4.58 (IU/L), whereas the titer in patients in the GO group was 13.07 (IU/L), indicating diagnostic significance for GO (Figures 2(c) and 3(b)). Our findings are consistent with those of several published studies [34, 35]. It has also been shown that levels of TSAb and TBII may be related to the incidence and severity of GO, with more significance for the former [34, 36]. At present, it is believed that a lack of CD28 on the surface of CD4+ and CD8+ T cell subsets in the peripheral blood is the reason that patients are prone to immune-mediated diseases [37, 38]. The amount of CD28 in the serum of patients with GD is higher than that of normal subjects, and levels of endogenous IFN-γ and IL-6 cytokines in patients with GD with GO are higher than those in GO-free patients, with a positive correlation with TRAb [3741].

TSH is the most sensitive indicator for thyroid function. As T3 and T4 levels are usually affected by binding proteins and albumin, FT3 and FT4 can well reflect thyroid function. During the 7-month follow-up, we found that basal FT3 levels and the difference between basal levels and endpoint values had diagnostic values for GO, based on Mann–Whitney test results (, ; , ). Nonetheless, there were no statistically significant differences between the endpoints of FT4 and TSH and their basal levels. Thus, a true representation of thyroid function should be evaluated using serum FT4 levels instead of serum FT3 levels. However, in our study, FT3 levels in the serum of GO patients showed a positive result that was significantly different from that for the serum of GO-free patients, which may be because FT3 is more sensitive than FT4 during the early stages. We collected patient serum during the early stages of GD, and FT3 might be more easily overactivated than FT4. Hence, patients with GD with significantly elevated FT3 during the early stages were more likely to develop GO, therefore, such an increase in FT3 was easy to detect.

Based on our study, we speculate that RhoA is a potential biomarker for GO and may play an important role in the differentiation and migration of T cells and promote the pathogenesis of GO. We also found that GD patients with a significant increase in TRAb and FT3 in the early stage were more likely to develop GO. TRAb is an autoantibody that is highly expressed in the adipose tissue of GO patients and correlates positively with the secretion of endogenous INF-γ and IL-6 cytokines [42]. Therefore, we speculate that RhoA and TRAb promote the development of GO at the immune level, but whether there is a direct relationship between them requires more detailed evidence

5. Limitation

Relatively few patients met the inclusion criteria, and some patients were excluded because of loss to follow-up, poor compliance, and the use of other oral drugs. These factors may have affected the results of the experiment. Therefore, the study was limited by the inadequate sample size, which we will aim to expand by continuing to collect relevant medical records. At the same time, the conclusions we reached based on the current sample size are consistent with those drawn by others, indicating that the results of this experiment are credible despite the limitations.

6. Conclusion

RhoA is a potential biomarker in the mild stage of GO. If elevation of RhoA in the serum of GD patients is accompanied by elevation of TRAb or FT3, GO is highly likely to occur, even in the absence of obvious ocular symptoms. Our study provides a new biomarker for early diagnosis of GO, and we hope that it will allow patients with GO to be diagnosed and treated in a timelier manner in the future.

Data Availability

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

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

Sidi Zhao, Shuangshuang Shi, and Wanchen Yang contributed equally to this work and share first authorship.

Acknowledgments

This study was supported by the Science and Technology Project of Tianjin Municipal Education Commission (NO. 2020KJ176), Science and Technology Project of Tianjin Binhai New Area Health Care Committee (NO. 2019BWKZ009), and Chengdu Medical Scientific Research Project of Chengdu Municipal Health Commission.

Supplementary Materials

Supplementary Table 1: demographic and ophthalmological data of the study participants, i.e., GO-free and GO groups. Supplementary Table 2: Mauchly’s sphericity test of thyroid hormone levels in the two groups. Supplementary Table 3: analysis of variance of repeated-measurement data for the two groups. (Supplementary Materials)