Clinical and Developmental Immunology

Volume 2013 (2013), Article ID 891306, 6 pages

http://dx.doi.org/10.1155/2013/891306

Research Article

Replication of British Rheumatoid Arthritis Susceptibility Loci in Two Unrelated Chinese Population Groups

Hua Li,¹ Yonghe Hu,² Tao Zhang,¹ Yang Liu,³ Yantang Wang,³ Tai Yang,³ Minhui Li,³ Qiaoli Luo,¹ Yu Cheng,¹ and Qiang Zou³

¹Department of Oncology, Chengdu Military General Hospital, Chengdu 610083, China

²Department of Traditional Chinese Medicine, Chengdu Military General Hospital, Chengdu 610083, China

³Department of Immunology, Chengdu Medical College, Chengdu 610083, China

Received 29 April 2013; Accepted 2 August 2013

Academic Editor: Jianying Zhang

Copyright © 2013 Hua Li 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.

Previous genome-wide association study by WTCCC identified many susceptibility loci of common autoimmune diseases in British, including rheumatoid arthritis (RA). Because of the genetic heterogeneity of RA, it is necessary to replicate these susceptibility loci in other populations. Here, three SNPs with strong RA association signal in the British were analyzed in Han Chinese, and two SNPs (rs6457617 and rs11761231) were genotyped in the test cohort firstly. The rs6457617 was significantly associated with RA in the test cohort. The individuals bearing the homozygous genotype CC had 0.39-fold risk than these bearing the wild-type genotype TT (, OR 0.39, [95% CI 0.21–0.74]). And the protective effect of allele C was confirmed in another validation cohort with 1514 samples (   ×  10⁻¹⁰, OR 0.34, [95% CI 0.24–0.48]). The rs6457617 can be used as a tagSNP of HLA-DQA1*03 which encoded MHC-II α chain. Since MHC restriction is important for primary T-cells in positive selection and negative selection stages, MHC protein polymorphisms may be implicated in shaping the T-cell repertoire, including the emergence of a T-cell clone involved in the inflammatory arthritis.

1. Introduction

Rheumatoid arthritis (RA (MIM 180300)) is a common autoimmune disease characterized by chronic inflammatory, destructive, and debilitating arthritis. The etiology of RA, like that of other autoimmune disorders, is complex and remains elusive. The occurrence of RA is relatively constant with a prevalence of between 0.5 and 1.0% in several European populations [1, 2], North-American populations [3], Japan [4], and China [5]. But some native American-Indian populations have high prevalence of RA such as the Pima Indians (5.3%) [6] and Chippewa Indians (6.8%) [7]. Although it is difficult to elucidate whether the environmental or genetic effect influence the differences between populations in different countries, it is thought to have both a genetic and an environment basis [8–11]. The heritability of RA has been estimated to be about 60% [12]. The highly polymorphic HLA region is estimated to account for about one-third of the total genetic component of susceptibility [13]. Many genes outside the HLA region also account for the RA risk genetic factor.

There are two strategies used commonly to detect the RA risk genetic loci. One is the candidate gene association study, and another is the genome-wide association study (GWAS). Candidate gene association studies rely on prior knowledge of the biology of the disease and the highly subjective selection of potential genes. It is hard to find new RA susceptibility loci outside of HLA. Genome-wide association studies (GWAS) have an advantage that it entails a systematically search of the entire genome for susceptibility variants without any clue about pathogenesis. With the advent of GWAS, relatively large number of new potential susceptibility loci for RA in some populations has been identified [14, 15]. One of the most impressive GWASs was the one by the WTCCC, which scanned 500,568 SNPs (Affymetrix Chip) in 14,000 cases of seven major autoimmune diseases, including 2000 United Kingdom RA cases and 3000 controls [16]. Among the loci showing the strongest association signals, there were 3 RA susceptibility SNPs (rs6679677, rs6457617, and rs11761231) with a value of less than .

Meanwhile, much evidence shows that there is the genetic heterogeneity of RA across the major racial groups. For example, the haplotype in STAT4 gene associated with RA in Caucasians is not associated in the Han Chinese population, but with the presence of rheumatoid factor [17]. The rs247661 in PTPN22 gene, which was associated with RA in Europeans [18], does not exist in Han Chinese. The PADI4 gene, which was originally identified as RA susceptibility gene in large Japanese and Korean cohorts [19], was not associated with RA and unlikely to be responsible for the presence of anti-CCP autoantibodies in Caucasian populations [20]. Therefore, it is necessary to test RA susceptibility loci of GWAS in different populations. Here, the three SNPs identified by the WTCCC with a value of less than were taken into our consideration to investigate whether they are associated with RA in Han Chinese population.

