Role of Nicotine Dependence on the Relationship between Variants in the Nicotinic Receptor Genes and Risk of Lung Adenocarcinoma

Several variations in the nicotinic receptor genes have been identified to be associated with both lung cancer risk and smoking in the genome-wide association (GWA) studies. However, the relationships among these three factors (genetic variants, nicotine dependence, and lung cancer) remain unclear. In an attempt to elucidate these relationships, we applied mediation analysis to quantify the impact of nicotine dependence on the association between the nicotinic receptor genetic variants and lung adenocarcinoma risk. We evaluated 23 single nucleotide polymorphisms (SNPs) in the five nicotinic receptor related genes (CHRNB3, CHRNA6, and CHRNA5/A3/B4) previously reported to be associated with lung cancer risk and smoking behavior and 14 SNPs in the four ‘control’ genes (TERT, CLPTM1L, CYP1A1, and TP53), which were not reported in the smoking GWA studies. A total of 661 lung adenocarcinoma cases and 1,347 controls with a smoking history, obtained from the Environment and Genetics in Lung Cancer Etiology case-control study, were included in the study. Results show that nicotine dependence is a mediator of the association between lung adenocarcinoma and gene variations in the regions of CHRNA5/A3/B4 and accounts for approximately 15% of this relationship. The top two CHRNA3 SNPs associated with the risk for lung adenocarcinoma were rs1051730 and rs12914385 (p-value = 1.9×10⁻¹⁰ and 1.1×10⁻¹⁰, respectively). Also, these two SNPs had significant indirect effects on lung adenocarcinoma risk through nicotine dependence (p = 0.003 and 0.007). Gene variations rs2736100 and rs2853676 in TERT and rs401681 and rs31489 in CLPTM1L had significant direct associations on lung adenocarcinoma without indirect effects through nicotine dependence. Our findings suggest that nicotine dependence plays an important role between genetic variants in the CHRNA5/A3/B4 region, especially CHRNA3, and lung adenocarcinoma. This may provide valuable information for understanding the pathogenesis of lung adenocarcinoma and for conducting personalized smoking cessation interventions.     

Genome-Wide Investigation of DNA Methylation Marks Associated with FV Leiden Mutation

In order to investigate whether DNA methylation marks could contribute to the incomplete penetrance of the FV Leiden mutation, a major genetic risk factor for venous thrombosis (VT), we measured genome-wide DNA methylation levels in peripheral blood samples of 98 VT patients carrying the mutation and 251 VT patients without the mutation using the dedicated Illumina HumanMethylation450 array. The genome-wide analysis of 388,120 CpG probes identified three sites mapping to the SLC19A2 locus whose DNA methylation levels differed significantly (p<3 10⁻⁸) between carriers and non-carriers. The three sites replicated (p<2 10⁻⁷) in an independent sample of 214 individuals from five large families ascertained on VT and FV Leiden mutation among which 53 were carriers and 161 were non-carriers of the mutation. In both studies, these three CpG sites were also associated (2.33 10⁻¹¹<p<3.02 10⁻⁴) with biomarkers of the Protein C pathway known to be influenced by the FV Leiden mutation. A comprehensive linkage disequilibrium (LD) analysis of the whole locus revealed that the original associations were due to LD between the FV Leiden mutation and a block of single nucleotide polymorphisms (SNP) located in SLC19A2. After adjusting for this block of SNPs, the FV Leiden mutation was no longer associated with any CpG site (p>0.05). In conclusion, our work clearly illustrates some promises and pitfalls of DNA methylation investigations on peripheral blood DNA in large epidemiological cohorts. DNA methylation levels at SLC19A2 are influenced by SNPs in LD with FV Leiden, but these DNA methylation marks do not explain the incomplete penetrance of the FV Leiden mutation.     

Identification of Novel Genetic Risk Loci in Maltese Dogs with Necrotizing Meningoencephalitis and Evidence of a Shared Genetic Risk across Toy Dog Breeds

