Identification of novel genetic markers associated with clinical phenotypes of systemic sclerosis through a genome-wide association strategy

The aim of this study was to determine, through a genome-wide association study (GWAS), the genetic components contributing to different clinical sub-phenotypes of systemic sclerosis (SSc). We considered limited (lcSSc) and diffuse (dcSSc) cutaneous involvement, and the relationships with presence of the SSc-specific auto-antibodies, anti-centromere (ACA), and anti-topoisomerase I (ATA). Four GWAS cohorts, comprising 2,296 SSc patients and 5,171 healthy controls, were meta-analyzed looking for associations in the selected subgroups. Eighteen polymorphisms were further tested in nine independent cohorts comprising an additional 3,175 SSc patients and 4,971 controls. Conditional analysis for associated SNPs in the HLA region was performed to explore their independent association in antibody subgroups. Overall analysis showed that non-HLA polymorphism rs11642873 in IRF8 gene to be associated at GWAS level with lcSSc (P = 2.32×10⁻¹², OR = 0.75). Also, rs12540874 in GRB10 gene (P = 1.27 × 10⁻⁶, OR = 1.15) and rs11047102 in SOX5 gene (P = 1.39×10⁻⁷, OR = 1.36) showed a suggestive association with lcSSc and ACA subgroups respectively. In the HLA region, we observed highly associated allelic combinations in the HLA-DQB1 locus with ACA (P = 1.79×10⁻⁶¹, OR = 2.48), in the HLA-DPA1/B1 loci with ATA (P = 4.57×10⁻⁷⁶, OR = 8.84), and in NOTCH4 with ACA P = 8.84×10⁻²¹, OR = 0.55) and ATA (P = 1.14×10⁻⁸, OR = 0.54). We have identified three new non-HLA genes (IRF8, GRB10, and SOX5) associated with SSc clinical and auto-antibody subgroups. Within the HLA region, HLA-DQB1, HLA-DPA1/B1, and NOTCH4 associations with SSc are likely confined to specific auto-antibodies. These data emphasize the differential genetic components of subphenotypes of SSc.     

��-Catenin Is Genetically and Biologically Associated with Cortical Cataract and Future Alzheimer-Related Structural and Functional Brain Changes

Multiple lines of evidence suggest that specific subtypes of age-related cataract (ARC) and Alzheimer disease (AD) are related etiologically. To identify shared genetic factors for ARC and AD, we estimated co-heritability of quantitative measures of cataract subtypes with AD-related brain MRI traits among 1,249 members of the Framingham Eye Study who had a brain MRI scan approximately ten years after the eye exam. Cortical cataract (CC) was found to be co-heritable with future development of AD and with several MRI traits, especially temporal horn volume (THV, ρ = 0.24, P<10⁻⁴). A genome-wide association study using 187,657 single nucleotide polymorphisms (SNPs) for the bivariate outcome of CC and THV identified genome-wide significant association with CTNND2 SNPs rs17183619, rs13155993 and rs13170756 (P<2.6×10⁻⁷). These SNPs were also significantly associated with bivariate outcomes of CC and scores on several highly heritable neuropsychological tests (5.7×10⁻⁹≤P<3.7×10⁻⁶). Statistical interaction was demonstrated between rs17183619 and APP SNP rs2096488 on CC (P = 0.0015) and CC-THV (P = 0.038). A rare CTNND2 missense mutation (G810R) 249 base pairs from rs17183619 altered δ-catenin localization and increased secreted amyloid-β_1–42 in neuronal cell culture. Immunohistopathological analysis of lens tissue obtained from two autopsy-confirmed AD subjects and two non-AD controls revealed elevated expression of δ-catenin in epithelial and cortical regions of lenses from AD subjects compared to controls. Our findings suggest that genetic variation in delta catenin may underlie both cortical lens opacities in mid-life and subsequent MRI and cognitive changes that presage the development of AD.     

Thyroid stimulating hormone receptor (TSHR) intron 1 variants are major risk factors for Graves' disease in three European Caucasian cohorts

Background

The thyroid stimulating hormone receptor (TSHR) gene is an established susceptibility locus for Graves'' disease (GD), with recent studies refining association to two single nucleotide polymorphisms (SNPs), rs179247 and rs12101255, within TSHR intron 1.

Methodology and Principal Findings

