MANOCCA: a robust and computationally efficient test of covariance in high-dimension multivariate omics data.

Multivariate analysis is becoming central in studies investigating high-throughput molecular data, yet, some important features of these data are seldom explored. Here, we present MANOCCA (Multivariate Analysis of Conditional CovAriance), a powerful method to test for the effect of a predictor on the covariance matrix of a multivariate outcome. The proposed test is by construction orthogonal to tests based on the mean and variance and is able to capture effects that are missed by both approaches. We first compare the performances of MANOCCA with existing correlation-based methods and show that MANOCCA is the only test correctly calibrated in simulation mimicking omics data. We then investigate the impact of reducing the dimensionality of the data using principal component analysis when the sample size is smaller than the number of pairwise covariance terms analysed. We show that, in many realistic scenarios, the maximum power can be achieved with a limited number of components. Finally, we apply MANOCCA to 1000 healthy individuals from the Milieu Interieur cohort, to assess the effect of health, lifestyle and genetic factors on the covariance of two sets of phenotypes, blood biomarkers and flow cytometry-based immune phenotypes. Our analyses identify significant associations between multiple factors and the covariance of both omics data.

Portability of 245 polygenic scores when derived from the UK Biobank and applied to 9 ancestry groups from the same cohort.

The low portability of polygenic scores (PGSs) across global populations is a major concern that must be addressed before PGSs can be used for everyone in the clinic. Indeed, prediction accuracy has been shown to decay as a function of the genetic distance between the training and test cohorts. However, such cohorts differ not only in their genetic distance but also in their geographical distance and their data collection and assaying, conflating multiple factors. In this study, we examine the extent to which PGSs are transferable between ancestries by deriving polygenic scores for 245 curated traits from the UK Biobank data and applying them in nine ancestry groups from the same cohort. By restricting both training and testing to the UK Biobank data, we reduce the risk of environmental and genotyping confounding from using different cohorts. We define the nine ancestry groups at a sub-continental level, based on a simple, robust, and effective method that we introduce here. We then apply two different predictive methods to derive polygenic scores for all 245 phenotypes and show a systematic and dramatic reduction in portability of PGSs trained using Northwestern European individuals and applied to nine ancestry groups. These analyses demonstrate that prediction already drops off within European ancestries and reduces globally in proportion to genetic distance. Altogether, our study provides unique and robust insights into the PGS portability problem.

Multi-trait GWAS for diverse ancestries: mapping the knowledge gap.

BACKGROUND: Approximately 95% of samples analyzed in univariate genome-wide association studies (GWAS) are of European ancestry. This bias toward European ancestry populations in association screening also exists for other analyses and methods that are often developed and tested on European ancestry only. However, existing data in non-European populations, which are often of modest sample size, could benefit from innovative approaches as recently illustrated in the context of polygenic risk scores. METHODS: Here, we extend and assess the potential limitations and gains of our multi-trait GWAS pipeline, JASS (Joint Analysis of Summary Statistics), for the analysis of non-European ancestries. To this end, we conducted the joint GWAS of 19 hematological traits and glycemic traits across five ancestries (European (EUR), admixed American (AMR), African (AFR), East Asian (EAS), and South-East Asian (SAS)). RESULTS: We detected 367 new genome-wide significant associations in non-European populations (15 in Admixed American (AMR), 72 in African (AFR) and 280 in East Asian (EAS)). New associations detected represent 5%, 17% and 13% of associations in the AFR, AMR and EAS populations, respectively. Overall, multi-trait testing increases the replication of European associated loci in non-European ancestry by 15%. Pleiotropic effects were highly similar at significant loci across ancestries (e.g. the mean correlation between multi-trait genetic effects of EUR and EAS ancestries was 0.88). For hematological traits, strong discrepancies in multi-trait genetic effects are tied to known evolutionary divergences: the ARKC1 loci, which is adaptive to overcome p.vivax induced malaria. CONCLUSIONS: Multi-trait GWAS can be a valuable tool to narrow the genetic knowledge gap between European and non-European populations.

Powerful gene-based testing by integrating long-range chromatin interactions and knockoff genotypes.

