International Journal of Genomics

Volume 2014 (2014), Article ID 537969, 7 pages

http://dx.doi.org/10.1155/2014/537969

Research Article

Associations of CTLA4 Gene Polymorphisms with Graves’ Ophthalmopathy: A Meta-Analysis

Pengfei Du,¹ Xiaojie Ma,¹ and Changjiang Wang²

¹Department of Endocrinology, The Second Affiliated Hospital of Jiaxing University, Jiaxing 314000, China

²Department of Clinical Pharmacy, The Second Affiliated Hospital of Jiaxing University, Jiaxing 314000, China

Received 5 May 2014; Accepted 18 June 2014; Published 9 July 2014

Academic Editor: Giulia Piaggio

Copyright © 2014 Pengfei Du et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Many studies have established that T-lymphocyte antigen-4 (CTLA4) is a susceptible gene for Graves’ disease (GD). Also many studies showed the association between the CTLA4 exon-1 49A/G polymorphism and the risk of developing Graves’ ophthalmopathy (GO) in GD patients. But those results were inconsistent. In recent years many new studies were published which helped to shed light on the relationship of CTLA4 SNP49 with GO. So we performed the meta-analysis to explore the association between the SNP49 and GO susceptibility in GD patients. Studies up to February 29, 2012, were searched by using PubMed. The odds ratio was used to evaluate the strength of the association. Altogether 12 case-control studies involving 2,505 participants were included in the meta-analysis. Results showed that the G allele was related to the increased risk of GO compared with the A allele under allelic genetic model (OR = 1.14, 95% CI: 1.14–1.72, ) in European subgroup. No publication bias was detected. Our results showed that the SNP49 polymorphism of CTLA4 gene was related to increased risk of GO.

1. Introduction

Graves’ disease is a thyroid autoimmune disorder with 25–50% individuals having ophthalmopathy which is called thyroid-associated ophthalmopathy (TAO) or GO. The clinical manifestations of GO can mostly be explained by the increased volume of extraocular muscles and orbital connective tissues [1]. The symptoms of GO progress in some patients which eventually will lead to blindness. Though many studies have established the role of T cells in GO development [2–4], the understanding of the mechanism of GO is still poor. More effective therapy may develop on a better understanding of the T cell-fibroblast interaction.

Human CTLA-4 gene, located on 2q33, encodes a molecule which plays an important role in the downregulation of CD28 interaction with the ligands on the surface of antigen-presenting cells (APCs). The important inhibitory role of CTLA-4 in T-cell function has made it be one candidate gene when exploring autoimmune diseases.

Many reports have established the association of CTLA-4 with GD [5–9]. The CTLA-4 exon-1 +49A/G polymorphism has been found to be associated with the development of GO in GD patients in some studies but with inconsistent results [10–24].

One meta-analysis exploring the association of CTLA-4 gene polymorphism with GO was performed in 2006 which found no definite result [21]. Many new researches have been done since 2006 which will help to clarify the association. So we performed the meta-analysis with more case-control study results.

2. Methods

2.1. Search Strategy

We searched PubMed up to February 29, 2012, for all possible publications on the association between the CTLA-4 49A/G polymorphism and GO. Search terms were used as follows: “Cytotoxic T-lymphocyte associated antigen-4 or CTLA4,” “Graves’ disease,” and “ophthalmopathy.” We also manually searched all reference lists of the relevant articles for additional papers. Language was limited to English. All eligible studies met the following criteria: (1) case-control design investigating the associations of CTLA4 49A/G polymorphism and GO; (2) providing genotype distribution information in both cases and controls or provided odds ratio (OR) with 95% confidence interval (95% CI) (or sufficient data that allowed us to calculate these). When a study reported results on different subpopulations, we treated each subpopulation as a separate comparison in the meta-analysis. For overlapping data, the most complete or recent study was included.

2.2. Data Extraction

Two authors (Pengfei Du and Changjiang Wang) independently extracted data and reached consensus on all items. The following information was sought from each report: authors, journal and year of publication, country of origin, selection and characteristics of GO cases and controls, demographics, ethnic group of the study population, eligible and genotyped cases and controls, and number of cases and controls for the genotype. For studies including subjects of different ethnic groups, data were extracted separately for each ethnicity whenever possible.

2.3. Statistical Analysis

