Search

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.

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.

MotivationCurrent 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.ResultsWe 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 implementationSoftware, documentation and scripts to reproduce our results are available from https://github.com/OpenMendel/MendelImpute.jl.Supplementary informationSupplementary data are available at Bioinformatics online.

BackgroundStatistical 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.ResultsWe 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.ConclusionThe 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.

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.

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.

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.

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.

BackgroundDisparities 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.MethodsWith 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.ResultsIn 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.ConclusionsOur 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.

BackgroundImmune-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.MethodsWe 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.ResultsIn 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 (

Abstract: 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,170 women, including 537 who developed breast cancer, from the Women’s Health Initiative Database for Genotypes and Phenotypes Study, which consists of 5 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 pathways linking to breast cancer risk, and suggest potential clinical use of CRP to predict the specific cancer subtypes and CRP–inflammatory marker targeting interventions to reduce postmenopausal breast cancer risk.

BackgroundConsecutive 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.ResultsWe 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.ConclusionsOur 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.

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.

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.

BackgroundImmune-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.MethodsWe 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.ResultsIn 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 (

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.

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.

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 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.

BackgroundImpaired 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.MethodUsing 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.ResultsSeveral 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.ConclusionsOur 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.

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.

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.

PurposeImpaired glucose metabolism-related genetic variants likely interact with obesity-modifiable factors in response to glucose intolerance, yet their interconnected pathways have not been fully characterized.MethodsWith data from 1,027 postmenopausal participants of the Genomics and Randomized Trials Network study and 15 single-nucleotide polymorphisms (SNPs) associated with glucose homeostasis, we assessed whether obesity, physical activity, and high dietary fat intake interact with the SNP-glucose variations. We used regression analysis plus stratification and graphic approaches.ResultsAcross carriers of the 15 SNPs, fasting levels of glucose, insulin, and homeostatic model assessment-insulin resistance (HOMA-IR) were higher in obese, inactive, and high fat-diet women than in their respective counterparts. Carriers within subgroups differently demonstrated the direction and/or magnitude of the variants' effect on glucose-relevant traits. Variants in GCKR, GCK, DGKB/TMEM195 (P for interactions = 0.02, 0.02, and 0.01), especially, showed interactions with obesity: obese, inactive, and high fat-diet women had greater increases in fasting glucose, insulin, and HOMA-IR levels. Obese carriers at TCF7L2 variant had greater increases in fasting glucose levels than nonobese carriers (P for interaction = 0.04), whereas active women had greater decreases in insulin and HOMA-IR levels than inactive women (P for interaction = 0.02 in both levels).ConclusionsOur data support the important role of obesity in modifying glucose homeostasis in response to glucose metabolism-relevant variants. These findings may inform research on the role of glucose homeostasis in the etiology of chronic disease and the development of intervention strategies to reduce risk in postmenopausal women.

Background/aimsMaternal and offspring cell contact at the site of placentation presents a plausible setting for maternal-fetal genotype (MFG) interactions affecting fetal growth. We test hypotheses regarding killer cell immunoglobulin-like receptor (KIR) and HLA-C MFG effects on human birth weight by extending the quantitative MFG (QMFG) test.MethodsUntil recently, association testing for MFG interactions had limited applications. To improve the ability to test for these interactions, we developed the extended QMFG test, a linear mixed-effect model that can use multi-locus genotype data from families.ResultsWe demonstrate the extended QMFG test's statistical properties. We also show that if an offspring-only model is fit when MFG effects exist, associations can be missed or misattributed. Furthermore, imprecisely modeling the effects of both KIR and HLA-C could result in a failure to replicate if these loci's allele frequencies differ among populations. To further illustrate the extended QMFG test's advantages, we apply the extended QMFG test to a UK cohort study and the Norwegian Mother and Child Cohort (MoBa) study.ConclusionWe find a significant KIR-HLA-C interaction effect on birth weight. More generally, the QMFG test can detect genetic associations that may be missed by standard genome-wide association studies for quantitative traits.

Maternal-offspring gene interactions, aka maternal-fetal genotype (MFG) incompatibilities, are neglected in complex diseases and quantitative trait studies. They are implicated in birth to adult onset diseases but there are limited ways to investigate their influence on quantitative traits. We present the quantitative-MFG (QMFG) test, a linear mixed model where maternal and offspring genotypes are fixed effects and residual correlations between family members are random effects. The QMFG handles families of any size, common or general scenarios of MFG incompatibility, and additional covariates. We develop likelihood ratio tests (LRTs) and rapid score tests and show they provide correct inference. In addition, the LRT's alternative model provides unbiased parameter estimates. We show that testing the association of SNPs by fitting a standard model, which only considers the offspring genotypes, has very low power or can lead to incorrect conclusions. We also show that offspring genetic effects are missed if the MFG modeling assumptions are too restrictive. With genome-wide association study data from the San Antonio Family Heart Study, we demonstrate that the QMFG score test is an effective and rapid screening tool. The QMFG test therefore has important potential to identify pathways of complex diseases for which the genetic etiology remains to be discovered.

BackgroundAlthough cardiovascular disease (CVD) and type 2 diabetes mellitus (T2D) share many common risk factors, potential molecular mechanisms that may also be shared for these 2 disorders remain unknown.Methods and resultsUsing an integrative pathway and network analysis, we performed genome-wide association studies in 8155 blacks, 3494 Hispanic American, and 3697 Caucasian American women who participated in the national Women's Health Initiative single-nucleotide polymorphism (SNP) Health Association Resource and the Genomics and Randomized Trials Network. Eight top pathways and gene networks related to cardiomyopathy, calcium signaling, axon guidance, cell adhesion, and extracellular matrix seemed to be commonly shared between CVD and T2D across all 3 ethnic groups. We also identified ethnicity-specific pathways, such as cell cycle (specific for Hispanic American and Caucasian American) and tight junction (CVD and combined CVD and T2D in Hispanic American). In network analysis of gene-gene or protein-protein interactions, we identified key drivers that included COL1A1, COL3A1, and ELN in the shared pathways for both CVD and T2D. These key driver genes were cross-validated in multiple mouse models of diabetes mellitus and atherosclerosis.ConclusionsOur integrative analysis of American women of 3 ethnicities identified multiple shared biological pathways and key regulatory genes for the development of CVD and T2D. These prospective findings also support the notion that ethnicity-specific susceptibility genes and process are involved in the pathogenesis of CVD and T2D.

The linkage era left a rich legacy of pedigree samples that can be used for modern genome-wide association sequencing (GWAS) or next-generation sequencing (NGS) studies. Family designs are naturally equipped to detect rare variants, control for population stratification, and facilitate the study of parent-of-origin effects. Unfortunately, pedigree likelihoods are notoriously hard to compute, and current software for association mapping in pedigrees is prohibitively slow in processing dense marker maps. In a recent release of the comprehensive genetic analysis software MENDEL, we implemented an ultra-fast score test for association mapping with pedigree-based GWAS or NGS study data. Our implementation (a) works for random sample data, pedigree data, or a mix of both;(b) allows for covariate adjustment, including correction for population stratification;(c) accommodates both univariate and multivariate quantitative traits; and (d) allows missing values in multivariate traits. In this paper, we assess the capabilities of MENDEL on the Genetic Analysis Workshop 18 sequencing data. For instance, when jointly testing the 4 longitudinally measured diastolic blood pressure traits, it takes MENDEL less than 51 minutes on a standard laptop computer to read, quality check, and analyze a data set with 959 individuals and 8.3 million single-nucleotide polymorphisms (SNPs). Our analysis reveals association of one SNP in the q32.2 region of chromosome 1. MENDEL is freely available on http://www.genetics.ucla.edu/software.

Statistical genetics is undergoing the same transition to big data that all branches of applied statistics are experiencing. With the advent of inexpensive DNA sequencing, the transition is only accelerating. This brief review highlights some modern techniques with recent successes in statistical genetics. These include: (a) lasso penalized regression and association mapping, (b) ethnic admixture estimation, (c) matrix completion for genotype and sequence data, (d) the fused lasso and copy number variation, (e) haplotyping, (f) estimation of relatedness, (g) variance components models, and (h) rare variant testing. For more than a century, genetics has been both a driver and beneficiary of statistical theory and practice. This symbiotic relationship will persist for the foreseeable future.

The Systems Genetics Resource (SGR) (http://systems.genetics.ucla.edu) is a new open-access web application and database that contains genotypes and clinical and intermediate phenotypes from both human and mouse studies. The mouse data include studies using crosses between specific inbred strains and studies using the Hybrid Mouse Diversity Panel. SGR is designed to assist researchers studying genes and pathways contributing to complex disease traits, including obesity, diabetes, atherosclerosis, heart failure, osteoporosis, and lipoprotein metabolism. Over the next few years, we hope to add data relevant to deafness, addiction, hepatic steatosis, toxin responses, and vascular injury. The intermediate phenotypes include expression array data for a variety of tissues and cultured cells, metabolite levels, and protein levels. Pre-computed tables of genetic loci controlling intermediate and clinical phenotypes, as well as phenotype correlations, are accessed via a user-friendly web interface. The web site includes detailed protocols for all of the studies. Data from published studies are freely available; unpublished studies have restricted access during their embargo period.

Although genetics is the most significant known determinant of human intelligence, specific gene contributions remain largely unknown. To accelerate understanding in this area, we have taken a new approach by studying the relationship between quantitative gene expression and intelligence in a cohort of 65 patients with Williams Syndrome (WS), a neurodevelopmental disorder caused by a 1.5 Mb deletion on chromosome 7q11.23. We find that variation in the transcript levels of the brain gene STX1A correlates significantly with intelligence in WS patients measured by principal component analysis (PCA) of standardized WAIS-R subtests, r = 0.40 (Pearson correlation, Bonferroni corrected p-value = 0.007), accounting for 15.6% of the cognitive variation. These results suggest that syntaxin 1A, a neuronal regulator of presynaptic vesicle release, may play a role in WS and be a component of the cellular pathway determining human intelligence.

BackgroundSystems biologic approaches such as Weighted Gene Co-expression Network Analysis (WGCNA) can effectively integrate gene expression and trait data to identify pathways and candidate biomarkers. Here we show that the additional inclusion of genetic marker data allows one to characterize network relationships as causal or reactive in a chronic fatigue syndrome (CFS) data set.ResultsWe combine WGCNA with genetic marker data to identify a disease-related pathway and its causal drivers, an analysis which we refer to as "Integrated WGCNA" or IWGCNA. Specifically, we present the following IWGCNA approach: 1) construct a co-expression network, 2) identify trait-related modules within the network, 3) use a trait-related genetic marker to prioritize genes within the module, 4) apply an integrated gene screening strategy to identify candidate genes and 5) carry out causality testing to verify and/or prioritize results. By applying this strategy to a CFS data set consisting of microarray, SNP and clinical trait data, we identify a module of 299 highly correlated genes that is associated with CFS severity. Our integrated gene screening strategy results in 20 candidate genes. We show that our approach yields biologically interesting genes that function in the same pathway and are causal drivers for their parent module. We use a separate data set to replicate findings and use Ingenuity Pathways Analysis software to functionally annotate the candidate gene pathways.ConclusionWe show how WGCNA can be combined with genetic marker data to identify disease-related pathways and the causal drivers within them. The systems genetics approach described here can easily be used to generate testable genetic hypotheses in other complex disease studies.

