Rare genetic variants explain missing heritability in smoking.

Common genetic variants explain less variation in complex phenotypes than inferred from family-based studies, and there is a debate on the source of this 'missing heritability'. We investigated the contribution of rare genetic variants to tobacco use with whole-genome sequences from up to 26,257 unrelated individuals of European ancestries and 11,743 individuals of African ancestries. Across four smoking traits, single-nucleotide-polymorphism-based heritability ([Formula: see text]) was estimated from 0.13 to 0.28 (s.e., 0.10-0.13) in European ancestries, with 35-74% of it attributable to rare variants with minor allele frequencies between 0.01% and 1%. These heritability estimates are 1.5-4 times higher than past estimates based on common variants alone and accounted for 60% to 100% of our pedigree-based estimates of narrow-sense heritability ([Formula: see text], 0.18-0.34). In the African ancestry samples, [Formula: see text] was estimated from 0.03 to 0.33 (s.e., 0.09-0.14) across the four smoking traits. These results suggest that rare variants are important contributors to the heritability of smoking.

Associations between DNA methylation and BMI vary by metabolic health status: a potential link to disparate cardiovascular outcomes.

BACKGROUND: Body mass index (BMI), a well-known risk factor for poor cardiovascular outcomes, is associated with differential DNA methylation (DNAm). Similarly, metabolic health has also been associated with changes in DNAm. It is unclear how overall metabolic health outside of BMI may modify the relationship between BMI and methylation profiles, and what consequences this may have on downstream cardiovascular disease. The purpose of this study was to identify cytosine-phosphate-guanine (CpG) sites at which the association between BMI and DNAm could be modified by overall metabolic health. RESULTS: The discovery study population was derived from three Women's Health Initiative (WHI) ancillary studies (n = 3977) and two Atherosclerosis Risk in Communities (ARIC) ancillary studies (n = 3520). Findings were validated in the Multi-Ethnic Study of Atherosclerosis (MESA) cohort (n = 1200). Generalized linear models regressed methylation ß values on the interaction between BMI and metabolic health Z score (BMI × MHZ) adjusted for BMI, MHZ, cell composition, chip number and location, study characteristics, top three ancestry principal components, smoking, age, ethnicity (WHI), and sex (ARIC). Among the 429,566 sites examined, differential associations between BMI × MHZ and DNAm were identified at 22 CpG sites (FDR q < 0.05), with one site replicated in MESA (cg18989722, in the TRAPPC9 gene). Three of the 22 sites were associated with incident coronary heart disease (CHD) in WHI. For each 0.01 unit increase in DNAm ß value, the risk of incident CHD increased by 9% in one site and decreased by 6-10% in two sites over 25 years. CONCLUSIONS: Differential associations between DNAm and BMI by MHZ were identified at 22 sites, one of which was validated (cg18989722) and three of which were predictive of incident CHD. These sites are located in several genes related to NF-kappa-B signaling, suggesting a potential role for inflammation between DNA methylation and BMI-associated metabolic health.

Genome-wide association study identifying novel variant for fasting insulin and allelic heterogeneity in known glycemic loci in Chilean adolescents: The Santiago Longitudinal Study.

BACKGROUND: The genetic underpinnings of glycemic traits have been understudied in adolescent and Hispanic/Latino (H/L) populations in comparison to adults and populations of European ancestry. OBJECTIVE: To identify genetic factors underlying glycemic traits in an adolescent H/L population. METHODS: We conducted a genome-wide association study (GWAS) of fasting glucose (FG) and fasting insulin (FI) in H/L adolescents from the Santiago Longitudinal Study. RESULTS: We identified one novel variant positioned in the CSMD1 gene on chromosome 8 (rs77465890, effect allele frequency = 0.10) that was associated with FI (ß = -0.299, SE = 0.054, p = 2.72×10-8 ) and was only slightly attenuated after adjusting for body mass index z-scores (ß = -0.252, SE = 0.047, p = 1.03×10-7 ). We demonstrated directionally consistent, but not statistically significant results in African and Hispanic adults of the Population Architecture Using Genomics and Epidemiology Consortium. We also identified secondary signals for two FG loci after conditioning on known variants, which demonstrate allelic heterogeneity in well-known glucose loci. CONCLUSION: Our results exemplify the importance of including populations with diverse ancestral origin and adolescent participants in GWAS of glycemic traits to uncover novel risk loci and expand our understanding of disease aetiology.

A genome-wide association study of polycystic ovary syndrome identified from electronic health records.

