Genetic determinants of cardiovascular events among women with migraine: a genome-wide association study

Background

Migraine is associated with an increased risk for cardiovascular disease (CVD). Both migraine and CVD are highly heritable. However, the genetic liability for CVD among migraineurs is unclear.

Methods

We performed a genome-wide association study for incident CVD events during 12 years of follow-up among 5,122 migraineurs participating in the population-based Women''s Genome Health Study. Migraine was self-reported and CVD events were confirmed after medical records review. We calculated odds ratios (OR) and 95% confidence intervals (CI) and considered a genome-wide p-value <5×10⁻⁸ as significant.

Results

Among the 5,122 women with migraine 164 incident CVD events occurred during follow-up. No SNP was associated with major CVD, ischemic stroke, myocardial infarction, or CVD death at the genome-wide level; however, five SNPs showed association with p<5×10⁻⁶. Among migraineurs with aura rs7698623 in MEPE (OR = 6.37; 95% CI 3.15–12.90; p = 2.7×10⁻⁷) and rs4975709 in IRX4 (OR = 5.06; 95% CI 2.66–9.62; p = 7.7×10⁻⁷) appeared to be associated with ischemic stroke, rs2143678 located close to MDF1 with major CVD (OR = 3.05; 95% CI 1.98–4.69; p = 4.3×10⁻⁷), and the intergenic rs1406961 with CVD death (OR = 12.33; 95% CI 4.62–32.87; p = 5.2×10⁻⁷). Further, rs1047964 in BACE1 appeared to be associated with CVD death among women with any migraine (OR = 4.67; 95% CI 2.53–8.62; p = 8.0×10⁻⁷).

Conclusion

Our results provide some suggestion for an association of five SNPs with CVD events among women with migraine; none of the results was genome-wide significant. Four associations appeared among migraineurs with aura, two of those with ischemic stroke. Although our population is among the largest with migraine and incident CVD information, these results must be treated with caution, given the limited number of CVD events among women with migraine and the low minor allele frequencies for three of the SNPs. Our results await independent replication and should be considered hypothesis generating for future research.     

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.     

A genome-wide scan of Ashkenazi Jewish Crohn's disease suggests novel susceptibility loci

Crohn''s disease (CD) is a complex disorder resulting from the interaction of intestinal microbiota with the host immune system in genetically susceptible individuals. The largest meta-analysis of genome-wide association to date identified 71 CD–susceptibility loci in individuals of European ancestry. An important epidemiological feature of CD is that it is 2–4 times more prevalent among individuals of Ashkenazi Jewish (AJ) descent compared to non-Jewish Europeans (NJ). To explore genetic variation associated with CD in AJs, we conducted a genome-wide association study (GWAS) by combining raw genotype data across 10 AJ cohorts consisting of 907 cases and 2,345 controls in the discovery stage, followed up by a replication study in 971 cases and 2,124 controls. We confirmed genome-wide significant associations of 9 known CD loci in AJs and replicated 3 additional loci with strong signal (p<5×10⁻⁶). Novel signals detected among AJs were mapped to chromosomes 5q21.1 (rs7705924, combined p = 2×10⁻⁸; combined odds ratio OR = 1.48), 2p15 (rs6545946, p = 7×10⁻⁹; OR = 1.16), 8q21.11 (rs12677663, p = 2×10⁻⁸; OR = 1.15), 10q26.3 (rs10734105, p = 3×10⁻⁸; OR = 1.27), and 11q12.1 (rs11229030, p = 8×10⁻⁹; OR = 1.15), implicating biologically plausible candidate genes, including RPL7, CPAMD8, PRG2, and PRG3. In all, the 16 replicated and newly discovered loci, in addition to the three coding NOD2 variants, accounted for 11.2% of the total genetic variance for CD risk in the AJ population. This study demonstrates the complementary value of genetic studies in the Ashkenazim.     

Common variants at 9p21 and 8q22 are associated with increased susceptibility to optic nerve degeneration in glaucoma