Gene-based tests are valuable techniques for identifying genetic factors in complex traits. Here, we propose a gene-based testing framework that incorporates data on long-range chromatin interactions, several recent technical advances for region-based tests, and leverages the knockoff framework for synthetic genotype generation for improved gene discovery. Through simulations and applications to genome-wide association studies (GWAS) and whole-genome sequencing data for multiple diseases and traits, we show that the proposed test increases the power over state-of-the-art gene-based tests in the literature, identifies genes that replicate in larger studies, and can provide a more narrow focus on the possible causal genes at a locus by reducing the confounding effect of linkage disequilibrium. Furthermore, our results show that incorporating genetic variation in distal regulatory elements tends to improve power over conventional tests. Results for UK Biobank and BioBank Japan traits are also available in a publicly accessible database that allows researchers to query gene-based results in an easy fashion.

Multitrait GWAS to connect disease variants and biological mechanisms.

Genome-wide association studies (GWASs) have uncovered a wealth of associations between common variants and human phenotypes. Here, we present an integrative analysis of GWAS summary statistics from 36 phenotypes to decipher multitrait genetic architecture and its link with biological mechanisms. Our framework incorporates multitrait association mapping along with an investigation of the breakdown of genetic associations into clusters of variants harboring similar multitrait association profiles. Focusing on two subsets of immunity and metabolism phenotypes, we then demonstrate how genetic variants within clusters can be mapped to biological pathways and disease mechanisms. Finally, for the metabolism set, we investigate the link between gene cluster assignment and the success of drug targets in randomized controlled trials.

Polygenic risk scores identify heterogeneity in asthma and chronic obstructive pulmonary disease.

BACKGROUND: Asthma and chronic obstructive pulmonary disease (COPD) have distinct and overlapping genetic and clinical features. OBJECTIVE: We sought to test the hypothesis that polygenic risk scores (PRSs) for asthma (PRSAsthma) and spirometry (FEV1 and FEV1/forced vital capacity; PRSspiro) would demonstrate differential associations with asthma, COPD, and asthma-COPD overlap (ACO). METHODS: We developed and tested 2 asthma PRSs and applied the higher performing PRSAsthma and a previously published PRSspiro to research (Genetic Epidemiology of COPD study and Childhood Asthma Management Program, with spirometry) and electronic health record-based (Mass General Brigham Biobank and Genetic Epidemiology Research on Adult Health and Aging [GERA]) studies. We assessed the association of PRSs with COPD and asthma using modified random-effects and binary-effects meta-analyses, and ACO and asthma exacerbations in specific cohorts. Models were adjusted for confounders and genetic ancestry. RESULTS: In meta-analyses of 102,477 participants, the PRSAsthma (odds ratio [OR] per SD, 1.16 [95% CI, 1.14-1.19]) and PRSspiro (OR per SD, 1.19 [95% CI, 1.17-1.22]) both predicted asthma, whereas the PRSspiro predicted COPD (OR per SD, 1.25 [95% CI, 1.21-1.30]). However, results differed by cohort. The PRSspiro was not associated with COPD in GERA and Mass General Brigham Biobank. In the Genetic Epidemiology of COPD study, the PRSAsthma (OR per SD: Whites, 1.3; African Americans, 1.2) and PRSspiro (OR per SD: Whites, 2.2; African Americans, 1.6) were both associated with ACO. In GERA, the PRSAsthma was associated with asthma exacerbations (OR, 1.18) in Whites; the PRSspiro was associated with asthma exacerbations in White, LatinX, and East Asian participants. CONCLUSIONS: PRSs for asthma and spirometry are both associated with ACO and asthma exacerbations. Genetic prediction performance differs in research versus electronic health record-based cohorts.

The Association of Physical Activity with Glaucoma and Related Traits in the UK Biobank.

