Unsupervised discovery of ancestry-informative markers and genetic admixture proportions in biobank-scale datasets.

Admixture estimation plays a crucial role in ancestry inference and genome-wide association studies (GWASs). Computer programs such as ADMIXTURE and STRUCTURE are commonly employed to estimate the admixture proportions of sample individuals. However, these programs can be overwhelmed by the computational burdens imposed by the 105 to 106 samples and millions of markers commonly found in modern biobanks. An attractive strategy is to run these programs on a set of ancestry-informative SNP markers (AIMs) that exhibit substantially different frequencies across populations. Unfortunately, existing methods for identifying AIMs require knowing ancestry labels for a subset of the sample. This supervised learning approach creates a chicken and the egg scenario. In this paper, we present an unsupervised, scalable framework that seamlessly carries out AIM selection and likelihood-based estimation of admixture proportions. Our simulated and real data examples show that this approach is scalable to modern biobank datasets. OpenADMIXTURE, our Julia implementation of the method, is open source and available for free.

Genetic variants of glucose metabolism and exposure to smoking in African American breast cancer.

Insulin resistance (IR) is a well-established risk factor for breast cancer (BC) development in African American (AA) postmenopausal women. While obesity and IR are more prevalent in AA than in white women, they are under-represented in genome-wide studies for systemic regulation of IR. By examining 780 genome-wide IR single-nucleotide polymorphisms (SNPs) available in our data, we tested 4689 AA women in a Random Survival Forest framework. With 37 BC-associated lifestyle factors, we conducted a gene-environment interaction analysis to estimate risk prediction for BC with the most influential genetic and behavioral factors and evaluated their combined and joint effects on BC risk. By accounting for variations of individual SNPs in BC in the prediction model, we detected four fasting glucose-associated SNPs in PCSK1, SPC25, ADCY5, and MTNR1B and three lifestyle factors (smoking, oral contraceptive use, and age at menopause) as the most predictive markers for BC risk. Our joint analysis of risk genotypes and lifestyle with smoking revealed a synergistic effect on the increased risk of BC, particularly estrogen/progesterone positive (ER/PR+) BC, in a gene-lifestyle dose-dependent manner. The joint effect of smoking was more substantial in women with prolonged exposure to cigarette smoking and female hormones. The top genome-wide association-SNPs associated with metabolic biomarkers in combination with lifestyles synergistically increase the predictability of invasive ER/PR+ BC risk among AA women. Our findings highlight generically targeted preventive interventions for women who carry particular risk genotypes and lifestyles.

Synergistic Effects of Genetic Variants of Glucose Homeostasis and Lifelong Exposures to Cigarette Smoking, Female Hormones, and Dietary Fat Intake on Primary Colorectal Cancer Development in African and Hispanic/Latino American Women.

BACKGROUND: Disparities in cancer genomic science exist among racial/ethnic minorities. Particularly, African American (AA) and Hispanic/Latino American (HA) women, the 2 largest minorities, are underrepresented in genetic/genome-wide studies for cancers and their risk factors. We conducted on AA and HA postmenopausal women a genomic study for insulin resistance (IR), the main biologic mechanism underlying colorectal cancer (CRC) carcinogenesis owing to obesity. METHODS: With 780 genome-wide IR-specific single-nucleotide polymorphisms (SNPs) among 4,692 AA and 1,986 HA women, we constructed a CRC-risk prediction model. Along with these SNPs, we incorporated CRC-associated lifestyles in the model of each group and detected the topmost influential genetic and lifestyle factors. Further, we estimated the attributable risk of the topmost risk factors shared by the groups to explore potential factors that differentiate CRC risk between these groups. RESULTS: In both groups, we detected IR-SNPs in PCSK1 (in AA) and IFT172, GCKR, and NRBP1 (in HA) and risk lifestyles, including long lifetime exposures to cigarette smoking and endogenous female hormones and daily intake of polyunsaturated fatty acids (PFA), as the topmost predictive variables for CRC risk. Combinations of those top genetic- and lifestyle-markers synergistically increased CRC risk. Of those risk factors, dietary PFA intake and long lifetime exposure to female hormones may play a key role in mediating racial disparity of CRC incidence between AA and HA women. CONCLUSIONS: Our results may improve CRC risk prediction performance in those medically/scientifically underrepresented groups and lead to the development of genetically informed interventions for cancer prevention and therapeutic effort, thus contributing to reduced cancer disparities in those minority subpopulations.