2. Materials and Methods

2.1. Subjects Evaluated in the Association Study

Total 1894 RA patients and healthy individuals were recruited for this study. The test cohort was 380 Chengdu residents recruited from the Chengdu Medical College, and the validation cohort was 1514 Chongqing residents recruited from the Southwest Hospital, as reported in our previous study [21]. All the subjects were unrelated ethnic Han Chinese, according to self-reported ancestry. Patients have an established diagnosis of RA according to the 1987 revised criteria of the ACR for the classification of the disease. In each cohort, healthy controls were individually matched to RA cases on the basis of sex, age, ethnicity, and local residential region. There was no significant difference between the 2 cohorts in terms of the distribution of age, age at RA onset, or percentage of female subjects (Table 2). Written informed consent was obtained from all subjects, and the study was performed with the approval of the ethical committee of the Chengdu Military General Hospital.

2.2. SNP Selection and Genotyping Methods

Three SNPs with strong RA association signal in British population were rs6457617, rs11761231, and rs6679677. The rs6679677 has no diversity in Han Chinese population based on the HCB data of HapMap and did not genotyped in our study. The rs6457617 is in the MHC region, and the minor allele frequency (MAF) in HCB is 0.465. The rs11761231 locates on chromosome 7q32, and the MAF is 0.233. The rs6457617 and rs11761231 were genotyped using the SNPStream Ultra High Throughput Genotyping system (Beckman Coulter, Fullerton, CA, USA) according to the manufacturer’s instructions. The primers and probes of each SNP were listed in Table 1. The genotype concordance rate was 100% as assessed by random retyping across different plates.

Table 1: PCR primers and tagged extension probes used for SNP detection.

SNPs
Locus
Gene
Primer
Sequence 5′-3′

rs6457617
6q
MHC
Forward
AAATGCAGTCAGTGGACTCAA

Reverse
AAAACAAAAAAAACCCTTCAATC

Probe
GACCTGGGTGTCGATACCTACTGTTTGTTGAGTCCATGAGCAGAT

rs11761231
7q32
Unknown
Forward
TGTCTTATCATGAGAACGTGCA

Reverse
TTTTGTGTTCAAAGAATTCTGTCTT

Probe
AGAGCGAGTGACGCATACTATGAAATCAAGAAGGTCTGAAAA

Table 2: Clinical characteristics of the test and validation cohorts.

Subjects
Test cohort
Validation cohort

Number of RA cases
190
757

Number of healthy controls
190
757

Study design
Unrelated cases with individually matched controls
Unrelated cases with individually matched controls

Ethnicity
Han Chinese
Han Chinese

City of residence
Chengdu
Chongqing

Female, number (%)
300 (78.9)
1253 (82.7)

Age, mean ± SD years
46.4 ± 8.0
46.8 ± 10.2

Age at RA onset, mean ± SD years*
39.4 ± 9.6
41.0 ± 10.1

Disease duration, mean ± SD years*
6.97 ± 6.2
5.7 ± 5.4

*Refers to rheumatoid arthritis (RA) patients only.

2.3. Statistical Analysis

Hardy-Weinberg equilibrium was assessed using a chi-square goodness-of-fit test. Case-control association analyses were performed by chi-square or Fisher’s exact test, as appropriate. values, odds ratios (ORs), and 95% confidence intervals (95% CIs) were calculated.

3. Results and Discussion

All SNPs were in Hardy-Weinberg equilibrium () in both patients and healthy controls, indicating that our subjects were random-mating population with no selection, mutation, or migration. As shown in Table 3, rs6457617 was significantly associated with RA in the test cohort. The frequency of the minor allele C was significantly lower in RA patients compared to healthy control (28.9% of the patients versus 40.8% of the controls; , OR 0.59, [95% CI 0.44–0.80]), implying that the minor allele C was protective against RA. When we compared the distribution of the homozygous genotype CC or the heterozygous genotype CT among the RA patients and controls to the one of the wild type genotype TT, the protective effect of allele C was obvious and presented a dose effect relationship. The individuals bearing the homozygous genotype CC had 0.39-fold decreased risk than these bearing the wild-type genotype TT (, OR 0.39, [95%  CI 0.21–0.74]). while individuals bearing the heterozygous CT had 0.59-fold decreased risk for RA than those bearing wild-type genotype (, OR 0.59, [95%  CI 0.38–0.92], Table 3).