Necrotizing meningoencephalitis (NME) affects toy and small breed dogs causing progressive, often fatal, inflammation and necrosis in the brain. Genetic risk loci for NME previously were identified in pug dogs, particularly associated with the dog leukocyte antigen (DLA) class II complex on chromosome 12, but have not been investigated in other susceptible breeds. We sought to evaluate Maltese and Chihuahua dogs, in addition to pug dogs, to identify novel or shared genetic risk factors for NME development. Genome-wide association testing of single nucleotide polymorphisms (SNPs) in Maltese dogs with NME identified 2 regions of genome-wide significance on chromosomes 4 (chr4:74522353T>A, p = 8.1×10⁻⁷) and 15 (chr15:53338796A>G, p = 1.5×10⁻⁷). Haplotype analysis and fine-mapping suggests that ILR7 and FBXW7, respectively, both important for regulation of immune system function, could be the underlying associated genes. Further evaluation of these regions and the previously identified DLA II locus across all three breeds, revealed an enrichment of nominal significant SNPs associated with chromosome 15 in pug dogs and DLA II in Maltese and Chihuahua dogs. Meta-analysis confirmed effect sizes the same direction in all three breeds for both the chromosome 15 and DLA II loci (p = 8.6×10^–11 and p = 2.5×10⁻⁷, respectively). This suggests a shared genetic background exists between all breeds and confers susceptibility to NME, but effect sizes might be different among breeds. In conclusion, we identified the first genetic risk factors for NME development in the Maltese, chromosome 4 and chromosome 15, and provide evidence for a shared genetic risk between breeds associated with chromosome 15 and DLA II. Last, DLA II and IL7R both have been implicated in human inflammatory diseases of the central nervous system such as multiple sclerosis, suggesting that similar pharmacotherapeutic targets across species should be investigated.     

Fetal DNA Methylation Associates with Early Spontaneous Preterm Birth and Gestational Age

Spontaneous preterm birth (PTB, <37 weeks gestation) is a major public health concern, and children born preterm have a higher risk of morbidity and mortality throughout their lives. Recent studies suggest that fetal DNA methylation of several genes varies across a range of gestational ages (GA), but it is not yet clear if fetal epigenetic changes associate with PTB. The objective of this study is to interrogate methylation patterns across the genome in fetal leukocyte DNA from African Americans with early PTB (24^1/7–34^0/7 weeks; N = 22) or term births (39^0/7–40^6/7weeks; N = 28) and to evaluate the association of each CpG site with PTB and GA. DNA methylation was assessed across the genome with the HumanMethylation450 BeadChip. For each individual sample and CpG site, the proportion of DNA methylation was estimated. The associations between methylation and PTB or GA were evaluated by fitting a separate linear model for each CpG site, adjusting for relevant covariates. Overall, 29 CpG sites associated with PTB (FDR<.05; 5.7×10⁻¹⁰<p<2.9×10⁻⁶) independent of GA. Also, 9637 sites associated with GA (FDR<.05; 9.5×10⁻¹⁶<p<1.0×10⁻³), with 61.8% decreasing in methylation with shorter GA. GA-associated CpG sites were depleted in the CpG islands of their respective genes (p<2.2×10⁻¹⁶). Gene set enrichment analysis (GSEA) supported enrichment of GA-associated CpG sites in genes that play a role in embryonic development as well as the extracellular matrix. Additionally, this study replicated the association of several CpG sites associated with gestational age in other studies (CRHBP, PIK3CD and AVP). Dramatic differences in fetal DNA methylation are evident in fetuses born preterm versus at term, and the patterns established at birth may provide insight into the long-term consequences associated with PTB.     

Genome-Wide Association of Body Fat Distribution in African Ancestry Populations Suggests New Loci

Central obesity, measured by waist circumference (WC) or waist-hip ratio (WHR), is a marker of body fat distribution. Although obesity disproportionately affects minority populations, few studies have conducted genome-wide association study (GWAS) of fat distribution among those of predominantly African ancestry (AA). We performed GWAS of WC and WHR, adjusted and unadjusted for BMI, in up to 33,591 and 27,350 AA individuals, respectively. We identified loci associated with fat distribution in AA individuals using meta-analyses of GWA results for WC and WHR (stage 1). Overall, 25 SNPs with single genomic control (GC)-corrected p-values<5.0×10⁻⁶ were followed-up (stage 2) in AA with WC and with WHR. Additionally, we interrogated genomic regions of previously identified European ancestry (EA) WHR loci among AA. In joint analysis of association results including both Stage 1 and 2 cohorts, 2 SNPs demonstrated association, rs2075064 at LHX2, p = 2.24×10⁻⁸ for WC-adjusted-for-BMI, and rs6931262 at RREB1, p = 2.48×10⁻⁸ for WHR-adjusted-for-BMI. However, neither signal was genome-wide significant after double GC-correction (LHX2: p = 6.5×10⁻⁸; RREB1: p = 5.7×10⁻⁸). Six of fourteen previously reported loci for waist in EA populations were significant (p<0.05 divided by the number of independent SNPs within the region) in AA studied here (TBX15-WARS2, GRB14, ADAMTS9, LY86, RSPO3, ITPR2-SSPN). Further, we observed associations with metabolic traits: rs13389219 at GRB14 associated with HDL-cholesterol, triglycerides, and fasting insulin, and rs13060013 at ADAMTS9 with HDL-cholesterol and fasting insulin. Finally, we observed nominal evidence for sexual dimorphism, with stronger results in AA women at the GRB14 locus (p for interaction = 0.02). In conclusion, we identified two suggestive loci associated with fat distribution in AA populations in addition to confirming 6 loci previously identified in populations of EA. These findings reinforce the concept that there are fat distribution loci that are independent of generalized adiposity.     

