The landscape of regional missense mutational intolerance quantified from 125,748 exomes.

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.

Digenic inheritance involving a muscle-specific protein kinase and the giant titin protein causes a skeletal muscle myopathy.

In digenic inheritance, pathogenic variants in two genes must be inherited together to cause disease. Only very few examples of digenic inheritance have been described in the neuromuscular disease field. Here we show that predicted deleterious variants in SRPK3, encoding the X-linked serine/argenine protein kinase 3, lead to a progressive early onset skeletal muscle myopathy only when in combination with heterozygous variants in the TTN gene. The co-occurrence of predicted deleterious SRPK3/TTN variants was not seen among 76,702 healthy male individuals, and statistical modeling strongly supported digenic inheritance as the best-fitting model. Furthermore, double-mutant zebrafish (srpk3-/-; ttn.1+/-) replicated the myopathic phenotype and showed myofibrillar disorganization. Transcriptome data suggest that the interaction of srpk3 and ttn.1 in zebrafish occurs at a post-transcriptional level. We propose that digenic inheritance of deleterious changes impacting both the protein kinase SRPK3 and the giant muscle protein titin causes a skeletal myopathy and might serve as a model for other genetic diseases.

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory 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.

Mitochondrial DNA variation across 56,434 individuals in gnomAD.

Genomic databases of allele frequency are extremely helpful for evaluating clinical variants of unknown significance; however, until now, databases such as the Genome Aggregation Database (gnomAD) have focused on nuclear DNA and have ignored the mitochondrial genome (mtDNA). Here, we present a pipeline to call mtDNA variants that addresses three technical challenges: (1) detecting homoplasmic and heteroplasmic variants, present, respectively, in all or a fraction of mtDNA molecules; (2) circular mtDNA genome; and (3) misalignment of nuclear sequences of mitochondrial origin (NUMTs). We observed that mtDNA copy number per cell varied across gnomAD cohorts and influenced the fraction of NUMT-derived false-positive variant calls, which can account for the majority of putative heteroplasmies. To avoid false positives, we excluded contaminated samples, cell lines, and samples prone to NUMT misalignment due to few mtDNA copies. Furthermore, we report variants with heteroplasmy ≥10%. We applied this pipeline to 56,434 whole-genome sequences in the gnomAD v3.1 database that includes individuals of European (58%), African (25%), Latino (10%), and Asian (5%) ancestry. Our gnomAD v3.1 release contains population frequencies for 10,850 unique mtDNA variants at more than half of all mtDNA bases. Importantly, we report frequencies within each nuclear ancestral population and mitochondrial haplogroup. Homoplasmic variants account for most variant calls (98%) and unique variants (85%). We observed that 1/250 individuals carry a pathogenic mtDNA variant with heteroplasmy above 10%. These mtDNA population allele frequencies are freely accessible and will aid in diagnostic interpretation and research studies.

Cystin genetic variants cause autosomal recessive polycystic kidney disease associated with altered Myc expression.

Mutation of the Cys1 gene underlies the renal cystic disease in the Cys1cpk/cpk (cpk) mouse that phenocopies human autosomal recessive polycystic kidney disease (ARPKD). Cystin, the protein product of Cys1, is expressed in the primary apical cilia of renal ductal epithelial cells. In previous studies, we showed that cystin regulates Myc expression via interaction with the tumor suppressor, necdin. Here, we demonstrate rescue of the cpk renal phenotype by kidney-specific expression of a cystin-GFP fusion protein encoded by a transgene integrated into the Rosa26 locus. In addition, we show that expression of the cystin-GFP fusion protein in collecting duct cells down-regulates expression of Myc in cpk kidneys. Finally, we report the first human patient with an ARPKD phenotype due to homozygosity for a deleterious splicing variant in CYS1. These findings suggest that mutations in Cys1/CYS1 cause an ARPKD phenotype in mouse and human, respectively, and that the renal cystic phenotype in the mouse is driven by overexpression of the Myc proto-oncogene.