Molecular Biology Networks and Key Gene Regulators for Inflammatory Biomarkers Shared by Breast Cancer Development: Multi-Omics Systems Analysis.

As key inflammatory biomarkers C-reactive protein (CRP) and interleukin-6 (IL6) play an important role in the pathogenesis of non-inflammatory diseases, including specific cancers, such as breast cancer (BC). Previous genome-wide association studies (GWASs) have neither explained the large proportion of genetic heritability nor provided comprehensive understanding of the underlying regulatory mechanisms. We adopted an integrative genomic network approach by incorporating our previous GWAS data for CRP and IL6 with multi-omics datasets, such as whole-blood expression quantitative loci, molecular biologic pathways, and gene regulatory networks to capture the full range of genetic functionalities associated with CRP/IL6 and tissue-specific key drivers (KDs) in gene subnetworks. We applied another systematic genomics approach for BC development to detect shared gene sets in enriched subnetworks across BC and CRP/IL6. We detected the topmost significant common pathways across CRP/IL6 (e.g., immune regulatory; chemokines and their receptors; interferon γ, JAK-STAT, and ERBB4 signaling), several of which overlapped with BC pathways. Further, in gene-gene interaction networks enriched by those topmost pathways, we identified KDs-both well-established (e.g., JAK1/2/3, STAT3) and novel (e.g., CXCR3, CD3D, CD3G, STAT6)-in a tissue-specific manner, for mechanisms shared in regulating CRP/IL6 and BC risk. Our study may provide robust, comprehensive insights into the mechanisms of CRP/IL6 regulation and highlight potential novel genetic targets as preventive and therapeutic strategies for associated disorders, such as BC.

A fast data-driven method for genotype imputation, phasing and local ancestry inference: MendelImpute.jl.

MOTIVATION: Current methods for genotype imputation and phasing exploit the volume of data in haplotype reference panels and rely on hidden Markov models (HMMs). Existing programs all have essentially the same imputation accuracy, are computationally intensive and generally require prephasing the typed markers. RESULTS: We introduce a novel data-mining method for genotype imputation and phasing that substitutes highly efficient linear algebra routines for HMM calculations. This strategy, embodied in our Julia program MendelImpute.jl, avoids explicit assumptions about recombination and population structure while delivering similar prediction accuracy, better memory usage and an order of magnitude or better run-times compared to the fastest competing method. MendelImpute operates on both dosage data and unphased genotype data and simultaneously imputes missing genotypes and phase at both the typed and untyped SNPs (single nucleotide polymorphisms). Finally, MendelImpute naturally extends to global and local ancestry estimation and lends itself to new strategies for data compression and hence faster data transport and sharing. AVAILABILITY AND IMPLEMENTATION: Software, documentation and scripts to reproduce our results are available from https://github.com/OpenMendel/MendelImpute.jl. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genetically determined elevated C-reactive protein associated with primary colorectal cancer risk: Mendelian randomization with lifestyle interactions.

Systemic inflammation-related etiologic pathways via inflammatory cytokines in the development of colorectal cancer (CRC) have not been convincingly determined and may be confounded by lifestyle factors or reverse causality. We investigated the genetically predicted C-reactive protein (CRP) phenotype in the potential causal pathway of primary CRC risk in postmenopausal women in a Mendelian randomization (MR) framework. We employed individual-level data of the Women's Health Initiative Database for Genotypes and Phenotypes Study, which consists of 5 genome-wide association (GWA) studies, including 10,142 women, 737 of whom developed primary CRC. We examined 61 GWA single-nucleotide polymorphisms (SNPs) associated with CRP by using weighted/penalized MR weighted-medians and MR gene-environment interactions that allow some relaxation of the strict variable requirements and attenuate the heterogeneous estimates of outlying SNPs. In lifestyle-stratification analyses, genetically determined CRP exhibited its effects on the decreased CRC risk in non-viscerally obese and high-fat diet subgroups. In contrast, genetically driven CRP was associated with an increased risk for CRC in women who smoked ≥ 15 cigarettes/day, with significant interaction of the gene-smoking relationship. Further, a substantially increased risk of CRC induced by CRP was observed in relatively short-term users (< 5 years) of estrogen (E)-only and also longer-term users (5 to > 10 years) of E plus progestin. Our findings may provide novel evidence on immune-related etiologic pathways connected to CRC risk and suggest the possible use of CRP as a CRC-predictive biomarker in women with particular behaviors and CRP marker-informed interventions to reduce CRC risk.

