Association of AIRE Polymorphisms with Genetic Susceptibility to Rheumatoid Arthritis in a Chinese Population

Recently, genetic polymorphisms within the autoimmune regulator (AIRE) have been implicated in the genetic susceptibility to rheumatoid arthritis (RA) in Japanese and Spanish. The aim of this case–control study involving 232 patients with RA and 313 ethnically matched control subjects was to investigate the association of AIRE rs2075876 and rs760426 polymorphisms with genetic predisposition to RA in a Chinese population. The genotypes of AIRE rs2075876 and rs760426 polymorphisms were determined by SNaPshot assay. A significant difference in the allele frequency of AIRE rs2075876 polymorphism between cases and controls was detected (A versus G, OR 1.33, 95 %CI 1.04–1.69, P?=?0.02, P corrected (Bonferroni correction) Pc?=?0.04). Significant evidence was found for the association between the minor allele A of AIRE rs2075876 polymorphism and the risk of RA under the recessive model (AA versus AG?+?GG, P?=?7.15?×?10^?3, Pc?=?1.43?×?10^?2). The frequency of the minor allele G of AIRE rs760426 polymorphism was higher in patients compared with controls (47.8 % versus 42.1 %), and this deviation showed a trend towards significant level (P?=?0.06, Pc?=?0.12). The association between the minor allele G of AIRE rs760426 polymorphism with RA risk under the dominant model and the recessive model revealed that significant evidence was detected under the recessive model (GG versus GA?+?AA, P?=?0.02, Pc?=?0.04). Our results indicated that AIRE rs2075876 and rs760426 polymorphisms were involved in the genetic background of RA in the Chinese population.     

Fine mapping the TAGAP risk locus in rheumatoid arthritis

A common allele at the TAGAP gene locus demonstrates a suggestive, but not conclusive association with risk of rheumatoid arthritis (RA). To fine map the locus, we conducted comprehensive imputation of CEU HapMap single-nucleotide polymorphisms (SNPs) in a genome-wide association study (GWAS) of 5,500 RA cases and 22,621 controls (all of European ancestry). After controlling for population stratification with principal components analysis, the strongest signal of association was to an imputed SNP, rs212389 (P=3.9 × 10(-8), odds ratio=0.87). This SNP remained highly significant upon conditioning on the previous RA risk variant (rs394581, P=2.2 × 10(-5)) or on a SNP previously associated with celiac disease and type I diabetes (rs1738074, P=1.7 × 10(-4)). Our study has refined the TAGAP signal of association to a single haplotype in RA, and in doing so provides conclusive statistical evidence that the TAGAP locus is associated with RA risk. Our study also underscores the utility of comprehensive imputation in large GWAS data sets to fine map disease risk alleles.     

Genetic analysis of autoimmune regulator haplotypes in alopecia areata

Alopecia areata is an immune-mediated disorder, occurring with the highest observed frequency in the rare recessive autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy (APECED) syndrome caused by mutations of the autoimmune regulator (AIRE) gene on chromosome 21q22.3. We have previously detected association between alopecia areata and a single nucleotide polymorphism (SNP) in the AIRE gene in patients without APECED, and we now report the findings of an extended examination of the association of alopecia areata with haplotype analysis including six SNPs in the AIRE gene: C-103T, C4144G, T5238C, G6528A, T7215C and T11787C. In Caucasian groups of 295 patients and 363 controls, we found strong association between the AIRE 7215C allele and AA [P = 3.8 x 10(-8), OR (95% CI): 2.69 (1.8-4.0)]. The previously reported association between AA and the AIRE 4144G allele was no longer significant on correction for multiple testing. The AIRE haplotypes CCTGCT and CGTGCC showed a highly significant association with AA [P = 6.05 x 10(-6), 9.47 (2.91-30.8) and P = 0.001, 3.51 (1.55-7.95), respectively]. To select the haplotypes most informative for analysis, we tagged the polymorphisms using SNPTag software. Employing AIRE C-103T, G6528A, T7215C and T11787C as tag SNPs, two haplotypes were associated with AA; AIRE CGCT and AIRE CGCC [P = 3.84 x 10(-7), 11.40 (3.53-36.9) and P = 3.94 x 10(-4), 2.13 (1.39-3.24) respectively]. The AIRE risk haplotypes identified in this study potentially account for a major component of the genetic risk of developing alopecia areata.     