Identification of Allelic Heterogeneity at Type-2 Diabetes Loci and Impact on Prediction

Although over 60 single nucleotide polymorphisms (SNPs) have been identified by meta-analysis of genome-wide association studies for type-2 diabetes (T2D) among individuals of European descent, much of the genetic variation remains unexplained. There are likely many more SNPs that contribute to variation in T2D risk, some of which may lie in the regions surrounding established SNPs - a phenomenon often referred to as allelic heterogeneity. Here, we use the summary statistics from the DIAGRAM consortium meta-analysis of T2D genome-wide association studies along with linkage disequilibrium patterns inferred from a large reference sample to identify novel SNPs associated with T2D surrounding each of the previously established risk loci. We then examine the extent to which the use of these additional SNPs improves prediction of T2D risk in an independent validation dataset. Our results suggest that multiple SNPs at each of 3 loci contribute to T2D susceptibility (TCF7L2, CDKN2A/B, and KCNQ1; p<5×10⁻⁸). Using a less stringent threshold (p<5×10⁻⁴), we identify 34 additional loci with multiple associated SNPs. The addition of these SNPs slightly improves T2D prediction compared to the use of only the respective lead SNPs, when assessed using an independent validation cohort. Our findings suggest that some currently established T2D risk loci likely harbor multiple polymorphisms which contribute independently and collectively to T2D risk. This opens a promising avenue for improving prediction of T2D, and for a better understanding of the genetic architecture of T2D.     

Results of Genome-Wide Analyses on Neurodevelopmental Phenotypes at Four-Year Follow-Up following Cardiac Surgery in Infancy

Background

Adverse neurodevelopmental sequelae are reported among children who undergo early cardiac surgery to repair congenital heart defects (CHD). APOE genotype has previously been determined to contribute to the prediction of these outcomes. Understanding further genetic causes for the development of poor neurobehavioral outcomes should enhance patient risk stratification and improve both prevention and treatment strategies.

Methods

We performed a prospective observational study of children who underwent cardiac surgery before six months of age; this included a neurodevelopmental evaluation between their fourth and fifth birthdays. Attention and behavioral skills were assessed through parental report utilizing the Attention Deficit-Hyperactivity Disorder-IV scale preschool edition (ADHD-IV), and Child Behavior Checklist (CBCL/1.5-5), respectively. Of the seven investigated, three neurodevelopmental phenotypes met genomic quality control criteria. Linear regression was performed to determine the effect of genome-wide genetic variation on these three neurodevelopmental measures in 316 subjects.

Results

This genome-wide association study identified single nucleotide polymorphisms (SNPs) associated with three neurobehavioral phenotypes in the postoperative children ADHD-IV Impulsivity/Hyperactivity, CBCL/1.5-5 PDPs, and CBCL/1.5-5 Total Problems. The most predictive SNPs for each phenotype were: a LGALS8 intronic SNP, rs4659682, associated with ADHD-IV Impulsivity (P = 1.03×10⁻⁶); a PCSK5 intronic SNP, rs2261722, associated with CBCL/1.5-5 PDPs (P = 1.11×10⁻⁶); and an intergenic SNP, rs11617488, 50 kb from FGF9, associated with CBCL/1.5-5 Total Problems (P = 3.47×10⁻⁷). 10 SNPs (3 for ADHD-IV Impulsivity, 5 for CBCL/1.5-5 PDPs, and 2 for CBCL/1.5-5 Total Problems) had p<10⁻⁵.

Conclusions

No SNPs met genome-wide significance for our three neurobehavioral phenotypes; however, 10 SNPs reached a threshold for suggestive significance (p<10⁻⁵). Given the unique nature of this cohort, larger studies and/or replication are not possible. Studies to further investigate the mechanisms through which these newly identified genes may influence neurodevelopment dysfunction are warranted.     