Modern simulation utilities for genetic analysis.

BACKGROUND: Statistical geneticists employ simulation to estimate the power of proposed studies, test new analysis tools, and evaluate properties of causal models. Although there are existing trait simulators, there is ample room for modernization. For example, most phenotype simulators are limited to Gaussian traits or traits transformable to normality, while ignoring qualitative traits and realistic, non-normal trait distributions. Also, modern computer languages, such as Julia, that accommodate parallelization and cloud-based computing are now mainstream but rarely used in older applications. To meet the challenges of contemporary big studies, it is important for geneticists to adopt new computational tools. RESULTS: We present TraitSimulation, an open-source Julia package that makes it trivial to quickly simulate phenotypes under a variety of genetic architectures. This package is integrated into our OpenMendel suite for easy downstream analyses. Julia was purpose-built for scientific programming and provides tremendous speed and memory efficiency, easy access to multi-CPU and GPU hardware, and to distributed and cloud-based parallelization. TraitSimulation is designed to encourage flexible trait simulation, including via the standard devices of applied statistics, generalized linear models (GLMs) and generalized linear mixed models (GLMMs). TraitSimulation also accommodates many study designs: unrelateds, sibships, pedigrees, or a mixture of all three. (Of course, for data with pedigrees or cryptic relationships, the simulation process must include the genetic dependencies among the individuals.) We consider an assortment of trait models and study designs to illustrate integrated simulation and analysis pipelines. Step-by-step instructions for these analyses are available in our electronic Jupyter notebooks on Github. These interactive notebooks are ideal for reproducible research. CONCLUSION: The TraitSimulation package has three main advantages. (1) It leverages the computational efficiency and ease of use of Julia to provide extremely fast, straightforward simulation of even the most complex genetic models, including GLMs and GLMMs. (2) It can be operated entirely within, but is not limited to, the integrated analysis pipeline of OpenMendel. And finally (3), by allowing a wider range of more realistic phenotype models, TraitSimulation brings power calculations and diagnostic tools closer to what investigators might see in real-world analyses.

Pro-inflammatory cytokine polymorphisms and interactions with dietary alcohol and estrogen, risk factors for invasive breast cancer using a post genome-wide analysis for gene-gene and gene-lifestyle interaction.

Molecular and genetic immune-related pathways connected to breast cancer and lifestyles in postmenopausal women are not fully characterized. In this study, we explored the role of pro-inflammatory cytokines such as C-reactive protein (CRP) and interleukin-6 (IL-6) in those pathways at the genome-wide level. With single-nucleotide polymorphisms (SNPs) in the biomarkers and lifestyles together, we further constructed risk profiles to improve predictability for breast cancer. Our earlier genome-wide association gene-environment interaction study used large cohort data from the Women's Health Initiative Database for Genotypes and Phenotypes Study and identified 88 SNPs associated with CRP and IL-6. For this study, we added an additional 68 SNPs from previous GWA studies, and together with 48 selected lifestyles, evaluated for the association with breast cancer risk via a 2-stage multimodal random survival forest and generalized multifactor dimensionality reduction methods. Overall and in obesity strata (by body mass index, waist, waist-to-hip ratio, exercise, and dietary fat intake), we identified the most predictive genetic and lifestyle variables. Two SNPs (SALL1 rs10521222 and HLA-DQA1 rs9271608) and lifestyles, including alcohol intake, lifetime cumulative exposure to estrogen, and overall and visceral obesity, are the most common and strongest predictive markers for breast cancer across the analyses. The risk profile that combined those variables presented their synergistic effect on the increased breast cancer risk in a gene-lifestyle dose-dependent manner. Our study may contribute to improved predictability for breast cancer and suggest potential interventions for the women with the risk genotypes and lifestyles to reduce their breast cancer risk.