Sequential targeted exome sequencing of 1001 patients affected by unexplained limb-girdle weakness.

PURPOSE: Several hundred genetic muscle diseases have been described, all of which are rare. Their clinical and genetic heterogeneity means that a genetic diagnosis is challenging. We established an international consortium, MYO-SEQ, to aid the work-ups of muscle disease patients and to better understand disease etiology. METHODS: Exome sequencing was applied to 1001 undiagnosed patients recruited from more than 40 neuromuscular disease referral centers; standardized phenotypic information was collected for each patient. Exomes were examined for variants in 429 genes associated with muscle conditions. RESULTS: We identified suspected pathogenic variants in 52% of patients across 87 genes. We detected 401 novel variants, 116 of which were recurrent. Variants in CAPN3, DYSF, ANO5, DMD, RYR1, TTN, COL6A2, and SGCA collectively accounted for over half of the solved cases; while variants in newer disease genes, such as BVES and POGLUT1, were also found. The remaining well-characterized unsolved patients (48%) need further investigation. CONCLUSION: Using our unique infrastructure, we developed a pathway to expedite muscle disease diagnoses. Our data suggest that exome sequencing should be used for pathogenic variant detection in patients with suspected genetic muscle diseases, focusing first on the most common disease genes described here, and subsequently in rarer and newly characterized disease genes.

A structural variation reference for medical and population genetics.

Structural variants (SVs) rearrange large segments of DNA1 and can have profound consequences in evolution and human disease2,3. As national biobanks, disease-association studies, and clinical genetic testing have grown increasingly reliant on genome sequencing, population references such as the Genome Aggregation Database (gnomAD)4 have become integral in the interpretation of single-nucleotide variants (SNVs)5. However, there are no reference maps of SVs from high-coverage genome sequencing comparable to those for SNVs. Here we present a reference of sequence-resolved SVs constructed from 14,891 genomes across diverse global populations (54% non-European) in gnomAD. We discovered a rich and complex landscape of 433,371 SVs, from which we estimate that SVs are responsible for 25-29% of all rare protein-truncating events per genome. We found strong correlations between natural selection against damaging SNVs and rare SVs that disrupt or duplicate protein-coding sequence, which suggests that genes that are highly intolerant to loss-of-function are also sensitive to increased dosage6. We also uncovered modest selection against noncoding SVs in cis-regulatory elements, although selection against protein-truncating SVs was stronger than all noncoding effects. Finally, we identified very large (over one megabase), rare SVs in 3.9% of samples, and estimate that 0.13% of individuals may carry an SV that meets the existing criteria for clinically important incidental findings7. This SV resource is freely distributed via the gnomAD browser8 and will have broad utility in population genetics, disease-association studies, and diagnostic screening.

The mutational constraint spectrum quantified from variation in 141,456 humans.

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.

TBC1D8B Mutations Implicate RAB11-Dependent Vesicular Trafficking in the Pathogenesis of Nephrotic Syndrome.

