Search

Mitochondria carry their own circular genome and disruption of the mitochondrial genome is associated with various aging-related diseases. Unlike the nuclear genome, mitochondrial DNA (mtDNA) can be present at 1000 s to 10,000 s copies in somatic cells and variants may exist in a state of heteroplasmy, where only a fraction of the DNA molecules harbors a particular variant. We quantify mtDNA heteroplasmy in 194,871 participants in the UK Biobank and find that heteroplasmy is associated with a 1.5-fold increased risk of all-cause mortality. Additionally, we functionally characterize mtDNA single nucleotide variants (SNVs) using a constraint-based score, mitochondrial local constraint score sum (MSS) and find it associated with all-cause mortality, and with the prevalence and incidence of cancer and cancer-related mortality, particularly leukemia. These results indicate that mitochondria may have a functional role in certain cancers, and mitochondrial heteroplasmic SNVs may serve as a prognostic marker for cancer, especially for leukemia., (© 2023. Springer Nature Limited.)

Inter-individual variation in the number of copies of the mitochondrial genome, called mitochondrial DNA copy number (mtDNA-CN), reflects mitochondrial function and has been associated with various aging-related diseases. We examined 415,422 exomes of self-reported White ancestry individuals from the UK Biobank and tested the impact of rare variants, at the level of single variants and through aggregate variant-set tests, on mtDNA-CN. A survey across nine variant sets tested enrichment of putatively causal variants and identified 14 genes at experiment-wide significance and three genes at marginal significance. These included associations at known mtDNA depletion syndrome genes (mtDNA helicase TWNK , p = 1.1 × 10 -30 ; mitochondrial transcription factor TFAM , p = 4.3 × 10 -15 ; mtDNA maintenance exonuclease MGME1 , p = 2.0 × 10 -6 ) and the V617F dominant gain-of-function mutation in the tyrosine kinase JAK2 (p = 2.7 × 10 -17 ), associated with myeloproliferative disease. Novel genes included the ATP-dependent protease CLPX (p = 8.4 × 10 -9 ), involved in mitochondrial proteome quality, and the mitochondrial adenylate kinase AK2 (p = 4.7 × 10 -8 ), involved in hematopoiesis. The most significant association was a missense variant in SAMHD1 (p = 4.2 × 10 -28 ), found on a rare, 1.2-Mb shared ancestral haplotype on chromosome 20. SAMHD1 encodes a cytoplasmic host restriction factor involved in viral defense response and the mitochondrial nucleotide salvage pathway, and is associated with Aicardi-Goutières syndrome 5, a childhood encephalopathy and chronic inflammatory response disorder. Rare variants were enriched in Mendelian mtDNA depletion syndrome loci, and these variants implicated core processes in mtDNA replication, nucleoid structure formation, and maintenance. These data indicate that strong-effect mutations from the nuclear genome contribute to the genetic architecture of mtDNA-CN., Competing Interests: The authors declare no competing interests., (© 2022 The Author(s).)

A challenge of next generation sequencing is read contamination. We use Genotype-Tissue Expression (GTEx) datasets and technical metadata along with RNA-seq datasets from other studies to understand factors that contribute to contamination. Here we report, of 48 analyzed tissues in GTEx, 26 have variant co-expression clusters of four highly expressed and pancreas-enriched genes (PRSS1, PNLIP, CLPS, and/or CELA3A). Fourteen additional highly expressed genes from other tissues also indicate contamination. Sample contamination is strongly associated with a sample being sequenced on the same day as a tissue that natively expresses those genes. Discrepant SNPs across four contaminating genes validate the contamination. Low-level contamination affects ~40% of samples and leads to numerous eQTL assignments in inappropriate tissues among these 18 genes. This type of contamination occurs widely, impacting bulk and single cell (scRNA-seq) data set analysis. In conclusion, highly expressed, tissue-enriched genes basally contaminate GTEx and other datasets impacting analyses.

Mycobacterium tuberculosis infects over 10 million people annually and kills more people each year than any other human pathogen. The current tuberculosis (TB) vaccine is only partially effective in preventing infection, while current TB treatment is problematic in terms of length, complexity and patient compliance. There is an urgent need for new drugs to combat the burden of TB disease and the natural environment has re-emerged as a rich source of bioactive molecules for development of lead compounds. In this study, one species of marine sponge from the Tedania genus was found to yield samples with exceptionally potent activity against M. tuberculosis. Bioassay-guided fractionation identified bengamide B as the active component, which displayed activity in the nanomolar range against both drug-sensitive and drug-resistant M. tuberculosis. The active compound inhibited in vitro activity of M. tuberculosis MetAP1c protein, suggesting the potent inhibitory action may be due to interference with methionine aminopeptidase activity. Tedania-derived bengamide B was non-toxic against human cell lines, synergised with rifampicin for in vitro inhibition of bacterial growth and reduced intracellular replication of M. tuberculosis. Thus, bengamides isolated from Tedania sp. show significant potential as a new class of compounds for the treatment of drug-resistant M. tuberculosis.