Genome-wide Association Analysis of Proinflammatory Cytokines and Gene-lifestyle Interaction for Invasive Breast Cancer Risk: The WHI dbGaP Study.

Immune-related etiologic pathways to influence invasive breast cancer risk may interact with lifestyle factors, but the interrelated molecular genetic pathways are incompletely characterized. We used data from the Women's Health Initiative Database for Genotypes and Phenotypes Study including 16,088 postmenopausal women, a population highly susceptible to inflammation, obesity, and increased risk for breast cancer. With 21,784,812 common autosomal single-nucleotide polymorphisms (SNP), we conducted a genome-wide association (GWA) gene-environment interaction (G × E) analysis in six independent GWA Studies for proinflammatory cytokines [IL6 and C-reactive protein (CRP)] and their gene-lifestyle interactions. Subsequently, we tested for the association of the GWA SNPs with breast cancer risk. In women overall and stratified by obesity status (body mass index, waist circumference, and waist-to-hip ratio) and obesity-related lifestyle factors (exercise and high-fat diet), 88 GWA SNPs in 10 loci were associated with proinflammatory cytokines: 3 associated with IL6 (1 index SNP in MAPK1 and 1 independent SNP in DEC1); 85 with CRP (3 index SNPs in CRPP1, CRP, RP11-419N10.5, HNF1A-AS1, HNF1A, and C1q2orf43; and two independent SNPs in APOE and APOC1). Of those, 27 in HNF1A-AS1, HNF1A, and C1q2orf43 displayed significantly increased risk for breast cancer. We found a number of novel top markers for CRP and IL6, which interacted with obesity factors. A substantial proportion of those SNPs' susceptibility influenced breast cancer risk. Our findings may contribute to better understanding of genetic associations between pro-inflammation and cancer and suggest intervention strategies for women who carry the risk genotypes, reducing breast cancer risk. PREVENTION RELEVANCE: The top GWA-SNPs associated with pro-inflammatory biomarkers have implications for breast carcinogenesis by interacting with obesity factors. Our findings may suggest interventions for women who carry the inflammatory-risk genotypes to reduce breast cancer risk.

Pro-inflammatory cytokine polymorphisms in ONECUT2 and HNF4A and primary colorectal carcinoma: a post genome-wide gene-lifestyle interaction study.

Immune-related molecular and genetic pathways that are connected to colorectal cancer (CRC) and lifestyles in postmenopausal women are incompletely characterized. In this study, we examined the role of pro-inflammatory biomarkers such as C-reactive protein (CRP) and interleukin-6 (IL-6) in those pathways. Through selection of the best predictive single-nucleotide polymorphisms (SNPs) and lifestyles, our goal was to improve the prediction accuracy and ability for CRC risk. Using large cohort data of postmenopausal women from the Women's Health Initiative Database for Genotypes and Phenotypes Study, we previously conducted a genome-wide association (GWA) for a CRP and IL-6 gene-behavioral interaction study. For the present study, we added GWA-SNPs from outside GWA studies, resulting in a total of 152 SNPs. Together with 41 selected lifestyles, we performed a 2-stage multimodal random survival forest analysis with generalized multifactor dimensionality reduction approach to construct CRC risk profiles. Overall and in obesity strata (by body mass index, waist circumference, waist-to-hip ratio, exercise, and dietary fat intake), we identified the best predictive genetic markers in inflammatory cytokines and lifestyles. Across the strata, 2 SNPs (ONECUT2 rs4092465 and HNF4A rs1800961) and 1 lifestyle factor (relatively short-term past use of oral contraceptives) were the most common and strongest predictive markers for CRC risk. The risk profile that combined those variables exhibited synergistically increased risk for CRC; this pattern appeared more strongly in obese and inactive subgroups. Our results may contribute to improved predictability for CRC and suggest genetically targeted lifestyle interventions for women carrying the inflammatory-risk genotypes, reducing CRC risk.