t-10, c-12 CLA Dietary Supplementation Inhibits Atherosclerotic Lesion Development Despite Adverse Cardiovascular and Hepatic Metabolic Marker Profiles

Animal and human studies have indicated that fatty acids such as the conjugated linoleic acids (CLA) found in milk could potentially alter the risk of developing metabolic disorders including diabetes and cardiovascular disease (CVD). Using susceptible rodent models (apoE^−/− and LDLr^−/− mice) we investigated the interrelationship between mouse strain, dietary conjugated linoleic acids and metabolic markers of CVD. Despite an adverse metabolic risk profile, atherosclerosis (measured directly by lesion area), was significantly reduced with t-10, c-12 CLA and mixed isomer CLA (Mix) supplementation in both apoE^−/− (p<0.05, n = 11) and LDLr^−/− mice (p<0.01, n = 10). Principal component analysis was utilized to delineate the influence of multiple plasma and tissue metabolites on the development of atherosclerosis. Group clustering by dietary supplementation was evident, with the t-10, c-12 CLA supplemented animals having distinct patterns, suggestive of hepatic insulin resistance, regardless of mouse strain. The effect of CLA supplementation on hepatic lipid and fatty acid composition was explored in the LDLr^−/− strain. Dietary supplementation with t-10, c-12 CLA significantly increased liver weight (p<0.05, n = 10), triglyceride (p<0.01, n = 10) and cholesterol ester content (p<0.01, n = 10). Furthermore, t-10, c-12 CLA also increased the ratio of 18∶1 to 18∶0 fatty acid in the liver suggesting an increase in the activity of stearoyl-CoA desaturase. Changes in plasma adiponectin and liver weight with t-10, c-12 CLA supplementation were evident within 3 weeks of initiation of the diet. These observations provide evidence that the individual CLA isomers have divergent mechanisms of action and that t-10, c-12 CLA rapidly changes plasma and liver markers of metabolic syndrome, despite evidence of reduction in atherosclerosis.     

Genome Wide Meta-analysis Highlights the Role of Genetic Variation in RARRES2 in the Regulation of Circulating Serum Chemerin

Chemerin is an adipokine proposed to link obesity and chronic inflammation of adipose tissue. Genetic factors determining chemerin release from adipose tissue are yet unknown. We conducted a meta-analysis of genome-wide association studies (GWAS) for serum chemerin in three independent cohorts from Europe: Sorbs and KORA from Germany and PPP-Botnia from Finland (total N = 2,791). In addition, we measured mRNA expression of genes within the associated loci in peripheral mononuclear cells by micro-arrays, and within adipose tissue by quantitative RT-PCR and performed mRNA expression quantitative trait and expression-chemerin association studies to functionally substantiate our loci. Heritability estimate of circulating chemerin levels was 16.2% in the Sorbs cohort. Thirty single nucleotide polymorphisms (SNPs) at chromosome 7 within the retinoic acid receptor responder 2 (RARRES2)/Leucine Rich Repeat Containing (LRRC61) locus reached genome-wide significance (p<5.0×10⁻⁸) in the meta-analysis (the strongest evidence for association at rs7806429 with p = 7.8×10⁻¹⁴, beta = −0.067, explained variance 2.0%). All other SNPs within the cluster were in linkage disequilibrium with rs7806429 (minimum r² = 0.43 in the Sorbs cohort). The results of the subgroup analyses of males and females were consistent with the results found in the total cohort. No significant SNP-sex interaction was observed. rs7806429 was associated with mRNA expression of RARRES2 in visceral adipose tissue in women (p<0.05 after adjusting for age and body mass index). In conclusion, the present meta-GWAS combined with mRNA expression studies highlights the role of genetic variation in the RARRES2 locus in the regulation of circulating chemerin concentrations.     

Genome-Wide Association Study Identifies GPC5 as a Novel Genetic Locus Protective against Sudden Cardiac Arrest

Background

Existing studies indicate a significant genetic component for sudden cardiac arrest (SCA) and genome-wide association studies (GWAS) provide an unbiased approach for identification of novel genes. We performed a GWAS to identify genetic determinants of SCA.

Methodology/Principal Findings