The primary analysis compared GO cases with controls for the frequency of G versus A alleles. This analysis aims to detect overall differences. Pooled OR with 95% CI was used to evaluate the strength of the associations of CTLA-4 49A/G polymorphism and GO. When zero events occurred, we treated this problem by adding 0.5 to all the 2 2 cells of the contingency table to calculate the study-specific OR just as implemented in RevMan software. We examined the association of 49A/G polymorphism and GO for allelic contrast (G versus A). Furthermore, we conducted subgroup analyses by stratifying ethnicity into Europeans and Asians separately. A test was used to determine the significance of the pooled OR and value was considered significant. The between-study heterogeneity across all eligible comparisons was estimated using the chi-square based Q statistic [25]. Heterogeneity was considered significant for . Data were combined using both fixed-effects (Mantel-Haenszel) and random effects (DerSimonian and Laird) models [26]. Random effects incorporate an estimate of the between-study variance and tend to provide wider confidence intervals, when the results of the constituent studies differ among themselves. In the absence of between-study heterogeneity, the two methods provide identical results. Random effects are more appropriate when heterogeneity is present [26]. Inverted funnel plots and the Begg-Mazumdar publication bias diagnostic (nonparametric correlation coefficient) [27] evaluated whether the magnitude of the observed association was related to the variance of each study; that is, whether large studies gave different results compared with smaller ones. All the statistical analyses were performed using RevMan5.0 software (the Cochrane Collaboration, Oxford, England). All values were two-sided.

3. Results

3.1. Data Summary

The primary literature search yielded 17 papers [10–24, 28, 29]. After further scrutiny, 12 potential relevant papers were identified. Of these papers, five papers were excluded for the following reasons: two papers include overlapping data [10, 29] with another two studies, and only the detailed one was selected [11, 21]; three studies [12, 14, 28] were not included because their genotype or allele distribution information was not available. One of the eligible studies [19] comprised subjects of two different racial descents, so it was evaluated as two studies. Finally, a total of 12 papers and 13 separate comparisons were included in our meta-analysis [11, 13, 15–24]. The general characteristics of the selected studies are shown in Table 1. Controls were GD patients without GO. These 12 papers all evaluated the association of 49A/G polymorphism with GO. Five studies were conducted in Asian populations [15, 19, 21–23] and the others were in European population. The sample size of individual studies ranged from 33 to 222. Polymerase chain reaction-restriction fragment-length polymorphism (PCR-RFLP), DCFH (dual-color fluorescence hybridization), and minisequencing were used as genotyping methods. Table 2 shows the genotype frequencies in included studies.

Table 1: General characteristics of the selected studies in the meta-analysis.

Ethnic
Country
SNP genotyping
Grouping method
Eligible subjects
First author (year) [reference]

GO
Non-GO
GO
Non-GO

European
USA
RFLP
NOSPEC3–6
No clinical feature
85
52
Villanueva (2000) [16]

UK
RFLP
NOSPEC3–6
NOSPEC < 3
161
323
Allahabadia (2001) [17]

Slovenia
RFLP
Clinical features
No clinical feature
33
34
Zaletel (2002) [18]

Iran
RFLP
NOSPEC3–6
NOSPEC < 3
105
100
Esteghamati (2009) [24]

Poland
RFLP
NOSPEC3–6
No clinical feature
95
169
Bednarczuk (2003) [19]

UK
RFLP
NOSPEC3–6
NOSPEC < 3
124
168
Vaidya (2003) [11]

Poland
Minisequencing
NOSPEC3–6
NOSPEC < 3
50
49
Frydecka (2004) [13]

Italy
RFLP
NOSPEC2–6
NOSPEC < 3
90
60
Petrone (2005) [20]

Asian
China
DCFH
NOSPEC3–6
NOSPEC < 3
142
119
Han (2006) [21]

China
RFLP
NOSPEC3–6
No clinical feature
33
56
Zhang (2006) [22]

China
RFLP
NOSPEC2
NOSPEC3, 4
34
22
Chong (2008) [23]

Japan
RFLP
NOSPEC3–6
No clinical feature
99
220
Bednarczuk (2003) [19]

Korea
RFLP
Clinical feature
No clinical feature
26
47
Park (2000) [15]

Table 2: Distribution of CTLA4 49A/G alleles among GO and controls in the included studies.

First author
Country
A/A ()
A/G ()
G/G ()