We aimed to validate association of rs179247 and rs12101255 in Polish and UK Caucasian GD case-control subjects, determine the mode of inheritance and to see if association correlates with specific GD clinical manifestations. We investigated three case-control populations; 558 GD patients and 520 controls from Warsaw, Poland, 196 GD patients and 198 controls from Gliwice, Poland and 2504 GD patients from the UK National collection and 2784 controls from the 1958 British Birth cohort. Both rs179247 (P = 1.2×10⁻²–6.2×10⁻¹⁵, OR = 1.38–1.45) and rs12101255 (P = 1.0×10⁻⁴–3.68×10⁻²¹, OR = 1.47–1.87) exhibited strong association with GD in all three cohorts. Logistic regression suggested association of rs179247 is secondary to rs12101255 in all cohorts. Inheritance modeling suggested a co-dominant mode of inheritance in all cohorts. Genotype-phenotype correlations provided no clear evidence of association with any specific clinical characteristics.

Conclusions

We have validated association of TSHR intron 1 SNPs with GD in three independent European cohorts and have demonstrated that the aetiological variant within the TSHR is likely to be in strong linkage disequilibrium with rs12101255. Fine mapping is now required to determine the exact location of the aetiological DNA variants within the TSHR.     

Genome-wide gene-environment study identifies glutamate receptor gene GRIN2A as a Parkinson's disease modifier gene via interaction with coffee

Our aim was to identify genes that influence the inverse association of coffee with the risk of developing Parkinson''s disease (PD). We used genome-wide genotype data and lifetime caffeinated-coffee-consumption data on 1,458 persons with PD and 931 without PD from the NeuroGenetics Research Consortium (NGRC), and we performed a genome-wide association and interaction study (GWAIS), testing each SNP''s main-effect plus its interaction with coffee, adjusting for sex, age, and two principal components. We then stratified subjects as heavy or light coffee-drinkers and performed genome-wide association study (GWAS) in each group. We replicated the most significant SNP. Finally, we imputed the NGRC dataset, increasing genomic coverage to examine the region of interest in detail. The primary analyses (GWAIS, GWAS, Replication) were performed using genotyped data. In GWAIS, the most significant signal came from rs4998386 and the neighboring SNPs in GRIN2A. GRIN2A encodes an NMDA-glutamate-receptor subunit and regulates excitatory neurotransmission in the brain. Achieving P_2df = 10⁻⁶, GRIN2A surpassed all known PD susceptibility genes in significance in the GWAIS. In stratified GWAS, the GRIN2A signal was present in heavy coffee-drinkers (OR = 0.43; P = 6×10⁻⁷) but not in light coffee-drinkers. The a priori Replication hypothesis that “Among heavy coffee-drinkers, rs4998386_T carriers have lower PD risk than rs4998386_CC carriers” was confirmed: OR_Replication = 0.59, P_Replication = 10⁻³; OR_Pooled = 0.51, P_Pooled = 7×10⁻⁸. Compared to light coffee-drinkers with rs4998386_CC genotype, heavy coffee-drinkers with rs4998386_CC genotype had 18% lower risk (P = 3×10⁻³), whereas heavy coffee-drinkers with rs4998386_TC genotype had 59% lower risk (P = 6×10⁻¹³). Imputation revealed a block of SNPs that achieved P_2df<5×10⁻⁸ in GWAIS, and OR = 0.41, P = 3×10⁻⁸ in heavy coffee-drinkers. This study is proof of concept that inclusion of environmental factors can help identify genes that are missed in GWAS. Both adenosine antagonists (caffeine-like) and glutamate antagonists (GRIN2A-related) are being tested in clinical trials for treatment of PD. GRIN2A may be a useful pharmacogenetic marker for subdividing individuals in clinical trials to determine which medications might work best for which patients.     

Quantile-specific penetrance of genes affecting lipoproteins, adiposity and height