We used a case-control design within the ongoing Oregon Sudden Unexpected Death Study (Oregon-SUDS). Cases (n = 424) were SCAs with coronary artery disease (CAD) among residents of Portland, OR (2002–07, population ∼1,000,000) and controls (n = 226) were residents with CAD, but no history of SCA. All subjects were of White-European ancestry and GWAS was performed using Affymetrix 500K/5.0 and 6.0 arrays. High signal markers were genotyped in SCA cases (n = 521) identified from the Atherosclerosis Risk in Communities Study (ARIC) and the Cardiovascular Health Study (CHS) (combined n = 19,611). No SNPs reached genome-wide significance (p<5×10⁻⁸). SNPs at 6 loci were prioritized for follow-up primarily based on significance of p<10⁻⁴ and proximity to a known gene (CSMD2, GPR37L1, LIN9, B4GALNT3, GPC5, and ZNF592). The minor allele of GPC5 (GLYPICAN 5, rs3864180) was associated with a lower risk of SCA in Oregon-SUDS, an effect that was also observed in ARIC/CHS whites (p<0.05) and blacks (p<0.04). In a combined Cox proportional hazards model analysis that adjusted for race, the minor allele exhibited a hazard ratio of 0.85 (95% CI 0.74 to 0.98; p<0.01).

Conclusions/Significance

A novel genetic locus for SCA, GPC5, was identified from Oregon-SUDS and successfully validated in the ARIC and CHS cohorts. Three other members of the Glypican family have been previously implicated in human disease, including cardiac conditions. The mechanism of this specific association requires further study.     

Genetically Dependent ERBB3 Expression Modulates Antigen Presenting Cell Function and Type 1 Diabetes Risk

Type 1 diabetes (T1D) is an autoimmune disease resulting from the complex interaction between multiple susceptibility genes, environmental factors and the immune system. Over 40 T1D susceptibility regions have been suggested by recent genome-wide association studies; however, the specific genes and their role in the disease remain elusive. The objective of this study is to identify the susceptibility gene(s) in the 12q13 region and investigate the functional link to the disease pathogenesis. A total of 19 SNPs in the 12q13 region were analyzed by the TaqMan assay for 1,434 T1D patients and 1,865 controls. Thirteen of the SNPs are associated with T1D (best p = 4×10⁻¹¹), thus providing confirmatory evidence for at least one susceptibility gene in this region. To identify candidate genes, expression of six genes in the region was analyzed by real-time RT-PCR for PBMCs from 192 T1D patients and 192 controls. SNP genotypes in the 12q13 region are the main factors that determine ERBB3 mRNA levels in PBMCs. The protective genotypes for T1D are associated with higher ERBB3 mRNA level (p<10⁻¹⁰). Furthermore, ERBB3 protein is expressed on the surface of CD11c⁺ cells (dendritic cells and monocytes) in peripheral blood after stimulation with LPS, polyI:C or CpG. Subjects with protective genotypes have significantly higher percentages of ERBB3⁺ monocytes and dendritic cells (p = 1.1×10⁻⁹); and the percentages of ERBB3⁺ cells positively correlate with the ability of APC to stimulate T cell proliferation (R² = 0.90, p<0.0001). Our results indicate that ERBB3 plays a critical role in determining APC function and potentially T1D pathogenesis.     

Identification of Non-HLA Genes Associated with Celiac Disease and Country-Specific Differences in a Large,International Pediatric Cohort

Objectives

There are significant geographical differences in the prevalence and incidence of celiac disease that cannot be explained by HLA alone. More than 40 loci outside of the HLA region have been associated with celiac disease. We investigated the roles of these non-HLA genes in the development of tissue transglutaminase autoantibodies (tTGA) and celiac disease in a large international prospective cohort study.

Methods

A total of 424,788 newborns from the US and European general populations and first-degree relatives with type 1 diabetes were screened for specific HLA genotypes. Of these, 21,589 carried 1 of the 9 HLA genotypes associated with increased risk for type 1 diabetes and celiac disease; we followed 8676 of the children in a 15 y prospective follow-up study. Genotype analyses were performed on 6010 children using the Illumina ImmunoChip. Levels of tTGA were measured in serum samples using radio-ligand binding assays; diagnoses of celiac disease were made based on persistent detection of tTGA and biopsy analysis. Data were analyzed using Cox proportional hazards analyses.

Results