Table 3: Allele and genotype frequencies of SNPs in test cohort.

SNP
Chromosome
Gene
Allele
Allele analysis
Genotype
Genotype analysis

Controls number (%)
Cases number (%)
OR (95% CI)

Controls number (%)
Cases number (%)
OR (95% CI)

rs6457617
6q
MHC
T
218 (59.2)
263 (71.1)
1 (reference)
—
TT
68 (37.0)
97 (52.4)
1 (reference)
—

C
150 (40.8)
107 (28.9)
0.59 (0.44–0.80)
0.001
CT
82 (44.6)
69 (37.3)
0.59 (0.38–0.92)
0.02

CC
34 (18.5)
19 (10.3)
0.39 (0.21–0.74)
0.004

rs11761231
7q32
Unknown
T
287 (75.9)
290 (76.7)
1 (reference)
—
TT
108 (57.1)
112 (59.3)
1 (reference)
—

C
91 (24.1)
88 (23.3)
0.96 (0.68–1.34)
0.80
CT
71 (37.6)
66 (34.9)
0.90 (0.59–1.37)
0.62

CC
10 (5.3)
11 (5.8)
1.06 (0.43–2.60)
0.90

rs11761231 in female group
7q32
Unknown
T
218 (73.2)
220 (73.8)
1 (reference)
—
TT
79 (53.0)
82 (55.0)
1 (reference)
—

C
80 (26.8)
78 (26.2)
0.97 (0.67–1.39)
0.85
CT
60 (40.3)
56 (37.6)
0.90 (0.56–1.45)
0.66

CC
10 (6.7)
11 (7.4)
1.06 (0.43–2.63)
0.90

rs11761231 in male group
7q32
Unknown
T
69 (86.3)
70 (87.5)
1 (reference)
—
TT
29 (72.5)
30 (75.0)
1 (reference)
—

C
11 (13.7)
10 (12.5)
0.90 (0.36–2.25)
0.82
CT
11 (27.5)
10 (25.0)
1.14 (0.42–3.08)
1.0

CC
0
0

Data shown in (%). OR = odds ratio; 95% CI = 95% confidence interval; MHC: major histocompatibility complex; allele analyses for association were used by Pearson chi-square test, and genotype analyses for association were used by logistic regression.

The strong association signal of rs6457617 was replicated in the validation cohort with 1514 subjects. Among 1470 samples successfully genotyped, each genotype was assessed with the use of codominant, dominant, and recessive genetic models (Table 4). The significant association was found in all genetic models. The individuals with the homozygous genotype CC had the lowest risk than those with wild-type genotype (, OR 0.34, [95% CI 0.24–0.48]). Because the allele frequencies were similar across the controls tested in both cohorts, we combined all of the data and obtained stronger statistical evidence of an association between this SNP and RA susceptibility. Compared to individuals with genotype TT, the Han Chinese bearing genotype CC had a lower decreased risk (, OR = 0.35, 95%  CI: 0.26–0.47) than heterozygous carriers had (, OR = 0.68, 95%  CI: 0.56–0.82). The allele T was RA risk allele both in the British and the Han Chinese. But the Han Chinese bearing genotype TT had a 2.86-fold (1/0.35) increased risk compared to those with genotype CC, much lower than it in the British ().

Table 4: Validation of SNPs in test cohort, validation cohort, and combined analyses (combined validation cohort with test cohort).

Genotype
Subjects number (%)
Codominant genetic model
Dominant genetic model
Recessive genetic model

Controls
Cases
OR (95% CI)

OR (95% CI)

OR (95% CI)

Test cohort rs6457617
TT
68 (37.0)
97 (52.4)
1 (ref)
—
0.53 (0.35–0.81)
0.003
0.51 (0.28–0.92)
0.025

CT
82 (44.6)
69 (37.3)
0.59 (0.38–0.92)
0.02

CC
34 (18.5)
19 (10.3)
0.39 (0.21–0.74)
0.004

Validation cohort rs6457617
TT
263 (35.7)
353 (48.2)
1 (ref)
—
0.60 (0.49–0.74)

0.41 (0.30–0.56)

CT
340 (46.1)
319 (43.5)
0.70 (0.56–0.87)

CC
134 (18.2)
61 (8.3)
0.34 (0.24–0.48)