Optic nerve degeneration caused by glaucoma is a leading cause of blindness worldwide. Patients affected by the normal-pressure form of glaucoma are more likely to harbor risk alleles for glaucoma-related optic nerve disease. We have performed a meta-analysis of two independent genome-wide association studies for primary open angle glaucoma (POAG) followed by a normal-pressure glaucoma (NPG, defined by intraocular pressure (IOP) less than 22 mmHg) subgroup analysis. The single-nucleotide polymorphisms that showed the most significant associations were tested for association with a second form of glaucoma, exfoliation-syndrome glaucoma. The overall meta-analysis of the GLAUGEN and NEIGHBOR dataset results (3,146 cases and 3,487 controls) identified significant associations between two loci and POAG: the CDKN2BAS region on 9p21 (rs2157719 [G], OR = 0.69 [95%CI 0.63–0.75], p = 1.86×10⁻¹⁸), and the SIX1/SIX6 region on chromosome 14q23 (rs10483727 [A], OR = 1.32 [95%CI 1.21–1.43], p = 3.87×10⁻¹¹). In sub-group analysis two loci were significantly associated with NPG: 9p21 containing the CDKN2BAS gene (rs2157719 [G], OR = 0.58 [95% CI 0.50–0.67], p = 1.17×10⁻¹²) and a probable regulatory region on 8q22 (rs284489 [G], OR = 0.62 [95% CI 0.53–0.72], p = 8.88×10⁻¹⁰). Both NPG loci were also nominally associated with a second type of glaucoma, exfoliation syndrome glaucoma (rs2157719 [G], OR = 0.59 [95% CI 0.41–0.87], p = 0.004 and rs284489 [G], OR = 0.76 [95% CI 0.54–1.06], p = 0.021), suggesting that these loci might contribute more generally to optic nerve degeneration in glaucoma. Because both loci influence transforming growth factor beta (TGF-beta) signaling, we performed a genomic pathway analysis that showed an association between the TGF-beta pathway and NPG (permuted p = 0.009). These results suggest that neuro-protective therapies targeting TGF-beta signaling could be effective for multiple forms of glaucoma.     

Association of variants at 1q32 and STAT3 with ankylosing spondylitis suggests genetic overlap with Crohn's disease

Ankylosing spondylitis (AS) is a common inflammatory arthritic condition. Overt inflammatory bowel disease (IBD) occurs in about 10% of AS patients, and in addition 70% of AS cases may have subclinical terminal ileitis. Spondyloarthritis is also common in IBD patients. We therefore tested Crohn''s disease susceptibility genes for association with AS, aiming to identify pleiotropic genetic associations with both diseases. Genotyping was carried out using Sequenom and Applied Biosystems TaqMan and OpenArray technologies on 53 markers selected from 30 Crohn''s disease associated genomic regions. We tested genotypes in a population of unrelated individual cases (n = 2,773) and controls (n = 2,215) of white European ancestry for association with AS. Statistical analysis was carried out using a Cochran-Armitage test for trend in PLINK. Strong association was detected at chr1q32 near KIF21B (rs11584383, P = 1.6×10⁻¹⁰, odds ratio (OR) = 0.74, 95% CI:0.68–0.82). Association with disease was also detected for 2 variants within STAT3 (rs6503695, P = 4.6×10⁻⁴. OR = 0.86 (95% CI:0.79–0.93); rs744166, P = 2.6×10⁻⁵, OR = 0.84 (95% CI:0.77–0.91)). Association was confirmed for IL23R (rs11465804, P = 1.2×10⁻⁵, OR = 0.65 (95% CI:0.54–0.79)), and further associations were detected for IL12B (rs10045431, P = 5.2×10⁻⁵, OR = 0.83 (95% CI:0.76–0.91)), CDKAL1 (rs6908425, P = 1.1×10⁻⁴, OR = 0.82 (95% CI:0.74–0.91)), LRRK2/MUC19 (rs11175593, P = 9.9×10⁻⁵, OR = 1.92 (95% CI: 1.38–2.67)), and chr13q14 (rs3764147, P = 5.9×10⁻⁴, OR = 1.19 (95% CI: 1.08–1.31)). Excluding cases with clinical IBD did not significantly affect these findings. This study identifies chr1q32 and STAT3 as ankylosing spondylitis susceptibility loci. It also further confirms association for IL23R and detects suggestive association with another 4 loci. STAT3 is a key signaling molecule within the Th17 lymphocyte differentiation pathway and further enhances the case for a major role of this T-lymphocyte subset in ankylosing spondylitis. Finally these findings suggest common aetiopathogenic pathways for AS and Crohn''s disease and further highlight the involvement of common risk variants across multiple diseases.     