PURPOSE: To examine the association of physical activity (PA) with glaucoma and related traits, to assess whether genetic predisposition to glaucoma modified these associations, and to probe causal relationships using Mendelian randomization (MR). DESIGN: Cross-sectional observational and gene-environment interaction analyses in the UK Biobank. Two-sample MR experiments using summary statistics from large genetic consortia. PARTICIPANTS: UK Biobank participants with data on self-reported or accelerometer-derived PA and intraocular pressure (IOP; n = 94 206 and n = 27 777, respectively), macular inner retinal OCT measurements (n = 36 274 and n = 9991, respectively), and glaucoma status (n = 86 803 and n = 23 556, respectively). METHODS: We evaluated multivariable-adjusted associations of self-reported (International Physical Activity Questionnaire) and accelerometer-derived PA with IOP and macular inner retinal OCT parameters using linear regression and with glaucoma status using logistic regression. For all outcomes, we examined gene-PA interactions using a polygenic risk score (PRS) that combined the effects of 2673 genetic variants associated with glaucoma. MAIN OUTCOME MEASURES: Intraocular pressure, macular retinal nerve fiber layer (mRNFL) thickness, macular ganglion cell-inner plexiform layer (mGCIPL) thickness, and glaucoma status. RESULTS: In multivariable-adjusted regression models, we found no association of PA level or time spent in PA with glaucoma status. Higher overall levels and greater time spent in higher levels of both self-reported and accelerometer-derived PA were associated positively with thicker mGCIPL (P < 0.001 for trend for each). Compared with the lowest quartile of PA, participants in the highest quartiles of accelerometer-derived moderate- and vigorous-intensity PA showed a thicker mGCIPL by +0.57 µm (P < 0.001) and +0.42 µm (P = 0.005). No association was found with mRNFL thickness. High overall level of self-reported PA was associated with a modestly higher IOP of +0.08 mmHg (P = 0.01), but this was not replicated in the accelerometry data. No associations were modified by a glaucoma PRS, and MR analyses did not support a causal relationship between PA and any glaucoma-related outcome. CONCLUSIONS: Higher overall PA level and greater time spent in moderate and vigorous PA were not associated with glaucoma status but were associated with thicker mGCIPL. Associations with IOP were modest and inconsistent. Despite the well-documented acute reduction in IOP after PA, we found no evidence that high levels of habitual PA are associated with glaucoma status or IOP in the general population. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Ultrarare heterozygous pathogenic variants of genes causing dominant forms of early-onset deafness underlie severe presbycusis.

Presbycusis, or age-related hearing loss (ARHL), is a major public health issue. About half the phenotypic variance has been attributed to genetic factors. Here, we assessed the contribution to presbycusis of ultrarare pathogenic variants, considered indicative of Mendelian forms. We focused on severe presbycusis without environmental or comorbidity risk factors and studied multiplex family age-related hearing loss (mARHL) and simplex/sporadic age-related hearing loss (sARHL) cases and controls with normal hearing by whole-exome sequencing. Ultrarare variants (allele frequency [AF] < 0.0001) of 35 genes responsible for autosomal dominant early-onset forms of deafness, predicted to be pathogenic, were detected in 25.7% of mARHL and 22.7% of sARHL cases vs. 7.5% of controls (P = 0.001); half were previously unknown (AF < 0.000002). MYO6, MYO7A, PTPRQ, and TECTA variants were present in 8.9% of ARHL cases but less than 1% of controls. Evidence for a causal role of variants in presbycusis was provided by pathogenicity prediction programs, documented haploinsufficiency, three-dimensional structure/function analyses, cell biology experiments, and reported early effects. We also established Tmc1N321I/+ mice, carrying the TMC1:p.(Asn327Ile) variant detected in an mARHL case, as a mouse model for a monogenic form of presbycusis. Deafness gene variants can thus result in a continuum of auditory phenotypes. Our findings demonstrate that the genetics of presbycusis is shaped by not only well-studied polygenic risk factors of small effect size revealed by common variants but also, ultrarare variants likely resulting in monogenic forms, thereby paving the way for treatment with emerging inner ear gene therapy.

Fitting Gaussian mixture models on incomplete data.

BACKGROUND: Bioinformatics investigators often gain insights by combining information across multiple and disparate data sets. Merging data from multiple sources frequently results in data sets that are incomplete or contain missing values. Although missing data are ubiquitous, existing implementations of Gaussian mixture models (GMMs) either cannot accommodate missing data, or do so by imposing simplifying assumptions that limit the applicability of the model. In the presence of missing data, a standard ad hoc practice is to perform complete case analysis or imputation prior to model fitting. Both approaches have serious drawbacks, potentially resulting in biased and unstable parameter estimates. RESULTS: Here we present missingness-aware Gaussian mixture models (MGMM), an R package for fitting GMMs in the presence of missing data. Unlike existing GMM implementations that can accommodate missing data, MGMM places no restrictions on the form of the covariance matrix. Using three case studies on real and simulated 'omics data sets, we demonstrate that, when the underlying data distribution is near-to a GMM, MGMM is more effective at recovering the true cluster assignments than either the existing GMM implementations that accommodate missing data, or fitting a standard GMM after state of the art imputation. Moreover, MGMM provides an accurate assessment of cluster assignment uncertainty, even when the generative distribution is not a GMM. CONCLUSION: Compared to state-of-the-art competitors, MGMM demonstrates a better ability to recover the true cluster assignments for a wide variety of data sets and a large range of missingness rates. MGMM provides the bioinformatics community with a powerful, easy-to-use, and statistically sound tool for performing clustering and density estimation in the presence of missing data. MGMM is publicly available as an R package on CRAN: https://CRAN.R-project.org/package=MGMM .