Large-scale fine mapping of the HNF1B locus and prostate cancer risk

Previous genome-wide association studies have identified two independent variants in HNF1B as susceptibility loci for prostate cancer risk. To fine-map common genetic variation in this region, we genotyped 79 single nucleotide polymorphisms (SNPs) in the 17q12 region harboring HNF1B in 10 272 prostate cancer cases and 9123 controls of European ancestry from 10 case-control studies as part of the Cancer Genetic Markers of Susceptibility (CGEMS) initiative. Ten SNPs were significantly related to prostate cancer risk at a genome-wide significance level of P < 5 × 10(-8) with the most significant association with rs4430796 (P = 1.62 × 10(-24)). However, risk within this first locus was not entirely explained by rs4430796. Although modestly correlated (r(2)= 0.64), rs7405696 was also associated with risk (P = 9.35 × 10(-23)) even after adjustment for rs4430769 (P = 0.007). As expected, rs11649743 was related to prostate cancer risk (P = 3.54 × 10(-8)); however, the association within this second locus was stronger for rs4794758 (P = 4.95 × 10(-10)), which explained all of the risk observed with rs11649743 when both SNPs were included in the same model (P = 0.32 for rs11649743; P = 0.002 for rs4794758). Sequential conditional analyses indicated that five SNPs (rs4430796, rs7405696, rs4794758, rs1016990 and rs3094509) together comprise the best model for risk in this region. This study demonstrates a complex relationship between variants in the HNF1B region and prostate cancer risk. Further studies are needed to investigate the biological basis of the association of variants in 17q12 with prostate cancer.     

A genome-wide association study in the Japanese population confirms 9p21 and 14q23 as susceptibility loci for primary open angle glaucoma

Primary open angle glaucoma (POAG) is one of leading causes of adult blindness worldwide. To identify genetic variants associated with susceptibility to POAG, we conducted a genome-wide association study (GWAS) using 1394 cases and 6599 controls. Subsequently, we analyzed 33 single nucleotide polymorphisms (SNPs) which showed suggestive association (P < 1 × 10(-4)) by GWAS, using an additional set of 1802 cases and 7212 controls. In addition to confirmation of the association of the chromosome 9p21 locus [rs1063192, P= 5.2 × 10(-11), odds ratio (OR) = 0.75], and 14q23 (rs10483727, P = 9.49 × 10(-8), OR = 0.79) with POAG in Caucasians reported recently, we identified a suggestive-associated locus on 2q21 (rs7588567, P = 3.89 × 10(-7), OR = 0.85). For these described SNPs, minor alleles are suspected to have a protective effect from the disease. An linkage disequilibrium block containing rs10483727 includes the SIX6 gene that was implicated to have a critical role in eye development, and genes in both represented loci (SIX6 on chromosome 14q23, and CDKN2A-CDKN2B on chromosome 9p21) are known to be expressed in human ocular tissues, including the retina. Our GWAS results should contribute to better insight into the genetic basis of POAG.     

The chromosome 2p21 region harbors a complex genetic architecture for association with risk for renal cell carcinoma