Genetic association study identifies HSPB7 as a risk gene for idiopathic dilated cardiomyopathy

Dilated cardiomyopathy (DCM) is a structural heart disease with strong genetic background. Monogenic forms of DCM are observed in families with mutations located mostly in genes encoding structural and sarcomeric proteins. However, strong evidence suggests that genetic factors also affect the susceptibility to idiopathic DCM. To identify risk alleles for non-familial forms of DCM, we carried out a case-control association study, genotyping 664 DCM cases and 1,874 population-based healthy controls from Germany using a 50K human cardiovascular disease bead chip covering more than 2,000 genes pre-selected for cardiovascular relevance. After quality control, 30,920 single nucleotide polymorphisms (SNP) were tested for association with the disease by logistic regression adjusted for gender, and results were genomic-control corrected. The analysis revealed a significant association between a SNP in HSPB7 gene (rs1739843, minor allele frequency 39%) and idiopathic DCM (p = 1.06×10⁻⁶, OR = 0.67 [95% CI 0.57–0.79] for the minor allele T). Three more SNPs showed p < 2.21×10⁻⁵. De novo genotyping of these four SNPs was done in three independent case-control studies of idiopathic DCM. Association between SNP rs1739843 and DCM was significant in all replication samples: Germany (n = 564, n = 981 controls, p = 2.07×10⁻³, OR = 0.79 [95% CI 0.67–0.92]), France 1 (n = 433 cases, n = 395 controls, p = 3.73×10⁻³, OR = 0.74 [95% CI 0.60–0.91]), and France 2 (n = 249 cases, n = 380 controls, p = 2.26×10⁻⁴, OR = 0.63 [95% CI 0.50–0.81]). The combined analysis of all four studies including a total of n = 1,910 cases and n = 3,630 controls showed highly significant evidence for association between rs1739843 and idiopathic DCM (p = 5.28×10⁻¹³, OR = 0.72 [95% CI 0.65–0.78]). None of the other three SNPs showed significant results in the replication stage.This finding of the HSPB7 gene from a genetic search for idiopathic DCM using a large SNP panel underscores the influence of common polymorphisms on DCM susceptibility.     

A genome-wide association study of nephrolithiasis in the Japanese population identifies novel susceptible Loci at 5q35.3, 7p14.3, and 13q14.1

Nephrolithiasis is a common nephrologic disorder with complex etiology. To identify the genetic factor(s) for nephrolithiasis, we conducted a three-stage genome-wide association study (GWAS) using a total of 5,892 nephrolithiasis cases and 17,809 controls of Japanese origin. Here we found three novel loci for nephrolithiasis: RGS14-SLC34A1-PFN3-F12 on 5q35.3 (rs11746443; P = 8.51×10⁻¹², odds ratio (OR) = 1.19), INMT-FAM188B-AQP1 on 7p14.3 (rs1000597; P = 2.16×10⁻¹⁴, OR = 1.22), and DGKH on 13q14.1 (rs4142110; P = 4.62×10⁻⁹, OR = 1.14). Subsequent analyses in 21,842 Japanese subjects revealed the association of SNP rs11746443 with the reduction of estimated glomerular filtration rate (eGFR) (P = 6.54×10⁻⁸), suggesting a crucial role for this variation in renal function. Our findings elucidated the significance of genetic variations for the pathogenesis of nephrolithiasis.     

Variation in MSRA modifies risk of neonatal intestinal obstruction in cystic fibrosis