Quantile-dependent penetrance is proposed to occur when the phenotypic expression of a SNP depends upon the population percentile of the phenotype. To illustrate the phenomenon, quantiles of height, body mass index (BMI), and plasma lipids and lipoproteins were compared to genetic risk scores (GRS) derived from single nucleotide polymorphisms (SNP)s having established genome-wide significance: 180 SNPs for height, 32 for BMI, 37 for low-density lipoprotein (LDL)-cholesterol, 47 for high-density lipoprotein (HDL)-cholesterol, 52 for total cholesterol, and 31 for triglycerides in 1930 subjects. Both phenotypes and GRSs were adjusted for sex, age, study, and smoking status. Quantile regression showed that the slope of the genotype-phenotype relationships increased with the percentile of BMI (P = 0.002), LDL-cholesterol (P = 3×10⁻⁸), HDL-cholesterol (P = 5×10⁻⁶), total cholesterol (P = 2.5×10⁻⁶), and triglyceride distribution (P = 7.5×10⁻⁶), but not height (P = 0.09). Compared to a GRS''s phenotypic effect at the 10^th population percentile, its effect at the 90^th percentile was 4.2-fold greater for BMI, 4.9-fold greater for LDL-cholesterol, 1.9-fold greater for HDL-cholesterol, 3.1-fold greater for total cholesterol, and 3.3-fold greater for triglycerides. Moreover, the effect of the rs1558902 (FTO) risk allele was 6.7-fold greater at the 90^th than the 10^th percentile of the BMI distribution, and that of the rs3764261 (CETP) risk allele was 2.4-fold greater at the 90^th than the 10^th percentile of the HDL-cholesterol distribution. Conceptually, it maybe useful to distinguish environmental effects on the phenotype that in turn alters a gene''s phenotypic expression (quantile-dependent penetrance) from environmental effects affecting the gene''s phenotypic expression directly (gene-environment interaction).     

Evidence for more than one Parkinson's disease-associated variant within the HLA region

Parkinson''s disease (PD) was recently found to be associated with HLA in a genome-wide association study (GWAS). Follow-up GWAS''s replicated the PD-HLA association but their top hits differ. Do the different hits tag the same locus or is there more than one PD-associated variant within HLA? We show that the top GWAS hits are not correlated with each other (0.00≤r²≤0.15). Using our GWAS (2000 cases, 1986 controls) we conducted step-wise conditional analysis on 107 SNPs with P<10⁻³ for PD-association; 103 dropped-out, four remained significant. Each SNP, when conditioned on the other three, yielded P_SNP1 = 5×10⁻⁴, P_SNP2 = 5×10⁻⁴, P_SNP3 = 4×10⁻³ and P_SNP4 = 0.025. The four SNPs were not correlated (0.01≤r²≤0.20). Haplotype analysis (excluding rare SNP2) revealed increasing PD risk with increasing risk alleles from OR = 1.27, P = 5×10⁻³ for one risk allele to OR = 1.65, P = 4×10⁻⁸ for three. Using additional 843 cases and 856 controls we replicated the independent effects of SNP1 (P_{conditioned-on-SNP4} = 0.04) and SNP4 (P_{conditioned-on-SNP1} = 0.04); SNP2 and SNP3 could not be replicated. In pooled GWAS and replication, SNP1 had OR_{conditioned-on-SNP4} = 1.23, P_{conditioned-on-SNP4} = 6×10⁻⁷; SNP4 had OR_{conditioned-on-SNP1} = 1.18, P_{conditioned-on-SNP1} = 3×10⁻³; and the haplotype with both risk alleles had OR = 1.48, P = 2×10⁻¹². Genotypic OR increased with the number of risk alleles an individual possessed up to OR = 1.94, P = 2×10⁻¹¹ for individuals who were homozygous for the risk allele at both SNP1 and SNP4. SNP1 is a variant in HLA-DRA and is associated with HLA-DRA, DRB5 and DQA2 gene expression. SNP4 is correlated (r² = 0.95) with variants that are associated with HLA-DQA2 expression, and with the top HLA SNP from the IPDGC GWAS (r² = 0.60). Our findings suggest more than one PD-HLA association; either different alleles of the same gene, or separate loci.     

A genome-wide association scan on the levels of markers of inflammation in Sardinians reveals associations that underpin its complex regulation