Mendelian Randomization Study: The Association Between Metabolic Pathways and Colorectal Cancer Risk.

Background: The roles of obesity-related biomarkers and their molecular pathways in the development of postmenopausal colorectal cancer (CRC) have been inconclusive. We examined insulin resistance (IR) as a major hormonal pathway mediating the association between obesity and CRC risk in a Mendelian randomization (MR) framework. Methods: We performed MR analysis using individual-level data of 11,078 non-Hispanic white postmenopausal women from our earlier genome-wide association study. We identified four independent single-nucleotide polymorphisms associated with fasting glucose (FG), three with fasting insulin (FI), and six with homeostatic model assessment-IR (HOMA-IR), which were not associated with obesity. We estimated hazard ratios (HRs) for CRC by adjusting for potential confounding factors plus genetic principal components. Results: Overall, we observed no direct association between combined 13 IR genetic instruments and CRC risk (HR = 0.96, 95% confidence interval [CI]: 0.78-1.17). In phenotypic analysis, genetically raised HOMA-IR exhibited its effects on the increased risk and FG and FI on the reduced risk for CRC, but with a lack of statistical power. Subgroup analyses by physical activity level and dietary fat intake with combined phenotypes showed that genetically determined IR was associated with reduced CRC risk in both physical activity-stratified (single contributor: MTRR rs722025; HR = 0.12, 95% CI: 0.02-0.62) and high-fat diet subgroups (main contributor: G6PC2 rs560887; HR = 0.59, 95% CI: 0.37-0.94). Conclusions: Complex evidence was observed for a potential causal association between IR and CRC risk. Our findings may provide an additional value of intervention trials to lower IR and reduce CRC risk.

Iterative hard thresholding in genome-wide association studies: Generalized linear models, prior weights, and double sparsity.

BACKGROUND: Consecutive testing of single nucleotide polymorphisms (SNPs) is usually employed to identify genetic variants associated with complex traits. Ideally one should model all covariates in unison, but most existing analysis methods for genome-wide association studies (GWAS) perform only univariate regression. RESULTS: We extend and efficiently implement iterative hard thresholding (IHT) for multiple regression, treating all SNPs simultaneously. Our extensions accommodate generalized linear models, prior information on genetic variants, and grouping of variants. In our simulations, IHT recovers up to 30% more true predictors than SNP-by-SNP association testing and exhibits a 2-3 orders of magnitude decrease in false-positive rates compared with lasso regression. We also test IHT on the UK Biobank hypertension phenotypes and the Northern Finland Birth Cohort of 1966 cardiovascular phenotypes. We find that IHT scales to the large datasets of contemporary human genetics and recovers the plausible genetic variants identified by previous studies. CONCLUSIONS: Our real data analysis and simulation studies suggest that IHT can (i) recover highly correlated predictors, (ii) avoid over-fitting, (iii) deliver better true-positive and false-positive rates than either marginal testing or lasso regression, (iv) recover unbiased regression coefficients, (v) exploit prior information and group-sparsity, and (vi) be used with biobank-sized datasets. Although these advances are studied for genome-wide association studies inference, our extensions are pertinent to other regression problems with large numbers of predictors.

Genetically Predicted C-Reactive Protein Associated With Postmenopausal Breast Cancer Risk: Interrelation With Estrogen and Cancer Molecular Subtypes Using Mendelian Randomization.