Meconium ileus (MI), a life-threatening intestinal obstruction due to meconium with abnormal protein content, occurs in approximately 15 percent of neonates with cystic fibrosis (CF). Analysis of twins with CF demonstrates that MI is a highly heritable trait, indicating that genetic modifiers are largely responsible for this complication. Here, we performed regional family-based association analysis of a locus that had previously been linked to MI and found that SNP haplotypes 5′ to and within the MSRA gene were associated with MI (P = 1.99×10⁻⁵ to 1.08×10⁻⁶; Bonferroni P = 0.057 to 3.1×10⁻³). The haplotype with the lowest P value showed association with MI in an independent sample of 1,335 unrelated CF patients (OR = 0.72, 95% CI [0.53–0.98], P = 0.04). Intestinal obstruction at the time of weaning was decreased in CF mice with Msra null alleles compared to those with wild-type Msra resulting in significant improvement in survival (P = 1.2×10⁻⁴). Similar levels of goblet cell hyperplasia were observed in the ilea of the Cftr ^−/− and Cftr ^−/− Msra ^−/− mice. Modulation of MSRA, an antioxidant shown to preserve the activity of enzymes, may influence proteolysis in the developing intestine of the CF fetus, thereby altering the incidence of obstruction in the newborn period. Identification of MSRA as a modifier of MI provides new insight into the biologic mechanism of neonatal intestinal obstruction caused by loss of CFTR function.     

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.     

Genome-wide association study identifies HLA-DP as a susceptibility gene for pediatric asthma in Asian populations

Asthma is a complex phenotype influenced by genetic and environmental factors. We conducted a genome-wide association study (GWAS) with 938 Japanese pediatric asthma patients and 2,376 controls. Single-nucleotide polymorphisms (SNPs) showing strong associations (P<1×10⁻⁸) in GWAS were further genotyped in an independent Japanese samples (818 cases and 1,032 controls) and in Korean samples (835 cases and 421 controls). SNP rs987870, located between HLA-DPA1 and HLA-DPB1, was consistently associated with pediatric asthma in 3 independent populations (P _combined = 2.3×10⁻¹⁰, odds ratio [OR] = 1.40). HLA-DP allele analysis showed that DPA1*0201 and DPB1*0901, which were in strong linkage disequilibrium, were strongly associated with pediatric asthma (DPA1*0201: P = 5.5×10⁻¹⁰, OR = 1.52, and DPB1*0901: P = 2.0×10⁻⁷, OR = 1.49). Our findings show that genetic variants in the HLA-DP locus are associated with the risk of pediatric asthma in Asian populations.     

ABCC5, a Gene That Influences the Anterior Chamber Depth,Is Associated with Primary Angle Closure Glaucoma

Anterior chamber depth (ACD) is a key anatomical risk factor for primary angle closure glaucoma (PACG). We conducted a genome-wide association study (GWAS) on ACD to discover novel genes for PACG on a total of 5,308 population-based individuals of Asian descent. Genome-wide significant association was observed at a sequence variant within ABCC5 (rs1401999; per-allele effect size = −0.045 mm, P = 8.17×10⁻⁹). This locus was associated with an increase in risk of PACG in a separate case-control study of 4,276 PACG cases and 18,801 controls (per-allele OR = 1.13 [95% CI: 1.06–1.22], P = 0.00046). The association was strengthened when a sub-group of controls with open angles were included in the analysis (per-allele OR = 1.30, P = 7.45×10⁻⁹; 3,458 cases vs. 3,831 controls). Our findings suggest that the increase in PACG risk could in part be mediated by genetic sequence variants influencing anterior chamber dimensions.     

Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals

Circulating levels of adiponectin, a hormone produced predominantly by adipocytes, are highly heritable and are inversely associated with type 2 diabetes mellitus (T2D) and other metabolic traits. We conducted a meta-analysis of genome-wide association studies in 39,883 individuals of European ancestry to identify genes associated with metabolic disease. We identified 8 novel loci associated with adiponectin levels and confirmed 2 previously reported loci (P = 4.5×10⁻⁸–1.2×10⁻⁴³). Using a novel method to combine data across ethnicities (N = 4,232 African Americans, N = 1,776 Asians, and N = 29,347 Europeans), we identified two additional novel loci. Expression analyses of 436 human adipocyte samples revealed that mRNA levels of 18 genes at candidate regions were associated with adiponectin concentrations after accounting for multiple testing (p<3×10⁻⁴). We next developed a multi-SNP genotypic risk score to test the association of adiponectin decreasing risk alleles on metabolic traits and diseases using consortia-level meta-analytic data. This risk score was associated with increased risk of T2D (p = 4.3×10⁻³, n = 22,044), increased triglycerides (p = 2.6×10⁻¹⁴, n = 93,440), increased waist-to-hip ratio (p = 1.8×10⁻⁵, n = 77,167), increased glucose two hours post oral glucose tolerance testing (p = 4.4×10⁻³, n = 15,234), increased fasting insulin (p = 0.015, n = 48,238), but with lower in HDL-cholesterol concentrations (p = 4.5×10⁻¹³, n = 96,748) and decreased BMI (p = 1.4×10⁻⁴, n = 121,335). These findings identify novel genetic determinants of adiponectin levels, which, taken together, influence risk of T2D and markers of insulin resistance.     