Identifying the genes that influence levels of pro-inflammatory molecules can help to elucidate the mechanisms underlying this process. We first conducted a two-stage genome-wide association scan (GWAS) for the key inflammatory biomarkers Interleukin-6 (IL-6), the general measure of inflammation erythrocyte sedimentation rate (ESR), monocyte chemotactic protein-1 (MCP-1), and high-sensitivity C-reactive protein (hsCRP) in a large cohort of individuals from the founder population of Sardinia. By analysing 731,213 autosomal or X chromosome SNPs and an additional ∼1.9 million imputed variants in 4,694 individuals, we identified several SNPs associated with the selected quantitative trait loci (QTLs) and replicated all the top signals in an independent sample of 1,392 individuals from the same population. Next, to increase power to detect and resolve associations, we further genotyped the whole cohort (6,145 individuals) for 293,875 variants included on the ImmunoChip and MetaboChip custom arrays. Overall, our combined approach led to the identification of 9 genome-wide significant novel independent signals—5 of which were identified only with the custom arrays—and provided confirmatory evidence for an additional 7. Novel signals include: for IL-6, in the ABO gene (rs657152, p = 2.13×10⁻²⁹); for ESR, at the HBB (rs4910472, p = 2.31×10⁻¹¹) and UCN119B/SPPL3 (rs11829037, p = 8.91×10⁻¹⁰) loci; for MCP-1, near its receptor CCR2 (rs17141006, p = 7.53×10⁻¹³) and in CADM3 (rs3026968, p = 7.63×10⁻¹³); for hsCRP, within the CRP gene (rs3093077, p = 5.73×10⁻²¹), near DARC (rs3845624, p = 1.43×10⁻¹⁰), UNC119B/SPPL3 (rs11829037, p = 1.50×10⁻¹⁴), and ICOSLG/AIRE (rs113459440, p = 1.54×10⁻⁰⁸) loci. Confirmatory evidence was found for IL-6 in the IL-6R gene (rs4129267); for ESR at CR1 (rs12567990) and TMEM57 (rs10903129); for MCP-1 at DARC (rs12075); and for hsCRP at CRP (rs1205), HNF1A (rs225918), and APOC-I (rs4420638). Our results improve the current knowledge of genetic variants underlying inflammation and provide novel clues for the understanding of the molecular mechanisms regulating this complex process.     

Genetic determinants for body iron store and type 2 diabetes risk in US men and women

Background

High body iron store has been associated with an increased risk of type 2 diabetes (T2D); it remains unknown whether the genetic variants related to body iron status affect T2D risk. We aimed at comprehensively investigating the associations between the genetic variants related to body iron status and the T2D risk.

Methodology/Principal Findings

Six common SNPs related to body iron status from recent genome-wide association (GWA) studies were determined in the Nurses’ Health Study (NHS; 1,467 diabetic cases and 1,754 controls) and the Health Professionals Follow-up Study (HPFS; 1,124, diabetic cases and 1,298 controls). Plasma levels of ferritin, soluble transferrin receptor (sTfR), and transferrin were measured in NHS. Significant associations were observed for loci in TPMRSS6 with sTfR (P = 3.47×10⁻⁶), TF with transferrin (P = 0.0002 to 1.72×10⁻¹⁰); and HFE with ferritin (P = 0.017 to 1.6×10⁻⁸), sTfR (P = 0.007 to 7.9×10⁻⁶), and transferrin (P = 0.006 to 0.0007). The six SNPs together explained 5.7%, 2.7%, and 13.3% of the variation in plasma levels of ferritin, sTfR, and transferrin. After adjustment for the conventional risk factors, the T allele of SNP rs855791 in the TPMRSS6 gene was significantly associated with a 19% decreased risk of T2D (OR = 0.81; 95% CI = 0.66–0.98; P = 0.03) in men. Multiple tests attenuated this significant association to null. No associations were observed in women. SNPs at HFE and TF were not associated with diabetes risk in either sex. Dietary iron intake did not modify the associations of the newly identified loci with diabetes risk.

Conclusions/Significance

The newly identified iron-related SNP rs855791 in TPMRSS6 was nominally associated with a decreased risk of T2D in men but not in women. The apparent differences by gender warrant further study.     

Phenotype restricted genome-wide association study using a gene-centric approach identifies three low-risk neuroblastoma susceptibility Loci

Neuroblastoma is a malignant neoplasm of the developing sympathetic nervous system that is notable for its phenotypic diversity. High-risk patients typically have widely disseminated disease at diagnosis and a poor survival probability, but low-risk patients frequently have localized tumors that are almost always cured with little or no chemotherapy. Our genome-wide association study (GWAS) has identified common variants within FLJ22536, BARD1, and LMO1 as significantly associated with neuroblastoma and more robustly associated with high-risk disease. Here we show that a GWAS focused on low-risk cases identified SNPs within DUSP12 at 1q23.3 (P = 2.07×10⁻⁶), DDX4 and IL31RA both at 5q11.2 (P = 2.94×10⁻⁶ and 6.54×10⁻⁷ respectively), and HSD17B12 at 11p11.2 (P = 4.20×10⁻⁷) as being associated with the less aggressive form of the disease. These data demonstrate the importance of robust phenotypic data in GWAS analyses and identify additional susceptibility variants for neuroblastoma.     