In this exciting era of "next-gen cytogenetics," integrating genomic sequencing into the prenatal diagnostic setting is possible within an actionable time frame and can provide precise delineation of balanced chromosomal rearrangements at the nucleotide level. Given the increased risk of congenital abnormalities in newborns with de novo balanced chromosomal rearrangements, comprehensive interpretation of breakpoints could substantially improve prediction of phenotypic outcomes and support perinatal medical care. Herein, we present and evaluate sequencing results of balanced chromosomal rearrangements in ten prenatal subjects with respect to the location of regulatory chromatin domains (topologically associated domains [TADs]). The genomic material from all subjects was interpreted to be "normal" by microarray analyses, and their rearrangements would not have been detected by cell-free DNA (cfDNA) screening. The findings of our systematic approach correlate with phenotypes of both pregnancies with untoward outcomes (5/10) and with healthy newborns (3/10). Two pregnancies, one with a chromosomal aberration predicted to be of unknown clinical significance and another one predicted to be likely benign, were terminated prior to phenotype-genotype correlation (2/10). We demonstrate that the clinical interpretation of structural rearrangements should not be limited to interruption, deletion, or duplication of specific genes and should also incorporate regulatory domains of the human genome with critical ramifications for the control of gene expression. As detailed in this study, our molecular approach to both detecting and interpreting the breakpoints of structural rearrangements yields unparalleled information in comparison to other commonly used first-tier diagnostic methods, such as non-invasive cfDNA screening and microarray analysis, to provide improved genetic counseling for phenotypic outcome in the prenatal setting., (Copyright © 2016 American Society of Human Genetics. All rights reserved.)

Analysis of chromosomal rearrangements has been highly successful in identifying genes involved in many congenital abnormalities including hearing loss. Herein, we report a subject, designated DGAP242, with congenital hearing loss (HL) and a de novo balanced translocation 46,XX,t(1;5)(q32;q15)dn. Using multiple next-generation sequencing techniques, we obtained high resolution of the breakpoints. This revealed disruption of the orphan receptor ESRRG on chromosome 1, which is differentially expressed in inner ear hair cells and has previously been implicated in HL, and disruption of KIAA0825 on chromosome 5. Given the translocation breakpoints and supporting literature, disruption of ESRRG is the most likely cause for DGAP242's phenotype and implicates ESRRG in a monogenic form of congenital HL, although a putative contributory role for KIAA0825 in the subject's disorder cannot be excluded.