BACKGROUND: Mutations in about 50 genes have been identified as monogenic causes of nephrotic syndrome, a frequent cause of CKD. These genes delineated the pathogenetic pathways and rendered significant insight into podocyte biology. METHODS: We used whole-exome sequencing to identify novel monogenic causes of steroid-resistant nephrotic syndrome (SRNS). We analyzed the functional significance of an SRNS-associated gene in vitro and in podocyte-like Drosophila nephrocytes. RESULTS: We identified hemizygous missense mutations in the gene TBC1D8B in five families with nephrotic syndrome. Coimmunoprecipitation assays indicated interactions between TBC1D8B and active forms of RAB11. Silencing TBC1D8B in HEK293T cells increased basal autophagy and exocytosis, two cellular functions that are independently regulated by RAB11. This suggests that TBC1D8B plays a regulatory role by inhibiting endogenous RAB11. Coimmunoprecipitation assays showed TBC1D8B also interacts with the slit diaphragm protein nephrin, and colocalizes with it in immortalized cell lines. Overexpressed murine Tbc1d8b with patient-derived mutations had lower affinity for endogenous RAB11 and nephrin compared with wild-type Tbc1d8b protein. Knockdown of Tbc1d8b in Drosophila impaired function of the podocyte-like nephrocytes, and caused mistrafficking of Sns, the Drosophila ortholog of nephrin. Expression of Rab11 RNAi in nephrocytes entailed defective delivery of slit diaphragm protein to the membrane, whereas RAB11 overexpression revealed a partial phenotypic overlap to Tbc1d8b loss of function. CONCLUSIONS: Novel mutations in TBC1D8B are monogenic causes of SRNS. This gene inhibits RAB11. Our findings suggest that RAB11-dependent vesicular nephrin trafficking plays a role in the pathogenesis of nephrotic syndrome.

Mutations in WDR4 as a new cause of Galloway-Mowat syndrome.

Galloway-Mowat syndrome (GAMOS) is a phenotypically heterogeneous disorder characterized by neurodevelopmental defects combined with renal-glomerular disease, manifesting with proteinuria. To identify additional monogenic disease causes, we here performed whole exome sequencing (WES), linkage analysis, and homozygosity mapping in three affected siblings of an Indian family with GAMOS. Applying established criteria for variant filtering, we identify a novel homozygous splice site mutation in the gene WDR4 as the likely disease-causing mutation in this family. In line with previous reports, we observe growth deficiency, microcephaly, developmental delay, and intellectual disability as phenotypic features resulting from WDR4 mutations. However, the newly identified allele additionally gives rise to proteinuria and nephrotic syndrome, a phenotype that was never reported in patients with WDR4 mutations. Our data thus expand the phenotypic spectrum of WDR4 mutations by demonstrating that, depending on the specific mutated allele, a renal phenotype may be present. This finding suggests that GAMOS may occupy a phenotypic spectrum with other microcephalic diseases. Furthermore, WDR4 is an additional example of a gene that encodes a tRNA modifying enzyme and gives rise to GAMOS, if mutated. Our findings thereby support the recent observation that, like neurons, podocytes of the renal glomerulus are particularly vulnerable to cellular defects resulting from altered tRNA modifications.

GAPVD1 and ANKFY1 Mutations Implicate RAB5 Regulation in Nephrotic Syndrome.

BACKGROUND: Steroid-resistant nephrotic syndrome (SRNS) is a frequent cause of CKD. The discovery of monogenic causes of SRNS has revealed specific pathogenetic pathways, but these monogenic causes do not explain all cases of SRNS. METHODS: To identify novel monogenic causes of SRNS, we screened 665 patients by whole-exome sequencing. We then evaluated the in vitro functional significance of two genes and the mutations therein that we discovered through this sequencing and conducted complementary studies in podocyte-like Drosophila nephrocytes. RESULTS: We identified conserved, homozygous missense mutations of GAPVD1 in two families with early-onset NS and a homozygous missense mutation of ANKFY1 in two siblings with SRNS. GAPVD1 and ANKFY1 interact with the endosomal regulator RAB5. Coimmunoprecipitation assays indicated interaction between GAPVD1 and ANKFY1 proteins, which also colocalized when expressed in HEK293T cells. Silencing either protein diminished the podocyte migration rate. Compared with wild-type GAPVD1 and ANKFY1, the mutated proteins produced upon ectopic expression of GAPVD1 or ANKFY1 bearing the patient-derived mutations exhibited altered binding affinity for active RAB5 and reduced ability to rescue the knockout-induced defect in podocyte migration. Coimmunoprecipitation assays further demonstrated a physical interaction between nephrin and GAPVD1, and immunofluorescence revealed partial colocalization of these proteins in rat glomeruli. The patient-derived GAPVD1 mutations reduced nephrin-GAPVD1 binding affinity. In Drosophila, silencing Gapvd1 impaired endocytosis and caused mistrafficking of the nephrin ortholog. CONCLUSIONS: Mutations in GAPVD1 and probably in ANKFY1 are novel monogenic causes of NS. The discovery of these genes implicates RAB5 regulation in the pathogenesis of human NS.