Making the Most of Clumping and Thresholding for Polygenic Scores.

Polygenic prediction has the potential to contribute to precision medicine. Clumping and thresholding (C+T) is a widely used method to derive polygenic scores. When using C+T, several p value thresholds are tested to maximize predictive ability of the derived polygenic scores. Along with this p value threshold, we propose to tune three other hyper-parameters for C+T. We implement an efficient way to derive thousands of different C+T scores corresponding to a grid over four hyper-parameters. For example, it takes a few hours to derive 123K different C+T scores for 300K individuals and 1M variants using 16 physical cores. We find that optimizing over these four hyper-parameters improves the predictive performance of C+T in both simulations and real data applications as compared to tuning only the p value threshold. A particularly large increase can be noted when predicting depression status, from an AUC of 0.557 (95% CI: [0.544-0.569]) when tuning only the p value threshold to an AUC of 0.592 (95% CI: [0.580-0.604]) when tuning all four hyper-parameters we propose for C+T. We further propose stacked clumping and thresholding (SCT), a polygenic score that results from stacking all derived C+T scores. Instead of choosing one set of hyper-parameters that maximizes prediction in some training set, SCT learns an optimal linear combination of all C+T scores by using an efficient penalized regression. We apply SCT to eight different case-control diseases in the UK biobank data and find that SCT substantially improves prediction accuracy with an average AUC increase of 0.035 over standard C+T.

The Association between Serum Lipids and Intraocular Pressure in 2 Large United Kingdom Cohorts.

PURPOSE: Serum lipids are modifiable, routinely collected blood test features associated with cardiovascular health. We examined the association of commonly collected serum lipid measures (total cholesterol [TC], high-density lipoprotein cholesterol [HDL-C], low-density lipoprotein cholesterol [LDL-C], and triglycerides) with intraocular pressure (IOP). DESIGN: Cross-sectional study in the UK Biobank and European Prospective Investigation into Cancer and Nutrition (EPIC)-Norfolk cohorts. PARTICIPANTS: We included 94 323 participants from the UK Biobank (mean age, 57 years) and 6230 participants from the EPIC-Norfolk (mean age, 68 years) cohorts with data on TC, HDL-C, LDL-C, and triglycerides collected between 2006 and 2009. METHODS: Multivariate linear regression adjusting for demographic, lifestyle, anthropometric, medical, and ophthalmic covariables was used to examine the associations of serum lipids with corneal-compensated IOP (IOPcc). MAIN OUTCOME MEASURES: Corneal-compensated IOP. RESULTS: Higher levels of TC, HDL-C, and LDL-C were associated independently with higher IOPcc in both cohorts after adjustment for key demographic, medical, and lifestyle factors. For each 1-standard deviation increase in TC, HDL-C, and LDL-C, IOPcc was higher by 0.09 mmHg (95% confidence interval [CI], 0.06-0.11 mmHg; P < 0.001), 0.11 mmHg (95% CI, 0.08-0.13 mmHg; P < 0.001), and 0.07 mmHg (95% CI, 0.05-0.09 mmHg; P < 0.001), respectively, in the UK Biobank cohort. In the EPIC-Norfolk cohort, each 1-standard deviation increase in TC, HDL-C, and LDL-C was associated with a higher IOPcc by 0.19 mmHg (95% CI, 0.07-0.31 mmHg; P = 0.001), 0.14 mmHg (95% CI, 0.03-0.25 mmHg; P = 0.016), and 0.17 mmHg (95% CI, 0.06-0.29 mmHg; P = 0.003). An inverse association between triglyceride levels and IOP in the UK Biobank (-0.05 mmHg; 95% CI, -0.08 to -0.03; P < 0.001) was not replicated in the EPIC-Norfolk cohort (P = 0.30). CONCLUSIONS: Our findings suggest that serum TC, HDL-C, and LDL-C are associated positively with IOP in 2 United Kingdom cohorts and that triglyceride levels may be associated negatively. Future research is required to assess whether these associations are causal in nature.