CAPZB is an actin-capping protein that caps the growing end of F-actin and modulates the cytoskeleton and tethers actin filaments to the Z-line of the sarcomere in muscles. Whole-genome sequencing was performed on a subject with micrognathia, cleft palate and hypotonia that harbored a de novo, balanced chromosomal translocation that disrupts the CAPZB gene. The function of capzb was analyzed in the zebrafish model. capzb(-/-) mutants exhibit both craniofacial and muscle defects that recapitulate the phenotypes observed in the human subject. Loss of capzb affects cell morphology, differentiation and neural crest migration. Differentiation of both myogenic stem cells and neural crest cells requires capzb. During palate morphogenesis, defective cranial neural crest cell migration in capzb(-/-) mutants results in loss of the median cell population, creating a cleft phenotype. capzb is also required for trunk neural crest migration, as evident from melanophores disorganization in capzb(-/-) mutants. In addition, capzb over-expression results in embryonic lethality. Therefore, proper capzb dosage is important during embryogenesis, and regulates both cell behavior and tissue morphogenesis., (© The Author 2016. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Copy-number variants (CNVs) have been the predominant focus of genetic studies of structural variation, and chromosomal microarray (CMA) for genome-wide CNV detection is the recommended first-tier genetic diagnostic screen in neurodevelopmental disorders. We compared CNVs observed by CMA to the structural variation detected by whole-genome large-insert sequencing in 259 individuals diagnosed with autism spectrum disorder (ASD) from the Simons Simplex Collection. These analyses revealed a diverse landscape of complex duplications in the human genome. One remarkably common class of complex rearrangement, which we term dupINVdup, involves two closely located duplications ("paired duplications") that flank the breakpoints of an inversion. This complex variant class is cryptic to CMA, but we observed it in 8.1% of all subjects. We also detected other paired-duplication signatures and duplication-mediated complex rearrangements in 15.8% of all ASD subjects. Breakpoint analysis showed that the predominant mechanism of formation of these complex duplication-associated variants was microhomology-mediated repair. On the basis of the striking prevalence of dupINVdups in this cohort, we explored the landscape of all inversion variation among the 235 highest-quality libraries and found abundant complexity among these variants: only 39.3% of inversions were canonical, or simple, inversions without additional rearrangement. Collectively, these findings indicate that dupINVdups, as well as other complex duplication-associated rearrangements, represent relatively common sources of genomic variation that is cryptic to population-based microarray and low-depth whole-genome sequencing. They also suggest that paired-duplication signatures detected by CMA warrant further scrutiny in genetic diagnostic testing given that they might mark complex rearrangements of potential clinical relevance., (Copyright © 2015 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Genomic rearrangements are a common cause of human congenital abnormalities. However, their origin and consequences are poorly understood. We performed molecular analysis of two patients with congenital disease who carried de novo genomic rearrangements. We found that the rearrangements in both patients hit genes that are recurrently rearranged in cancer (ETV1, FOXP1, and microRNA cluster C19MC) and drive formation of fusion genes similar to those described in cancer. Subsequent analysis of a large set of 552 de novo germline genomic rearrangements underlying congenital disorders revealed enrichment for genes rearranged in cancer and overlap with somatic cancer breakpoints. Breakpoints of common (inherited) germline structural variations also overlap with cancer breakpoints but are depleted for cancer genes. We propose that the same genomic positions are prone to genomic rearrangements in germline and soma but that timing and context of breakage determines whether developmental defects or cancer are promoted., (Copyright © 2014 The Authors. Published by Elsevier Inc. All rights reserved.)

Structural variation (SV) is a significant component of the genetic etiology of both neurodevelopmental and psychiatric disorders; however, routine guidelines for clinical genetic screening have been established only in the former category. Genome-wide chromosomal microarray (CMA) can detect genomic imbalances such as copy-number variants (CNVs), but balanced chromosomal abnormalities (BCAs) still require karyotyping for clinical detection. Moreover, submicroscopic BCAs and subarray threshold CNVs are intractable, or cryptic, to both CMA and karyotyping. Here, we performed whole-genome sequencing using large-insert jumping libraries to delineate both cytogenetically visible and cryptic SVs in a single test among 30 clinically referred youth representing a range of severe neuropsychiatric conditions. We detected 96 SVs per person on average that passed filtering criteria above our highest-confidence resolution (6,305 bp) and an additional 111 SVs per genome below this resolution. These SVs rearranged 3.8 Mb of genomic sequence and resulted in 42 putative loss-of-function (LoF) or gain-of-function mutations per person. We estimate that 80% of the LoF variants were cryptic to clinical CMA. We found myriad complex and cryptic rearrangements, including a "paired" duplication (360 kb, 169 kb) that flanks a 5.25 Mb inversion that appears in 7 additional cases from clinical CNV data among 47,562 individuals. Following convergent genomic profiling of these independent clinical CNV data, we interpreted three SVs to be of potential clinical significance. These data indicate that sequence-based delineation of the full SV mutational spectrum warrants exploration in youth referred for neuropsychiatric evaluation and clinical diagnostic SV screening more broadly., (Copyright © 2014 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

NRXN1 microdeletions occur at a relatively high frequency and confer increased risk for neurodevelopmental and neurobehavioral abnormalities. The mechanism that makes NRXN1 a deletion hotspot is unknown. Here, we identified deletions of the NRXN1 region in affected cohorts, confirming a strong association with the autism spectrum and other neurodevelopmental disorders. Interestingly, deletions in both affected and control individuals were clustered in the 5' portion of NRXN1 and its immediate upstream region. To explore the mechanism of deletion, we mapped and analyzed the breakpoints of 32 deletions. At the deletion breakpoints, frequent microhomology (68.8%, 2-19 bp) suggested predominant mechanisms of DNA replication error and/or microhomology-mediated end-joining. Long terminal repeat (LTR) elements, unique non-B-DNA structures, and MEME-defined sequence motifs were significantly enriched, but Alu and LINE sequences were not. Importantly, small-size inverted repeats (minus self chains, minus sequence motifs, and partial complementary sequences) were significantly overrepresented in the vicinity of NRXN1 region deletion breakpoints, suggesting that, although they are not interrupted by the deletion process, such inverted repeats can predispose a region to genomic instability by mediating single-strand DNA looping via the annealing of partially reverse complementary strands and the promoting of DNA replication fork stalling and DNA replication error. Our observations highlight the potential importance of inverted repeats of variable sizes in generating a rearrangement hotspot in which individual breakpoints are not recurrent. Mechanisms that involve short inverted repeats in initiating deletion may also apply to other deletion hotspots in the human genome., (Copyright © 2013 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Balanced chromosomal abnormalities (BCAs) represent a relatively untapped reservoir of single-gene disruptions in neurodevelopmental disorders (NDDs). We sequenced BCAs in patients with autism or related NDDs, revealing disruption of 33 loci in four general categories: (1) genes previously associated with abnormal neurodevelopment (e.g., AUTS2, FOXP1, and CDKL5), (2) single-gene contributors to microdeletion syndromes (MBD5, SATB2, EHMT1, and SNURF-SNRPN), (3) novel risk loci (e.g., CHD8, KIRREL3, and ZNF507), and (4) genes associated with later-onset psychiatric disorders (e.g., TCF4, ZNF804A, PDE10A, GRIN2B, and ANK3). We also discovered among neurodevelopmental cases a profoundly increased burden of copy-number variants from these 33 loci and a significant enrichment of polygenic risk alleles from genome-wide association studies of autism and schizophrenia. Our findings suggest a polygenic risk model of autism and reveal that some neurodevelopmental genes are sensitive to perturbation by multiple mutational mechanisms, leading to variable phenotypic outcomes that manifest at different life stages., (Copyright © 2012 Elsevier Inc. All rights reserved.)