In follow-up of a recent genome-wide association study (GWAS) that identified a locus in chromosome 2p21 associated with risk for renal cell carcinoma (RCC), we conducted a fine mapping analysis of a 120 kb region that includes EPAS1. We genotyped 59 tagged common single-nucleotide polymorphisms (SNPs) in 2278 RCC and 3719 controls of European background and observed a novel signal for rs9679290 [P = 5.75 × 10(-8), per-allele odds ratio (OR) = 1.27, 95% confidence interval (CI): 1.17-1.39]. Imputation of common SNPs surrounding rs9679290 using HapMap 3 and 1000 Genomes data yielded two additional signals, rs4953346 (P = 4.09 × 10(-14)) and rs12617313 (P = 7.48 × 10(-12)), both highly correlated with rs9679290 (r(2) > 0.95), but interestingly not correlated with the two SNPs reported in the GWAS: rs11894252 and rs7579899 (r(2) < 0.1 with rs9679290). Genotype analysis of rs12617313 confirmed an association with RCC risk (P = 1.72 × 10(-9), per-allele OR = 1.28, 95% CI: 1.18-1.39) In conclusion, we report that chromosome 2p21 harbors a complex genetic architecture for common RCC risk variants.     

A genome-wide association study identifies locus at 10q22 associated with clinical outcomes of adjuvant tamoxifen therapy for breast cancer patients in Japanese

Although many association studies of polymorphisms in candidate genes with the clinical outcomes of breast cancer patients receiving adjuvant tamoxifen therapy have been reported, genetic factors determining individual response to tamoxifen are not fully understood. To identify genetic polymorphisms associated with clinical outcomes of patients with tamoxifen treatment, we conducted a genome-wide association study (GWAS). We studied 462 Japanese patients with hormone receptor-positive, invasive breast cancer receiving adjuvant tamoxifen therapy. Of them, 240 patients were analyzed by genome-wide genotyping using the Illumina Human610-Quad BeadChips, and two independent sets of 105 and 117 cases were used for replication studies. In the GWAS, we detected significant associations with recurrence-free survival at 15 single-nucleotide polymorphisms (SNPs) on nine chromosomal loci (1p31, 1q41, 5q33, 7p11, 10q22, 12q13, 13q22, 18q12 and 19p13) that satisfied a genome-wide significant threshold (log-rank P= 2.87 × 10(-9)-9.41 × 10(-8)). Among them, rs10509373 in C10orf11 gene on 10q22 was significantly associated with recurrence-free survival in the replication study (log-rank P= 2.02 × 10(-4)) and a combined analysis indicated a strong association of this SNP with recurrence-free survival in breast cancer patients treated with tamoxifen (log-rank P= 1.26 × 10(-10)). Hazard ratio per C allele of rs10509373 was 4.51 [95% confidence interval (CI), 2.72-7.51; P= 6.29 × 10(-9)]. In a combined analysis of rs10509373 genotype with previously identified genetic makers, CYP2D6 and ABCC2, the number of risk alleles of these three genes had cumulative effects on recurrence-free survival among 345 patients receiving tamoxifen monotherapy (log-rank P= 2.28 × 10(-12)). In conclusion, we identified a novel locus associated with recurrence-free survival in Japanese breast cancer patients receiving adjuvant tamoxifen therapy.     

MYO9B gene polymorphisms are associated with autoimmune diseases in Spanish population

The aim of the study was to test MYO9B gene polymorphisms for association with three autoimmune diseases, systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), and celiac disease (CD), in a Spanish population. We analyzed three SNPs (rs2305767, rs1457092, and rs2305764) in a case-control cohort composed of 349 SLE patients, 356 RA patients, 90 CD patients, and 345 healthy controls. All three SNPs showed a consistent increased frequency of the A allele in SLE, RA, and CD patients compared with healthy controls. An association was observed between CD and rs2305764 (p=0.01, OR=2.3), between SLE and rs1457092 (p=0.002, OR=1.4), and between RA and rs1457092 (p=0.02, OR=1.3). The three autoimmune diseases combined showed significant association with rs1457092 and rs2305764 and with the AAA haplotype (p haplotype=0.005, OR=1.3). Our data demonstrate consistent association with the A allele and AAA haplotype of three SNPs in the MYO9B gene, which were previously reported to be associated with CD in the Dutch population. This suggests that genetic variation in MYO9B is associated with CD, SLE, and RA and that MYO9B is a general risk factor for autoimmunity.     

Genome-wide association study confirms extant PD risk loci among the Dutch

