Search

Genetic variants that inactivate protein-coding genes are a powerful source of information about the phenotypic consequences of gene disruption: genes that are crucial for the function of an organism will be depleted of such variants in natural populations, whereas non-essential genes will tolerate their accumulation. However, predicted loss-of-function variants are enriched for annotation errors, and tend to be found at extremely low frequencies, so their analysis requires careful variant annotation and very large sample sizes1. Here we describe the aggregation of 125,748 exomes and 15,708 genomes from human sequencing studies into the Genome Aggregation Database (gnomAD). We identify 443,769 high-confidence predicted loss-of-function variants in this cohort after filtering for artefacts caused by sequencing and annotation errors. Using an improved model of human mutation rates, we classify human protein-coding genes along a spectrum that represents tolerance to inactivation, validate this classification using data from model organisms and engineered human cells, and show that it can be used to improve the power of gene discovery for both common and rare diseases., A catalogue of predicted loss-of-function variants in 125,748 whole-exome and 15,708 whole-genome sequencing datasets from the Genome Aggregation Database (gnomAD) reveals the spectrum of mutational constraints that affect these human protein-coding genes.

Human genetic variants predicted to cause loss-of-function of protein-coding genes (pLoF variants) provide natural in vivo models of human gene inactivation, and can be valuable indicators of gene function and the potential toxicity of therapeutic inhibitors targeting these genes1,2. Gain-of-kinase-function variants in LRRK2 are known to significantly increase the risk of Parkinson’s disease3,4, suggesting that inhibition of LRRK2 kinase activity is a promising therapeutic strategy. While preclinical studies in model organisms have raised some on-target toxicity concerns5–8, the biological consequences of LRRK2 inhibition have not been well-characterized in humans. Here we systematically analyse pLoF variants in LRRK2 observed across 141,456 individuals sequenced in the Genome Aggregation Database (gnomAD)9, 49,960 exome sequenced individuals from the UK Biobank, and over 4 million participants in the 23andMe genotyped dataset. After stringent variant curation, we identify 1,455 individuals with high-confidence pLoF variants in LRRK2. Experimental validation of three variants, combined with prior work10, confirmed reduced protein levels in 82.5% of our cohort. We show that heterozygous pLoF variants in LRRK2 reduce LRRK2 protein levels but are not strongly associated with any specific phenotype or disease state. Our results demonstrate the value of large-scale genomic databases and phenotyping of human LoF carriers for target validation in drug discovery.

Naturally occurring human genetic variants predicted to inactivate protein-coding genes provide an in vivo model of human gene inactivation that complements cell and model organism knockout studies. Here we report three key findings regarding assessment of candidate drug targets using human loss-of-function variants. First, even essential genes, where loss-of-function variants are not tolerated, can be highly successful as targets of inhibitory drugs. Second, in most genes, loss-of-function variants are sufficiently rare that genotype-based ascertainment of homozygous or compound heterozygous “knockout” humans will await sample sizes ~1,000 times those presently available, unless recruitment focuses on consanguineous individuals. Third, automated variant annotation and filtering are powerful, but manual curation remains critical for removing artifacts, and is a prerequisite for recall-by-genotype efforts. Our results provide a roadmap for human “knockout” studies and should guide interpretation of loss-of-function variants in drug development.

Cranial growth and development is a complex process which affects the closely related traits of head circumference (HC) and intracranial volume (ICV). The underlying genetic influences shaping these traits during the transition from childhood to adulthood are little understood, but might include both age-specific genetic factors and low-frequency genetic variation. Here, we model the developmental genetic architecture of HC, showing this is genetically stable and correlated with genetic determinants of ICV. Investigating up to 46,000 children and adults of European descent, we identify association with final HC and/or final ICV + HC at 9 novel common and low-frequency loci, illustrating that genetic variation from a wide allele frequency spectrum contributes to cranial growth. The largest effects are reported for low-frequency variants within TP53, with 0.5 cm wider heads in increaser-allele carriers versus non-carriers during mid-childhood, suggesting a previously unrecognized role of TP53 transcripts in human cranial development.