BACKGROUND: Immune-related etiologic pathways that influence breast cancer risk are incompletely understood and may be confounded by lifestyles or reverse causality. Using a Mendelian randomization (MR) approach, we investigated the potential causal relationship between genetically elevated C-reactive protein (CRP) concentrations and primary invasive breast cancer risk in postmenopausal women. METHODS: We used individual-level data obtained from 10,179 women, including 537 who developed breast cancer, from the Women's Health Initiative Database for Genotypes and Phenotypes Study, which consists of five genome-wide association (GWA) studies. We examined 61 GWA single-nucleotide polymorphisms (SNPs) previously associated with CRP. We employed weighted/penalized weighted-medians and MR gene-environment interactions that allow instruments' invalidity to some extent and attenuate the heterogeneous estimates of outlying SNPs. RESULTS: In lifestyle-stratification analyses, genetically elevated CRP decreased risk for breast cancer in exogenous estrogen-only, estrogen + progestin, and past oral contraceptive (OC) users, but only among relatively short-term users (<5 years). Estrogen-only users for ≥5 years had more profound CRP-decreased breast cancer risk in dose-response fashion, whereas past OC users for ≥5 years had CRP-increased cancer risk. Also, genetically predicted CRP was strongly associated with increased risk for hormone-receptor positive or human epidermal growth factor receptor-2 negative breast cancer. CONCLUSIONS: Our findings may provide novel evidence on the immune-related molecular pathways linking to breast cancer risk and suggest potential clinical use of CRP to predict the specific cancer subtypes. Our findings suggest potential interventions targeting CRP-inflammatory markers to reduce breast cancer risk.

OPENMENDEL: a cooperative programming project for statistical genetics.