Multi-ancestry genome-wide gene-sleep interactions identify novel loci for blood pressure.

Long and short sleep duration are associated with elevated blood pressure (BP), possibly through effects on molecular pathways that influence neuroendocrine and vascular systems. To gain new insights into the genetic basis of sleep-related BP variation, we performed genome-wide gene by short or long sleep duration interaction analyses on four BP traits (systolic BP, diastolic BP, mean arterial pressure, and pulse pressure) across five ancestry groups in two stages using 2 degree of freedom (df) joint test followed by 1df test of interaction effects. Primary multi-ancestry analysis in 62,969 individuals in stage 1 identified three novel gene by sleep interactions that were replicated in an additional 59,296 individuals in stage 2 (stage 1 + 2 Pjoint < 5 × 10-8), including rs7955964 (FIGNL2/ANKRD33) that increases BP among long sleepers, and rs73493041 (SNORA26/C9orf170) and rs10406644 (KCTD15/LSM14A) that increase BP among short sleepers (Pint < 5 × 10-8). Secondary ancestry-specific analysis identified another novel gene by long sleep interaction at rs111887471 (TRPC3/KIAA1109) in individuals of African ancestry (Pint = 2 × 10-6). Combined stage 1 and 2 analyses additionally identified significant gene by long sleep interactions at 10 loci including MKLN1 and RGL3/ELAVL3 previously associated with BP, and significant gene by short sleep interactions at 10 loci including C2orf43 previously associated with BP (Pint < 10-3). 2df test also identified novel loci for BP after modeling sleep that has known functions in sleep-wake regulation, nervous and cardiometabolic systems. This study indicates that sleep and primary mechanisms regulating BP may interact to elevate BP level, suggesting novel insights into sleep-related BP regulation.

Gene-educational attainment interactions in a multi-ancestry genome-wide meta-analysis identify novel blood pressure loci.

Educational attainment is widely used as a surrogate for socioeconomic status (SES). Low SES is a risk factor for hypertension and high blood pressure (BP). To identify novel BP loci, we performed multi-ancestry meta-analyses accounting for gene-educational attainment interactions using two variables, "Some College" (yes/no) and "Graduated College" (yes/no). Interactions were evaluated using both a 1 degree of freedom (DF) interaction term and a 2DF joint test of genetic and interaction effects. Analyses were performed for systolic BP, diastolic BP, mean arterial pressure, and pulse pressure. We pursued genome-wide interrogation in Stage 1 studies (N = 117 438) and follow-up on promising variants in Stage 2 studies (N = 293 787) in five ancestry groups. Through combined meta-analyses of Stages 1 and 2, we identified 84 known and 18 novel BP loci at genome-wide significance level (P < 5 × 10-8). Two novel loci were identified based on the 1DF test of interaction with educational attainment, while the remaining 16 loci were identified through the 2DF joint test of genetic and interaction effects. Ten novel loci were identified in individuals of African ancestry. Several novel loci show strong biological plausibility since they involve physiologic systems implicated in BP regulation. They include genes involved in the central nervous system-adrenal signaling axis (ZDHHC17, CADPS, PIK3C2G), vascular structure and function (GNB3, CDON), and renal function (HAS2 and HAS2-AS1, SLIT3). Collectively, these findings suggest a role of educational attainment or SES in further dissection of the genetic architecture of BP.

Genetic effects on the commensal microbiota in inflammatory bowel disease patients.

Several bacteria in the gut microbiota have been shown to be associated with inflammatory bowel disease (IBD), and dozens of IBD genetic variants have been identified in genome-wide association studies. However, the role of the microbiota in the etiology of IBD in terms of host genetic susceptibility remains unclear. Here, we studied the association between four major genetic variants associated with an increased risk of IBD and bacterial taxa in up to 633 IBD cases. We performed systematic screening for associations, identifying and replicating associations between NOD2 variants and two taxa: the Roseburia genus and the Faecalibacterium prausnitzii species. By exploring the overall association patterns between genes and bacteria, we found that IBD risk alleles were significantly enriched for associations concordant with bacteria-IBD associations. To understand the significance of this pattern in terms of the study design and known effects from the literature, we used counterfactual principles to assess the fitness of a few parsimonious gene-bacteria-IBD causal models. Our analyses showed evidence that the disease risk of these genetic variants were likely to be partially mediated by the microbiome. We confirmed these results in extensive simulation studies and sensitivity analyses using the association between NOD2 and F. prausnitzii as a case study.