GO
Non-GO
GO
Non-GO
GO
Non-GO

Villanueva (2000) [16]
USA
26
16
42
25
17
11

Allahabadia (2001) [17]
UK
53
83
82
180
26
60

Zaletel (2002) [18]
Slovenia
11
14
17
17
5
3

Esteghamati (2009) [24]
Iran
48
66
43
28
14
6

Bednarczuk (2003) [19]
Poland
27
48
41
82
27
39

Vaidya (2003) [11]
UK
27
61
62
77
40
34

Frydecka (2004) [13]
Poland
11
21
29
21
10
7

Petrone (2005) [20]
Italy
NA
NA
NA
NA
NA
NA

Han (2006) [21]
China
18
14
55
39
69
66

Zhang (2006) [22]
China
1
1
11
18
21
37

Chong (2008) [23]
China
1
0
9
8
24
14

Bednarczuk (2003) [19]
Japan
12
16
41
99
46
105

Park (2000) [15]
Korea
1
3
10
15
15
29

NA: not available.

3.2. Association of the 49A/G Polymorphism with GO

A total of 12 studies concerning the 49A/G polymorphism included 2,505 individuals (1,082 cases and 1423 controls). The pooled frequency of G allele was 68.1% and 51.8% for cases and controls, respectively. Two studies suggested an at-risk effect of G allele [11, 24] and the others did not produce significant results. The overall analysis for evaluating the association of 49A/G polymorphism with GO under allelic model showed that there was large heterogeneity between these 13 studies ( heterogeneity = 0.01), and the random effects pooled OR was not significant: OR = 1.14, 95% CI: 0.94–1.38, (Figure 1). One study was excluded from the meta-analysis because of the contribution to the heterogeneity in the European subgroup (random OR = 1.40, 95% CI 1.14–1.72, , heterogeneity = 0.19) [17]. Sensitivity analyses evaluating the influence of single studies on overall risk estimate by omitting one study in each turn did not materially alter the overall effect size, with a range from OR = 1.32 (95% CI 1.09–1.60) to OR = 1.48 (95% CI 1.19–1.85). The GG genotype is associated with increased risk of GO in European subgroup (random OR = 1.89, 95% CI 1.26–2.83, , heterogeneity = 0.32) (Figure 2). The G carrier was more likely to have GO than A carrier in European subgroup (random OR = 1.31, 95% CI 1.07–1.60, , heterogeneity = 0.01) (Figure 3).

(a)

(b)

Figure 1: Meta-analysis for the effect of the G allele versus the A allele on the risk of GO in GD patients. Each comparison is presented by the name of the first author and the year of publication. The point estimate of the odds ratio and the accompanying 95% confidence interval (CI) are shown. “Total” represents the summary random effects estimation for the comparison along with the respective 95% confidence interval. Values above 1 denote an increased risk for GO with the G allele.

Figure 2: Meta-analysis for the effect of the GG versus the AA genotype on the risk of GO in GD patients in European subgroup. Each comparison is presented by the name of the first author and the year of publication. The point estimate of the odds ratio and the accompanying 95% confidence interval (CI) are shown. “Total” represents the summary random effects estimation for the comparison along with the respective 95% confidence interval. Values above 1 denote an increased risk for GO with the GG genotype.

Figure 3: Meta-analysis for the effect of the CTLA4 49A/G allele in G carriers and in A carriers on the risk of GO in GD patients. Each comparison is presented by the name of the first author and the year of publication. The point estimate of the odds ratio and the accompanying 95% confidence interval (CI) are shown. “Total” represents the summary random effects estimation for the comparison along with the respective 95% confidence interval. Values above 1 denote an increased risk for GO with the G allele.

4. Discussion

The current meta-analysis including 12 case-control studies was in an effort to clarify the relationship between CTLA-4 gene polymorphism and GO susceptibility. The overall results indicated that the CTLA-4 49A/G polymorphism was associated with susceptibility of GO which was more significant in European population. Although obvious heterogeneity was detected for 49A/G associations, sensitivity analyses did not materially alter the overall and subgroup results under different genetic models, indicating that the results were stable and reliable.