Statistical methods for genome-wide association studies (GWAS) continue to improve. However, the increasing volume and variety of genetic and genomic data make computational speed and ease of data manipulation mandatory in future software. In our view, a collaborative effort of statistical geneticists is required to develop open source software targeted to genetic epidemiology. Our attempt to meet this need is called the OPENMENDEL project (https://openmendel.github.io). It aims to (1) enable interactive and reproducible analyses with informative intermediate results, (2) scale to big data analytics, (3) embrace parallel and distributed computing, (4) adapt to rapid hardware evolution, (5) allow cloud computing, (6) allow integration of varied genetic data types, and (7) foster easy communication between clinicians, geneticists, statisticians, and computer scientists. This article reviews and makes recommendations to the genetic epidemiology community in the context of the OPENMENDEL project.

Post Genome-Wide Gene-Environment Interaction Study Using Random Survival Forest: Insulin Resistance, Lifestyle Factors, and Colorectal Cancer Risk.

Molecular and genetic pathways of insulin resistance (IR) connecting colorectal cancer and obesity factors in postmenopausal women remain inconclusive. We examined the IR pathways on both genetic and phenotypic perspectives at the genome-wide level. We further constructed colorectal cancer risk profiles with the most predictive IR SNPs and lifestyle factors. In our earlier genome-wide association gene-environmental interaction study, we used data from a large cohort of postmenopausal women in the Women's Health Initiative Database for Genotypes and Phenotypes Study and identified 58 SNPs in relation to IR phenotypes. In this study, we evaluated the identified IR SNPs and selected 34 lifestyles for their association with colorectal cancer risk in a total of 11,078 women (including 736 women with colorectal cancer) using a 2-stage multimodal random survival forest analysis. In overall and subgroup (defined via body mass index, exercise, and dietary-fat intake) analyses, we identified 2 SNPs (LINC00460 rs1725459 and MTRR rs722025) and lifetime cumulative exposure to estrogen (oral contraceptive use) and cigarette smoking as the most common and strongest predictive markers for colorectal cancer risk across the analyses. The combinations of genetic and lifestyle factors had much greater impact on colorectal cancer risk than any individual risk factors, and a possible synergism existed to increase colorectal cancer risk in a gene-behavior dose-dependent manner. Our findings may inform research on the role of IR in the etiology of colorectal cancer and contribute to more accurate prediction of colorectal cancer risk, suggesting potential intervention strategies for women with specific genotypes and lifestyles to reduce their colorectal cancer risk.

Breast Cancer Risk and Insulin Resistance: Post Genome-Wide Gene-Environment Interaction Study Using a Random Survival Forest.

Obesity-insulin connections have been considered potential risk factors for postmenopausal breast cancer, and the association between insulin resistance (IR) genotypes and phenotypes can be modified by obesity-lifestyle factors, affecting breast cancer risk. In this study, we explored the role of IR in those pathways at the genome-wide level. We identified IR-genetic factors and selected lifestyles to generate risk profiles for postmenopausal breast cancer. Using large-scale cohort data from postmenopausal women in the Women's Health Initiative Database for Genotypes and Phenotypes Study, our previous genome-wide association gene-behavior interaction study identified 58 loci for associations with IR phenotypes (homeostatic model assessment-IR, hyperglycemia, and hyperinsulinemia). We evaluated those single-nucleotide polymorphisms (SNP) and additional 31 lifestyles in relation to breast cancer risk by conducting a two-stage multimodal random survival forest analysis. We identified the most predictive genetic and lifestyle variables in overall and subgroup analyses [stratified by body mass index (BMI), exercise, and dietary fat intake]. Two SNPs (LINC00460 rs17254590 and MKLN1 rs117911989), exogenous factors related to lifetime cumulative exposure to estrogen, BMI, and dietary alcohol consumption were the most common influential factors across the analyses. Individual SNPs did not have significant associations with breast cancer, but SNPs and lifestyles combined synergistically increased the risk of breast cancer in a gene-behavior, dose-dependent manner. These findings may contribute to more accurate predictions of breast cancer and suggest potential intervention strategies for women with specific genetic and lifestyle factors to reduce their breast cancer risk. SIGNIFICANCE: These findings identify insulin resistance SNPs in combination with lifestyle as synergistic factors for breast cancer risk, suggesting lifestyle changes can prevent breast cancer in women who carry the risk genotypes.

Genetic Variants in Metabolic Signaling Pathways and Their Interaction with Lifestyle Factors on Breast Cancer Risk: A Random Survival Forest Analysis.

Genetic variants in the insulin-like growth factor-I (IGF-I)/insulin resistance axis may interact with lifestyle factors, influencing postmenopausal breast cancer risk, but these interrelated pathways are not fully understood. In this study, we examined 54 single-nucleotide polymorphisms (SNP) in genes related to IGF-I/insulin phenotypes and signaling pathways and lifestyle factors in relation to postmenopausal breast cancer, using data from 6,567 postmenopausal women in the Women's Health Initiative Harmonized and Imputed Genome-Wide Association Studies. We used a machine-learning method, two-stage random survival forest analysis. We identified three genetic variants (AKT1 rs2494740, AKT1 rs2494744, and AKT1 rs2498789) and two lifestyle factors [body mass index (BMI) and dietary alcohol intake] as the top five most influential predictors for breast cancer risk. The combination of the three SNPs, BMI, and alcohol consumption (≥1 g/day) significantly increased the risk of breast cancer in a gene and lifestyle dose-dependent manner. Our findings provide insight into gene-lifestyle interactions and will enable researchers to focus on individuals with risk genotypes to promote intervention strategies. These data also suggest potential genetic targets in future intervention/clinical trials for cancer prevention in order to reduce the risk for breast cancer in postmenopausal women. Cancer Prev Res; 11(1); 44-51. ©2017 AACR.

Effect of genetic variants and traits related to glucose metabolism and their interaction with obesity on breast and colorectal cancer risk among postmenopausal women.

BACKGROUND: Impaired glucose metabolism-related genetic variants and traits likely interact with obesity and related lifestyle factors, influencing postmenopausal breast and colorectal cancer (CRC), but their interconnected pathways are not fully understood. By stratifying via obesity and lifestyles, we partitioned the total effect of glucose metabolism genetic variants on cancer risk into two putative mechanisms: 1) indirect (risk-associated glucose metabolism genetic variants mediated by glucose metabolism traits) and 2) direct (risk-associated glucose metabolism genetic variants through pathways other than glucose metabolism traits) effects. METHOD: Using 16 single-nucleotide polymorphisms (SNPs) associated with glucose metabolism and data from 5379 postmenopausal women in the Women's Health Initiative Harmonized and Imputed Genome-Wide Association Studies, we retrospectively assessed the indirect and direct effects of glucose metabolism-traits (fasting glucose, insulin, and homeostatic model assessment-insulin resistance [HOMA-IR]) using two quantitative tests. RESULTS: Several SNPs were associated with breast cancer and CRC risk, and these SNP-cancer associations differed between non-obese and obese women. In both strata, the direct effect of cancer risk associated with the SNP accounted for the majority of the total effect for most SNPs, with roughly 10% of cancer risk due to the SNP that was from an indirect effect mediated by glucose metabolism traits. No apparent differences in the indirect (glucose metabolism-mediated) effects were seen between non-obese and obese women. It is notable that among obese women, 50% of cancer risk was mediated via glucose metabolism trait, owing to two SNPs: in breast cancer, in relation to GCKR through glucose, and in CRC, in relation to DGKB/TMEM195 through HOMA-IR. CONCLUSIONS: Our findings suggest that glucose metabolism genetic variants interact with obesity, resulting in altered cancer risk through pathways other than those mediated by glucose metabolism traits.