Combined samples rs6457617
TT
331 (35.9)
450 (49.0)
1 (ref)
—
0.58 (0.48–0.70)

0.43 (0.32–0.57)

CT
422 (45.8)
388 (42.3)
0.68 (0.56–0.82)

CC
168 (18.3)
80 (8.7)
0.35 (0.26–0.47)

Data shown in (%). OR = odds ratio; 95% CI = 95% confidence interval; allele analysis and genotype analysis in dominant genetic model or in recessive genetic model for association were used by Pearson chi-square test, and genotype analyses in codominant genetic model for association were used by logistic regression.

This rs6457617 was identified in another GWAS in 800 Spanish subjects and replicated in an independent cohort of 794 Spanish subjects [22]. For the Asian population, this SNP was also replicated in the Korean [23]. The minor allele was C in Korean and Chinese but was T in the WTCCC British samples. With the data of each genotype frequency in RA cases and controls in Korean in their Table 3 of the previous report [23], we found that the genotype TT in Korean had a 4.77-fold increased risk compared to the Korean with genotype CC. The OR in Korean is between the one (OR 2.86) in Chinese and another (OR 5.21) in British. Recently, this SNP was replicated in another replication study in north Indians () [24]. But this report has not shown more detail data about rs6457617. Although there are different frequencies of the alleles C and T in different populations, the association between rs6457617 and RA susceptibility exists in all these different populations, and the risk allele for RA is the same T allele, strongly suggesting that rs6457617 in HLA-DQ region is the real RA susceptibility locus.

The rs6457617 locates at 32,663,850 on 6p21.3 between the HLA-DQB1 (32,627,241-32,634,466) encoded chain of MHC-II protein and HLA-DQA2 (32,709,163-32,714,664) encoded chain of MHC-II protein. The MHC protein is an heterodimer protein receptor that is typically expressed on the surface of antigen-presenting cells. MHC loci are the most genetically variable loci in the human population. HLA DQ is highly variable, the subunit more so than the chain. rs6457617 can be used as a tagSNP for HLA-DQA1*03 [25, 26], which encodes MHC-II chain. As we know, DRB1*0401, *0404, *0408, and *01, DQA*01, *0201, and *04 are the shared epitope (SE) associated with RA [27], while HLA-DQA1*03 has been reported to be associated with chronic HCV infection [28], type 1 diabetes [29, 30], and childhood-onset ocular myasthenia gravis [31]. Up to now, whether HLA-DQA1*03 is associated with RA has not been identified. Since rs6457617 is the RA susceptibility locus and can be a tagSNP for HLA-DQA1*03, HLA-DQA1*03 may be a SE of RA. Highly variable MHC-II molecules recognize and present different antigens to T-cells. Since MHC restriction is particularly important for primary T-cells in positive selection and negative selection stages, MHC protein polymorphisms may play an important role in shaping the T cell repertoire, including the emergence of an unusual T cell clone characterized by the potential of inflammatory arthritis. Whether the rs6457617 is a functional SNP and influences the expression or the structure of chain is needed to identify.

As for rs11761231, no association evidence was found in the allele and genotype analyses. We noticed that in the male Han Chinese, including RA patients and healthy controls, none of them bears genotype CC. The stratification of patients and controls for gender was performed. However, there is no association between rs11761231 and RA disease either in female group or in male group (Table 3). This SNP locates within the eukaryotic translation elongation factor 1 beta 2 pseudogene. Although it had very strong association signal of RA in British population, it did not show any evidence of association in Han Chinese populations. The minor allele C frequency (MAF) in Han Chinese is only 0.278, quite different from the MAF in CEU (0.407). Moreover, the prevalence of RA in China is lower than it in European. Therefore, the test cohort with 380 subjects might have no enough statistical power to detect the genetic effect of rs11761231. So, to investigate the association between rs11761231 and RA in Han Chinese, the case-control study with sufficiently large sample size is needed.

4. Conclusions

In this study, rs6457617 in the MHC region, which was identified as RA susceptibility locus in previous GWAS study in British, was associated with RA in the Han Chinese populations. The minor allele C has dependent protective effect of RA risk. The rs6457617 is a real RA susceptibility locus and could be a good predictor of RA risk.

Conflict of Interests

The authors declare that there is no conflict of interests in this study.

Acknowledgments