In view of the population-specific heterogeneity in reported genetic risk factors for Parkinson's disease (PD), we conducted a genome-wide association study (GWAS) in a large sample of PD cases and controls from the Netherlands. After quality control (QC), a total of 514,799 SNPs genotyped in 772 PD cases and 2024 controls were included in our analyses. Direct replication of SNPs within SNCA and BST1 confirmed these two genes to be associated with PD in the Netherlands (SNCA, rs2736990: P = 1.63 × 10(-5), OR = 1.325 and BST1, rs12502586: P = 1.63 × 10(-3), OR = 1.337). Within SNCA, two independent signals in two different linkage disequilibrium (LD) blocks in the 3' and 5' ends of the gene were detected. Besides, post-hoc analysis confirmed GAK/DGKQ, HLA and MAPT as PD risk loci among the Dutch (GAK/DGKQ, rs2242235: P = 1.22 × 10(-4), OR = 1.51; HLA, rs4248166: P = 4.39 × 10(-5), OR = 1.36; and MAPT, rs3785880: P = 1.9 × 10(-3), OR = 1.19).     

A genome-wide association study identifies a novel susceptibility locus for renal cell carcinoma on 12p11.23

Renal cell carcinoma (RCC) is the most lethal urologic cancer. Only two common susceptibility loci for RCC have been confirmed to date. To identify additional RCC common susceptibility loci, we conducted an independent genome-wide association study (GWAS). We analyzed 533 191 single nucleotide polymorphisms (SNPs) for association with RCC in 894 cases and 1516 controls of European descent recruited from MD Anderson Cancer Center in the primary scan, and validated the top 500 SNPs in silico in 3772 cases and 8505 controls of European descent involved in the only published GWAS of RCC. We identified two common variants in linkage disequilibrium, rs718314 and rs1049380 (r(2) = 0.64, D?' = 0.84), in the inositol 1,4,5-triphosphate receptor, type 2 (ITPR2) gene on 12p11.23 as novel susceptibility loci for RCC (P = 8.89 × 10(-10) and P = 6.07 × 10(-9), respectively, in meta-analysis) with an allelic odds ratio of 1.19 [95% confidence interval (CI): 1.13-1.26] for rs718314 and 1.18 (95% CI: 1.12-1.25) for rs1049380. It has been recently identified that rs718314 in ITPR2 is associated with waist-hip ratio (WHR) phenotype. To our knowledge, this is the first genetic locus associated with both cancer risk and WHR.     

Positive association of the <Emphasis Type="Italic">FTSJ1</Emphasis> gene polymorphisms with nonsyndromic X-linked mental retardation in young Chinese male subjects

To investigate the possible genetic association of nonsyndromic X-linked mental retardation (NS-XLMR) with FTSJ1 gene polymorphisms, a case-control association study was performed focusing on the Chinese Han population in the Qinba mountain region. Three common single nucleotide polymorphisms (SNPs) (rs2268954, rs2070991, rs5905692) in the gene were selected and genotyped using the polymerase chain reaction single-strand confirmation polymorphism (PCR-SSCP) method. Pairwise linkage disequilibrium (LD) analysis showed that the three SNPs were in strong LD (all D' > 0.8). There were significant differences between cases and controls in allele frequency distribution of rs2268954 (P = 0.036), rs2070991 (P = 0.043), and rs5905692 (P = 0.014) and in the distributions of common haplotypes combined by these SNPs (global P = 0.01236) in male subjects. In female subjects, however, no positive results were found. Our results suggest a positive association between the genetic variants of the FTSJ1 gene and NS-XLMR in young male subjects in the Chinese Han population in the Qinba region.     

Genetic variants on chromosome 6p21.1 and 6p22.3 are associated with type 2 diabetes risk: a case-control study in Han Chinese

Recent genome-wide association studies have identified several single nucleotide polymorphisms (SNPs) on chromosome 6p21.1 and 6p22.3 as type 2 diabetes (T2D) susceptibility loci in the European and Japanese populations. However, these SNPs have not been well evaluated in Chinese population. Here, we performed a case-control study with 2925 T2D cases and 3281 controls in a Chinese population. We used TaqMan OpenArray and Sequenom MassARRAY to genotype the four SNPs (rs4712523, rs7756992, rs4712524 and rs6931514) in CDKAL1 (cyclin-dependent kinase 5 regulatory subunit-associated protein 1-like 1) at 6p22.3 and one SNP (rs9472138) near vascular endothelial growth factor A (VEGFA) at 6p21.1. All the five SNPs were significantly associated with T2D risk with overall effects (odds ratio, OR) from 1.19 to 1.29 in the additive genetic model (rs6931514: OR=1.29, 95% confidence intervals (95% CI)=1.19-1.39, P=5.6 × 10(-10); rs7756992: OR=1.23, 95% CI=1.15-1.32, P=1.2 × 10(-8); rs4712523: OR=1.25, 95% CI=1.15-1.35, P=3.8 × 10(-8); rs4712524: OR=1.24, 95% CI=1.15-1.35, P=6.8 × 10(-8); rs9472138: OR=1.19, 95% CI=1.05-1.34, P=006). Conditional analysis identified two independent signals (rs6931514 at 6p22.3 and rs9472138 at 6p21.1) that were significantly associated with T2D. Compared with the wild homozygote of rs6931514 and rs9472138, subjects with variant alleles of the two SNPs had increased risk for T2D susceptibility in a dose-response manner (P(trend)=7.4 × 10(-12)). Our findings indicated that genetic variants of CDKAL1 and VEGFA on chromosome 6 may contribute to T2D risk in Chinese population, especially for rs9472138 at 6p21.1 identified for the first time to significantly increase the T2D risk in Chinese individuals.     

A genome-wide association study of bladder cancer identifies a new susceptibility locus within SLC14A1, a urea transporter gene on chromosome 18q12.3

Genome-wide and candidate-gene association studies of bladder cancer have identified 10 susceptibility loci thus far. We conducted a meta-analysis of two previously published genome-wide scans (4501 cases and 6076 controls of European background) and followed up the most significant association signals [17 single nucleotide polymorphisms (SNPs) in 10 genomic regions] in 1382 cases and 2201 controls from four studies. A combined analysis adjusted for study center, age, sex, and smoking status identified a novel susceptibility locus that mapped to a region of 18q12.3, marked by rs7238033 (P = 8.7 × 10(-9); allelic odds ratio 1.20 with 95% CI: 1.13-1.28) and two highly correlated SNPs, rs10775480/rs10853535 (r(2)= 1.00; P = 8.9 × 10(-9); allelic odds ratio 1.16 with 95% CI: 1.10-1.22). The signal localizes to the solute carrier family 14 member 1 gene, SLC14A1, a urea transporter that regulates cellular osmotic pressure. In the kidney, SLC14A1 regulates urine volume and concentration whereas in erythrocytes it determines the Kidd blood groups. Our findings suggest that genetic variation in SLC14A1 could provide new etiological insights into bladder carcinogenesis.     

ANKRD55 and DHCR7 are novel multiple sclerosis risk loci

Multiple sclerosis (MS) shares some risk genes with other disorders hallmarked by an autoimmune pathogenesis, most notably IL2RA and CLEC16A. We analyzed 10 single-nucleotide polymorphisms (SNPs) in nine risk genes, which recently emerged from a series of non-MS genome-wide association studies (GWAS), in a Spanish cohort comprising 2895 MS patients and 2942 controls. We identified two SNPs associated with MS. The first SNP, rs6859219, located in ANKRD55 (Chr5), was recently discovered in a meta-analysis of GWAS on rheumatoid arthritis (RA), and emerged from this study with genome-wide significance (odds ratio (OR) = 1.35; P = 2.3 × 10(-9)). The second SNP, rs12785878, is located near DHCR7 (Chr11), a genetic determinant of vitamin D insufficiency, and showed a size effect in MS similar to that recently observed in Type 1 diabetes (T1D; OR = 1.10; P = 0.009). ANKRD55 is a gene of unknown function, and is flanked proximally by the IL6ST-IL31RA gene cluster. However, rs6859219 did not show correlation with a series of haplotype-tagging SNPs covering IL6ST-IL31RA, analyzed in a subset of our dataset (D'< 0.31; r(2)< 0.011). Our results expand the number of risk genes shared between MS, RA and T1D.     

Variations in/nearby genes coding for JAZF1, TSPAN8/LGR5 and HHEX-IDE and risk of type 2 diabetes in Han Chinese

Several genetic loci (JAZF1, CDC123/CAMK1D, TSPAN8/LGR5, ADAMTS9, VEGFA and HHEX-IDE) were identified to be significantly related to the risk of type 2 diabetes and quantitative metabolic traits in European populations. Here, we aimed to evaluate the impacts of these novel loci on type 2 diabetes risk in a population-based case-control study of Han Chinese (1912 cases and 2041 controls). We genotyped 13 single-nucleotide polymorphisms (SNPs) in/near these genes and examined the differences in allele/genotype frequency between cases and controls. We found that both IDE rs11187007 and HHEX rs1111875 were associated with type 2 diabetes risk (for both variants: odds ratio (OR)=1.15, 95% confidence interval (CI) 1.04-1.28, P=0.009). In a meta-analysis where we pooled our data with the three previous studies conducted in East Asians, we found that the variants of JAZF1 rs864745 (1.09 (1.03-1.16); P=3.49 × 10(-3)) and TSPAN8/LGR5 rs7961581 (1.11(1.05-1.17); P=1.89 × 10(-4)) were significantly associated with type 2 diabetes risk. In addition, the meta-analysis (7207 cases and 8260 controls) also showed that HHEX rs1111875 did have effects on type 2 diabetes in Chinese population (OR=1.15(1.10-1.21); P=1.93 × 10(-8)). This large population-based study and meta-analysis further confirmed the modest effects of the JAZF1, TSPAN8/LGR5 and HHEX-IDE loci on type 2 diabetes in Chinese and other East Asians.     

Fine mapping of a region of chromosome 11q13 reveals multiple independent loci associated with risk of prostate cancer

Genome-wide association studies have identified prostate cancer susceptibility alleles on chromosome 11q13. As part of the Cancer Genetic Markers of Susceptibility (CGEMS) Initiative, the region flanking the most significant marker, rs10896449, was fine mapped in 10 272 cases and 9123 controls of European origin (10 studies) using 120 common single nucleotide polymorphisms (SNPs) selected by a two-staged tagging strategy using HapMap SNPs. Single-locus analysis identified 18 SNPs below genome-wide significance (P< 10(-8)) with rs10896449 the most significant (P= 7.94 × 10(-19)). Multi-locus models that included significant SNPs sequentially identified a second association at rs12793759 [odds ratio (OR) = 1.14, P= 4.76 × 10(-5), adjusted P= 0.004] that is independent of rs10896449 and remained significant after adjustment for multiple testing within the region. rs10896438, a proxy of previously reported rs12418451 (r(2)= 0.96), independent of both rs10896449 and rs12793759 was detected (OR = 1.07, P= 5.92 × 10(-3), adjusted P= 0.054). Our observation of a recombination hotspot that separates rs10896438 from rs10896449 and rs12793759, and low linkage disequilibrium (rs10896449-rs12793759, r(2)= 0.17; rs10896449-rs10896438, r(2)= 0.10; rs12793759-rs10896438, r(2)= 0.12) corroborate our finding of three independent signals. By analysis of tagged SNPs across ～123 kb using next generation sequencing of 63 controls of European origin, 1000 Genome and HapMap data, we observed multiple surrogates for the three independent signals marked by rs10896449 (n= 31), rs10896438 (n= 24) and rs12793759 (n= 8). Our results indicate that a complex architecture underlying the common variants contributing to prostate cancer risk at 11q13. We estimate that at least 63 common variants should be considered in future studies designed to investigate the biological basis of the multiple association signals.