Fast Genome-Wide QTL Association Mapping on Pedigree and Population Data.

Since most analysis software for genome-wide association studies (GWAS) currently exploit only unrelated individuals, there is a need for efficient applications that can handle general pedigree data or mixtures of both population and pedigree data. Even datasets thought to consist of only unrelated individuals may include cryptic relationships that can lead to false positives if not discovered and controlled for. In addition, family designs possess compelling advantages. They are better equipped to detect rare variants, control for population stratification, and facilitate the study of parent-of-origin effects. Pedigrees selected for extreme trait values often segregate a single gene with strong effect. Finally, many pedigrees are available as an important legacy from the era of linkage analysis. Unfortunately, pedigree likelihoods are notoriously hard to compute. In this paper, we reexamine the computational bottlenecks and implement ultra-fast pedigree-based GWAS analysis. Kinship coefficients can either be based on explicitly provided pedigrees or automatically estimated from dense markers. Our strategy (a) works for random sample data, pedigree data, or a mix of both; (b) entails no loss of power; (c) allows for any number of covariate adjustments, including correction for population stratification; (d) allows for testing SNPs under additive, dominant, and recessive models; and (e) accommodates both univariate and multivariate quantitative traits. On a typical personal computer (six CPU cores at 2.67 GHz), analyzing a univariate HDL (high-density lipoprotein) trait from the San Antonio Family Heart Study (935,392 SNPs on 1,388 individuals in 124 pedigrees) takes less than 2 min and 1.5 GB of memory. Complete multivariate QTL analysis of the three time-points of the longitudinal HDL multivariate trait takes less than 5 min and 1.5 GB of memory. The algorithm is implemented as the Ped-GWAS Analysis (Option 29) in the Mendel statistical genetics package, which is freely available for Macintosh, Linux, and Windows platforms from http://genetics.ucla.edu/software/mendel.

Genome-wide QTL and eQTL analyses using Mendel.

Pedigree genome-wide association studies (GWAS) (Option 29) in the current version of the Mendel software is an optimized subroutine for performing large-scale genome-wide quantitative trait locus (QTL) analysis. This analysis (a) works for random sample data, pedigree data, or a mix of both; (b) is highly efficient in both run time and memory requirement; (c) accommodates both univariate and multivariate traits; (d) works for autosomal and x-linked loci; (e) correctly deals with missing data in traits, covariates, and genotypes; (f) allows for covariate adjustment and constraints among parameters; (g) uses either theoretical or single nucleotide polymorphism (SNP)-based empirical kinship matrix for additive polygenic effects; (h) allows extra variance components such as dominant polygenic effects and household effects; (i) detects and reports outlier individuals and pedigrees; and (j) allows for robust estimation via the t-distribution. This paper assesses these capabilities on the genetics analysis workshop 19 (GAW19) sequencing data. We analyzed simulated and real phenotypes for both family and random sample data sets. For instance, when jointly testing the 8 longitudinally measured systolic blood pressure and diastolic blood pressure traits, it takes Mendel 78 min on a standard laptop computer to read, quality check, and analyze a data set with 849 individuals and 8.3 million SNPs. Genome-wide expression QTL analysis of 20,643 expression traits on 641 individuals with 8.3 million SNPs takes 30 h using 20 parallel runs on a cluster. Mendel is freely available at http://www.genetics.ucla.edu/software.