Stratifying type 2 diabetes cases by BMI identifies genetic risk variants in LAMA1 and enrichment for risk variants in lean compared to obese cases

Common diseases such as type 2 diabetes are phenotypically heterogeneous. Obesity is a major risk factor for type 2 diabetes, but patients vary appreciably in body mass index. We hypothesized that the genetic predisposition to the disease may be different in lean (BMI<25 Kg/m²) compared to obese cases (BMI≥30 Kg/m²). We performed two case-control genome-wide studies using two accepted cut-offs for defining individuals as overweight or obese. We used 2,112 lean type 2 diabetes cases (BMI<25 kg/m²) or 4,123 obese cases (BMI≥30 kg/m²), and 54,412 un-stratified controls. Replication was performed in 2,881 lean cases or 8,702 obese cases, and 18,957 un-stratified controls. To assess the effects of known signals, we tested the individual and combined effects of SNPs representing 36 type 2 diabetes loci. After combining data from discovery and replication datasets, we identified two signals not previously reported in Europeans. A variant (rs8090011) in the LAMA1 gene was associated with type 2 diabetes in lean cases (P = 8.4×10⁻⁹, OR = 1.13 [95% CI 1.09–1.18]), and this association was stronger than that in obese cases (P = 0.04, OR = 1.03 [95% CI 1.00–1.06]). A variant in HMG20A—previously identified in South Asians but not Europeans—was associated with type 2 diabetes in obese cases (P = 1.3×10⁻⁸, OR = 1.11 [95% CI 1.07–1.15]), although this association was not significantly stronger than that in lean cases (P = 0.02, OR = 1.09 [95% CI 1.02–1.17]). For 36 known type 2 diabetes loci, 29 had a larger odds ratio in the lean compared to obese (binomial P = 0.0002). In the lean analysis, we observed a weighted per-risk allele OR = 1.13 [95% CI 1.10–1.17], P = 3.2×10⁻¹⁴. This was larger than the same model fitted in the obese analysis where the OR = 1.06 [95% CI 1.05–1.08], P = 2.2×10⁻¹⁶. This study provides evidence that stratification of type 2 diabetes cases by BMI may help identify additional risk variants and that lean cases may have a stronger genetic predisposition to type 2 diabetes.     

A genome-wide association study of female sexual dysfunction

Background

Female sexual dysfunction (FSD) is an important but controversial problem with serious negative impact on women’s quality of life. Data from twin studies have shown a genetic contribution to the development and maintenance of FSD.

Methodology/Principal Findings

We performed a genome-wide association study (GWAS) on 2.5 million single-nucleotide polymorphisms (SNPs) in 1,104 female twins (25–81 years of age) in a population-based register and phenotypic data on lifelong sexual functioning. Although none reached conventional genome-wide level of significance (10×-8), we found strongly suggestive associations with the phenotypic dimension of arousal (rs13202860, P = 1.2×10⁻⁷; rs1876525, P = 1.2×10⁻⁷; and rs13209281 P = 8.3×10⁻⁷) on chromosome 6, around 500kb upstream of the locus HTR1E (5-hydroxytryptamine receptor 1E) locus, related to the serotonin brain pathways. We could not replicate previously reported candidate SNPs associated with FSD in the DRD4, 5HT2A and IL-1B loci.

Conclusions/Significance

We report the first GWAS of FSD symptoms in humans. This has pointed to several “risk alleles” and the implication of the serotonin and GABA pathways. Ultimately, understanding key mechanisms via this research may lead to new FSD treatments and inform clinical practice and developments in psychiatric nosology.     

Associations of HLA-DP variants with hepatitis B virus infection in southern and northern Han Chinese populations: a multicenter case-control study

Background

Human leukocyte antigen DP (HLA-DP) locus has been reported to be associated with hepatitis B virus (HBV) infection in populations of Japan and Thailand. We aimed to examine whether the association can be replicated in Han Chinese populations.

Methodology/Principal Findings

