@article {99, title = {Dynamic incorporation of multiple in silico functional annotations empowers rare variant association analysis of large whole-genome sequencing studies at scale.}, journal = {Nat Genet}, volume = {52}, year = {2020}, month = {2020 09}, pages = {969-983}, abstract = {

Large-scale whole-genome sequencing studies have enabled the analysis of rare variants (RVs) associated with complex phenotypes. Commonly used RV association tests have limited scope to leverage variant functions. We propose STAAR (variant-set test for association using annotation information), a scalable and powerful RV association test method that effectively incorporates both variant categories and multiple complementary annotations using a dynamic weighting scheme. For the latter, we introduce {\textquoteright}annotation principal components{\textquoteright}, multidimensional summaries of in silico variant annotations. STAAR accounts for population structure and relatedness and is scalable for analyzing very large cohort and biobank whole-genome sequencing studies of continuous and dichotomous traits. We applied STAAR to identify RVs associated with four lipid traits in 12,316 discovery and 17,822 replication samples from the Trans-Omics for Precision Medicine Program. We discovered and replicated new RV associations, including disruptive missense RVs of NPC1L1 and an intergenic region near APOC1P1 associated with low-density lipoprotein cholesterol.

}, keywords = {Cholesterol, LDL, Computer Simulation, Genetic Predisposition to Disease, Genetic Variation, Genome, Genome-Wide Association Study, Humans, Models, Genetic, Molecular Sequence Annotation, Phenotype, Whole Genome Sequencing}, issn = {1546-1718}, doi = {10.1038/s41588-020-0676-4}, author = {Li, Xihao and Li, Zilin and Zhou, Hufeng and Gaynor, Sheila M and Liu, Yaowu and Chen, Han and Sun, Ryan and Dey, Rounak and Arnett, Donna K and Aslibekyan, Stella and Ballantyne, Christie M and Bielak, Lawrence F and Blangero, John and Boerwinkle, Eric and Bowden, Donald W and Broome, Jai G and Conomos, Matthew P and Correa, Adolfo and Cupples, L Adrienne and Curran, Joanne E and Freedman, Barry I and Guo, Xiuqing and Hindy, George and Irvin, Marguerite R and Kardia, Sharon L R and Kathiresan, Sekar and Khan, Alyna T and Kooperberg, Charles L and Laurie, Cathy C and Liu, X Shirley and Mahaney, Michael C and Manichaikul, Ani W and Martin, Lisa W and Mathias, Rasika A and McGarvey, Stephen T and Mitchell, Braxton D and Montasser, May E and Moore, Jill E and Morrison, Alanna C and O{\textquoteright}Connell, Jeffrey R and Palmer, Nicholette D and Pampana, Akhil and Peralta, Juan M and Peyser, Patricia A and Psaty, Bruce M and Redline, Susan and Rice, Kenneth M and Rich, Stephen S and Smith, Jennifer A and Tiwari, Hemant K and Tsai, Michael Y and Vasan, Ramachandran S and Wang, Fei Fei and Weeks, Daniel E and Weng, Zhiping and Wilson, James G and Yanek, Lisa R and Neale, Benjamin M and Sunyaev, Shamil R and Abecasis, Gon{\c c}alo R and Rotter, Jerome I and Willer, Cristen J and Peloso, Gina M and Natarajan, Pradeep and Lin, Xihong} } @article {120, title = {Evidence for secondary-variant genetic burden and non-random distribution across biological modules in a recessive ciliopathy.}, journal = {Nat Genet}, volume = {52}, year = {2020}, month = {2020 11}, pages = {1145-1150}, abstract = {

The influence of genetic background on driver mutations is well established; however, the mechanisms by which the background interacts with Mendelian loci remain unclear. We performed a systematic secondary-variant burden analysis of two independent cohorts of patients with Bardet-Biedl syndrome (BBS) with known recessive biallelic pathogenic mutations in one of 17 BBS genes for each individual. We observed a significant enrichment of trans-acting rare nonsynonymous secondary variants in patients with BBS compared with either population controls or a cohort of individuals with a non-BBS diagnosis and recessive variants in the same gene set. Strikingly, we found a significant over-representation of secondary alleles in chaperonin-encoding genes-a finding corroborated by the observation of epistatic interactions involving this complex in vivo. These data indicate a complex genetic architecture for BBS that informs the biological properties of disease modules and presents a model for secondary-variant burden analysis in recessive disorders.

}, keywords = {Alleles, Bardet-Biedl Syndrome, Cohort Studies, Exome, Genetic Variation, Humans}, issn = {1546-1718}, doi = {10.1038/s41588-020-0707-1}, author = {Kousi, Maria and S{\"o}ylemez, Onuralp and Ozanturk, Ayseg{\"u}l and Mourtzi, Niki and Akle, Sebastian and Jungreis, Irwin and Muller, Jean and Cassa, Christopher A and Brand, Harrison and Mokry, Jill Anne and Wolf, Maxim Y and Sadeghpour, Azita and McFadden, Kelsey and Lewis, Richard A and Talkowski, Michael E and Dollfus, H{\'e}l{\`e}ne and Kellis, Manolis and Davis, Erica E and Sunyaev, Shamil R and Katsanis, Nicholas} } @article {81, title = {Identification of cancer driver genes based on nucleotide context.}, journal = {Nat Genet}, volume = {52}, year = {2020}, month = {2020 02}, pages = {208-218}, abstract = {

Cancer genomes contain large numbers of somatic mutations but few of these mutations drive tumor development. Current approaches either identify driver genes on the basis of mutational recurrence or approximate the functional consequences of nonsynonymous mutations by using bioinformatic scores. Passenger mutations are enriched in characteristic nucleotide contexts, whereas driver mutations occur in functional positions, which are not necessarily surrounded by a particular nucleotide context. We observed that mutations in contexts that deviate from the characteristic contexts around passenger mutations provide a signal in favor of driver genes. We therefore developed a method that combines this feature with the signals traditionally used for driver-gene identification. We applied our method to whole-exome sequencing data from 11,873 tumor-normal pairs and identified 460 driver genes that clustered into 21 cancer-related pathways. Our study provides a resource of driver genes across 28 tumor types with additional driver genes identified according to mutations in unusual nucleotide contexts.

}, keywords = {Cluster Analysis, Computational Biology, Humans, Mutation, Neoplasms, Nucleotides, Proteins, Whole Exome Sequencing}, issn = {1546-1718}, doi = {10.1038/s41588-019-0572-y}, author = {Dietlein, Felix and Weghorn, Donate and Taylor-Weiner, Amaro and Richters, Andr{\'e} and Reardon, Brendan and Liu, David and Lander, Eric S and Van Allen, Eliezer M and Sunyaev, Shamil R} } @article {70, title = {Applicability of the Mutation-Selection Balance Model to Population Genetics of Heterozygous Protein-Truncating Variants in Humans.}, journal = {Mol Biol Evol}, volume = {36}, year = {2019}, month = {2019 08 01}, pages = {1701-1710}, abstract = {

The fate of alleles in the human population is believed to be highly affected by the stochastic force of genetic drift. Estimation of the strength of natural selection in humans generally necessitates a careful modeling of drift including complex effects of the population history and structure. Protein-truncating variants (PTVs) are expected to evolve under strong purifying selection and to have a relatively high per-gene mutation rate. Thus, it is appealing to model the population genetics of PTVs under a simple deterministic mutation-selection balance, as has been proposed earlier (Cassa et~al. 2017). Here, we investigated the limits of this approximation using both computer simulations and data-driven approaches. Our simulations rely on a model of demographic history estimated from 33,370 individual exomes of the Non-Finnish European subset of the ExAC data set (Lek et~al. 2016). Additionally, we compared the African and European subset of the ExAC study and analyzed de novo PTVs. We show that the mutation-selection balance model is applicable to the majority of human genes, but not to genes under the weakest selection.

}, issn = {1537-1719}, doi = {10.1093/molbev/msz092}, author = {Weghorn, Donate and Balick, Daniel J and Cassa, Christopher and Kosmicki, Jack A and Daly, Mark J and Beier, David R and Sunyaev, Shamil R} } @article {42, title = {Error-prone bypass of DNA lesions during lagging-strand replication is a common source of germline and cancer mutations.}, journal = {Nat Genet}, volume = {51}, year = {2019}, month = {2019 Jan}, pages = {36-41}, abstract = {

Studies in experimental systems have identified a multitude of mutational mechanisms including DNA replication infidelity and DNA damage followed by inefficient repair or replicative bypass. However, the relative contributions of these mechanisms to human germline mutation remain unknown. Here, we show that error-prone damage bypass on the lagging strand plays a major role in human mutagenesis. Transcription-coupled DNA repair removes lesions on the transcribed strand; lesions on the non-transcribed strand are preferentially converted into mutations. In human polymorphism we detect a striking similarity between mutation types predominant on the non-transcribed strand and on the strand lagging during replication. Moreover, damage-induced mutations in cancers accumulate asymmetrically with respect to the direction of replication, suggesting that DNA lesions are resolved asymmetrically. We experimentally demonstrate that replication delay greatly attenuates the mutagenic effect of ultraviolet irradiation, confirming that replication converts DNA damage into mutations. We estimate that at least 10\% of human mutations arise due to DNA damage.

}, keywords = {Cells, Cultured, DNA, DNA Damage, DNA Repair, DNA Replication, Germ-Line Mutation, Humans, Mutagenesis, Neoplasms, Polymorphism, Single Nucleotide, Transcription, Genetic}, issn = {1546-1718}, doi = {10.1038/s41588-018-0285-7}, author = {Seplyarskiy, Vladimir B and Akkuratov, Evgeny E and Akkuratova, Natalia and Andrianova, Maria A and Nikolaev, Sergey I and Bazykin, Georgii A and Adameyko, Igor and Sunyaev, Shamil R} } @article {53, title = {Polygenic adaptation on height is overestimated due to uncorrected stratification in genome-wide association studies.}, journal = {Elife}, volume = {8}, year = {2019}, month = {2019 Mar 21}, abstract = {

Genetic predictions of height differ among human populations and these differences have been interpreted as evidence of polygenic adaptation. These differences were first detected using SNPs genome-wide significantly associated with height, and shown to grow stronger when large numbers of sub-significant SNPs were included, leading to excitement about the prospect of analyzing large fractions of the genome to detect polygenic adaptation for multiple traits. Previous studies of height have been based on SNP effect size measurements in the GIANT Consortium meta-analysis. Here we repeat the analyses in the UK Biobank, a much more homogeneously designed study. We show that polygenic adaptation signals based on large numbers of SNPs below genome-wide significance are extremely sensitive to biases due to uncorrected population stratification. More generally, our results imply that typical constructions of polygenic scores are sensitive to population stratification and that population-level differences should be interpreted with caution.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor{\textquoteright}s assessment is that all the issues have been addressed (see decision letter).

}, issn = {2050-084X}, doi = {10.7554/eLife.39702}, author = {Sohail, Mashaal and Maier, Robert M and Ganna, Andrea and Bloemendal, Alex and Martin, Alicia R and Turchin, Michael C and Chiang, Charleston Wk and Hirschhorn, Joel and Daly, Mark J and Patterson, Nick and Neale, Benjamin and Mathieson, Iain and Reich, David and Sunyaev, Shamil R} } @article {48, title = {Quantification of frequency-dependent genetic architectures in 25 UK Biobank traits reveals action of negative selection.}, journal = {Nat Commun}, volume = {10}, year = {2019}, month = {2019 Feb 15}, pages = {790}, abstract = {

Understanding the role of rare variants is important in elucidating the genetic basis of human disease. Negative selection can cause rare variants to have larger per-allele effect sizes than common variants. Here, we develop a method to estimate the minor allele frequency (MAF) dependence of SNP effect sizes. We use a model in which per-allele effect sizes have variance proportional to [p(1 - p)], where p is the MAF and negative values of α imply larger effect sizes for rare variants. We estimate α for 25 UK Biobank diseases and complex traits. All traits produce negative α estimates, with best-fit mean of -0.38 (s.e. 0.02) across traits. Despite larger rare variant effect sizes, rare variants (MAF < 1\%) explain less than 10\% of total SNP-heritability for most traits analyzed. Using evolutionary modeling and forward simulations, we validate the α model of MAF-dependent trait effects and assess plausible values of relevant evolutionary parameters.

}, keywords = {Algorithms, Alleles, Biological Specimen Banks, Gene Frequency, Genome-Wide Association Study, Genotype, Humans, Models, Genetic, Polymorphism, Single Nucleotide, Quantitative Trait, Heritable, Selection, Genetic, United Kingdom}, issn = {2041-1723}, doi = {10.1038/s41467-019-08424-6}, author = {Schoech, Armin P and Jordan, Daniel M and Loh, Po-Ru and Gazal, Steven and O{\textquoteright}Connor, Luke J and Balick, Daniel J and Palamara, Pier F and Finucane, Hilary K and Sunyaev, Shamil R and Price, Alkes L} } @article {11, title = {Estimating the selective effects of heterozygous protein-truncating variants from human exome data.}, journal = {Nat Genet}, volume = {49}, year = {2017}, month = {2017 May}, pages = {806-810}, abstract = {

The evolutionary cost of gene loss is a central question in genetics and has been investigated in model organisms and human cell lines. In humans, tolerance of the loss of one or both functional copies of a gene is related to the gene{\textquoteright}s causal role in disease. However, estimates of the selection and dominance coefficients in humans have been elusive. Here we analyze exome sequence data from 60,706 individuals to make genome-wide estimates of selection against heterozygous loss of gene function. Using this distribution of selection coefficients for heterozygous protein-truncating variants (PTVs), we provide corresponding Bayesian estimates for individual genes. We find that genes under the strongest selection are enriched in embryonic lethal mouse knockouts, Mendelian disease-associated genes, and regulators of transcription. Screening by essentiality, we find a large set of genes under strong selection that are likely to have crucial functions but have not yet been thoroughly characterized.

}, keywords = {Algorithms, Animals, Bayes Theorem, Exome, Gene Frequency, Genetic Predisposition to Disease, Genetic Variation, Genome-Wide Association Study, Genotype, Heterozygote, Humans, Mice, Knockout, Models, Genetic, Mutation, Selection, Genetic, Sequence Analysis, DNA}, issn = {1546-1718}, doi = {10.1038/ng.3831}, author = {Cassa, Christopher A and Weghorn, Donate and Balick, Daniel J and Jordan, Daniel M and Nusinow, David and Samocha, Kaitlin E and O{\textquoteright}Donnell-Luria, Anne and MacArthur, Daniel G and Daly, Mark J and Beier, David R and Sunyaev, Shamil R} } @article {8, title = {Limited statistical evidence for shared genetic effects of eQTLs and autoimmune-disease-associated loci in three major immune-cell types.}, journal = {Nat Genet}, volume = {49}, year = {2017}, month = {2017 Apr}, pages = {600-605}, abstract = {

Most autoimmune-disease-risk effects identified by genome-wide association studies (GWAS) localize to open chromatin with gene-regulatory activity. GWAS loci are also enriched in expression quantitative trait loci (eQTLs), thus suggesting that most risk variants alter gene expression. However, because causal variants are difficult to identify, and cis-eQTLs occur frequently, it remains challenging to identify specific instances of disease-relevant changes to gene regulation. Here, we used a novel joint likelihood framework with higher resolution than that of previous methods to identify loci where autoimmune-disease risk and an eQTL are driven by a single shared genetic effect. Using eQTLs from three major immune subpopulations, we found shared effects in only \~{}25\% of the loci examined. Thus, we show that a fraction of gene-regulatory changes suggest strong mechanistic hypotheses for disease risk, but we conclude that most risk mechanisms are not likely to involve changes in basal gene expression.

}, keywords = {Autoimmune Diseases, Gene Expression, Gene Regulatory Networks, Genetic Predisposition to Disease, Genome-Wide Association Study, Humans, Immunity, Polymorphism, Single Nucleotide, Quantitative Trait Loci}, issn = {1546-1718}, doi = {10.1038/ng.3795}, author = {Chun, Sung and Casparino, Alexandra and Patsopoulos, Nikolaos A and Croteau-Chonka, Damien C and Raby, Benjamin A and De Jager, Philip L and Sunyaev, Shamil R and Cotsapas, Chris} } @article {13, title = {Negative selection in humans and fruit flies involves synergistic epistasis.}, journal = {Science}, volume = {356}, year = {2017}, month = {2017 May 05}, pages = {539-542}, abstract = {

Negative selection against deleterious alleles produced by mutation influences within-population variation as the most pervasive form of natural selection. However, it is not known whether deleterious alleles affect fitness independently, so that cumulative fitness loss depends exponentially on the number of deleterious alleles, or synergistically, so that each additional deleterious allele results in a larger decrease in relative fitness. Negative selection with synergistic epistasis should produce negative linkage disequilibrium between deleterious alleles and, therefore, an underdispersed distribution of the number of deleterious alleles in the genome. Indeed, we detected underdispersion of the number of rare loss-of-function alleles in eight independent data sets from human and fly populations. Thus, selection against rare protein-disrupting alleles is characterized by synergistic epistasis, which may explain how human and fly populations persist despite high genomic mutation rates.

}, keywords = {Alleles, Animals, Drosophila melanogaster, Epistasis, Genetic, Genetic Fitness, Genome, Human, Genome, Insect, Humans, Linkage Disequilibrium, Mutation Rate, Mutation, Missense, Selection, Genetic}, issn = {1095-9203}, doi = {10.1126/science.aah5238}, author = {Sohail, Mashaal and Vakhrusheva, Olga A and Sul, Jae Hoon and Pulit, Sara L and Francioli, Laurent C and van den Berg, Leonard H and Veldink, Jan H and de Bakker, Paul I W and Bazykin, Georgii A and Kondrashov, Alexey S and Sunyaev, Shamil R} }