Search

Mitochondrial DNA (mtDNA) has an important yet often overlooked role in health and disease. Constraint models quantify the removal of deleterious variation from the population by selection and represent powerful tools for identifying genetic variation that underlies human phenotypes 1-4 . However, nuclear constraint models are not applicable to mtDNA, owing to its distinct features. Here we describe the development of a mitochondrial genome constraint model and its application to the Genome Aggregation Database (gnomAD), a large-scale population dataset that reports mtDNA variation across 56,434 human participants 5 . Specifically, we analyse constraint by comparing the observed variation in gnomAD to that expected under neutrality, which was calculated using a mtDNA mutational model and observed maximum heteroplasmy-level data. Our results highlight strong depletion of expected variation, which suggests that many deleterious mtDNA variants remain undetected. To aid their discovery, we compute constraint metrics for every mitochondrial protein, tRNA and rRNA gene, which revealed a range of intolerance to variation. We further characterize the most constrained regions within genes through regional constraint and identify the most constrained sites within the entire mitochondrial genome through local constraint, which showed enrichment of pathogenic variation. Constraint also clustered in three-dimensional structures, which provided insight into functionally important domains and their disease relevance. Notably, we identify constraint at often overlooked sites, including in rRNA and noncoding regions. Last, we demonstrate that these metrics can improve the discovery of deleterious variation that underlies rare and common phenotypes., Competing Interests: Competing interests The authors declare no competing interests., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

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 generate 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. We evaluate enhancer activity for 59 elements using an in vivo transgenic assay and validate 44 (75%), demonstrating that single cell accessibility can be a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieve significant reduction in our variant search space and nominate candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work delivers 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., (© 2024. The Author(s).)

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., (© 2024. The Author(s).)

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., (© 2022 Laricchia et al.; Published by Cold Spring Harbor Laboratory Press.)

Mutation of the Cys1 gene underlies the renal cystic disease in the Cys1 cpk/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., (© 2021. The Author(s).)

Congenital anomalies of the kidney and urinary tract (CAKUT) constitute one of the most frequent birth defects and represent the most common cause of chronic kidney disease in the first three decades of life. Despite the discovery of dozens of monogenic causes of CAKUT, most pathogenic pathways remain elusive. We performed whole-exome sequencing (WES) in 551 individuals with CAKUT and identified a heterozygous de novo stop-gain variant in ZMYM2 in two different families with CAKUT. Through collaboration, we identified in total 14 different heterozygous loss-of-function mutations in ZMYM2 in 15 unrelated families. Most mutations occurred de novo, indicating possible interference with reproductive function. Human disease features are replicated in X. tropicalis larvae with morpholino knockdowns, in which expression of truncated ZMYM2 proteins, based on individual mutations, failed to rescue renal and craniofacial defects. Moreover, heterozygous Zmym2-deficient mice recapitulated features of CAKUT with high penetrance. The ZMYM2 protein is a component of a transcriptional corepressor complex recently linked to the silencing of developmentally regulated endogenous retrovirus elements. Using protein-protein interaction assays, we show that ZMYM2 interacts with additional epigenetic silencing complexes, as well as confirming that it binds to FOXP1, a transcription factor that has also been linked to CAKUT. In summary, our findings establish that loss-of-function mutations of ZMYM2, and potentially that of other proteins in its interactome, as causes of human CAKUT, offering new routes for studying the pathogenesis of the disorder., (Copyright © 2020 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

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.

Structural variants (SVs) rearrange large segments of DNA 1 and can have profound consequences in evolution and human disease 2,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 dosage 6 . 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 findings 7 . This SV resource is freely distributed via the gnomAD browser 8 and will have broad utility in population genetics, disease-association studies, and diagnostic screening.

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 sizes 1 . 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.

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., (Copyright © 2019 by the American Society of Nephrology.)

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., (© 2018 Wiley Periodicals, Inc.)

Background: Congenital anomalies of the kidney and urinary tract (CAKUT) are the most prevalent cause of kidney disease in the first three decades of life. Previous gene panel studies showed monogenic causation in up to 12% of patients with CAKUT., Methods: We applied whole-exome sequencing to analyze the genotypes of individuals from 232 families with CAKUT, evaluating for mutations in single genes known to cause human CAKUT and genes known to cause CAKUT in mice. In consanguineous or multiplex families, we additionally performed a search for novel monogenic causes of CAKUT., Results: In 29 families (13%), we detected a causative mutation in a known gene for isolated or syndromic CAKUT that sufficiently explained the patient's CAKUT phenotype. In three families (1%), we detected a mutation in a gene reported to cause a phenocopy of CAKUT. In 15 of 155 families with isolated CAKUT, we detected deleterious mutations in syndromic CAKUT genes. Our additional search for novel monogenic causes of CAKUT in consanguineous and multiplex families revealed a potential single, novel monogenic CAKUT gene in 19 of 232 families (8%)., Conclusions: We identified monogenic mutations in a known human CAKUT gene or CAKUT phenocopy gene as the cause of disease in 14% of the CAKUT families in this study. Whole-exome sequencing provides an etiologic diagnosis in a high fraction of patients with CAKUT and will provide a new basis for the mechanistic understanding of CAKUT., (Copyright © 2018 by the American Society of Nephrology.)

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., (Copyright © 2018 by the American Society of Nephrology.)

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., (© 2018 The Authors. American Journal of Botany is published by Wiley Periodicals, Inc. on behalf of the Botanical Society of America.)

Transposons are mobile DNA elements that generate both adaptive and deleterious phenotypic variation thereby driving genome evolution. For these reasons, genomes have mechanisms to regulate transposable element (TE) activity. Approximately 12-16% of the Caenorhabditis elegans genome is composed of TEs, of which the majority are likely inactive. However, most studies of TE activity have been conducted in the laboratory strain N2, which limits our knowledge of the effects of these mobile elements across natural populations. We analyzed the distribution and abundance of TEs in 208 wild C. elegans strains to better understand how transposons contribute to variation in natural populations. We identified 3,397 TEs as compared with the reference strain, of which 2,771 are novel insertions and 241 are TEs that have been excised in at least one wild strain. Likely because of their hypothesized deleterious effects, we find that TEs are found at low allele frequencies throughout the population, and we predict functional effects of TE insertions. The abundances of TEs reflect their activities, and these data allowed us to perform both genome-wide association mappings and rare variant correlations to reveal several candidate genes that impact TE regulation, including small regulatory piwi-interacting RNAs and chromatin factors. Because TE variation in natural populations could underlie phenotypic variation for organismal and behavioral traits, the transposons that we identified and their regulatory mechanisms can be used in future studies to explore the genomics of complex traits and evolutionary changes., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

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.

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.