Two HLA-DP variants rs2395309 and rs9277535 (the most strongly associated SNPs from each HLA-DP locus) were genotyped in three independent Han cohorts consisting of 2 805 cases and 1 796 controls. By using logistic regression analysis, these two SNPs in the HLA-DPA1 and HLA-DPB1 genes were significantly associated with HBV infection in Han Chinese populations (P = 0.021∼3.36×10⁻⁸ at rs2395309; P = 8.37×10⁻³∼2.68×10⁻¹⁰ at rs9277535). In addition, the genotype distributions of both sites (rs2395309 and rs9277535) were clearly different between southern and northern Chinese population (P = 8.95×10⁻⁵ at rs2395309; P = 1.64×10⁻⁹ at rs9277535). By using asymptomatic HBV carrier as control group, our study showed that there were no associations of two HLA-DP variants with HBV progression (P = 0.305∼0.822 and 0.163∼0.881 in southern Chinese population, respectively; P = 0.097∼0.697 and 0.198∼0.615 in northern Chinese population, respectively).

Conclusions

Our results confirmed that two SNPs (rs2395309 and rs9277535) in the HLA-DP loci were strongly associated with HBV infection in southern and northern Han Chinese populations, but not with HBV progression.     

Identification of novel susceptibility Loci for kawasaki disease in a Han chinese population by a genome-wide association study

Kawasaki disease (KD) is an acute systemic vasculitis syndrome that primarily affects infants and young children. Its etiology is unknown; however, epidemiological findings suggest that genetic predisposition underlies disease susceptibility. Taiwan has the third-highest incidence of KD in the world, after Japan and Korea. To investigate novel mechanisms that might predispose individuals to KD, we conducted a genome-wide association study (GWAS) in 250 KD patients and 446 controls in a Han Chinese population residing in Taiwan, and further validated our findings in an independent Han Chinese cohort of 208 cases and 366 controls. The most strongly associated single-nucleotide polymorphisms (SNPs) detected in the joint analysis corresponded to three novel loci. Among these KD-associated SNPs three were close to the COPB2 (coatomer protein complex beta-2 subunit) gene: rs1873668 (p = 9.52×10⁻⁵), rs4243399 (p = 9.93×10⁻⁵), and rs16849083 (p = 9.93×10⁻⁵). We also identified a SNP in the intronic region of the ERAP1 (endoplasmic reticulum amino peptidase 1) gene (rs149481, p_best = 4.61×10⁻⁵). Six SNPs (rs17113284, rs8005468, rs10129255, rs2007467, rs10150241, and rs12590667) clustered in an area containing immunoglobulin heavy chain variable regions genes, with p_best-values between 2.08×10⁻⁵ and 8.93×10⁻⁶, were also identified. This is the first KD GWAS performed in a Han Chinese population. The novel KD candidates we identified have been implicated in T cell receptor signaling, regulation of proinflammatory cytokines, as well as antibody-mediated immune responses. These findings may lead to a better understanding of the underlying molecular pathogenesis of KD.     

The role of osteopontin (OPN/SPP1) haplotypes in the susceptibility to Crohn's disease

Background

Osteopontin represents a multifunctional molecule playing a pivotal role in chronic inflammatory and autoimmune diseases. Its expression is increased in inflammatory bowel disease (IBD). The aim of our study was to analyze the association of osteopontin (OPN/SPP1) gene variants in a large cohort of IBD patients.

Methodology/Principal Findings

Genomic DNA from 2819 Caucasian individuals (n = 841 patients with Crohn''s disease (CD), n = 473 patients with ulcerative colitis (UC), and n = 1505 healthy unrelated controls) was analyzed for nine OPN SNPs (rs2728127, rs2853744, rs11730582, rs11739060, rs28357094, rs4754 = p.Asp80Asp, rs1126616 = p.Ala236Ala, rs1126772 and rs9138). Considering the important role of osteopontin in Th17-mediated diseases, we performed analysis for epistasis with IBD-associated IL23R variants and analyzed serum levels of the Th17 cytokine IL-22. For four OPN SNPs (rs4754, rs1126616, rs1126772 and rs9138), we observed significantly different distributions between male and female CD patients. rs4754 was protective in male CD patients (p = 0.0004, OR = 0.69). None of the other investigated OPN SNPs was associated with CD or UC susceptibility. However, several OPN haplotypes showed significant associations with CD susceptibility. The strongest association was found for a haplotype consisting of the 8 OPN SNPs rs2728127-rs2853744-rs11730582-rs11439060-rs28357094-rs112661-rs1126772-rs9138 (omnibus p-value = 2.07×10⁻⁸). Overall, the mean IL-22 secretion in the combined group of OPN minor allele carriers with CD was significantly lower than that of CD patients with OPN wildtype alleles (p = 3.66×10⁻⁵). There was evidence for weak epistasis between the OPN SNP rs28357094 with the IL23R SNP rs10489629 (p = 4.18×10⁻²) and between OPN SNP rs1126616 and IL23R SNP rs2201841 (p = 4.18×10⁻²) but none of these associations remained significant after Bonferroni correction.