The polymorphism of CTLA-4 had been suggested in many diseases which included Addison’s disease, autoimmune hypothyroidism, and rheumatoid arthritis. The CTLA-4 49A/G SNP in exon-1 leads to the substitution of Ala with Thr in the signal peptide part which was reported to cause misprocessing of CTLA-4 in the ER resulting in less efficient glycosylation and diminished surface expression of CTLA-4 protein [30]. Studies have demonstrated that the longer repeats of the 3′UTR microsatellite are associated with reduced CTLA-4 inhibitory function [31]. Some workers have shown association between the G allele of 49A/G and reduced control of T-cell proliferation [9, 32]. One study found an association between the protective allele A of the 49A/G SNP, which is associated with increased CTLA-4 function/expression, and interferon induced thyroiditis [33]. Hence, we postulated that, under most cases of spontaneous autoimmunity, alleles associated with reduced CTLA-4 function and/or expression could lead to increased activation of T cells thereby triggering autoimmunity.

The association of CTLA-4 49A/G SNP with GO was more significant in European population but not in Asian population. The CTLA4 49A/G SNP shows the same effect on the GO development with different significance. The results are confusing but not conflicting. This can be explained by the following reasons. Firstly, both GD and GO are multigenic conditions. It is unlikely to have one major susceptibility gene contributing to the development of GO. In the studies of Xu et al., they transiently transfected a T-cell line with a CTLA-4 construct harboring either the G or the A allele of the 49A/G SNP and found no difference in CTLA-4 expression and/or function harboring the A or the G allele [34]. It is likely that the 49A/G SNP did not directly affect the T-cell activation but may in LD with some major effect alleles. Secondly, the impact of smoking on GO was not specifically considered in the association studies. As was reported, the odds ratio for both current smokers and exsmokers was high and the GO was more severe in smokers than in nonsmokers [35–37].

In conclusion, our analysis first intensified the association between 49A/G and GO susceptibility, which showed that the G allele is a risk factor for GO. More risk factors should be considered in future association studies of the CTLA-4 SNP with GO. The gene-gene and gene-environment interaction should also be taken into consideration in future research. The identification of the candidate gene and understanding of the mechanisms by which they cause disease will identify those individuals at risk for developing GO in the future. Treatment will evolve towards more preventive strategies. The treatment will be personalized by targeting the specific pathways that are most contributory to GO development. We believe that the future will witness the engagement of novel molecules, such as altered peptide ligands [38] and specific monoclonal antibodies, to modulate specific autoimmune pathways in affected individuals, thus treating the cause of the disease.

Conflict of Interests

The authors declare that they have no financial and personal relationships with other people or organizations that can inappropriately influence their work.

References

R. S. Bahn and A. E. Heufelder, “Pathogenesis of Graves' ophthalmopathy,” The New England Journal of Medicine, vol. 329, no. 20, pp. 1468–1475, 1993.

T. F. Davies, A. Martin, E. S. Concepcion, P. Graves, L. Cohen, and A. Ben-Nun, “Evidence of limited variability of antigen receptors on intrathyroidal T cells in autoimmune thyroid disease,” The New England Journal of Medicine, vol. 325, no. 4, pp. 238–244, 1991.

T. F. Davies, A. Martin, E. S. Conception et al., “Widespread Vβ and preferential VαT cell receptor gene utilization by intrathyroidal T cells in human autoimmune thyroid disease,” The Journal of Clinical Investigation, vol. 89, no. 1, pp. 157–162, 1992.

A. E. Heufelder, S. Herterich, G. Ernst, R. S. Bahn, and P. C. Scriba, “Analysis of retroorbital T cell antigen receptor variable region gene usage in patients with Graves' ophthalmopathy,” European Journal of Endocrinology, vol. 132, no. 3, pp. 266–277, 1995.

T. Yanagawa, Y. Hidaka, V. Guimaraes, M. Soliman, and L. J. DeGroot, “CTLA-4 gene polymorphism associated with Graves' disease in a Caucasian population,” Journal of Clinical Endocrinology and Metabolism, vol. 80, no. 1, pp. 41–45, 1995.

T. Yanagawa, M. Taniyama, S. Enomoto et al., “CTLA4 gene polymorphism confers susceptibility to Graves' disease in Japanese,” Thyroid, vol. 7, no. 6, pp. 843–846, 1997.

J. M. Heward, A. Allahabadia, M. Armitage et al., “The development of Graves' disease and the CTLA-4 gene on chromosome 2q33,” Journal of Clinical Endocrinology and Metabolism, vol. 84, no. 7, pp. 2398–2401, 1999.