This work was supported by grants from the New Century Excellent Talents in University (No. NCET-09-0888), National Natural Science Foundation of China (no. 81072455), The Basic Research Program of Sichuan Province of China (no. 2012JY0031), Sichuan Youth Science and Technology Fund (no. 2010JQ0035), Research Fund of Chengdu Military General Hospital (no. 2011YG-B02), the Key Programs of the Military Medical Sciences during the 12th five-year plan period (no. BWS11J067), and Key Project of Chengdu Military Region during the 12th five-year plan period (no. A12003). The authors are grateful to the RA patients and the healthy subjects who made this study possible. The aurhors would also like to thank the Director Xiaolan Zhao, Dr. Haiyan Fu, other clinicians at the Health Care Center, Director Yongfei Fang, Dr. Bin Zhong, Dr. Jinyi Li, and other clinicians in the Department of Rheumatology for contributing to the study and collecting the samples and thank Dr. Ying Wang for genotyping assistance.

References

L. Carmona, V. Villaverde, C. Hernández-García, J. Ballina, R. Gabriel, and A. Laffon, “The prevalence of rheumatoid arthritis in the general population of Spain,” Rheumatology, vol. 41, no. 1, pp. 88–95, 2002.

T. Riise, B. K. Jacobsen, and J. T. Gran, “Incidence and prevalence of rheumatoid arthritis in the county of Troms, northern Norway,” Journal of Rheumatology, vol. 27, no. 6, pp. 1386–1389, 2000.

S. E. Gabriel, “The epidemiology of rheumatoid arthritis,” Rheumatic Disease Clinics of North America, vol. 27, no. 2, pp. 269–281, 2001.

K. Shichikawa, K. Inoue, S. Hirota et al., “Changes in the incidence and prevalence of rheumatoid arthritis in Kamitonda, Wakayama, Japan, 1965–1996,” Annals of the Rheumatic Diseases, vol. 58, no. 12, pp. 751–756, 1999.

Q. Zeng, S. Huang, and R. Chen, “10-year epidemiological study on rheumatic diseases in Shantou area,” Zhonghua Nei Ke Za Zhi, vol. 36, no. 3, pp. 193–197, 1997.

A. Del Puente, W. C. Knowler, D. J. Pettitt, and P. H. Bennett, “High incidence and prevalence of rheumatoid arthritis in Pima Indians,” American Journal of Epidemiology, vol. 129, no. 6, pp. 1170–1178, 1989.

J. Harvey, M. Lotze, and M. B. Stevens, “Rheumatoid arthritis in a Chippewa band. I. Pilot screening study of disease prevalence,” Arthritis & Rheumatism, vol. 24, no. 5, pp. 717–721, 1981.

J. Bowes and A. Barton, “Recent advances in the genetics of RA susceptibility,” Rheumatology, vol. 47, no. 4, pp. 399–402, 2008.

H. Källberg, L. Padyukov, R. M. Plenge et al., “Gene-gene and gene-environment interactions involving HLA-BRB1, PTPN22, and smoking in two subsets of rheumatoid arthritis,” American Journal of Human Genetics, vol. 80, no. 5, pp. 867–875, 2007.

L. Padyukov, C. Suva, P. Stolt, L. Alfredsson, and L. Klareskog, “A gene-environment interaction between smoking and shared epitope genes in HLA-DR provides a high risk of seropositive rheumatoid arthritis,” Arthritis & Rheumatism, vol. 50, no. 10, pp. 3085–3092, 2004.

S. Viatte, D. Plant, and S. Raychaudhuri, “Genetics and epigenetics of rheumatoid arthritis,” Nature Reviews Rheumatology, vol. 9, no. 3, pp. 141–153, 2013.

A. J. MacGregor, H. Snieder, A. S. Rigby et al., “Characterizing the quantitative genetic contribution to rheumatoid arthritis using data from twins,” Arthritis & Rheumatism, vol. 43, no. 1, pp. 30–37, 2000.

C. M. Deighton, D. J. Walker, I. D. Griffiths, and D. F. Roberts, “The contribution of HLA to rheumatoid arthritis,” Clinical Genetics, vol. 36, no. 3, pp. 178–182, 1989.

W. Thomson, A. Barton, X. Ke et al., “Rheumatoid arthritis association at 6q23,” Nature Genetics, vol. 39, no. 12, pp. 1431–1433, 2007.

Z. Wei and M. Li, “Genome-wide linkage and association analysis of rheumatoid arthritis in a Canadian population,” BMC Proceedings, vol. 1, supplement 1, article S19, 2007.