BACKGROUND: Polycystic ovary syndrome is the most common endocrine disorder affecting women of reproductive age. A number of criteria have been developed for clinical diagnosis of polycystic ovary syndrome, with the Rotterdam criteria being the most inclusive. Evidence suggests that polycystic ovary syndrome is significantly heritable, and previous studies have identified genetic variants associated with polycystic ovary syndrome diagnosed using different criteria. The widely adopted electronic health record system provides an opportunity to identify patients with polycystic ovary syndrome using the Rotterdam criteria for genetic studies. OBJECTIVE: To identify novel associated genetic variants under the same phenotype definition, we extracted polycystic ovary syndrome cases and unaffected controls based on the Rotterdam criteria from the electronic health records and performed a discovery-validation genome-wide association study. STUDY DESIGN: We developed a polycystic ovary syndrome phenotyping algorithm on the basis of the Rotterdam criteria and applied it to 3 electronic health record-linked biobanks to identify cases and controls for genetic study. In the discovery phase, we performed an individual genome-wide association study using the Geisinger MyCode and the Electronic Medical Records and Genomics cohorts, which were then meta-analyzed. We attempted validation of the significant association loci (P<1×10-6) in the BioVU cohort. All association analyses used logistic regression, assuming an additive genetic model, and adjusted for principal components to control for population stratification. An inverse-variance fixed-effect model was adopted for meta-analysis. In addition, we examined the top variants to evaluate their associations with each criterion in the phenotyping algorithm. We used the STRING database to characterize protein-protein interaction network. RESULTS: Using the same algorithm based on the Rotterdam criteria, we identified 2995 patients with polycystic ovary syndrome and 53,599 population controls in total (2742 cases and 51,438 controls from the discovery phase; 253 cases and 2161 controls in the validation phase). We identified 1 novel genome-wide significant variant rs17186366 (odds ratio [OR]=1.37 [1.23, 1.54], P=2.8×10-8) located near SOD2. In addition, 2 loci with suggestive association were also identified: rs113168128 (OR=1.72 [1.42, 2.10], P=5.2×10-8), an intronic variant of ERBB4 that is independent from the previously published variants, and rs144248326 (OR=2.13 [1.52, 2.86], P=8.45×10-7), a novel intronic variant in WWTR1. In the further association tests of the top 3 single-nucleotide polymorphisms with each criterion in the polycystic ovary syndrome algorithm, we found that rs17186366 (SOD2) was associated with polycystic ovaries and hyperandrogenism, whereas rs11316812 (ERBB4) and rs144248326 (WWTR1) were mainly associated with oligomenorrhea or infertility. We also validated the previously reported association with DENND1A1. Using the STRING database to characterize protein-protein interactions, we found both ERBB4 and WWTR1 can interact with YAP1, which has been previously associated with polycystic ovary syndrome. CONCLUSION: Through a discovery-validation genome-wide association study on polycystic ovary syndrome identified from electronic health records using an algorithm based on Rotterdam criteria, we identified and validated a novel genome-wide significant association with a variant near SOD2. We also identified a novel independent variant within ERBB4 and a suggestive association with WWTR1. With previously identified polycystic ovary syndrome gene YAP1, the ERBB4-YAP1-WWTR1 network suggests involvement of the epidermal growth factor receptor and the Hippo pathway in the multifactorial etiology of polycystic ovary syndrome.

Association of dietary folate and vitamin B-12 intake with genome-wide DNA methylation in blood: a large-scale epigenome-wide association analysis in 5841 individuals.

BACKGROUND: Folate and vitamin B-12 are essential micronutrients involved in the donation of methyl groups in cellular metabolism. However, associations between intake of these nutrients and genome-wide DNA methylation levels have not been studied comprehensively in humans. OBJECTIVE: The aim of this study was to assess whether folate and/or vitamin B-12 intake are asssociated with genome-wide changes in DNA methylation in leukocytes. METHODS: A large-scale epigenome-wide association study of folate and vitamin B-12 intake was performed on DNA from 5841 participants from 10 cohorts using Illumina 450k arrays. Folate and vitamin B-12 intakes were calculated from food-frequency questionnaires (FFQs). Continuous and categorical (low compared with high intake) linear regression mixed models were applied per cohort, controlling for confounders. A meta-analysis was performed to identify significant differentially methylated positions (DMPs) and regions (DMRs), and a pathway analysis was performed on the DMR annotated genes. RESULTS: The categorical model resulted in 6 DMPs, which are all negatively associated with folate intake, annotated to FAM64A, WRAP73, FRMD8, CUX1, and LCN8 genes, which have a role in cellular processes including centrosome localization, cell proliferation, and tumorigenesis. Regional analysis showed 74 folate-associated DMRs, of which 73 were negatively associated with folate intake. The most significant folate-associated DMR was a 400-base pair (bp) spanning region annotated to the LGALS3BP gene. In the categorical model, vitamin B-12 intake was associated with 29 DMRs annotated to 48 genes, of which the most significant was a 1100-bp spanning region annotated to the calcium-binding tyrosine phosphorylation-regulated gene (CABYR). Vitamin B-12 intake was not associated with DMPs. CONCLUSIONS: We identified novel epigenetic loci that are associated with folate and vitamin B-12 intake. Interestingly, we found a negative association between folate and DNA methylation. Replication of these methylation loci is necessary in future studies.