Genetic loci associated with plasma phospholipid n-3 fatty acids: a meta-analysis of genome-wide association studies from the CHARGE Consortium

Long-chain n-3 polyunsaturated fatty acids (PUFAs) can derive from diet or from α-linolenic acid (ALA) by elongation and desaturation. We investigated the association of common genetic variation with plasma phospholipid levels of the four major n-3 PUFAs by performing genome-wide association studies in five population-based cohorts comprising 8,866 subjects of European ancestry. Minor alleles of SNPs in FADS1 and FADS2 (desaturases) were associated with higher levels of ALA (p = 3×10⁻⁶⁴) and lower levels of eicosapentaenoic acid (EPA, p = 5×10⁻⁵⁸) and docosapentaenoic acid (DPA, p = 4×10⁻¹⁵⁴). Minor alleles of SNPs in ELOVL2 (elongase) were associated with higher EPA (p = 2×10⁻¹²) and DPA (p = 1×10⁻⁴³) and lower docosahexaenoic acid (DHA, p = 1×10⁻¹⁵). In addition to genes in the n-3 pathway, we identified a novel association of DPA with several SNPs in GCKR (glucokinase regulator, p = 1×10⁻⁸). We observed a weaker association between ALA and EPA among carriers of the minor allele of a representative SNP in FADS2 (rs1535), suggesting a lower rate of ALA-to-EPA conversion in these subjects. In samples of African, Chinese, and Hispanic ancestry, associations of n-3 PUFAs were similar with a representative SNP in FADS1 but less consistent with a representative SNP in ELOVL2. Our findings show that common variation in n-3 metabolic pathway genes and in GCKR influences plasma phospholipid levels of n-3 PUFAs in populations of European ancestry and, for FADS1, in other ancestries.     

Genome-wide scan identifies TNIP1, PSORS1C1, and RHOB as novel risk loci for systemic sclerosis

Systemic sclerosis (SSc) is an orphan, complex, inflammatory disease affecting the immune system and connective tissue. SSc stands out as a severely incapacitating and life-threatening inflammatory rheumatic disease, with a largely unknown pathogenesis. We have designed a two-stage genome-wide association study of SSc using case-control samples from France, Italy, Germany, and Northern Europe. The initial genome-wide scan was conducted in a French post quality-control sample of 564 cases and 1,776 controls, using almost 500 K SNPs. Two SNPs from the MHC region, together with the 6 loci outside MHC having at least one SNP with a P<10⁻⁵ were selected for follow-up analysis. These markers were genotyped in a post-QC replication sample of 1,682 SSc cases and 3,926 controls. The three top SNPs are in strong linkage disequilibrium and located on 6p21, in the HLA-DQB1 gene: rs9275224, P = 9.18×10⁻⁸, OR = 0.69, 95% CI [0.60–0.79]; rs6457617, P = 1.14×10⁻⁷ and rs9275245, P = 1.39×10⁻⁷. Within the MHC region, the next most associated SNP (rs3130573, P = 1.86×10⁻⁵, OR = 1.36 [1.18–1.56]) is located in the PSORS1C1 gene. Outside the MHC region, our GWAS analysis revealed 7 top SNPs (P<10⁻⁵) that spanned 6 independent genomic regions. Follow-up of the 17 top SNPs in an independent sample of 1,682 SSc and 3,926 controls showed associations at PSORS1C1 (overall P = 5.70×10⁻¹⁰, OR:1.25), TNIP1 (P = 4.68×10⁻⁹, OR:1.31), and RHOB loci (P = 3.17×10⁻⁶, OR:1.21). Because of its biological relevance, and previous reports of genetic association at this locus with connective tissue disorders, we investigated TNIP1 expression. A markedly reduced expression of the TNIP1 gene and also its protein product were observed both in lesional skin tissue and in cultured dermal fibroblasts from SSc patients. Furthermore, TNIP1 showed in vitro inhibitory effects on inflammatory cytokine-induced collagen production. The genetic signal of association with TNIP1 variants, together with tissular and cellular investigations, suggests that this pathway has a critical role in regulating autoimmunity and SSc pathogenesis.     

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.     

