ABSTRACT
Although both short and long sleep duration are associated with elevated hypertension risk, our understanding of their interplay with biological pathways governing blood pressure remains limited. To address this, we carried out genome-wide cross-population gene-by-short-sleep and long-sleep duration interaction analyses for three blood pressure traits (systolic, diastolic, and pulse pressure) in 811,405 individuals from diverse population groups. We discover 22 novel gene-sleep duration interaction loci for blood pressure, mapped to 23 genes. Investigating these genes' functional implications shed light on neurological, thyroidal, bone metabolism, and hematopoietic pathways that necessitate future investigation for blood pressure management that caters to sleep health lifestyle. Non-overlap between short sleep (12) and long sleep (10) interactions underscores the plausible nature of distinct influences of both sleep duration extremes in cardiovascular health. Several of our loci are specific towards a particular population background or sex, emphasizing the importance of addressing heterogeneity entangled in gene-environment interactions, when considering precision medicine design approaches for blood pressure management.
ABSTRACT
Genome-wide association studies (GWAS) for biomarkers important for clinical phenotypes can lead to clinically relevant discoveries. Conventional GWAS for quantitative traits are based on simplified regression models modeling the conditional mean of a phenotype as a linear function of genotype. We draw attention here to an alternative, lesser known approach, namely quantile regression that naturally extends linear regression to the analysis of the entire conditional distribution of a phenotype of interest. Quantile regression can be applied efficiently at biobank scale, while having some unique advantages such as (1) identifying variants with heterogeneous effects across quantiles of the phenotype distribution; (2) accommodating a wide range of phenotype distributions including non-normal distributions, with invariance of results to trait transformations; and (3) providing more detailed information about genotype-phenotype associations even for those associations identified by conventional GWAS. We show in simulations that quantile regression is powerful across both homogeneous and various heterogeneous models. Applications to 39 quantitative traits in the UK Biobank demonstrate that quantile regression can be a helpful complement to linear regression in GWAS and can identify variants with larger effects on high-risk subgroups of individuals but with lower or no contribution overall.
Subject(s)
Biological Specimen Banks , Biomarkers , Genome-Wide Association Study , Genome-Wide Association Study/methods , Humans , Phenotype , Polymorphism, Single Nucleotide , Quantitative Trait Loci , Linear Models , Regression Analysis , Models, Genetic , Genotype , Computer Simulation , United KingdomABSTRACT
Since the first genome-wide association studies (GWASs), thousands of variant-trait associations have been discovered. However, comprehensively mapping the genetic determinant of complex traits through univariate testing can require prohibitive sample sizes. Multi-trait GWAS can circumvent this issue and improve statistical power by leveraging the joint genetic architecture of human phenotypes. Although many methodological hurdles of multi-trait testing have been solved, the strategy to select traits has been overlooked. In this study, we conducted multi-trait GWAS on approximately 20,000 combinations of 72 traits using an omnibus test as implemented in the Joint Analysis of Summary Statistics. We assessed which genetic features of the sets of traits analyzed were associated with an increased detection of variants compared with univariate screening. Several features of the set of traits, including the heritability, the number of traits, and the genetic correlation, drive the multi-trait test gain. Using these features jointly in predictive models captures a large fraction of the power gain of the multi-trait test (Pearson's r between the observed and predicted gain equals 0.43, p < 1.6 × 10-60). Applying an alternative multi-trait approach (Multi-Trait Analysis of GWAS), we identified similar features of interest, but with an overall 70% lower number of new associations. Finally, selecting sets based on our data-driven models systematically outperformed the common strategy of selecting clinically similar traits. This work provides a unique picture of the determinant of multi-trait GWAS statistical power and outlines practical strategies for multi-trait testing.
Subject(s)
Genome-Wide Association Study , Phenotype , Polymorphism, Single Nucleotide , Genome-Wide Association Study/methods , Humans , Polymorphism, Single Nucleotide/genetics , Quantitative Trait Loci/genetics , Models, Genetic , Quantitative Trait, HeritableABSTRACT
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.
Subject(s)
Principal Component Analysis , Humans , Multivariate Analysis , Computational Biology/methods , Phenotype , Algorithms , Genomics/methods , Biomarkers/blood , Computer SimulationABSTRACT
PURPOSE: Excessive dietary sodium intake has known adverse effects on intravascular fluid volume and systemic blood pressure, which may influence intraocular pressure (IOP) and glaucoma risk. This study aimed to assess the association of urinary sodium excretion, a biomarker of dietary intake, with glaucoma and related traits, and determine whether this relationship is modified by genetic susceptibility to disease. DESIGN: Cross-sectional observational and gene-environment interaction analyses in the population-based UK Biobank study. PARTICIPANTS: Up to 103 634 individuals (mean age: 57 years; 51% women) with complete urinary, ocular, and covariable data. METHODS: Urine sodium:creatinine ratio (UNa:Cr; mmol:mmol) was calculated from a midstream urine sample. Ocular parameters were measured as part of a comprehensive eye examination, and glaucoma case ascertainment was through a combination of self-report and linked national hospital records. Genetic susceptibility to glaucoma was calculated based on a glaucoma polygenic risk score comprising 2673 common genetic variants. Multivariable linear and logistic regression, adjusted for key sociodemographic, medical, anthropometric, and lifestyle factors, were used to model associations and gene-environment interactions. MAIN OUTCOME MEASURES: Corneal-compensated IOP, OCT derived macular retinal nerve fiber layer and ganglion cell-inner plexiform layer (GCIPL) thickness, and prevalent glaucoma. RESULTS: In maximally adjusted regression models, a 1 standard deviation increase in UNa:Cr was associated with higher IOP (0.14 mmHg; 95% confidence interval [CI], 0.12-0.17; P < 0.001) and greater prevalence of glaucoma (odds ratio, 1.11; 95% CI, 1.07-1.14; P < 0.001) but not macular retinal nerve fiber layer or ganglion cell-inner plexiform layer thickness. Compared with those with UNa:Cr in the lowest quintile, those in the highest quintile had significantly higher IOP (0.45 mmHg; 95% CI, 0.36-0.53, P < 0.001) and prevalence of glaucoma (odds ratio, 1.30; 95% CI, 1.17-1.45; P < 0.001). Stronger associations with glaucoma (P interaction = 0.001) were noted in participants with a higher glaucoma polygenic risk score. CONCLUSIONS: Urinary sodium excretion, a biomarker of dietary intake, may represent an important modifiable risk factor for glaucoma, especially in individuals at high underlying genetic risk. These findings warrant further investigation because they may have important clinical and public health implications. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.
ABSTRACT
There is a long-standing debate about the magnitude of the contribution of gene-environment interactions to phenotypic variations of complex traits owing to the low statistical power and few reported interactions to date. To address this issue, the Gene-Lifestyle Interactions Working Group within the Cohorts for Heart and Aging Research in Genetic Epidemiology Consortium has been spearheading efforts to investigate G × E in large and diverse samples through meta-analysis. Here, we present a powerful new approach to screen for interactions across the genome, an approach that shares substantial similarity to the Mendelian randomization framework. We identify and confirm 5 loci (6 independent signals) interacted with either cigarette smoking or alcohol consumption for serum lipids, and empirically demonstrate that interaction and mediation are the major contributors to genetic effect size heterogeneity across populations. The estimated lower bound of the interaction and environmentally mediated heritability is significant (P < 0.02) for low-density lipoprotein cholesterol and triglycerides in Cross-Population data. Our study improves the understanding of the genetic architecture and environmental contributions to complex traits.
Subject(s)
Gene-Environment Interaction , Genome-Wide Association Study , Multifactorial Inheritance , Humans , Multifactorial Inheritance/genetics , Male , Triglycerides/blood , Female , Alcohol Drinking/genetics , Polymorphism, Single Nucleotide , Phenotype , Cholesterol, LDL/blood , Cholesterol, LDL/metabolism , Cigarette Smoking/genetics , Quantitative Trait Loci , Middle AgedABSTRACT
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.
Subject(s)
Asian People , Black People , Genome-Wide Association Study , Humans , Asian People/genetics , Black People/genetics , Genetic Predisposition to Disease , Genome-Wide Association Study/methods , Phenotype , Polymorphism, Single Nucleotide , European People/geneticsABSTRACT
Although both short and long sleep duration are associated with elevated hypertension risk, our understanding of their interplay with biological pathways governing blood pressure remains limited. To address this, we carried out genome-wide cross-population gene-by-short-sleep and long-sleep duration interaction analyses for three blood pressure traits (systolic, diastolic, and pulse pressure) in 811,405 individuals from diverse population groups. We discover 22 novel gene-sleep duration interaction loci for blood pressure, mapped to genes involved in neurological, thyroidal, bone metabolism, and hematopoietic pathways. Non-overlap between short sleep (12) and long sleep (10) interactions underscores the plausibility of distinct influences of both sleep duration extremes in cardiovascular health. With several of our loci reflecting specificity towards population background or sex, our discovery sheds light on the importance of embracing granularity when addressing heterogeneity entangled in gene-environment interactions, and in therapeutic design approaches for blood pressure management.
ABSTRACT
Large-scale gene-environment interaction (GxE) discovery efforts often involve compromises in the definition of outcomes and choice of covariates for the sake of data harmonization and statistical power. Consequently, refinement of exposures, covariates, outcomes, and population subsets may be helpful to establish often-elusive replication and evaluate potential clinical utility. Here, we used additional datasets, an expanded set of statistical models, and interrogation of lipoprotein metabolism via nuclear magnetic resonance (NMR)-based lipoprotein subfractions to refine a previously discovered GxE modifying the relationship between physical activity (PA) and HDL-cholesterol (HDL-C). This GxE was originally identified by Kilpeläinen et al., with the strongest cohort-specific signal coming from the Women's Genome Health Study (WGHS). We thus explored this GxE further in the WGHS (N = 23,294), with follow-up in the UK Biobank (UKB; N = 281,380), and the Multi-Ethnic Study of Atherosclerosis (MESA; N = 4,587). Self-reported PA (MET-hrs/wk), genotypes at rs295849 (nearest gene: LHX1), and NMR metabolomics data were available in all three cohorts. As originally reported, minor allele carriers of rs295849 in WGHS had a stronger positive association between PA and HDL-C (pint = 0.002). When testing a range of NMR metabolites (primarily lipoprotein and lipid subfractions) to refine the HDL-C outcome, we found a stronger interaction effect on medium-sized HDL particle concentrations (M-HDL-P; pint = 1.0×10-4) than HDL-C. Meta-regression revealed a systematically larger interaction effect in cohorts from the original meta-analysis with a greater fraction of women (p = 0.018). In the UKB, GxE effects were stronger both in women and using M-HDL-P as the outcome. In MESA, the primary interaction for HDL-C showed nominal significance (pint = 0.013), but without clear differences by sex and with a greater magnitude using large, rather than medium, HDL-P as an outcome. Towards reconciling these observations, further exploration leveraging NMR platform-specific HDL subfraction diameter annotations revealed modest agreement across all cohorts in the interaction affecting medium-to-large particles. Taken together, our work provides additional insights into a specific known gene-PA interaction while illustrating the importance of phenotype and model refinement towards understanding and replicating GxEs.
ABSTRACT
Motivation: Genome-wide association studies (GWAS) have identified thousands of genetic variants associated with common diseases. These results include a mix of causal and non-causal variants related through strong linkage disequilibrium (LD, i.e. highly correlated). Fine-mapping methods have been developed to decipher the causal from non-causal variants using GWAS results and LD information, assigning to each variant a probability of being causal. In this field, the PAINTOR program has become a standard, one of its advantages being its ability to take into account functional annotations. This approach requires many pre- and post-processing steps. Here, we developed a Nextflow pipeline called PaintorPipe that wraps all these steps and the fine-mapping itself together. PaintorPipe uses three independent sources of information: GWAS summary statistics, LD information and functional annotations, to rank the variants according to their susceptibility to be involved in the disease development. The PAINTOR framework is used to calculate the posterior probability of each variant (single nucleotide polymorphism) to be causal (a.k.a. Bayesian fine-mapping). The resulting credible sets of variants are annotated with their biological functions and visualized using CANVIS. This pipeline requires minimal input from users (a GWAS summary statistics file and a set of functional annotation files) and is designed to be modular and customizable, allowing for an easy integration of diverse functional annotations. Availability and implementation: PaintorPipe is implemented in the Nextflow pipeline specific language, can be run locally or on a slurm cluster and handles containerization using Singularity. PaintorPipe is freely available on GitHub (https://github.com/sdjebali/PaintorPipe).
ABSTRACT
Since the first Genome-Wide Association Studies (GWAS), thousands of variant-trait associations have been discovered. However, the sample size required to detect additional variants using standard univariate association screening is increasingly prohibitive. Multi-trait GWAS offers a relevant alternative: it can improve statistical power and lead to new insights about gene function and the joint genetic architecture of human phenotypes. Although many methodological hurdles of multi-trait testing have been discussed, the strategy to select trait, among overwhelming possibilities, has been overlooked. In this study, we conducted extensive multi-trait tests using JASS (Joint Analysis of Summary Statistics) and assessed which genetic features of the analysed sets were associated with an increased detection of variants as compared to univariate screening. Our analyses identified multiple factors associated with the gain in the association detection in multi-trait tests. Together, these factors of the analysed sets are predictive of the gain of the multi-trait test (Pearson's ρ equal to 0.43 between the observed and predicted gain, P < 1.6 × 10-60). Applying an alternative multi-trait approach (MTAG, multi-trait analysis of GWAS), we found that in most scenarios but particularly those with larger numbers of traits, JASS outperformed MTAG. Finally, we benchmark several strategies to select set of traits including the prevalent strategy of selecting clinically similar traits, which systematically underperformed selecting clinically heterogenous traits or selecting sets that issued from our data-driven models. This work provides a unique picture of the determinant of multi-trait GWAS statistical power and outline practical strategies for multi-trait testing.
ABSTRACT
There is a long-standing debate about the magnitude of the contribution of gene-environment interactions to phenotypic variations of complex traits owing to the low statistical power and few reported interactions to date. To address this issue, the CHARGE Gene-Lifestyle Interactions Working Group has been spearheading efforts to investigate G×E in large and diverse samples through meta-analysis. Here, we present a powerful new approach to screen for interactions across the genome, an approach that shares substantial similarity to the Mendelian randomization framework. We identified and confirmed 5 loci (6 independent signals) interacting with either cigarette smoking or alcohol consumption for serum lipids, and empirically demonstrated that interaction and mediation are the major contributors to genetic effect size heterogeneity across populations. The estimated lower bound of the interaction and environmentally mediated contribution ranges from 1.76% to 14.05% of SNP heritability of serum lipids in Cross-Population data. Our study improves the understanding of the genetic architecture and environmental contributions to complex traits.
ABSTRACT
The predictive performance of polygenic scores (PGS) is largely dependent on the number of samples available to train the PGS. Increasing the sample size for a specific phenotype is expensive and takes time, but this sample size can be effectively increased by using genetically correlated phenotypes. We propose a framework to generate multi-PGS from thousands of publicly available genome-wide association studies (GWAS) with no need to individually select the most relevant ones. In this study, the multi-PGS framework increases prediction accuracy over single PGS for all included psychiatric disorders and other available outcomes, with prediction R2 increases of up to 9-fold for attention-deficit/hyperactivity disorder compared to a single PGS. We also generate multi-PGS for phenotypes without an existing GWAS and for case-case predictions. We benchmark the multi-PGS framework against other methods and highlight its potential application to new emerging biobanks.
Subject(s)
Attention Deficit Disorder with Hyperactivity , Genome-Wide Association Study , Humans , Attention Deficit Disorder with Hyperactivity/genetics , Phenotype , Multifactorial Inheritance/geneticsABSTRACT
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.
Subject(s)
Asthma , Pulmonary Disease, Chronic Obstructive , Adult , Humans , Child , Pulmonary Disease, Chronic Obstructive/epidemiology , Pulmonary Disease, Chronic Obstructive/genetics , Asthma/epidemiology , Asthma/genetics , Vital Capacity , Respiratory Function Tests , Forced Expiratory VolumeABSTRACT
Genome-wide association studies (GWAS) for biomarkers important for clinical phenotypes can lead to clinically relevant discoveries. GWAS for quantitative traits are based on simplified regression models modeling the conditional mean of a phenotype as a linear function of genotype. An alternative and easy to apply approach is quantile regression that naturally extends linear regression to the analysis of the entire conditional distribution of a phenotype of interest by modeling conditional quantiles within a regression framework. Quantile regression can be applied efficiently at biobank scale using standard statistical packages in much the same way as linear regression, while having some unique advantages such as identifying variants with heterogeneous effects across different quantiles, including non-additive effects and variants involved in gene-environment interactions; accommodating a wide range of phenotype distributions with invariance to trait transformation; and overall providing more detailed information about the underlying genotype-phenotype associations. Here, we demonstrate the value of quantile regression in the context of GWAS by applying it to 39 quantitative traits in the UK Biobank (n>300,000 individuals). Across these 39 traits we identify 7,297 significant loci, including 259 loci only detected by quantile regression. We show that quantile regression can help uncover replicable but unmodelled gene-environment interactions, and can provide additional key insights into poorly understood genotype-phenotype correlations for clinically relevant biomarkers at minimal additional cost.
ABSTRACT
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.
Subject(s)
Glaucoma , Macula Lutea , Humans , Biological Specimen Banks , Cross-Sectional Studies , Glaucoma/genetics , Intraocular Pressure , Retinal Ganglion Cells , Tomography, Optical Coherence , United Kingdom/epidemiology , Mendelian Randomization AnalysisABSTRACT
Chronic Obstructive Pulmonary Disease (COPD) has a simple physiological diagnostic criterion but a wide range of clinical characteristics. The mechanisms underlying this variability in COPD phenotypes are unclear. To investigate the potential contribution of genetic variants to phenotypic heterogeneity, we examined the association of genome-wide associated lung function, COPD, and asthma variants with other phenotypes using phenome-wide association results derived in the UK Biobank. Our clustering analysis of the variants-phenotypes association matrix identified three clusters of genetic variants with different effects on white blood cell counts, height, and body mass index (BMI). To assess the potential clinical and molecular effects of these groups of variants, we investigated the association between cluster-specific genetic risk scores and phenotypes in the COPDGene cohort. We observed differences in steroid use, BMI, lymphocyte counts, chronic bronchitis, and differential gene and protein expression across the three genetic risk scores. Our results suggest that multi-phenotype analysis of obstructive lung disease-related risk variants may identify genetically driven phenotypic patterns in COPD.
ABSTRACT
Differences in cells' functions arise from differential activity of regulatory elements, including enhancers. Enhancers are cis-regulatory elements that cooperate with promoters through transcription factors to activate the expression of one or several genes by getting physically close to them in the 3D space of the nucleus. There is increasing evidence that genetic variants associated with common diseases are enriched in enhancers active in cell types relevant to these diseases. Identifying the enhancers associated with genes and conversely, the sets of genes activated by each enhancer (the so-called enhancer/gene or E/G relationships) across cell types, can help understanding the genetic mechanisms underlying human diseases. There are three broad approaches for the genome-wide identification of E/G relationships in a cell type: 1) genetic link methods or eQTL, 2) functional link methods based on 1D functional data such as open chromatin, histone mark or gene expression and 3) spatial link methods based on 3D data such as HiC. Since 1) and 3) are costly, the current strategy is to develop functional link methods and to use data from 1) and 3) as reference to evaluate them. However, there is still no consensus on the best functional link method to date, and method comparison remain seldom. Here, we compared the relative performances of three recent methods for the identification of enhancer-gene links, TargetFinder, Average-Rank, and the ABC model, using the three latest benchmarks from the field: a reference that combines 3D and eQTL data, called BENGI, and two genetic screening references, called CRiFF and CRiSPRi. Overall, none of the three methods performed best on the three references. CRiFF and CRISPRi reference sets are likely more reliable, but CRiFF is not genome-wide and CRiFF and CRISPRi are mostly available on the K562 cancer cell line. The BENGI reference set is genome-wide but likely contains many false positives. This study therefore calls for new reliable and genome-wide E/G reference data rather than new functional link E/G identification methods.
ABSTRACT
BACKGROUND: Worldwide, Escherichia coli is the leading cause of neonatal Gram-negative bacterial meningitis, but full understanding of the pathogenesis of this disease is not yet achieved. Moreover, to date, no vaccine is available against bacterial neonatal meningitis. METHODS: Here, we used Transposon Sequencing of saturated banks of mutants (TnSeq) to evaluate E. coli K1 genetic fitness in murine neonatal meningitis. We identified E. coli K1 genes encoding for factors important for systemic dissemination and brain infection, and focused on products with a likely outer-membrane or extra-cellular localization, as these are potential vaccine candidates. We used in vitro and in vivo models to study the efficacy of active and passive immunization. RESULTS: We selected for further study the conserved surface polysaccharide Poly-ß-(1-6)-N-Acetyl Glucosamine (PNAG), as a strong candidate for vaccine development. We found that PNAG was a virulence factor in our animal model. We showed that both passive and active immunization successfully prevented and/or treated meningitis caused by E. coli K1 in neonatal mice. We found an excellent opsonophagocytic killing activity of the antibodies to PNAG and in vitro these antibodies were also able to decrease binding, invasion and crossing of E. coli K1 through two blood brain barrier cell lines. Finally, to reinforce the potential of PNAG as a vaccine candidate in bacterial neonatal meningitis, we demonstrated that Group B Streptococcus, the main cause of neonatal meningitis in developed countries, also produced PNAG and that antibodies to PNAG could protect in vitro and in vivo against this major neonatal pathogen. INTERPRETATION: Altogether, these results indicate the utility of a high-throughput DNA sequencing method to identify potential immunotherapy targets for a pathogen, including in this study a potential broad-spectrum target for prevention of neonatal bacterial infections. FUNDINGS: ANR Seq-N-Vaq, Charles Hood Foundation, Hearst Foundation, and Groupe Pasteur Mutualité.
Subject(s)
Escherichia coli , Meningitis, Bacterial , Animals , Mice , Escherichia coli/genetics , Antibodies, Bacterial , Bacteria/genetics , Immunotherapy , High-Throughput Nucleotide SequencingABSTRACT
PURPOSE: To examine the associations of alcohol consumption with glaucoma and related traits, to assess whether a genetic predisposition to glaucoma modified these associations, and to perform Mendelian randomization (MR) experiments to probe causal effects. 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 intraocular pressure (IOP) (n = 109 097), OCT-derived macular inner retinal layer thickness measures (n = 46 236) and glaucoma status (n = 173 407). METHODS: Participants were categorized according to self-reported drinking behaviors. Quantitative estimates of alcohol intake were derived from touchscreen questionnaires and food composition tables. We performed a 2-step analysis, first comparing categories of alcohol consumption (never, infrequent, regular, and former drinkers) before assessing for a dose-response effect in regular drinkers only. Multivariable linear, logistic, and restricted cubic spline regression, adjusted for key sociodemographic, medical, anthropometric, and lifestyle factors, were used to examine associations. We assessed whether any association was modified by a multitrait glaucoma polygenic risk score. The inverse-variance weighted method was used for the main MR analyses. MAIN OUTCOME MEASURES: Intraocular pressure, macular retinal nerve fiber layer (mRNFL) thickness, macular ganglion cell-inner plexiform layer (mGCIPL) thickness, and prevalent glaucoma. RESULTS: Compared with infrequent drinkers, regular drinkers had higher IOP (+0.17 mmHg; P < 0.001) and thinner mGCIPL (-0.17 µm; P = 0.049), whereas former drinkers had a higher prevalence of glaucoma (odds ratio, 1.53; P = 0.002). In regular drinkers, alcohol intake was adversely associated with all outcomes in a dose-dependent manner (all P < 0.001). Restricted cubic spline regression analyses suggested nonlinear associations, with apparent threshold effects at approximately 50 g (â¼6 UK or 4 US alcoholic units)/week for mRNFL and mGCIPL thickness. Significantly stronger alcohol-IOP associations were observed in participants at higher genetic susceptibility to glaucoma (Pinteraction < 0.001). Mendelian randomization analyses provided evidence for a causal association with mGCIPL thickness. CONCLUSIONS: Alcohol intake was consistently and adversely associated with glaucoma and related traits, and at levels below current United Kingdom (< 112 g/week) and United States (women, < 98 g/week; men, < 196 g/week) guidelines. Although we cannot infer causality definitively, these results will be of interest to people with or at risk of glaucoma and their advising physicians. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found after the references.