B. Vaidya, H. Imrie, P. Perros et al., “The cytotoxic T lymphocyte antigen-4 is a major Graves' disease locus,” Human Molecular Genetics, vol. 8, no. 7, pp. 1195–1199, 1999.

T. Kouki, Y. Sawai, C. A. Gardine, M. E. Fisfalen, M. L. Alegre, and L. J. DeGroot, “CTLA-4 Gene polymorphism at position 49 in exon 1 reduces the inhibitory function of CTLA-4 and contributes to the pathogenesis of Graves' disease,” Journal of Immunology, vol. 165, no. 11, pp. 6606–6611, 2000.

B. Vaidya, H. Imrie, P. Perros et al., “Cytotoxic T lymphocyte antigen-4 (CTLA-4) gene polymorphism confers susceptibility to thyroid associated orbitopathy,” The Lancet, vol. 354, no. 9180, pp. 743–744, 1999.

B. Vaidya, E. J. C. Oakes, H. Imrie et al., “CTLA4 gene and Graves' disease: association of Graves' disease with the CTLA4 exon 1 and intron 1 polymorphisms, but not with the promoter polymorphism,” Clinical Endocrinology, vol. 58, no. 6, pp. 732–735, 2003.

R. Buzzetti, L. Nistico, A. Signore, and I. Cascino, “CTLA-4 and HLA gene susceptibility to thyroid-associated orbitopathy,” The Lancet, vol. 354, no. 9192, p. 1824, 1999.

I. Frydecka, J. Daroszewski, K. Suwalska et al., “CTLA-4 (CD152) gene polymorphism at position 49 in exon 1 in Graves' disease in a Polish population of the Lower Silesia region,” Archivum Immunologiae et Therapiae Experimentalis, vol. 52, no. 5, pp. 369–374, 2004.

K. Badenhoop, H. Donner, J. Braun, T. Siegmund, H. Rau, and K. H. Usadel, “Genetic markers in diagnosis and prediction of relapse in Graves' disease,” Experimental and Clinical Endocrinology and Diabetes, vol. 104, no. 4, pp. 98–100, 1996.

Y. J. Park, H. K. Chung, D. J. Park et al., “Polymorphism in the promoter and exon 1 of the cytotoxic T lymphocyte antigen-4 gene associated with autoimmune thyroid disease in Koreans,” Thyroid, vol. 10, no. 6, pp. 453–459, 2000.

R. Villanueva, A. M. Inzerillo, Y. Tomer et al., “Limited genetic susceptibility to severe graves' ophthalmopathy: no role for CTLA-4 but evidence for an environmental etiology,” Thyroid, vol. 10, no. 9, pp. 791–798, 2000.

A. Allahabadia, J. M. Heward, R. Nithiyananthan et al., “MHC class IIregion, CTLA4 gene, and ophthalmopathy in patients with Graves' disease,” The Lancet, vol. 358, no. 9286, pp. 984–985, 2001.

K. Zaletel, B. Krhin, S. Gaberšček, E. Pirnat, and S. Hojker, “The influence of the exon 1 polymorphism of the cytotoxic T lymphocyte antigen 4 gene on thyroid antibody production in patients with newly diagnosed Graves' disease,” Thyroid, vol. 12, no. 5, pp. 373–376, 2002.

T. Bednarczuk, Y. Hiromatsu, T. Fukutani et al., “Association of cytotoxic T-lymphocyte-associated antigen-4 (CTLA-4) gene polymorphism and non-genetic factors with Graves' ophthalmopathy in European and Japanese populations,” European Journal of Endocrinology, vol. 148, no. 1, pp. 13–18, 2003.

A. Petrone, G. Giorgi, A. Galgani et al., “CT60 single nucleotide polymorphisms of the cytotoxic T-lymphocyte- associated antigen-4 gene region is associated with graves' disease in an Italian population,” Thyroid, vol. 15, no. 3, pp. 232–238, 2005.

S. Han, S. Zhang, W. Zhang et al., “CTLA4 polymorphisms and ophthalmopathy in Graves' disease patients: association study and meta-analysis,” Human Immunology, vol. 67, no. 8, pp. 618–626, 2006.