A multi-ancestry genome-wide study incorporating gene-smoking interactions identifies multiple new loci for pulse pressure and mean arterial pressure.

Elevated blood pressure (BP), a leading cause of global morbidity and mortality, is influenced by both genetic and lifestyle factors. Cigarette smoking is one such lifestyle factor. Across five ancestries, we performed a genome-wide gene-smoking interaction study of mean arterial pressure (MAP) and pulse pressure (PP) in 129 913 individuals in stage 1 and follow-up analysis in 480 178 additional individuals in stage 2. We report here 136 loci significantly associated with MAP and/or PP. Of these, 61 were previously published through main-effect analysis of BP traits, 37 were recently reported by us for systolic BP and/or diastolic BP through gene-smoking interaction analysis and 38 were newly identified (P < 5 × 10-8, false discovery rate < 0.05). We also identified nine new signals near known loci. Of the 136 loci, 8 showed significant interaction with smoking status. They include CSMD1 previously reported for insulin resistance and BP in the spontaneously hypertensive rats. Many of the 38 new loci show biologic plausibility for a role in BP regulation. SLC26A7 encodes a chloride/bicarbonate exchanger expressed in the renal outer medullary collecting duct. AVPR1A is widely expressed, including in vascular smooth muscle cells, kidney, myocardium and brain. FHAD1 is a long non-coding RNA overexpressed in heart failure. TMEM51 was associated with contractile function in cardiomyocytes. CASP9 plays a central role in cardiomyocyte apoptosis. Identified only in African ancestry were 30 novel loci. Our findings highlight the value of multi-ancestry investigations, particularly in studies of interaction with lifestyle factors, where genomic and lifestyle differences may contribute to novel findings.

A Large-Scale Multi-ancestry Genome-wide Study Accounting for Smoking Behavior Identifies Multiple Significant Loci for Blood Pressure.

Genome-wide association analysis advanced understanding of blood pressure (BP), a major risk factor for vascular conditions such as coronary heart disease and stroke. Accounting for smoking behavior may help identify BP loci and extend our knowledge of its genetic architecture. We performed genome-wide association meta-analyses of systolic and diastolic BP incorporating gene-smoking interactions in 610,091 individuals. Stage 1 analysis examined ∼18.8 million SNPs and small insertion/deletion variants in 129,913 individuals from four ancestries (European, African, Asian, and Hispanic) with follow-up analysis of promising variants in 480,178 additional individuals from five ancestries. We identified 15 loci that were genome-wide significant (p < 5 × 10-8) in stage 1 and formally replicated in stage 2. A combined stage 1 and 2 meta-analysis identified 66 additional genome-wide significant loci (13, 35, and 18 loci in European, African, and trans-ancestry, respectively). A total of 56 known BP loci were also identified by our results (p < 5 × 10-8). Of the newly identified loci, ten showed significant interaction with smoking status, but none of them were replicated in stage 2. Several loci were identified in African ancestry, highlighting the importance of genetic studies in diverse populations. The identified loci show strong evidence for regulatory features and support shared pathophysiology with cardiometabolic and addiction traits. They also highlight a role in BP regulation for biological candidates such as modulators of vascular structure and function (CDKN1B, BCAR1-CFDP1, PXDN, EEA1), ciliopathies (SDCCAG8, RPGRIP1L), telomere maintenance (TNKS, PINX1, AKTIP), and central dopaminergic signaling (MSRA, EBF2).

Genetic meta-analysis of cancer diagnosis following statin use identifies new associations and implicates human leukocyte antigen (HLA) in women.

We sought to perform a genomic evaluation of the risk of incident cancer in statin users, free of cancer at study entry. Patients who previously participated in two phase IV trials (TNT and IDEAL) with genetic data were used (npooled = 11,196). A GWAS meta-analysis using Cox modeling for the prediction of incident cancer was conducted in the pooled cohort and sex-stratified. rs13210472 (near HLA-DOA gene) was associated with higher risk of incident cancer amongst women with prevalent coronary artery disease (CAD) taking statins (hazard ratio [HR]: 2.66, 95% confidence interval [CI]: 1.88-3.76, P = 3.5 × 10-8). Using the UK Biobank and focusing exclusively on women statin users with CAD (nfemale = 2952), rs13210472 remained significantly associated with incident cancer (HR: 1.71, 95% CI: 1.14-2.56, P = 9.0 × 10-3). The association was not observed in non-statin users. In this genetic meta-analysis, we have identified a variant in women statin users with prevalent CAD that was associated with incident cancer, possibly implicating the human leukocyte antigen pathway.