Direct and indirect genetic effects on triglycerides through omics and correlated phenotypes.

Even though there has been great success in identifying lipid-associated single-nucleotide polymorphisms (SNPs), the mechanisms through which the SNPs act on each trait are poorly understood. The emergence of large, complex biological data sets in well-characterized cohort studies offers an opportunity to investigate the genetic effects on trait variability as a way of informing the causal genes and biochemical pathways that are involved in lipoprotein metabolism. However, methods for simultaneously analyzing multiple omics, environmental exposures, and longitudinally measured, correlated phenotypes are lacking. The purpose of our study was to demonstrate the utility of the structural equation modeling (SEM) approach to inform our understanding of the pathways by which genetic variants lead to disease risk. With the SEM method, we examine multiple pathways directly and indirectly through previously identified triglyceride (TG)-associated SNPs, methylation, and high-density lipoprotein (HDL), including sex, age, and smoking behavior, while adding in biologically plausible direct and indirect pathways. We observed significant SNP effects (P < 0.05 and directionally consistent) on TGs at visit 4 (TG4) for five loci, including rs645040 (DOCK7), rs964184 (ZPR1/ZNF259), rs4765127 (ZNF664), rs1121980 (FTO), and rs10401969 (SUGP1). Across these loci, we identify three with strong evidence of an indirect genetic effect on TG4 through HDL, one with evidence of pleiotropic effect on HDL and TG4, and one variant that acts on TG4 indirectly through a nearby methylation site. Such information can be used to prioritize candidate genes in regions of interest, inform mechanisms of action of methylation effects, and highlight possible genes with pleiotropic effects.

Exome Sequencing-Based Screening for BRCA1/2 Expected Pathogenic Variants Among Adult Biobank Participants.

Importance: Detection of disease-associated variants in the BRCA1 and BRCA2 (BRCA1/2) genes allows for cancer prevention and early diagnosis in high-risk individuals. Objectives: To identify pathogenic and likely pathogenic (P/LP) BRCA1/2 variants in an unselected research cohort, and to characterize the features associated with P/LP variants. Design, Setting, and Participants: This is a cross-sectional study of adult volunteers (n = 50 726) who underwent exome sequencing at a single health care system (Geisinger Health System, Danville, Pennsylvania) from January 1, 2014, to March 1, 2016. Participants are part of the DiscovEHR cohort and were identified through the Geisinger MyCode Community Health Initiative. They consented to a research protocol that included sequencing and return of actionable test results. Clinical data from electronic health records and clinical visits were correlated with variants. Comparisons were made between those with (cases) and those without (controls) P/LP variants in BRCA1/2. Main Outcomes: Prevalence of P/LP BRCA1/2 variants in cohort, proportion of variant carriers not previously ascertained through clinical testing, and personal and family history of relevant cancers among BRCA1/2 variant carriers and noncarriers. Results: Of the 50 726 health system patients who underwent exome sequencing, 50 459 (99.5%) had no expected pathogenic BRCA1/2 variants and 267 (0.5%) were BRCA1/2 carriers. Of the 267 cases (148 [55.4%] were women and 119 [44.6%] were men with a mean [range] age of 58.9 [23-90] years), 183 (68.5%) received clinically confirmed results in their electronic health record. Among the 267 participants with P/LP BRCA1/2 variants, 219 (82.0%) had no prior clinical testing, 95 (35.6%) had BRCA1 variants, and 172 (64.4%) had BRCA2 variants. Syndromic cancer diagnoses were present in 11 (47.8%) of the 23 deceased BRCA1/2 carriers and in 56 (20.9%) of all 267 BRCA1/2 carriers. Among women, 31 (20.9%) of 148 variant carriers had a personal history of breast cancer, compared with 1554 (5.2%) of 29 880 noncarriers (odds ratio [OR], 5.95; 95% CI, 3.88-9.13; P < .001). Ovarian cancer history was present in 15 (10.1%) of 148 variant carriers and in 195 (0.6%) of 29 880 variant noncarriers (OR, 18.30; 95% CI, 10.48-31.4; P < .001). Among 89 BRCA1/2 carriers without prior testing but with comprehensive personal and family history data, 44 (49.4%) did not meet published guidelines for clinical testing. Conclusions and Relevance: This study found that compared with previous clinical care, exome sequencing-based screening identified 5 times as many individuals with P/LP BRCA1/2 variants. These findings suggest that genomic screening may identify BRCA1/2-associated cancer risk that might otherwise remain undetected within health care systems and may provide opportunities to reduce morbidity and mortality in patients.