Q. Zhang, Y. M. Yang, and X. Y. Lv, “Association of Graves' disease and Graves' ophthalmopathy with the polymorphisms in promoter and exon 1 of cytotoxic T lymphocyte associated antigen-4 gene,” Journal of Zhejiang University: Science B, vol. 7, no. 11, pp. 887–891, 2006.

K. K. Chong, S. W. Chiang, G. W. Wong et al., “Association of CTLA-4 and IL-13 gene polymorphisms with Graves' disease and ophthalmopathy in Chinese children,” Investigative Ophthalmology and Visual Science, vol. 49, no. 6, pp. 2409–2415, 2008.

A. Esteghamati, O. Khalilzadeh, Z. Mobarra et al., “Association of CTLA-4 gene polymorphism with Graves' disease and ophthalmopathy in Iranian patients,” European Journal of Internal Medicine, vol. 20, no. 4, pp. 424–428, 2009.

J. Lau, J. P. A. Ioannidis, and C. H. Schmid, “Quantitative synthesis in systematic reviews,” Annals of Internal Medicine, vol. 127, no. 9, pp. 820–826, 1997.

D. B. Petitti, Meta-Analysis, Decision Analysis and Cost-Effectiveness Analysis, Oxford University Press, New York, NY, USA, 1994.

C. B. Begg and M. Mazumdar, “Operating characteristics of a rank correlation test for publication bias,” Biometrics, vol. 50, no. 4, pp. 1088–1101, 1994.

J. Daroszewski, E. Pawlak, L. Karabon et al., “Soluble CTLA-4 receptor an immunological marker of Graves' disease and severity of ophthalmopathy is associated with CTLA-4 Jo31 and CT60 gene polymorphisms,” European Journal of Endocrinology, vol. 161, no. 5, pp. 787–793, 2009.

S. Z. Han, S. H. Zhang, R. Li, W. Y. Zhang, and Y. Li, “The common -318C/T polymorphism in the promoter region of CTLA4 gene is associated with reduced risk of ophthalmopathy in Chinese Graves' patients,” International Journal of Immunogenetics, vol. 33, no. 4, pp. 281–287, 2006.

S. Anjos, A. Nguyen, H. Ounissi-Benkalha, M. C. Tessier, and C. Polychronakos, “A common autoimmunity predisposing signal peptide variant of the cytotoxic T-lymphocyte antigen 4 results in inefficient glycosylation of the susceptibility allele,” The Journal of Biological Chemistry, vol. 277, no. 48, pp. 46478–46486, 2002.

M. Takara, T. Kouki, and L. J. DeGroot, “CTLA-4 AT-repeat polymorphism reduces the inhibitory function of CTLA-4 in Graves' disease,” Thyroid, vol. 13, no. 12, pp. 1083–1089, 2003.

Y. Ban, T. F. Davies, D. A. Greenberg et al., “Analysis of the CTLA-4, CD28, and inducible costimulator (ICOS) genes in autoimmune thyroid disease,” Genes and Immunity, vol. 4, no. 8, pp. 586–593, 2003.

E. M. Jacobson, S. Chaudhry, J. C. Mandac, E. Concepcion, and Y. Tomer, “Immune regulatory gene involvement in the etiology of interferon induced thyroiditis (IIT),” Thyroid, vol. 16, no. 9, article 926, 2006.

Y. Xu, P. N. Graves, Y. Tomer, and T. F. Davies, “CTLA-4 and autoimmune thyroid disease: lack of influence of the A49G signal peptide polymorphism on functional recombinant human CTLA-4,” Cellular Immunology, vol. 215, no. 2, pp. 133–140, 2002.

E. Hagg and K. Asplund, “Is endocrine ophthalmopathy related to smoking?” The British Medical Journal, vol. 295, no. 6599, pp. 634–635, 1987.

J. Pfeilschifter and R. Ziegler, “Smoking and endocrine ophthalmopathy: impact of smoking severity and current vs lifetime cigarette consumption,” Clinical Endocrinology, vol. 45, no. 4, pp. 477–481, 1996.

B. Shine, P. Fells, O. M. Edwards, and A. P. Weetman, “Association between Graves' ophthalmopathy and smoking,” The Lancet, vol. 335, no. 8700, pp. 1261–1263, 1990.

B. Bielekova and R. Martin, “Antigen-specific immunomodulation via altered peptide ligands,” Journal of Molecular Medicine, vol. 79, no. 10, pp. 552–565, 2001.