Search

Polygenic risk scores (PRSs) quantify the contribution of multiple genetic loci to an individual's likelihood of a complex trait or disease. However, existing PRSs estimate this likelihood with common genetic variants, excluding the impact of rare variants. Here, we report on a method to identify rare variants associated with outlier gene expression and integrate their impact into PRS predictions for body mass index (BMI), obesity, and bariatric surgery. Between the top and bottom 10%, we observed a 20.8% increase in risk for obesity (p = 3 × 10 -14 ), 62.3% increase in risk for severe obesity (p = 1 × 10 -6 ), and median 5.29 years earlier onset for bariatric surgery (p = 0.008), as a function of expression outlier-associated rare variant burden when controlling for common variant PRS. We show that these predictions were more significant than integrating the effects of rare protein-truncating variants (PTVs), observing a mean 19% increase in phenotypic variance explained with expression outlier-associated rare variants when compared with PTVs (p = 2 × 10 -15 ). We replicated these findings by using data from the Million Veteran Program and demonstrated that PRSs across multiple traits and diseases can benefit from the inclusion of expression outlier-associated rare variants identified through population-scale transcriptome sequencing., Competing Interests: Declaration of interests S.B.M. is a consultant for Myome Inc, Tenaya Therapeutics, and BioMarin. C.J.O. is an employee of Novartis Institute for Biomedical Research. S.B.M. and C.S. report a patent application related to this work., (Copyright © 2022 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Whole-genome sequencing studies applied to large populations or biobanks with extensive phenotyping raise new analytic challenges. The need to consider many variants at a locus or group of genes simultaneously and the potential to study many correlated phenotypes with shared genetic architecture provide opportunities for discovery not addressed by the traditional one variant, one phenotype association study. Here, we introduce a Bayesian model comparison approach called MRP (multiple rare variants and phenotypes) for rare-variant association studies that considers correlation, scale, and direction of genetic effects across a group of genetic variants, phenotypes, and studies, requiring only summary statistic data. We apply our method to exome sequencing data (n = 184,698) across 2,019 traits from the UK Biobank, aggregating signals in genes. MRP demonstrates an ability to recover signals such as associations between PCSK9 and LDL cholesterol levels. We additionally find MRP effective in conducting meta-analyses in exome data. Non-biomarker findings include associations between MC1R and red hair color and skin color, IL17RA and monocyte count, and IQGAP2 and mean platelet volume. Finally, we apply MRP in a multi-phenotype setting; after clustering the 35 biomarker phenotypes based on genetic correlation estimates, we find that joint analysis of these phenotypes results in substantial power gains for gene-trait associations, such as in TNFRSF13B in one of the clusters containing diabetes- and lipid-related traits. Overall, we show that the MRP model comparison approach improves upon useful features from widely used meta-analysis approaches for rare-variant association analyses and prioritizes protective modifiers of disease risk., Competing Interests: Declaration of interests M.A.R. is on the SAB of 54Gene and Related Sciences; is scientific founder of Broadwing Bio; and has advised BioMarin, Third Rock Ventures, and MazeTx. C.D.B. is the owner and president of C.D.B. Consulting, LTD, and also a director at EdenRoc Sciences, LLC and BigData Bio LLC (LLC) and Etalon DX; founder of Arc Bio LLC (formerly IdentifyGenomics LLC and BigData Bio LLC); and an SAB member of Imprimed, FaunaBio, Columbia Care, and Digitalis Ventures. He is also a venture partner at F-Prime Capital Partners. M.J.D. is a founder of MazeTx., (Copyright © 2021 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Precise interpretation of the effects of rare protein-truncating variants (PTVs) is important for accurate determination of variant impact. Current methods for assessing the ability of PTVs to induce nonsense-mediated decay (NMD) focus primarily on the position of the variant in the transcript. We used RNA sequencing of the Genotype Tissue Expression v.8 cohort to compute the efficiency of NMD using allelic imbalance for 2,320 rare (genome aggregation database minor allele frequency ≤ 1%) PTVs across 809 individuals in 49 tissues. We created an interpretable predictive model using penalized logistic regression in order to evaluate the comprehensive influence of variant annotation, tissue, and inter-individual variation on NMD. We found that variant position, allele frequency, the inclusion of ultra-rare and singleton variants, and conservation were predictive of allelic imbalance. Furthermore, we found that NMD effects were highly concordant across tissues and individuals. Due to this high consistency, we demonstrate in silico that utilizing peripheral tissues or cell lines provides accurate prediction of NMD for PTVs., (Copyright © 2021 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Follicular lymphomas (FLs) are slow-growing, indolent tumors containing extensive follicular dendritic cell (FDC) networks and recurrent EZH2 gain-of-function mutations. Paradoxically, FLs originate from highly proliferative germinal center (GC) B cells with proliferation strictly dependent on interactions with T follicular helper cells. Herein, we show that EZH2 mutations initiate FL by attenuating GC B cell requirement for T cell help and driving slow expansion of GC centrocytes that become enmeshed with and dependent on FDCs. By impairing T cell help, mutant EZH2 prevents induction of proliferative MYC programs. Thus, EZH2 mutation fosters malignant transformation by epigenetically reprograming B cells to form an aberrant immunological niche that reflects characteristic features of human FLs, explaining how indolent tumors arise from GC B cells., Competing Interests: Declaration of Interests A.M.M. is consulting for Epizyme and Constellation Pharmaceuticals, and receives research funding from Janssen Pharmaceuticals; S.H.K. is consulting for Northrop Grumman; C.E.M. is a co-founder and equity stake holder for Onegevity Health and Biotia, Inc.; C.S. has performed consultancy for Seattle Genetics, Curis Inc., Roche, AbbVie, Juno Therapeutics, and Bayer, and has received research funding from Bristol-Myers Squibb and Trillium Therapeutics, Inc. There are no competing interests., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Population-scale biobanks that combine genetic data and high-dimensional phenotyping for a large number of participants provide an exciting opportunity to perform genome-wide association studies (GWAS) to identify genetic variants associated with diverse quantitative traits and diseases. A major challenge for GWAS in population biobanks is ascertaining disease cases from heterogeneous data sources such as hospital records, digital questionnaire responses, or interviews. In this study, we use genetic parameters, including genetic correlation, to evaluate whether GWAS performed using cases in the UK Biobank ascertained from hospital records, questionnaire responses, and family history of disease implicate similar disease genetics across a range of effect sizes. We find that hospital record and questionnaire GWAS largely identify similar genetic effects for many complex phenotypes and that combining together both phenotyping methods improves power to detect genetic associations. We also show that family history GWAS using cases ascertained on family history of disease agrees with combined hospital record and questionnaire GWAS and that family history GWAS has better power to detect genetic associations for some phenotypes. Overall, this work demonstrates that digital phenotyping and unstructured phenotype data can be combined with structured data such as hospital records to identify cases for GWAS in biobanks and improve the ability of such studies to identify genetic associations., (Copyright © 2020 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Metastasis is the leading cause of cancer mortality. Chromatin remodeling provides the foundation for the cellular reprogramming necessary to drive metastasis. However, little is known about the nature of this remodeling and its regulation. Here, we show that metastasis-inducing pathways regulate histone chaperones to reduce canonical histone incorporation into chromatin, triggering deposition of H3.3 variant at the promoters of poor-prognosis genes and metastasis-inducing transcription factors. This specific incorporation of H3.3 into chromatin is both necessary and sufficient for the induction of aggressive traits that allow for metastasis formation. Together, our data clearly show incorporation of histone variant H3.3 into chromatin as a major regulator of cell fate during tumorigenesis, and histone chaperones as valuable therapeutic targets for invasive carcinomas., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Copy-number variations (CNVs) represent a significant proportion of the genetic differences between individuals and many CNVs associate causally with syndromic disease and clinical outcomes. Here, we characterize the landscape of copy-number variation and their phenome-wide effects in a sample of 472,228 array-genotyped individuals from the UK Biobank. In addition to population-level selection effects against genic loci conferring high mortality, we describe genetic burden from potentially pathogenic and previously uncharacterized CNV loci across more than 3,000 quantitative and dichotomous traits, with separate analyses for common and rare classes of variation. Specifically, we highlight the effects of CNVs at two well-known syndromic loci 16p11.2 and 22q11.2, previously uncharacterized variation at 9p23, and several genic associations in the context of acute coronary artery disease and high body mass index. Our data constitute a deeply contextualized portrait of population-wide burden of copy-number variation, as well as a series of dosage-mediated genic associations across the medical phenome., (Copyright © 2019 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