Ranking and characterization of established BMI and lipid associated loci as candidates for gene-environment interactions.

Phenotypic variance heterogeneity across genotypes at a single nucleotide polymorphism (SNP) may reflect underlying gene-environment (G×E) or gene-gene interactions. We modeled variance heterogeneity for blood lipids and BMI in up to 44,211 participants and investigated relationships between variance effects (Pv), G×E interaction effects (with smoking and physical activity), and marginal genetic effects (Pm). Correlations between Pv and Pm were stronger for SNPs with established marginal effects (Spearman's ρ = 0.401 for triglycerides, and ρ = 0.236 for BMI) compared to all SNPs. When Pv and Pm were compared for all pruned SNPs, only BMI was statistically significant (Spearman's ρ = 0.010). Overall, SNPs with established marginal effects were overrepresented in the nominally significant part of the Pv distribution (Pbinomial <0.05). SNPs from the top 1% of the Pm distribution for BMI had more significant Pv values (PMann-Whitney = 1.46×10-5), and the odds ratio of SNPs with nominally significant (<0.05) Pm and Pv was 1.33 (95% CI: 1.12, 1.57) for BMI. Moreover, BMI SNPs with nominally significant G×E interaction P-values (Pint<0.05) were enriched with nominally significant Pv values (Pbinomial = 8.63×10-9 and 8.52×10-7 for SNP × smoking and SNP × physical activity, respectively). We conclude that some loci with strong marginal effects may be good candidates for G×E, and variance-based prioritization can be used to identify them.

Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits.

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution.

Genome-wide association of trajectories of systolic blood pressure change.

BACKGROUND: There is great interindividual variation in systolic blood pressure (SBP) as a result of the influences of several factors, including sex, ancestry, smoking status, medication use, and, especially, age. The majority of genetic studies have examined SBP measured cross-sectionally; however, SBP changes over time, and not necessarily in a linear fashion. Therefore, this study conducted a genome-wide association (GWA) study of SBP change trajectories using data available through the Genetic Analysis Workshop 19 (GAW19) of 959 individuals from 20 extended Mexican American families from the San Antonio Family Studies with up to 4 measures of SBP. We performed structural equation modeling (SEM) while taking into account potential genetic effects to identify how, if at all, to include covariates in estimating the SBP change trajectories using a mixture model based latent class growth modeling (LCGM) approach for use in the GWA analyses. RESULTS: The semiparametric LCGM approach identified 5 trajectory classes that captured SBP changes across age. Each LCGM identified trajectory group was ranked based on the average number of cumulative years as hypertensive. Using a pairwise comparison of these classes the heritability estimates range from 12 to 94 % (SE = 17 to 40 %). CONCLUSION: These identified trajectories are significantly heritable, and we identified a total of 8 promising loci that influence one's trajectory in SBP change across age. Our results demonstrate the potential utility of capitalizing on extant genetic data and longitudinal SBP assessments available through GAW19 to explore novel analytical methods with promising results.

Mendelian Randomization Study of Body Mass Index and Colorectal Cancer Risk.

BACKGROUND: High body mass index (BMI) is consistently linked to increased risk of colorectal cancer for men, whereas the association is less clear for women. As risk estimates from observational studies may be biased and/or confounded, we conducted a Mendelian randomization study to estimate the causal association between BMI and colorectal cancer. METHODS: We used data from 10,226 colorectal cancer cases and 10,286 controls of European ancestry. The Mendelian randomization analysis used a weighted genetic risk score, derived from 77 genome-wide association study-identified variants associated with higher BMI, as an instrumental variable (IV). We compared the IV odds ratio (IV-OR) with the OR obtained using a conventional covariate-adjusted analysis. RESULTS: Individuals carrying greater numbers of BMI-increasing alleles had higher colorectal cancer risk [per weighted allele OR, 1.31; 95% confidence interval (CI), 1.10-1.57]. Our IV estimation results support the hypothesis that genetically influenced BMI is directly associated with risk for colorectal cancer (IV-OR per 5 kg/m(2), 1.50; 95% CI, 1.13-2.01). In the sex-specific IV analyses higher BMI was associated with higher risk of colorectal cancer among women (IV-OR per 5 kg/m(2), 1.82; 95% CI, 1.26-2.61). For men, genetically influenced BMI was not associated with colorectal cancer (IV-OR per 5 kg/m(2), 1.18; 95% CI, 0.73-1.92). CONCLUSIONS: High BMI was associated with increased colorectal cancer risk for women. Whether abdominal obesity, rather than overall obesity, is a more important risk factor for men requires further investigation. IMPACT: Overall, conventional epidemiologic and Mendelian randomization studies suggest a strong association between obesity and the risk of colorectal cancer.