JC Polyomavirus Infection Is Strongly Controlled by Human Leucocyte Antigen Class II Variants

JC polyomavirus (JCV) carriers with a compromised immune system, such as in HIV, or subjects on immune-modulating therapies, such as anti VLA-4 therapy may develop progressive multifocal leukoencephalopathy (PML) which is a lytic infection of oligodendrocytes in the brain. Serum antibodies to JCV mark infection occur only in 50–60% of infected individuals, and high JCV-antibody titers seem to increase the risk of developing PML. We here investigated the role of human leukocyte antigen (HLA), instrumental in immune defense in JCV antibody response. Anti-JCV antibody status, as a surrogate for JCV infection, were compared to HLA class I and II alleles in 1621 Scandinavian persons with MS and 1064 population-based Swedish controls and associations were replicated in 718 German persons with MS. HLA-alleles were determined by SNP imputation, sequence specific (SSP) kits and a reverse PCR sequence-specific oligonucleotide (PCR-SSO) method. An initial GWAS screen displayed a strong HLA class II region signal. The HLA-DRB1*15 haplotype was strongly negatively associated to JCV sero-status in Scandinavian MS cases (OR = 0.42, p = 7×10⁻¹⁵) and controls (OR = 0.53, p = 2×10⁻⁵). In contrast, the DQB1*06:03 haplotype was positively associated with JCV sero-status, in Scandinavian MS cases (OR = 1.63, p = 0.006), and controls (OR = 2.69, p = 1×10⁻⁵). The German dataset confirmed these findings (OR = 0.54, p = 1×10⁻⁴ and OR = 1.58, p = 0.03 respectively for these haplotypes). HLA class II restricted immune responses, and hence CD4+ T cell immunity is pivotal for JCV infection control. Alleles within the HLA-DR1*15 haplotype are associated with a protective effect on JCV infection. Alleles within the DQB1*06:03 haplotype show an opposite association. These associations between JC virus antibody response and human leucocyte antigens supports the notion that CD4+ T cells are crucial in the immune defence to JCV and lays the ground for risk stratification for PML and development of therapy and prevention.     

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.     

Replication of putative susceptibility loci from genome-wide association studies associated with coronary atherosclerosis in Chinese Han population

Background

Coronary atherosclerosis, the main cause of cardiovascular disease, is a progressive disease. Recent Genome Wide Association Studies (GWASs) discovered several novel loci associated with coronary artery disease (CAD) or its main complication myocardial infarction (MI). In this study, we investigated the associations between previously reported CAD- and MI-associated variants and coronary atherosclerosis in Chinese Han population.

Methodology/Principal Findings

We performed a case-control association study with 2,335 coronary atherosclerosis patients and 1,078 controls undergoing coronary angiography of Chinese Han from China. Fourteen single nucleotide polymorphisms (SNPs), located at 1p13.3, 1q41, 2q36.3, 6q25.1, 9p21.3, 10q11.21 and 15q22.33, were genotyped in our sample collection. Six SNPs at 9p21 were associated with coronary atherosclerosis susceptibility (P_trend<0.05) and rs10757274 showed the most significant association (P = 2.38×10⁻⁰⁸, OR = 1.34). These associations remained significant after adjustment for multiple comparisons. Rs17465637 at 1q41 (P_trend = 6.83×10⁻⁰³, OR = 0.86) also showed significant association with coronary atherosclerosis, but the association was not significant after multiple comparisons. Additionally, rs501120 (P = 8.36×10⁻⁰³, OR = 0.80) at 10q11.21 was associated with coronary atherosclerosis in females, but did not show association in males and all participants. Variants at 1p13.3, 2q36.3, 6q25.1 and 15q22.33 showed no associations with coronary atherosclerosis and main cardiovascular risk factors in our data.

Conclusions/Significance

Our findings indicated variants at 9p21 were significantly associated with coronary atherosclerosis in Han Chinese. Variants at 1q41 showed suggestive evidence of association and variants at 10q11.21 showed suggestive evidence of association in females, which warrant further study in a larger sample.     

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.