There is a limited understanding about the impact of rare protein-truncating variants across multiple phenotypes. We explore the impact of this class of variants on 13 quantitative traits and 10 diseases using whole-exome sequencing data from 100,296 individuals. Protein-truncating variants in genes intolerant to this class of mutations increased risk of autism, schizophrenia, bipolar disorder, intellectual disability, and ADHD. In individuals without these disorders, there was an association with shorter height, lower education, increased hospitalization, and reduced age at enrollment. Gene sets implicated from GWASs did not show a significant protein-truncating variants burden beyond what was captured by established Mendelian genes. In conclusion, we provide a thorough investigation of the impact of rare deleterious coding variants on complex traits, suggesting widespread pleiotropic risk., (Copyright © 2018 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Significant insights into disease pathogenesis have been gleaned from population-level genetic studies; however, many loci associated with complex genetic disease contain numerous genes, and phenotypic associations cannot be assigned unequivocally. In particular, a gene-dense locus on chromosome 11 (61.5-61.65 Mb) has been associated with inflammatory bowel disease, rheumatoid arthritis, and coronary artery disease. Here, we identify TMEM258 within this locus as a central regulator of intestinal inflammation. Strikingly, Tmem258 haploinsufficient mice exhibit severe intestinal inflammation in a model of colitis. At the mechanistic level, we demonstrate that TMEM258 is a required component of the oligosaccharyltransferase complex and is essential for N-linked protein glycosylation. Consequently, homozygous deficiency of Tmem258 in colonic organoids results in unresolved endoplasmic reticulum (ER) stress culminating in apoptosis. Collectively, our results demonstrate that TMEM258 is a central mediator of ER quality control and intestinal homeostasis., (Copyright © 2016 The Author(s). Published by Elsevier Inc. All rights reserved.)

Elevated low-density lipoprotein cholesterol (LDL-C) is a treatable, heritable risk factor for cardiovascular disease. Genome-wide association studies (GWASs) have identified 157 variants associated with lipid levels but are not well suited to assess the impact of rare and low-frequency variants. To determine whether rare or low-frequency coding variants are associated with LDL-C, we exome sequenced 2,005 individuals, including 554 individuals selected for extreme LDL-C (>98(th) or <2(nd) percentile). Follow-up analyses included sequencing of 1,302 additional individuals and genotype-based analysis of 52,221 individuals. We observed significant evidence of association between LDL-C and the burden of rare or low-frequency variants in PNPLA5, encoding a phospholipase-domain-containing protein, and both known and previously unidentified variants in PCSK9, LDLR and APOB, three known lipid-related genes. The effect sizes for the burden of rare variants for each associated gene were substantially higher than those observed for individual SNPs identified from GWASs. We replicated the PNPLA5 signal in an independent large-scale sequencing study of 2,084 individuals. In conclusion, this large whole-exome-sequencing study for LDL-C identified a gene not known to be implicated in LDL-C and provides unique insight into the design and analysis of similar experiments., (Copyright © 2014 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

The extent to which variants in the protein-coding sequence of genes contribute to risk of rheumatoid arthritis (RA) is unknown. In this study, we addressed this issue by deep exon sequencing and large-scale genotyping of 25 biological candidate genes located within RA risk loci discovered by genome-wide association studies (GWASs). First, we assessed the contribution of rare coding variants in the 25 genes to the risk of RA in a pooled sequencing study of 500 RA cases and 650 controls of European ancestry. We observed an accumulation of rare nonsynonymous variants exclusive to RA cases in IL2RA and IL2RB (burden test: p = 0.007 and p = 0.018, respectively). Next, we assessed the aggregate contribution of low-frequency and common coding variants to the risk of RA by dense genotyping of the 25 gene loci in 10,609 RA cases and 35,605 controls. We observed a strong enrichment of coding variants with a nominal signal of association with RA (p < 0.05) after adjusting for the best signal of association at the loci (p(enrichment) = 6.4 × 10(-4)). For one locus containing CD2, we found that a missense variant, rs699738 (c.798C>A [p.His266Gln]), and a noncoding variant, rs624988, reside on distinct haplotypes and independently contribute to the risk of RA (p = 4.6 × 10(-6)). Overall, our results indicate that variants (distributed across the allele-frequency spectrum) within the protein-coding portion of a subset of biological candidate genes identified by GWASs contribute to the risk of RA. Further, we have demonstrated that very large sample sizes will be required for comprehensively identifying the independent alleles contributing to the missing heritability of RA., (Copyright © 2013 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Omega-3 and omega-6 long-chain polyunsaturated fatty acids (LC-PUFAs) are essential for the development and function of the human brain. They can be obtained directly from food, e.g., fish, or synthesized from precursor molecules found in vegetable oils. To determine the importance of genetic variability to fatty-acid biosynthesis, we studied FADS1 and FADS2, which encode rate-limiting enzymes for fatty-acid conversion. We performed genome-wide genotyping (n = 5,652 individuals) and targeted resequencing (n = 960 individuals) of the FADS region in five European population cohorts. We also analyzed available genomic data from human populations, archaic hominins, and more distant primates. Our results show that present-day humans have two common FADS haplotypes-defined by 28 closely linked SNPs across 38.9 kb-that differ dramatically in their ability to generate LC-PUFAs. No independent effects on FADS activity were seen for rare SNPs detected by targeted resequencing. The more efficient, evolutionarily derived haplotype appeared after the lineage split leading to modern humans and Neanderthals and shows evidence of positive selection. This human-specific haplotype increases the efficiency of synthesizing essential long-chain fatty acids from precursors and thereby might have provided an advantage in environments with limited access to dietary LC-PUFAs. In the modern world, this haplotype has been associated with lifestyle-related diseases, such as coronary artery disease., (Copyright © 2012 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)