The Wellcome Trust Case Control Consortium, “Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls,” Nature, vol. 447, no. 7145, pp. 661–678, 2007.

H. Li, Q. Zou, Z. Xie et al., “A haplotype in STAT4 gene associated with rheumatoid arthritis in Caucasians is not associated in the Han Chinese population, but with the presence of rheumatoid factor,” Rheumatology, vol. 48, no. 11, pp. 1363–1368, 2009.

G. G. Song, S. C. Bae, J. H. Kim, and Y. H. Lee, “The PTPN22 C1858T polymorphism and rheumatoid arthritis: a meta-analysis,” Rheumatology International, vol. 33, no. 8, pp. 1991–1999, 2013.

A. Suzuki, R. Yamada, X. Chang et al., “Functional haplotypes of PADI4, encoding citrullinating enzyme peptidylarginine deiminase 4, are associated with rheumatoid arthritis,” Nature Genetics, vol. 34, no. 4, pp. 395–402, 2003.

G. Poór, Z. B. Nagy, Z. Schmidt, M. Brózik, K. Merétey, and P. Gergely Jr., “Genetic background of anticyclic citrullinated peptide autoantibody production in Hungarian patients with rheumatoid arthritis,” Annals of the New York Academy of Sciences, vol. 1110, pp. 23–32, 2007.

J. Sui, H. Li, Y. Fang et al., “NLRP1 gene polymorphism influences gene transcription and is a risk factor for rheumatoid arthritis in han Chinese,” Arthritis & Rheumatism, vol. 64, no. 3, pp. 647–654, 2012.

A. Julià, J. Ballina, J. D. Cañete et al., “Genome-wide association study of rheumatoid arthritis in the Spanish population: KLF12 as a risk locus for rheumatoid arthritis susceptibility,” Arthritis & Rheumatism, vol. 58, no. 8, pp. 2275–2286, 2008.

T. Han, S. Bang, C. Kang, and S. Bae, “TRAF1 polymorphisms associated with rheumatoid arthritis susceptibility in asians and in caucasians,” Arthritis & Rheumatism, vol. 60, no. 9, pp. 2577–2584, 2009.

P. Prasad, A. Kumar, R. Gupta, R. C. Juyal, and B. K. Thelma, “Caucasian and Asian specific rheumatoid arthritis risk loci reveal limited replication and apparent allelic heterogeneity in north Indians,” PLoS ONE, vol. 7, no. 2, Article ID e31584, 2012.

K. Nakanishi and Y. Shima, “Capture of type 1 diabetes-susceptible HLA DR-DQ haplotypes in Japanese subjects using a tag single nucleotide polymorphism,” Diabetes Care, vol. 33, no. 1, pp. 162–164, 2010.

P. I. W. de Bakker, G. McVean, P. C. Sabeti et al., “A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC,” Nature Genetics, vol. 38, no. 10, pp. 1166–1172, 2006.

S. Laivoranta-Nyman, T. Möttönen, R. Hermann et al., “HLA-DR-DQ haplotypes and genotypes in Finnish patients with rheumatoid arthritis,” Annals of the Rheumatic Diseases, vol. 63, no. 11, pp. 1406–1412, 2004.

C. Tibbs, P. Donaldson, J. Underhill, L. Thomson, K. Manabe, and R. Williams, “Evidence that the HLA DQA1*03 allele confers protection from chronic HCV-infection in northern European Caucasoids,” Hepatology, vol. 24, no. 6, pp. 1342–1345, 1996.

K. Nakanishi and H. Inoko, “Combination of HLA-A24, -DQA1*03, and -DR9 contributes to acute-onset and early complete β-cell destruction in type 1 diabetes: longitudinal study of residual β-cell function,” Diabetes, vol. 55, no. 6, pp. 1862–1868, 2006.

J. E. van Autreve, I. Weets, B. Gulbis, F. Vertongen, F. K. Gorus, and B. J. van der Auwera, “The rare HLA-DQA1*03-DQB1*02 haplotype confers susceptibility to type 1 diabetes in whites and is preferentially associated with early clinical disease onset in male subjects,” Human Immunology, vol. 65, no. 7, pp. 729–736, 2004.

W. Zhu, J. Lu, J. Lin et al., “HLA-DQA1*03:02/DQB1*03:03:02 is strongly associated with susceptibility to childhood-onset ocular myasthenia gravis in Southern Han Chinese,” Journal of Neuroimmunology, vol. 247, no. 1-2, pp. 81–85, 2012.