We found 54 single-nucleotide polymorphisms (SNPs) in 5 genes associated with celiac disease (TAGAP, IL18R1, RGS21, PLEK, and CCR9) in time to celiac disease analyses (10⁻⁴>P>5.8x10⁻⁶). The hazard ratios (HR) for the SNPs with the smallest P values in each region were 1.59, 1.45, 2.23, 2.64, and 1.40, respectively. Outside of regions previously associated with celiac disease, we identified 10 SNPs in 8 regions that could also be associated with the disease (P<10⁻⁴). A SNP near PKIA (rs117128341, P = 6.5x10⁻⁸, HR = 2.8) and a SNP near PFKFB3 (rs117139146, P<2.8x10⁻⁷, HR = 4.9) reached the genome-wide association threshold in subjects from Sweden. Analyses of time to detection of tTGA identified 29 SNPs in 2 regions previously associated with celiac disease (CTLA4, P = 1.3x10⁻⁶, HR = 0.76 and LPP, P = 2.8x10⁻⁵, HR = .80) and 6 SNPs in 5 regions not previously associated with celiac disease (P<10⁻⁴); non-HLA genes are therefore involved in development of tTGA.

Conclusions

In conclusion, using a genetic analysis of a large international cohort of children, we associated celiac disease development with 5 non-HLA regions previously associated with the disease and 8 regions not previously associated with celiac disease. We identified 5 regions associated with development of tTGA. Two loci associated with celiac disease progression reached a genome-wide association threshold in subjects from Sweden.     

Neutrophil-Derived Myeloperoxidase Aggravates Non-Alcoholic Steatohepatitis in Low-Density Lipoprotein Receptor-Deficient Mice

Background

Chronic inflammation and oxidative stress play fundamental roles in the pathogenesis of non-alcoholic steatohepatitis (NASH). Previously, we reported that myeloperoxidase (MPO), an aggressive oxidant-generating neutrophil enzyme, is associated with NASH severity in man. We now investigated the hypothesis that MPO contributes to the development and progression of NASH.

Methodology

Low-density lipoprotein receptor-deficient mice with an MPO-deficient hematopoietic system (LDLR^−/−/MPO^−/−tp mice) were generated and compared with LDLR^−/−/MPO^+/+tp mice after induction of NASH by high-fat feeding.

Results

High-fat feeding caused a ∼4-fold induction of liver MPO in LDLR^−/−/MPO^+/+ mice which was associated with hepatic sequestration of MPO-positive neutrophils and high levels of nitrotyrosine, a marker of MPO activity. Importantly, LDLR^−/−/MPO^−/−tp mice displayed markedly reduced hepatic neutrophil and T-lymphocyte infiltration (p<0.05), and strong down regulation of pro-inflammatory genes such as TNF-α and IL-6 (p<0.05, p<0.01) in comparison with LDLR^−/−/MPO^+/+tp mice. Next to the generalized reduction of inflammation, liver cholesterol accumulation was significantly diminished in LDLR^−/−/MPO^−/−tp mice (p = 0.01). Moreover, MPO deficiency appeared to attenuate the development of hepatic fibrosis as evident from reduced hydroxyproline levels (p<0.01). Interestingly, visceral adipose tissue inflammation was markedly reduced in LDLR^−/−/MPO^−/−tp mice, with a complete lack of macrophage crown-like structures. In conclusion, MPO deficiency attenuates the development of NASH and diminishes adipose tissue inflammation in response to a high fat diet, supporting an important role for neutrophils in the pathogenesis of metabolic disease.     

Targeting 160 Candidate Genes for Blood Pressure Regulation with a Genome-Wide Genotyping Array

The outcome of Genome-Wide Association Studies (GWAS) has challenged the field of blood pressure (BP) genetics as previous candidate genes have not been among the top loci in these scans. We used Affymetrix 500K genotyping data of KORA S3 cohort (n = 1,644; Southern-Germany) to address (i) SNP coverage in 160 BP candidate genes; (ii) the evidence for associations with BP traits in genome-wide and replication data, and haplotype analysis. In total, 160 gene regions (genic region±10 kb) covered 2,411 SNPs across 11.4 Mb. Marker densities in genes varied from 0 (n = 11) to 0.6 SNPs/kb. On average 52.5% of the HAPMAP SNPs per gene were captured. No evidence for association with BP was obtained for 1,449 tested SNPs. Considerable associations (P<10⁻³) were detected for the genes, where >50% of HAPMAP SNPs were tagged. In general, genes with higher marker density (>0.2 SNPs/kb) revealed a better chance to reach close to significance associations. Although, none of the detected P-values remained significant after Bonferroni correction (P<0.05/2319, P<2.15×10⁻⁵), the strength of some detected associations was close to this level: rs10889553 (LEPR) and systolic BP (SBP) (P = 4.5×10⁻⁵) as well as rs10954174 (LEP) and diastolic BP (DBP) (P = 5.20×10⁻⁵). In total, 12 markers in 7 genes (ADRA2A, LEP, LEPR, PTGER3, SLC2A1, SLC4A2, SLC8A1) revealed considerable association (P<10⁻³) either with SBP, DBP, and/or hypertension (HYP). None of these were confirmed in replication samples (KORA S4, HYPEST, BRIGHT). However, supportive evidence for the association of rs10889553 (LEPR) and rs11195419 (ADRA2A) with BP was obtained in meta-analysis across samples stratified either by body mass index, smoking or alcohol consumption. Haplotype analysis highlighted LEPR and PTGER3. In conclusion, the lack of associations in BP candidate genes may be attributed to inadequate marker coverage on the genome-wide arrays, small phenotypic effects of the loci and/or complex interaction with life-style and metabolic parameters.     