The clinical application of advanced next-generation sequencing technologies is increasingly uncovering novel classes of mutations that may serve as potential targets for precision medicine therapeutics. Here, we show that a deep intronic splice defect in the COL6A1 gene, originally discovered by applying muscle RNA sequencing in patients with clinical findings of collagen VI-related dystrophy (COL6-RD), inserts an in-frame pseudoexon into COL6A1 mRNA, encodes a mutant collagen a1(VI) protein that exerts a dominant-negative effect on collagen VI matrix assembly, and provides a unique opportunity for splice-correction approaches aimed at restoring normal gene expression. Using splice-modulating antisense oligomers, we efficiently skipped the pseudoexon in patient-derived fibroblast cultures and restored a wild-type matrix. Similarly, we used CRISPR/Cas9 to precisely delete an intronic sequence containing the pseudoexon and efficiently abolish its inclusion while preserving wild-type splicing. Considering that this splice defect is emerging as one of the single most frequent mutations in COL6-RD, the design of specific and effective splice-correction therapies offers a promising path for clinical translation.

Objective: Comprehensive clinical characterisation of congenital titinopathy to facilitate diagnosis and management of this important emerging disorder. Methods: Using massively parallel sequencing we identified 30 patients from 27 families with two pathogenic nonsense, frameshift and/or splice site TTN mutations in trans. We then undertook a detailed analysis of the clinical, histopathology and imaging features of these patients. Results: All patients had prenatal‐ or early‐onset hypotonia and/or congenital contractures. None had ophthalmoplegia. Scoliosis and respiratory insufficiency typically developed early and progressed rapidly, whereas limb weakness was often slowly progressive, and usually did not prevent independent walking. Cardiac involvement was present in 46% of patients. Relatives of two patients had dilated cardiomyopathy. Creatine kinase levels were normal to moderately elevated. Increased fibre size variation, internalised nuclei and cores were common histopathological abnormalities. Cap‐like structures, whorled or ring fibres, and mitochondrial accumulations were also observed. Muscle MRI showed gluteal, hamstring and calf muscle involvement. Western blot analysis showed a near‐normal sized titin protein in all samples. The presence of two mutations predicted to impact both N2BA and N2B cardiac isoforms appeared to be associated with greatest risk of cardiac involvement. One third of patients had one mutation predicted to impact exons present in fetal skeletal muscle, but not included within the mature skeletal muscle isoform transcript. This strongly suggests developmental isoforms are involved in the pathogenesis of this congenital/early‐onset disorder. Interpretation: This detailed clinical reference dataset will greatly facilitate diagnostic confirmation and management of patients and has provided important insights into disease pathogenesis.