BackgroundCorrectly merged data sets that have been independently genotyped can increase statistical power in linkage and association studies. However, alleles from microsatellite data sets genotyped with different experimental protocols or platforms cannot be accurately matched using base-pair size information alone. In a previous publication we introduced a statistical model for merging microsatellite data by matching allele frequencies between data sets. These methods are implemented in our software MicroMerge version 1 (v1). While MicroMerge v1 output can be analyzed by some genetic analysis programs, many programs can not analyze alignments that do not match alleles one-to-one between data sets. A consequence of such alignments is that codominant genotypes must often be analyzed as phenotypes. In this paper we describe several extensions that are implemented in MicroMerge version 2 (v2).ResultsNotably, MicroMerge v2 includes a new one-to-one alignment option that creates merged pedigree and locus files that can be handled by most genetic analysis software. Other features in MicroMerge v2 enhance the following aspects of control: 1) optimizing the algorithm for different merging scenarios, such as data sets with very different sample sizes or multiple data sets, 2) merging small data sets when a reliable set of allele frequencies are available, and 3) improving the quantity and 4) quality of merged data. We present results from simulated and real microsatellite genotype data sets, and conclude with an association analysis of three familial dyslipidemia (FD) study samples genotyped at different laboratories. Independent analysis of each FD data set did not yield consistent results, but analysis of the merged data sets identified strong association at locus D11S2002.ConclusionThe MicroMerge v2 features will enable merging for a variety of genotype data sets, which in turn will facilitate meta-analyses for powering association analysis.

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. [ABSTRACT FROM AUTHOR]