A transcriptome screen for positive selection in domesticated breadfruit and its wild relatives (Artocarpus spp.).

PREMISE OF THE STUDY: Underutilized crops, such as breadfruit (Artocarpus altilis, Moraceae) have the potential to improve global food security. Humans have artificially selected many cultivars of breadfruit since its domestication began approximately 3500 years ago. The goal of this research was to identify transcriptomic signals of positive selection and to develop genomic resources that may facilitate the development of improved breadfruit cultivars in the future. METHODS: A reference transcriptome of breadfruit was assembled de novo and annotated. Twenty-four transcriptomes of breadfruit and its wild relatives were generated and analyzed to reveal signals of positive selection that may have resulted from local adaptation or natural selection. Emphasis was placed on MADS-box genes, which are important because they often regulate fruiting timing and structures, and on carotenoid biosynthesis genes, which can impact the nutritional quality of the fruit. KEY RESULTS: Over 1000 genes showed signals of positive selection, and these genes were enriched for localization to plastids. Nucleotide sites and individuals under positive selection were discovered in MADS-box genes and carotenoid biosynthesis genes, with several sites located in cofactor or DNA-binding domains. A McDonald-Kreitman test comparing wild to cultivated samples revealed selection in one of the carotenoid biosynthesis genes, abscisic acid 8'-hydroxylase 3. CONCLUSIONS: This research highlights some of the many genes that may have been intentionally or unintentionally selected for during the human-mediated dispersal of breadfruit and stresses the importance of conserving a varied germplasm collection. It has revealed candidate genes for further study and produced new genomic resources for breadfruit.

Whole-Exome Sequencing Reveals FAT4 Mutations in a Clinically Unrecognizable Patient with Syndromic CAKUT: A Case Report.

We present the case of a patient of Macedonian origin with unilateral renal agenesis and ureterovesical junction obstruction in combination with further abnormalities including midface hypoplasia, scoliosis as well as camptodactyly of one toe. Whole-exome sequencing analysis revealed compound heterozygous variants in the FAT4 gene. Recessive variants in FAT4 are a known cause of van Maldergem syndrome (VMS) in which congenital anomalies of the kidney and urinary tract are a less characteristic but common feature. The initial presentation of our patient was not clinically recognizable. However, in view of the molecular findings, the most likely diagnosis is a mild manifestation of VMS. Only very few publications have reported patients with VMS and mutations in FAT4 to date. With this case, we hope to provide further insight into the phenotypic variability of this syndrome.

Chloroplast microsatellite markers for Artocarpus (Moraceae) developed from transcriptome sequences.

PREMISE OF THE STUDY: Chloroplast microsatellite loci were characterized from transcriptomes of Artocarpus altilis (breadfruit) and A. camansi (breadnut). They were tested in A. odoratissimus (terap) and A. altilis and evaluated in silico for two congeners. METHODS AND RESULTS: Fifteen simple sequence repeats (SSRs) were identified in chloroplast sequences from four Artocarpus transcriptome assemblies. The markers were evaluated using capillary electrophoresis in A. odoratissimus (105 accessions) and A. altilis (73). They were also evaluated in silico in A. altilis (10), A. camansi (6), and A. altilis × A. mariannensis (7) transcriptomes. All loci were polymorphic in at least one species, with all 15 polymorphic in A. camansi. Per species, average alleles per locus ranged between 2.2 and 2.5. Three loci had evidence of fragment-length homoplasy. CONCLUSIONS: These markers will complement existing nuclear markers by enabling confident identification of maternal and clone lines, which are often important in vegetatively propagated crops such as breadfruit.