Hand grip strength is a widely used proxy of muscular fitness, a marker of frailty, and predictor of a range of morbidities and all-cause mortality. To investigate the genetic determinants of variation in grip strength, we perform a large-scale genetic discovery analysis in a combined sample of 195,180 individuals and identify 16 loci associated with grip strength (P, Hand grip strength as a proxy of muscular fitness is a clinical predictor of mortality and morbidity. In a large-scale GWA study, the authors find 16 robustly associated genetic loci that highlight roles in muscle fibre structure and function, neuronal maintenance and nervous system signal transduction.

Glycosylphosphatidylinositol (GPI)-anchored proteins are ubiquitously expressed in the human body and are important for various functions at the cell surface. Mutations in many GPI biosynthesis genes have been described to date in patients with multi-system disease and together these constitute a subtype of congenital disorders of glycosylation. We used whole exome sequencing in two families to investigate the genetic basis of disease and used RNA and cellular studies to investigate the functional consequences of sequence variants in the PIGY gene. Two families with different phenotypes had homozygous recessive sequence variants in the GPI biosynthesis gene PIGY. Two sisters with c.137T>C (p.Leu46Pro) PIGY variants had multi-system disease including dysmorphism, seizures, severe developmental delay, cataracts and early death. There were significantly reduced levels of GPI-anchored proteins (CD55 and CD59) on the surface of patient-derived skin fibroblasts (∼20-50% compared with controls). In a second, consanguineous family, two siblings had moderate development delay and microcephaly. A homozygous PIGY promoter variant (c.-540G>A) was detected within a 7.7 Mb region of autozygosity. This variant was predicted to disrupt a SP1 consensus binding site and was shown to be associated with reduced gene expression. Mutations in PIGY can occur in coding and non-coding regions of the gene and cause variable phenotypes. This article contributes to understanding of the range of disease phenotypes and disease genes associated with deficiencies of the GPI-anchor biosynthesis pathway and also serves to highlight the potential importance of analysing variants detected in 5'-UTR regions despite their typically low coverage in exome data.

Interpreting variants, especially noncoding ones, in the increasing number of personal genomes is challenging. We used patterns of polymorphisms in functionally annotated regions in 1092 humans to identify deleterious variants; then we experimentally validated candidates. We analyzed both coding and noncoding regions, with the former corroborating the latter. We found regions particularly sensitive to mutations ("ultrasensitive") and variants that are disruptive because of mechanistic effects on transcription-factor binding (that is, "motif-breakers"). We also found variants in regions with higher network centrality tend to be deleterious. Insertions and deletions followed a similar pattern to single-nucleotide variants, with some notable exceptions (e.g., certain deletions and enhancers). On the basis of these patterns, we developed a computational tool (FunSeq), whose application to ~90 cancer genomes reveals nearly a hundred candidate noncoding drivers.

By characterizing the geographic and functional spectrum of human genetic variation, the 1000 Genomes Project aims to build a resource to help to understand the genetic contribution to disease. Here we describe the genomes of 1,092 individuals from 14 populations, constructed using a combination of low-coverage whole-genome and exome sequencing. By developing methods to integrate information across several algorithms and diverse data sources, we provide a validated haplotype map of 38 million single nucleotide polymorphisms, 1.4 million short insertions and deletions, and more than 14,000 larger deletions. We show that individuals from different populations carry different profiles of rare and common variants, and that low-frequency variants show substantial geographic differentiation, which is further increased by the action of purifying selection. We show that evolutionary conservation and coding consequence are key determinants of the strength of purifying selection, that rare-variant load varies substantially across biological pathways, and that each individual contains hundreds of rare non-coding variants at conserved sites, such as motif-disrupting changes in transcription-factor-binding sites. This resource, which captures up to 98% of accessible single nucleotide polymorphisms at a frequency of 1% in related populations, enables analysis of common and low-frequency variants in individuals from diverse, including admixed, populations., National Institutes of Health (U.S.) (Grant RC2HL102925), National Institutes of Health (U.S.) (Grant U54HG3067)

Hand grip strength is a widely used proxy of muscular fitness, a marker of frailty, and predictor of a range of morbidities and all-cause mortality. To investigate the genetic determinants of variation in grip strength, we perform a large-scale genetic discovery analysis in a combined sample of 195,180 individuals and identify 16 loci associated with grip strength (P

A major goal of biomedicine is to understand the function of every gene in the human genome. Loss-of-function mutations can disrupt both copies of a given gene in humans and phenotypic analysis of such 'human knockouts' can provide insight into gene function. Consanguineous unions are more likely to result in offspring carrying homozygous loss-of-function mutations. In Pakistan, consanguinity rates are notably high. Here we sequence the protein-coding regions of 10,503 adult participants in the Pakistan Risk of Myocardial Infarction Study (PROMIS), designed to understand the determinants of cardiometabolic diseases in individuals from South Asia. We identified individuals carrying homozygous predicted loss-of-function (pLoF) mutations, and performed phenotypic analysis involving more than 200 biochemical and disease traits. We enumerated 49,138 rare (

Understanding the genetic basis of gene expression can help us understand the molecular underpinnings of human traits and disease. Expression quantitative trait locus (eQTL) mapping can help in studying this relationship but have been shown to be very cell-type specific, motivating the use of single-cell RNA sequencing and single-cell eQTLs to obtain a more granular view of genetic regulation. Current methods for single-cell eQTL mapping either rely on the "pseudobulk" approach and traditional pipelines for bulk transcriptomics or do not scale well to large datasets. Here, we propose SAIGE-QTL, a robust and scalable tool that can directly map eQTLs using single-cell profiles without needing aggregation at the pseudobulk level. Additionally, SAIGE-QTL allows for testing the effects of less frequent/rare genetic variation through set-based tests, which is traditionally excluded from eQTL mapping studies. We evaluate the performance of SAIGE-QTL on both real and simulated data and demonstrate the improved power for eQTL mapping over existing pipelines.

Missense variants can have a range of functional impacts depending on factors such as the specific amino acid substitution and location within the gene. To interpret their deleteriousness, studies have sought to identify regions within genes that are specifically intolerant of missense variation 1-12 . Here, we leverage the patterns of rare missense variation in 125,748 individuals in the Genome Aggregation Database (gnomAD) 13 against a null mutational model to identify transcripts that display regional differences in missense constraint. Missense-depleted regions are enriched for ClinVar 14 pathogenic variants, de novo missense variants from individuals with neurodevelopmental disorders (NDDs) 15,16 , and complex trait heritability. Following ClinGen calibration recommendations for the ACMG/AMP guidelines, we establish that regions with less than 20% of their expected missense variation achieve moderate support for pathogenicity. We create a missense deleteriousness metric (MPC) that incorporates regional constraint and outperforms other deleteriousness scores at stratifying case and control de novo missense variation, with a strong enrichment in NDDs. These results provide additional tools to aid in missense variant interpretation.

Around 60% of individuals with neurodevelopmental disorders (NDD) remain undiagnosed after comprehensive genetic testing, primarily of protein-coding genes 1 . Increasingly, large genome-sequenced cohorts are improving our ability to discover new diagnoses in the non-coding genome. Here, we identify the non-coding RNA RNU4-2 as a novel syndromic NDD gene. RNU4-2 encodes the U4 small nuclear RNA (snRNA), which is a critical component of the U4/U6.U5 tri-snRNP complex of the major spliceosome 2 . We identify an 18 bp region of RNU4-2 mapping to two structural elements in the U4/U6 snRNA duplex (the T-loop and Stem III) that is severely depleted of variation in the general population, but in which we identify heterozygous variants in 119 individuals with NDD. The vast majority of individuals (77.3%) have the same highly recurrent single base-pair insertion (n.64_65insT). We estimate that variants in this region explain 0.41% of individuals with NDD. We demonstrate that RNU4-2 is highly expressed in the developing human brain, in contrast to its contiguous counterpart RNU4-1 and other U4 homologs, supporting RNU4-2 's role as the primary U4 transcript in the brain. Overall, this work underscores the importance of non-coding genes in rare disorders. It will provide a diagnosis to thousands of individuals with NDD worldwide and pave the way for the development of effective treatments for these individuals., Competing Interests: Competing interests NW receives research funding from Novo Nordisk and has consulted for ArgoBio studio. SJS receives research funding from BioMarin Pharmaceutical. AODL is on the scientific advisory board for Congenica, was a paid consultant for Tome Biosciences and Ono Pharma USA Inc., and received reagents from PacBio to support rare disease research. HLR has received support from Illumina and Microsoft to support rare disease gene discovery and diagnosis. MHW has consulted for Illumina and Sanofi and received speaking honoraria from Illumina and GeneDx. SBM is an advisor for BioMarin, Myome and Tenaya Therapeutics. SMS has received honoraria for educational events or advisory boards from Angelini Pharma, Biocodex, Eisai, Zogenix/UCB and institutional contributions for advisory boards, educational events or consultancy work from Eisai, Jazz/GW Pharma, Stoke Therapeutics, Takeda, UCB and Zogenix. The Department of Molecular and Human Genetics at Baylor College of Medicine receives revenue from clinical genetic testing completed at Baylor Genetics Laboratories. JMMH is a full-time employee of Novo Nordisk and holds shares in Novo Nordisk A/S. DGM is a paid consultant for GlaxoSmithKline, Insitro, and Overtone Therapeutics and receives research support from Microsoft.

Collagen VI-related dystrophies (COL6-RDs) manifest with a spectrum of clinical phenotypes, ranging from Ullrich congenital muscular dystrophy (UCMD), presenting with prominent congenital symptoms and characterised by progressive muscle weakness, joint contractures and respiratory insufficiency, to Bethlem muscular dystrophy, with milder symptoms typically recognised later and at times resembling a limb girdle muscular dystrophy, and intermediate phenotypes falling between UCMD and Bethlem muscular dystrophy. Despite clinical and immunohistochemical features highly suggestive of COL6-RD, some patients had remained without an identified causative variant in COL6A1 , COL6A2 or COL6A3 . With combined muscle RNA-sequencing and whole-genome sequencing we uncovered a recurrent, de novo deep intronic variant in intron 11 of COL6A1 (c.930+189C>T) that leads to a dominantly acting in-frame pseudoexon insertion. We subsequently identified and have characterised an international cohort of forty-four patients with this COL6A1 intron 11 causative variant, one of the most common recurrent causative variants in the collagen VI genes. Patients manifest a consistently severe phenotype characterised by a paucity of early symptoms followed by an accelerated progression to a severe form of UCMD, except for one patient with somatic mosaicism for this COL6A1 intron 11 variant who manifests a milder phenotype consistent with Bethlem muscular dystrophy. Characterisation of this individual provides a robust validation for the development of our pseudoexon skipping therapy. We have previously shown that splice-modulating antisense oligomers applied in vitro effectively decreased the abundance of the mutant pseudoexon-containing COL6A1 transcripts to levels comparable to the in vivo scenario of the somatic mosaicism shown here, indicating that this therapeutic approach carries significant translational promise for ameliorating the severe form of UCMD caused by this common recurrent COL6A1 causative variant to a Bethlem muscular dystrophy phenotype.

Genes encoding long non-coding RNAs (lncRNAs) comprise a large fraction of the human genome, yet haploinsufficiency of a lncRNA has not been shown to cause a Mendelian disease. CHASERR is a highly conserved human lncRNA adjacent to CHD2- a coding gene in which de novo loss-of-function variants cause developmental and epileptic encephalopathy. Here we report three unrelated individuals each harboring an ultra-rare heterozygous de novo deletion in the CHASERR locus. We report similarities in severe developmental delay, facial dysmorphisms, and cerebral dysmyelination in these individuals, distinguishing them from the phenotypic spectrum of CHD2 haploinsufficiency. We demonstrate reduced CHASERR mRNA expression and corresponding increased CHD2 mRNA and protein in whole blood and patient-derived cell lines-specifically increased expression of the CHD2 allele in cis with the CHASERR deletion, as predicted from a prior mouse model of Chaserr haploinsufficiency. We show for the first time that de novo structural variants facilitated by Alu-mediated non-allelic homologous recombination led to deletion of a non-coding element (the lncRNA CHASERR ) to cause a rare syndromic neurodevelopmental disorder. We also demonstrate that CHD2 has bidirectional dosage sensitivity in human disease. This work highlights the need to carefully evaluate other lncRNAs, particularly those upstream of genes associated with Mendelian disorders.

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generated single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. Seventy-five percent of elements (44 of 59) validated in an in vivo transgenic reporter assay, demonstrating that single cell accessibility is a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieved significant reduction in our variant search space and nominated candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1 , and EBF3 - as well as new candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work provides novel non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.