Genome-Wide Association Study of Celiac Disease in North America Confirms FRMD4B as New Celiac Locus

We performed a genome-wide association study (GWAS) of 1550 North American celiac disease cases and 3084 controls. Twelve SNPs, distributed across four regions (3p21.31, 4q27, 6q15, 6q25), were significantly associated with disease (p-value <1.0×10⁻⁷), and a further seven SNPs, across four additional regions (1q24.3, 10p15.1, 6q22.31, 17q21.32) had suggestive evidence (1.0×10⁻⁷ < p-value < 1.0×10⁻⁶). This study replicated a previous suggestive association within FRMD4B (3p14.1), confirming it as a celiac disease locus. All four regions with significant associations and two regions with suggestive results (1q24.3, 10p15.1) were known disease loci. The 6q22.31 and 10p11.23 regions were not replicated. A total of 410 SNPs distributed across the eight significant and suggestive regions were tested for association with dermatitis herpetiformis and microscopic colitis. Preliminary, suggestive statistical evidence for association with the two traits was found at chromosomes 3p21.31, 6q15, 6q25, 1q24.3 and 10p11.23, with future studies being required to validate the reported associations.     

Blood cis-eQTL Analysis Fails to Identify Novel Association Signals among Sub-Threshold Candidates from Genome-Wide Association Studies in Restless Legs Syndrome

Restless legs syndrome (RLS) is a common neurologic disorder characterized by nightly dysesthesias affecting the legs primarily during periods of rest and relieved by movement. RLS is a complex genetic disease and susceptibility factors in six genomic regions have been identified by means of genome-wide association studies (GWAS). For some complex genetic traits, expression quantitative trait loci (eQTLs) are enriched among trait-associated single nucleotide polymorphisms (SNPs). With the aim of identifying new genetic susceptibility factors for RLS, we assessed the 332 best-associated SNPs from the genome-wide phase of the to date largest RLS GWAS for cis-eQTL effects in peripheral blood from individuals of European descent. In 740 individuals belonging to the KORA general population cohort, 52 cis-eQTLs with p_nominal<10⁻³ were identified, while in 976 individuals belonging to the SHIP-TREND general population study 53 cis-eQTLs with p_nominal<10⁻³ were present. 23 of these cis-eQTLs overlapped between the two cohorts. Subsequently, the twelve of the 23 cis-eQTL SNPs, which were not located at an already published RLS-associated locus, were tested for association in 2449 RLS cases and 1462 controls. The top SNP, located in the DET1 gene, was nominally significant (p<0.05) but did not withstand correction for multiple testing (p = 0.42). Although a similar approach has been used successfully with regard to other complex diseases, we were unable to identify new genetic susceptibility factor for RLS by adding this novel level of functional assessment to RLS GWAS data.     

Genetic Determinants of Pelvic Organ Prolapse among African American and Hispanic Women in the Women’s Health Initiative