Conclusions/Significance

Our study identified OPN haplotypes as modifiers of CD susceptibility, while the combined effects of certain OPN variants may modulate IL-22 secretion.     

Caucasian and Asian specific rheumatoid arthritis risk loci reveal limited replication and apparent allelic heterogeneity in north Indians

Genome-wide association studies and meta-analysis indicate that several genes/loci are consistently associated with rheumatoid arthritis (RA) in European and Asian populations. To evaluate the transferability status of these findings to an ethnically diverse north Indian population, we performed a replication analysis. We investigated the association of 47 single-nucleotide polymorphisms (SNPs) at 43 of these genes/loci with RA in a north Indian cohort comprising 983 RA cases and 1007 age and gender matched controls. Genotyping was done using Infinium human 660w-quad. Association analysis by chi-square test implemented in plink was carried out in two steps. Firstly, association of the index or surrogate SNP (r2>0.8, calculated from reference GIH Hap-Map population) was tested. In the second step, evidence for allelic/locus heterogeneity at aforementioned genes/loci was assessed for by testing additional flanking SNPs in linkage equilibrium with index/surrogate marker.Of the 44 European specific index SNPs, neither index nor surrogate SNPs were present for nine SNPs in the genotyping array. Of the remaining 35, associations were replicated at seven genes namely PTPN22 (rs1217407, p = 3×10⁻³); IL2–21 (rs13119723, p = 0.008); HLA-DRB1 (rs660895, p = 2.56×10⁻⁵; rs6457617, p = 1.6×10⁻⁰⁹; rs13192471, p = 6.7×10⁻¹⁶); TNFA1P3 (rs9321637, p = 0.03); CCL21 (rs13293020, p = 0.01); IL2RA (rs2104286, p = 1.9×10⁻⁴) and ZEB1 (rs2793108, p = 0.006). Of the three Asian specific loci tested, rs2977227 in PADI4 showed modest association (p<0.02). Further, of the 140 SNPs (in LE with index/surrogate variant) tested, association was observed at 11 additional genes: PTPRC, AFF3, CD28, CTLA4, PXK, ANKRD55, TAGAP, CCR6, BLK, CD40 and IL2RB. This study indicates limited replication of European and Asian index SNPs and apparent allelic heterogeneity in RA etiology among north Indians warranting independent GWAS in this population. However, replicated associations of HLA-DRB1, PTPN22 (which confer ∼50% of the heritable risk to RA) and IL2RA suggest that cross-ethnicity fine mapping of such loci is apposite for identification of causal variants.     

Comparison of 6q25 Breast Cancer Hits from Asian and European Genome Wide Association Studies in the Breast Cancer Association Consortium (BCAC)

The 6q25.1 locus was first identified via a genome-wide association study (GWAS) in Chinese women and marked by single nucleotide polymorphism (SNP) rs2046210, approximately 180 Kb upstream of ESR1. There have been conflicting reports about the association of this locus with breast cancer in Europeans, and a GWAS in Europeans identified a different SNP, tagged here by rs12662670. We examined the associations of both SNPs in up to 61,689 cases and 58,822 controls from forty-four studies collaborating in the Breast Cancer Association Consortium, of which four studies were of Asian and 39 of European descent. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI). Case-only analyses were used to compare SNP effects in Estrogen Receptor positive (ER+) versus negative (ER−) tumours. Models including both SNPs were fitted to investigate whether the SNP effects were independent. Both SNPs are significantly associated with breast cancer risk in both ethnic groups. Per-allele ORs are higher in Asian than in European studies [rs2046210: OR (A/G) = 1.36 (95% CI 1.26–1.48), p = 7.6×10⁻¹⁴ in Asians and 1.09 (95% CI 1.07–1.11), p = 6.8×10⁻¹⁸ in Europeans. rs12662670: OR (G/T) = 1.29 (95% CI 1.19–1.41), p = 1.2×10⁻⁹ in Asians and 1.12 (95% CI 1.08–1.17), p = 3.8×10⁻⁹ in Europeans]. SNP rs2046210 is associated with a significantly greater risk of ER− than ER+ tumours in Europeans [OR (ER−) = 1.20 (95% CI 1.15–1.25), p = 1.8×10⁻¹⁷ versus OR (ER+) = 1.07 (95% CI 1.04–1.1), p = 1.3×10⁻⁷, p_{heterogeneity} = 5.1×10⁻⁶]. In these Asian studies, by contrast, there is no clear evidence of a differential association by tumour receptor status. Each SNP is associated with risk after adjustment for the other SNP. These results suggest the presence of two variants at 6q25.1 each independently associated with breast cancer risk in Asians and in Europeans. Of these two, the one tagged by rs2046210 is associated with a greater risk of ER− tumours.     