Intraocular Pressure, Glaucoma, and Dietary Caffeine Consumption: A Gene-Diet Interaction Study from the UK Biobank.

PURPOSE: We examined the association of habitual caffeine intake with intraocular pressure (IOP) and glaucoma and whether genetic predisposition to higher IOP modified these associations. We also assessed whether genetic predisposition to higher coffee consumption was related to IOP. DESIGN: Cross-sectional study in the UK Biobank. PARTICIPANTS: We included 121 374 participants (baseline ages, 39-73 years) with data on coffee and tea intake (collected 2006-2010) and corneal-compensated IOP measurements in 2009. In a subset of 77 906 participants with up to 5 web-based 24-hour-recall food frequency questionnaires (2009-2012), we evaluated total caffeine intake. We also assessed the same relationships with glaucoma (9286 cases and 189 763 controls). METHODS: We evaluated multivariable-adjusted associations with IOP using linear regression and with glaucoma using logistic regression. For both outcomes, we examined gene-diet interactions using a polygenic risk score (PRS) that combined the effects of 111 genetic variants associated with IOP. We also performed Mendelian randomization using 8 genetic variants associated with coffee intake to assess potential causal effects of coffee consumption on IOP. MAIN OUTCOME MEASURES: Intraocular pressure and glaucoma. RESULTS: Mendelian randomization analysis did not support a causal effect of coffee drinking on IOP (P > 0.1). Greater caffeine intake was associated weakly with lower IOP: the highest (≥232 mg/day) versus lowest (<87 mg/day) caffeine consumption was associated with a 0.10-mmHg lower IOP (Ptrend = 0.01). However, the IOP PRS modified this association: among those in the highest IOP PRS quartile, consuming > 480 mg/day versus < 80 mg/day was associated with a 0.35-mmHg higher IOP (Pinteraction = 0.01). The relationship between caffeine intake and glaucoma was null (P ≥ 0.1). However, the IOP PRS also modified this relationship: compared with those in the lowest IOP PRS quartile consuming no caffeine, those in the highest IOP PRS quartile consuming ≥ 321 mg/day showed a 3.90-fold higher glaucoma prevalence (Pinteraction = 0.0003). CONCLUSIONS: Habitual caffeine consumption was associated weakly with lower IOP, and the association between caffeine consumption and glaucoma was null. However, among participants with the strongest genetic predisposition to elevated IOP, greater caffeine consumption was associated with higher IOP and higher glaucoma prevalence.

Deriving stratified effects from joint models investigating gene-environment interactions.

BACKGROUND: Models including an interaction term and performing a joint test of SNP and/or interaction effect are often used to discover Gene-Environment (GxE) interactions. When the environmental exposure is a binary variable, analyses from exposure-stratified models which consist of estimating genetic effect in unexposed and exposed individuals separately can be of interest. In large-scale consortia focusing on GxE interactions in which only the joint test has been performed, it may be challenging to get summary statistics from both exposure-stratified and marginal (i.e not accounting for interaction) models. RESULTS: In this work, we developed a simple framework to estimate summary statistics in each stratum of a binary exposure and in the marginal model using summary statistics from the "joint" model. We performed simulation studies to assess our estimators' accuracy and examined potential sources of bias, such as correlation between genotype and exposure and differing phenotypic variances within exposure strata. Results from these simulations highlight the high theoretical accuracy of our estimators and yield insights into the impact of potential sources of bias. We then applied our methods to real data and demonstrate our estimators' retained accuracy after filtering SNPs by sample size to mitigate potential bias. CONCLUSIONS: These analyses demonstrated the accuracy of our method in estimating both stratified and marginal summary statistics from a joint model of gene-environment interaction. In addition to facilitating the interpretation of GxE screenings, this work could be used to guide further functional analyses. We provide a user-friendly Python script to apply this strategy to real datasets. The Python script and documentation are available at https://gitlab.pasteur.fr/statistical-genetics/j2s.