Current evidence suggests a multifactorial etiology to pelvic organ prolapse (POP), including genetic predisposition. We conducted a genome-wide association study of POP in African American (AA) and Hispanic (HP) women from the Women’s Health Initiative Hormone Therapy study. Cases were defined as any POP (grades 1–3) or moderate/severe POP (grades 2–3), while controls had grade 0 POP. We performed race-specific multiple logistic regression analyses between SNPs imputed to 1000 genomes in relation to POP (grade 0 vs 1–3; grade 0 vs 2–3) adjusting for age at diagnosis, body mass index, parity, and genetic ancestry. There were 1274 controls and 1427 cases of any POP and 317 cases of moderate/severe POP. Although none of the analyses reached genome-wide significance (p<5x10^-8), we noted variants in several loci that met p<10⁻⁶. In race-specific analysis of grade 0 vs 2–3, intronic SNPs in the CPE gene (rs28573326, OR:2.14; 95% CI 1.62–2.83; p = 1.0x10^-7) were associated with POP in AAs, and SNPs in the gene AL132709.5 (rs1950626, OR:2.96; 95% CI 1.96–4.48, p = 2.6x10^-7) were associated with POP in HPs. Inverse variance fixed-effect meta-analysis of the race-specific results showed suggestive signals for SNPs in the DPP6 gene (rs11243354, OR:1.36; p = 4.2x10^-7) in the grade 0 vs 1–3 analyses and for SNPs around PGBD5 (rs740494, OR:2.17; p = 8.6x10^-7) and SHC3 (rs2209875, OR:0.60; p = 9.3x10^-7) in the grade 0 vs 2–3 analyses. While we did not identify genome-wide significant findings, we document several SNPs reaching suggestive statistical significance. Further interrogation of POP in larger minority samples is warranted.     

Quantitative Assessment of the Association between rs2046210 at 6q25.1 and Breast Cancer Risk

Genome-wide association studies (GWAS) have identified several genetic susceptibility loci for breast cancer (BC). One of them, conducted among Chinese women, found an association of rs2046210 at 6q25.1 with the risk of BC recently. Since then, numerous association studies have been carried out to investigate the relationship between this polymorphism and BC risk in various populations. However, these have yielded contradictory results. We therefore performed a meta-analysis to clarify this inconsistency. Overall, a total of 235003 subjects based on 13 studies were included in our study. Significantly increased BC risk was detected in the pooled analysis [allele contrast: OR = 1.13, 95%CI = 1.10–1.17, P(Z) <10⁻⁵, P(Q) <10⁻⁴; dominant model: OR = 1.21, 95%CI = 1.14–1.27, P(Z) <10⁻⁵, P(Q) <10⁻⁴; recessive model: OR = 1.18, 95%CI = 1.12–1.24, P(Z) <10⁻⁵, P(Q) = 0.04]. In addition, our data revealed that rs2046210 conferred greater risk in estrogen receptor (ER)-negative tumors [OR = 1.27, 95%CI = 1.15–1.40, P(Z) <10⁻⁵, P(Q) <10⁻⁴] than in ER-positive ones [OR = 1.18, 95%CI = 1.09–1.28, P(Z) <10⁻⁴, P(Q) = 0.0003]. When stratified by ethnicity, significant associations were found in Caucasian and Asian populations, but not detected among Africans. There was evidence of heterogeneity (P<0.05), however, the heterogeneity largely disappeared after stratification by ethnicity. The present meta-analysis demonstrated that the rs2046210 polymorphism may be associated with increased BC susceptibility, but this association varies in different ethnicities.     

Genetic Association Analysis of Complex Diseases Incorporating Intermediate Phenotype Information

Genetic researchers often collect disease related quantitative traits in addition to disease status because they are interested in understanding the pathophysiology of disease processes. In genome-wide association (GWA) studies, these quantitative phenotypes may be relevant to disease development and serve as intermediate phenotypes or they could be behavioral or other risk factors that predict disease risk. Statistical tests combining both disease status and quantitative risk factors should be more powerful than case-control studies, as the former incorporates more information about the disease. In this paper, we proposed a modified inverse-variance weighted meta-analysis method to combine disease status and quantitative intermediate phenotype information. The simulation results showed that when an intermediate phenotype was available, the inverse-variance weighted method had more power than did a case-control study of complex diseases, especially in identifying susceptibility loci having minor effects. We further applied this modified meta-analysis to a study of imputed lung cancer genotypes with smoking data in 1154 cases and 1137 matched controls. The most significant SNPs came from the CHRNA3-CHRNA5-CHRNB4 region on chromosome 15q24–25.1, which has been replicated in many other studies. Our results confirm that this CHRNA region is associated with both lung cancer development and smoking behavior. We also detected three significant SNPs—rs1800469, rs1982072, and rs2241714—in the promoter region of the TGFB1 gene on chromosome 19 (p = 1.46×10⁻⁵, 1.18×10⁻⁵, and 6.57×10⁻⁶, respectively). The SNP rs1800469 is reported to be associated with chronic obstructive pulmonary disease and lung cancer in cigarette smokers. The present study is the first GWA study to replicate this result. Signals in the 3q26 region were also identified in the meta-analysis. We demonstrate the intermediate phenotype can potentially enhance the power of complex disease association analysis and the modified meta-analysis method is robust to incorporate intermediate phenotype or other quantitative risk factor in the analysis.