Effect of ABCB1 and ABCC3 polymorphisms on osteosarcoma survival after chemotherapy: a pharmacogenetic study

Background

Standard treatment for osteosarcoma patients consists of a combination of cisplatin, adriamycin, and methotrexate before surgical resection of the primary tumour, followed by postoperative chemotherapy including vincristine and cyclophosphamide. Unfortunately, many patients still relapse or suffer adverse events. We examined whether common germline polymorphisms in chemotherapeutic transporter and metabolic pathway genes of the drugs used in standard osteosarcoma treatment may predict treatment response.

Methodology/Principal Findings

In this study we screened 102 osteosarcoma patients for 346 Single Nucleotide Polymorphisms (SNPs) and 2 Copy Number Variants (CNVs) in 24 genes involved in the metabolism or transport of cisplatin, adriamycin, methotrexate, vincristine, and cyclophosphamide. We studied the association of the genotypes with tumour response and overall survival. We found that four SNPs in two ATP-binding cassette genes were significantly associated with overall survival: rs4148416 in ABCC3 (per-allele HR = 8.14, 95%CI = 2.73-20.2, p-value = 5.1×10⁻⁵), and three SNPs in ABCB1, rs4148737 (per-allele HR = 3.66, 95%CI = 1.85–6.11, p-value = 6.9×10⁻⁵), rs1128503 and rs10276036 (r² = 1, per-allele HR = 0.24, 95%CI = 0.11–0.47 p-value = 7.9×10⁻⁵). Associations with these SNPs remained statistically significant after correction for multiple testing (all corrected p-values [permutation test] ≤0.03).

Conclusions

Our findings suggest that these polymorphisms may affect osteosarcoma treatment efficacy. If these associations are independently validated, these variants could be used as genetic predictors of clinical outcome in the treatment of osteosarcoma, helping in the design of individualized therapy.     

GWAS identifies novel susceptibility loci on 6p21.32 and 21q21.3 for hepatocellular carcinoma in chronic hepatitis B virus carriers

Genome-wide association studies (GWAS) have recently identified KIF1B as susceptibility locus for hepatitis B virus (HBV)–related hepatocellular carcinoma (HCC). To further identify novel susceptibility loci associated with HBV–related HCC and replicate the previously reported association, we performed a large three-stage GWAS in the Han Chinese population. 523,663 autosomal SNPs in 1,538 HBV–positive HCC patients and 1,465 chronic HBV carriers were genotyped for the discovery stage. Top candidate SNPs were genotyped in the initial validation samples of 2,112 HBV–positive HCC cases and 2,208 HBV carriers and then in the second validation samples of 1,021 cases and 1,491 HBV carriers. We discovered two novel associations at rs9272105 (HLA-DQA1/DRB1) on 6p21.32 (OR = 1.30, P = 1.13×10⁻¹⁹) and rs455804 (GRIK1) on 21q21.3 (OR = 0.84, P = 1.86×10⁻⁸), which were further replicated in the fourth independent sample of 1,298 cases and 1,026 controls (rs9272105: OR = 1.25, P = 1.71×10⁻⁴; rs455804: OR = 0.84, P = 6.92×10⁻³). We also revealed the associations of HLA-DRB1*0405 and 0901*0602, which could partially account for the association at rs9272105. The association at rs455804 implicates GRIK1 as a novel susceptibility gene for HBV–related HCC, suggesting the involvement of glutamate signaling in the development of HBV–related HCC.