1
|
Weissensteiner H, Forer L, Kronenberg F, Schönherr S. mtDNA-Server 2: advancing mitochondrial DNA analysis through highly parallelized data processing and interactive analytics. Nucleic Acids Res 2024; 52:W102-W107. [PMID: 38709886 PMCID: PMC11223783 DOI: 10.1093/nar/gkae296] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2024] [Revised: 04/03/2024] [Accepted: 04/05/2024] [Indexed: 05/08/2024] Open
Abstract
Over the past decade, mtDNA-Server established itself as one of the most widely used variant calling web-services for human mitochondrial genomes. The service accepts sequencing data in BAM format and returns an annotated variant analysis report for both homoplasmic and heteroplasmic variants. In this work we present mtDNA-Server 2, which includes several new features highly requested by the community. Most importantly, it includes (a) the integration of a novel variant calling mode that accurately call insertions, deletions and single nucleotide variants at once, (b) the integration of additional quality control and input validation modules, (c) a method to estimate the required coverage to minimize false positives and (d) an interactive analytics dashboard. Furthermore, we migrated the complete analysis workflow to the Nextflow workflow manager for improved parallelization, reproducibility and local execution. Recognizing the importance of insertions and deletions as well as offering novel quality control, validation and reporting features, mtDNA-Server 2 provides researchers and clinicians a new state-of-the-art analysis platform for interpreting mitochondrial genomes. mtDNA-Server 2 is available via mitoverse, our analysis platform that offers a centralized place for mtDNA analysis in the cloud. The web-service, source code and its documentation are freely accessible at https://mitoverse.i-med.ac.at.
Collapse
Affiliation(s)
- Hansi Weissensteiner
- Institute of Genetic Epidemiology, Medical University of Innsbruck, Innsbruck, Austria
| | - Lukas Forer
- Institute of Genetic Epidemiology, Medical University of Innsbruck, Innsbruck, Austria
| | - Florian Kronenberg
- Institute of Genetic Epidemiology, Medical University of Innsbruck, Innsbruck, Austria
| | - Sebastian Schönherr
- Institute of Genetic Epidemiology, Medical University of Innsbruck, Innsbruck, Austria
| |
Collapse
|
2
|
Kadam PS, Yang Z, Lu Y, Zhu H, Atiyas Y, Shah N, Fisher S, Nordgren E, Kim J, Issadore D, Eberwine J. Single-Mitochondrion Sequencing Uncovers Distinct Mutational Patterns and Heteroplasmy Landscape in Mouse Astrocytes and Neurons. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.06.13.598906. [PMID: 38915628 PMCID: PMC11195285 DOI: 10.1101/2024.06.13.598906] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/26/2024]
Abstract
Background Mitochondrial (mt) heteroplasmy can cause adverse biological consequences when deleterious mtDNA mutations accumulate disrupting 'normal' mt-driven processes and cellular functions. To investigate the heteroplasmy of such mtDNA changes we developed a moderate throughput mt isolation procedure to quantify the mt single-nucleotide variant (SNV) landscape in individual mouse neurons and astrocytes In this study we amplified mt-genomes from 1,645 single mitochondria (mts) isolated from mouse single astrocytes and neurons to 1. determine the distribution and proportion of mt-SNVs as well as mutation pattern in specific target regions across the mt-genome, 2. assess differences in mtDNA SNVs between neurons and astrocytes, and 3. Study cosegregation of variants in the mouse mtDNA. Results 1. The data show that specific sites of the mt-genome are permissive to SNV presentation while others appear to be under stringent purifying selection. Nested hierarchical analysis at the levels of mitochondrion, cell, and mouse reveals distinct patterns of inter- and intra-cellular variation for mt-SNVs at different sites. 2. Further, differences in the SNV incidence were observed between mouse neurons and astrocytes for two mt-SNV 9027:G>A and 9419:C>T showing variation in the mutational propensity between these cell types. Purifying selection was observed in neurons as shown by the Ka/Ks statistic, suggesting that neurons are under stronger evolutionary constraint as compared to astrocytes. 3. Intriguingly, these data show strong linkage between the SNV sites at nucleotide positions 9027 and 9461. Conclusion This study suggests that segregation as well as clonal expansion of mt-SNVs is specific to individual genomic loci, which is important foundational data in understanding of heteroplasmy and disease thresholds for mutation of pathogenic variants.
Collapse
Affiliation(s)
- Parnika S Kadam
- Department of Pharmacology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Zijian Yang
- Department of Bioengineering, School of Engineering and Applied Science, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Youtao Lu
- Department of Biology, School of Arts and Sciences; University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Hua Zhu
- Department of Pharmacology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Yasemin Atiyas
- Department of Bioengineering, School of Engineering and Applied Science, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Nishal Shah
- Department of Bioengineering, School of Engineering and Applied Science, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Stephen Fisher
- Department of Biology, School of Arts and Sciences; University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Erik Nordgren
- Department of Biology, School of Arts and Sciences; University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Junhyong Kim
- Department of Biology, School of Arts and Sciences; University of Pennsylvania, Philadelphia, PA 19104, USA
| | - David Issadore
- Department of Bioengineering, School of Engineering and Applied Science, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - James Eberwine
- Department of Pharmacology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104, USA
| |
Collapse
|
3
|
Borbély N, Dudás D, Tapasztó A, Dudás-Boda E, Csáky V, Szeifert B, Mende BG, Egyed B, Szécsényi-Nagy A, Pamjav H. Phylogenetic insights into the genetic legacies of Hungarian-speaking communities in the Carpathian Basin. Sci Rep 2024; 14:11480. [PMID: 38769390 PMCID: PMC11106325 DOI: 10.1038/s41598-024-61978-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2024] [Accepted: 05/13/2024] [Indexed: 05/22/2024] Open
Abstract
This study focuses on exploring the uniparental genetic lineages of Hungarian-speaking minorities residing in rural villages of Baranja (Croatia) and the Zobor region (Slovakia). We aimed to identify ancestral lineages by examining genetic markers distributed across the entire mitogenome and on the Y-chromosome. This allowed us to discern disparities in regional genetic structures within these communities. By integrating our newly acquired genetic data from a total of 168 participants with pre-existing Eurasian and ancient DNA datasets, our goal was to enrich the understanding of the genetic history trajectories of Carpathian Basin populations. Our findings suggest that while population-based analyses may not be sufficiently robust to detect fine-scale uniparental genetic patterns with the sample sizes at hand, phylogenetic analysis of well-characterized Y-chromosomal Short Tandem Repeat (STR) data and entire mitogenome sequences did uncover multiple lineage ties to far-flung regions and eras. While the predominant portions of both paternal and maternal DNA align with the East-Central European spectrum, rarer subhaplogroups and lineages have unveiled ancient ties to both prehistoric and historic populations spanning Europe and Eastern Eurasia. This research augments the expansive field of phylogenetics, offering critical perspectives on the genetic constitution and heritage of the communities in East-Central Europe.
Collapse
Affiliation(s)
- Noémi Borbély
- Institute of Archaeogenomics, HUN-REN Research Centre for the Humanities, Tóth Kálmán utca 4, Budapest, 1097, Hungary
- Doctoral School of Biology, Institute of Biology, ELTE Eötvös Loránd University, Pázmány Péter sétány 1/C, Budapest, 1117, Hungary
| | - Dániel Dudás
- Department of Reference Sample Analysis, Institute of Forensic Genetics, Hungarian Institute for Forensic Sciences, Gyorskocsi u. 25, Budapest, 1027, Hungary
| | - Attila Tapasztó
- Department of Reference Sample Analysis, Institute of Forensic Genetics, Hungarian Institute for Forensic Sciences, Gyorskocsi u. 25, Budapest, 1027, Hungary
| | - Eszter Dudás-Boda
- Department of Reference Sample Analysis, Institute of Forensic Genetics, Hungarian Institute for Forensic Sciences, Gyorskocsi u. 25, Budapest, 1027, Hungary
| | - Veronika Csáky
- Institute of Archaeogenomics, HUN-REN Research Centre for the Humanities, Tóth Kálmán utca 4, Budapest, 1097, Hungary
| | - Bea Szeifert
- Institute of Archaeogenomics, HUN-REN Research Centre for the Humanities, Tóth Kálmán utca 4, Budapest, 1097, Hungary
| | - Balázs Gusztáv Mende
- Institute of Archaeogenomics, HUN-REN Research Centre for the Humanities, Tóth Kálmán utca 4, Budapest, 1097, Hungary
| | - Balázs Egyed
- Department of Genetics, ELTE Eötvös Loránd University, Pázmány Péter sétány 1/C, Budapest, 1117, Hungary
| | - Anna Szécsényi-Nagy
- Institute of Archaeogenomics, HUN-REN Research Centre for the Humanities, Tóth Kálmán utca 4, Budapest, 1097, Hungary.
| | - Horolma Pamjav
- Department of Reference Sample Analysis, Institute of Forensic Genetics, Hungarian Institute for Forensic Sciences, Gyorskocsi u. 25, Budapest, 1027, Hungary.
| |
Collapse
|
4
|
Zhuang X, Ye R, Zhou Y, Cheng MY, Cui H, Wang L, Zhang S, Wang S, Cui Y, Zhang W. Leveraging new methods for comprehensive characterization of mitochondrial DNA in esophageal squamous cell carcinoma. Genome Med 2024; 16:50. [PMID: 38566210 PMCID: PMC10985887 DOI: 10.1186/s13073-024-01319-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Accepted: 03/21/2024] [Indexed: 04/04/2024] Open
Abstract
BACKGROUND Mitochondria play essential roles in tumorigenesis; however, little is known about the contribution of mitochondrial DNA (mtDNA) to esophageal squamous cell carcinoma (ESCC). Whole-genome sequencing (WGS) is by far the most efficient technology to fully characterize the molecular features of mtDNA; however, due to the high redundancy and heterogeneity of mtDNA in regular WGS data, methods for mtDNA analysis are far from satisfactory. METHODS Here, we developed a likelihood-based method dMTLV to identify low-heteroplasmic mtDNA variants. In addition, we described fNUMT, which can simultaneously detect non-reference nuclear sequences of mitochondrial origin (non-ref NUMTs) and their derived artifacts. Using these new methods, we explored the contribution of mtDNA to ESCC utilizing the multi-omics data of 663 paired tumor-normal samples. RESULTS dMTLV outperformed the existing methods in sensitivity without sacrificing specificity. The verification using Nanopore long-read sequencing data showed that fNUMT has superior specificity and more accurate breakpoint identification than the current methods. Leveraging the new method, we identified a significant association between the ESCC overall survival and the ratio of mtDNA copy number of paired tumor-normal samples, which could be potentially explained by the differential expression of genes enriched in pathways related to metabolism, DNA damage repair, and cell cycle checkpoint. Additionally, we observed that the expression of CBWD1 was downregulated by the non-ref NUMTs inserted into its intron region, which might provide precursor conditions for the tumor cells to adapt to a hypoxic environment. Moreover, we identified a strong positive relationship between the number of mtDNA truncating mutations and the contribution of signatures linked to tumorigenesis and treatment response. CONCLUSIONS Our new frameworks promote the characterization of mtDNA features, which enables the elucidation of the landscapes and roles of mtDNA in ESCC essential for extending the current understanding of ESCC etiology. dMTLV and fNUMT are freely available from https://github.com/sunnyzxh/dMTLV and https://github.com/sunnyzxh/fNUMT , respectively.
Collapse
Affiliation(s)
- Xuehan Zhuang
- Cancer Institute, Department of Oncology, Peking University Shenzhen Hospital, Shenzhen Peking University-the Hong Kong University of Science and Technology (PKU-HKUST) Medical Center; Institute of Cancer Research, Shenzhen Bay Laboratory, Shenzhen, Guangdong, 518000, China
| | - Rui Ye
- Department of Psychiatry, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong, China
| | - Yong Zhou
- Cancer Institute, Department of Oncology, Peking University Shenzhen Hospital, Shenzhen Peking University-the Hong Kong University of Science and Technology (PKU-HKUST) Medical Center; Institute of Cancer Research, Shenzhen Bay Laboratory, Shenzhen, Guangdong, 518000, China
| | - Matthew Yibo Cheng
- Cancer Institute, Department of Oncology, Peking University Shenzhen Hospital, Shenzhen Peking University-the Hong Kong University of Science and Technology (PKU-HKUST) Medical Center; Institute of Cancer Research, Shenzhen Bay Laboratory, Shenzhen, Guangdong, 518000, China
| | - Heyang Cui
- Cancer Institute, Department of Oncology, Peking University Shenzhen Hospital, Shenzhen Peking University-the Hong Kong University of Science and Technology (PKU-HKUST) Medical Center; Institute of Cancer Research, Shenzhen Bay Laboratory, Shenzhen, Guangdong, 518000, China
| | - Longlong Wang
- Cancer Institute, Department of Oncology, Peking University Shenzhen Hospital, Shenzhen Peking University-the Hong Kong University of Science and Technology (PKU-HKUST) Medical Center; Institute of Cancer Research, Shenzhen Bay Laboratory, Shenzhen, Guangdong, 518000, China
| | - Shuangping Zhang
- The Department of Thoracic Surgery, Shanxi Cancer Hospital; Key Laboratory of Cellular Physiology of the Ministry of Education, Department of Pathology, Shanxi Medical University, Taiyuan, Shanxi, 030001, China
| | - Shubin Wang
- Cancer Institute, Department of Oncology, Peking University Shenzhen Hospital, Shenzhen Peking University-the Hong Kong University of Science and Technology (PKU-HKUST) Medical Center; Institute of Cancer Research, Shenzhen Bay Laboratory, Shenzhen, Guangdong, 518000, China
| | - Yongping Cui
- Cancer Institute, Department of Oncology, Peking University Shenzhen Hospital, Shenzhen Peking University-the Hong Kong University of Science and Technology (PKU-HKUST) Medical Center; Institute of Cancer Research, Shenzhen Bay Laboratory, Shenzhen, Guangdong, 518000, China.
- The Department of Thoracic Surgery, Shanxi Cancer Hospital; Key Laboratory of Cellular Physiology of the Ministry of Education, Department of Pathology, Shanxi Medical University, Taiyuan, Shanxi, 030001, China.
| | - Weimin Zhang
- Cancer Institute, Department of Oncology, Peking University Shenzhen Hospital, Shenzhen Peking University-the Hong Kong University of Science and Technology (PKU-HKUST) Medical Center; Institute of Cancer Research, Shenzhen Bay Laboratory, Shenzhen, Guangdong, 518000, China.
- State Key Laboratory of Molecular Oncology, Beijing Key Laboratory of Carcinogenesis and Translational Research, Laboratory of Molecular Oncology, Peking University Cancer Hospital & Institute; Research Unit of Molecular Cancer Research, Chinese Academy of Medical Sciences, Beijing, 100142, China.
| |
Collapse
|
5
|
Biró B, Gál Z, Fekete Z, Klecska E, Hoffmann OI. Mitochondrial genome plasticity of mammalian species. BMC Genomics 2024; 25:278. [PMID: 38486136 PMCID: PMC10941376 DOI: 10.1186/s12864-024-10201-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2023] [Accepted: 03/08/2024] [Indexed: 03/17/2024] Open
Abstract
There is an ongoing process in which mitochondrial sequences are being integrated into the nuclear genome. The importance of these sequences has already been revealed in cancer biology, forensic, phylogenetic studies and in the evolution of the eukaryotic genetic information. Human and numerous model organisms' genomes were described from those sequences point of view. Furthermore, recent studies were published on the patterns of these nuclear localised mitochondrial sequences in different taxa.However, the results of the previously released studies are difficult to compare due to the lack of standardised methods and/or using few numbers of genomes. Therefore, in this paper our primary goal is to establish a uniform mining pipeline to explore these nuclear localised mitochondrial sequences.Our results show that the frequency of several repetitive elements is higher in the flanking regions of these sequences than expected. A machine learning model reveals that the flanking regions' repetitive elements and different structural characteristics are highly influential during the integration process.In this paper, we introduce a general mining pipeline for all mammalian genomes. The workflow is publicly available and is believed to serve as a validated baseline for future research in this field. We confirm the widespread opinion, on - as to our current knowledge - the largest dataset, that structural circumstances and events corresponding to repetitive elements are highly significant. An accurate model has also been trained to predict these sequences and their corresponding flanking regions.
Collapse
Affiliation(s)
- Bálint Biró
- Agribiotechnology and Precision Breeding for Food Security National Laboratory, Department of Animal Biotechnology, Institute of Genetics and Biotechnology, Hungarian University of Agriculture and Life Sciences, Szent-Györgyi Albert str. 4, 2100, Gödöllő, Hungary.
- Group BM, Data Insights Team, _VOIS, Kerepesi str. 35, 1087, Budapest, Hungary.
| | - Zoltán Gál
- Agribiotechnology and Precision Breeding for Food Security National Laboratory, Department of Animal Biotechnology, Institute of Genetics and Biotechnology, Hungarian University of Agriculture and Life Sciences, Szent-Györgyi Albert str. 4, 2100, Gödöllő, Hungary
| | - Zsófia Fekete
- Department of Genetics and Genomics, Institute of Genetics and Biotechnology, Hungarian University of Agriculture and Life Sciences, Szent-Györgyi Albert str. 4, 2100, Gödöllő, Hungary
| | - Eszter Klecska
- FamiCord Group, Krio Institute, Kelemen László str, 1026, Budapest, Hungary
| | - Orsolya Ivett Hoffmann
- Agribiotechnology and Precision Breeding for Food Security National Laboratory, Department of Animal Biotechnology, Institute of Genetics and Biotechnology, Hungarian University of Agriculture and Life Sciences, Szent-Györgyi Albert str. 4, 2100, Gödöllő, Hungary.
| |
Collapse
|
6
|
Harutyunyan T. The known unknowns of mitochondrial carcinogenesis: de novo NUMTs and intercellular mitochondrial transfer. Mutagenesis 2024; 39:1-12. [PMID: 37804235 DOI: 10.1093/mutage/gead031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2023] [Accepted: 10/05/2023] [Indexed: 10/09/2023] Open
Abstract
The translocation of mitochondrial DNA (mtDNA) sequences into the nuclear genome, resulted in the occurrence of nuclear sequences of mitochondrial origin (NUMTs) which can be detected in nearly all sequenced eukaryotes. However, de novo mtDNA insertions can contribute to the development of pathological conditions including cancer. Recent data indicate that de novo mtDNA translocation into chromosomes can occur due to genotoxic influence of DNA double-strand break-inducing environmental mutagens. This confirms the hypothesis of the involvement of genome instability in the occurrence of mtDNA fragments in chromosomes. Mounting evidence indicates that mitochondria can be transferred from normal cells to cancer cells and recover cellular respiration. These exchanged mitochondria can facilitate cancer progression and metastasis. This review article provides a comprehensive overview of the potential carcinogenicity of mtDNA insertions, and the relevance of mtDNA escape in cancer progression, metastasis, and treatment resistance in humans. Potential molecular targets involved in mtDNA escape and exchange of mitochondria that can be of possible clinical benefits are presented and discussed. Understanding these processes could lead to improved diagnostic approaches, novel therapeutic strategies, and a deeper understanding of the intricate relationship between mitochondria, nuclear DNA, and cancer biology.
Collapse
Affiliation(s)
- Tigran Harutyunyan
- Department of Genetics and Cytology, Yerevan State University, 1 Alex Manoogian, 0025 Yerevan, Armenia
| |
Collapse
|
7
|
Uvizl M, Puechmaille SJ, Power S, Pippel M, Carthy S, Haerty W, Myers EW, Teeling EC, Huang Z. Comparative Genome Microsynteny Illuminates the Fast Evolution of Nuclear Mitochondrial Segments (NUMTs) in Mammals. Mol Biol Evol 2024; 41:msad278. [PMID: 38124445 PMCID: PMC10764098 DOI: 10.1093/molbev/msad278] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2023] [Revised: 11/16/2023] [Accepted: 12/12/2023] [Indexed: 12/23/2023] Open
Abstract
The escape of DNA from mitochondria into the nuclear genome (nuclear mitochondrial DNA, NUMT) is an ongoing process. Although pervasively observed in eukaryotic genomes, their evolutionary trajectories in a mammal-wide context are poorly understood. The main challenge lies in the orthology assignment of NUMTs across species due to their fast evolution and chromosomal rearrangements over the past 200 million years. To address this issue, we systematically investigated the characteristics of NUMT insertions in 45 mammalian genomes and established a novel, synteny-based method to accurately predict orthologous NUMTs and ascertain their evolution across mammals. With a series of comparative analyses across taxa, we revealed that NUMTs may originate from nonrandom regions in mtDNA, are likely found in transposon-rich and intergenic regions, and unlikely code for functional proteins. Using our synteny-based approach, we leveraged 630 pairwise comparisons of genome-wide microsynteny and predicted the NUMT orthology relationships across 36 mammals. With the phylogenetic patterns of NUMT presence-and-absence across taxa, we constructed the ancestral state of NUMTs given the mammal tree using a coalescent method. We found support on the ancestral node of Fereuungulata within Laurasiatheria, whose subordinal relationships are still controversial. This study broadens our knowledge on NUMT insertion and evolution in mammalian genomes and highlights the merit of NUMTs as alternative genetic markers in phylogenetic inference.
Collapse
Affiliation(s)
- Marek Uvizl
- Department of Zoology, National Museum, 19300 Prague, Czech Republic
- Department of Zoology, Faculty of Science, Charles University, 12844 Prague, Czech Republic
| | - Sebastien J Puechmaille
- Institut des Sciences de l’Evolution de Montpellier (ISEM), University of Montpellier, 34095 Montpellier, France
- Institut Universitaire de France, Paris, France
| | - Sarahjane Power
- School of Biology and Environmental Science, University College Dublin, Belfield, Dublin 4, Ireland
| | - Martin Pippel
- Max Planck Institute of Molecular Cell Biology and Genetics, 01307 Dresden, Germany
- National Bioinformatics Infrastructure Sweden, Uppsala, Sweden
| | - Samuel Carthy
- School of Biology and Environmental Science, University College Dublin, Belfield, Dublin 4, Ireland
| | - Wilfried Haerty
- Earlham Institute, Norwich Research Park, Colney Ln, NR4 7UZ Norwich, UK
- School of Biological Sciences, University of East Anglia, Norwich, UK
| | - Eugene W Myers
- Max Planck Institute of Molecular Cell Biology and Genetics, 01307 Dresden, Germany
| | - Emma C Teeling
- School of Biology and Environmental Science, University College Dublin, Belfield, Dublin 4, Ireland
| | - Zixia Huang
- School of Biology and Environmental Science, University College Dublin, Belfield, Dublin 4, Ireland
| |
Collapse
|
8
|
Kuprina K, Smorkatcheva A, Rudyk A, Galkina S. Numerous insertions of mitochondrial DNA in the genome of the northern mole vole, Ellobius talpinus. Mol Biol Rep 2023; 51:36. [PMID: 38157080 PMCID: PMC10756869 DOI: 10.1007/s11033-023-08913-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2023] [Accepted: 10/23/2023] [Indexed: 01/03/2024]
Abstract
BACKGROUND Ellobius talpinus is a subterranean rodent representing an attractive model in population ecology studies due to its highly special lifestyle and sociality. In such studies, mitochondrial DNA (mtDNA) is widely used. However, if nuclear copies of mtDNA, aka NUMTs, are present, they may co-amplify with the target mtDNA fragment, generating misleading results. The aim of this study was to determine whether NUMTs are present in E. talpinus. METHODS AND RESULTS PCR amplification of the putative mtDNA CytB-D-loop fragment using 'universal' primers from 56 E. talpinus samples produced multiple double peaks in 90% of the sequencing chromatograms. To reveal NUMTs, molecular cloning and sequencing of PCR products of three specimens was conducted, followed by phylogenetic analysis. The pseudogene nature of three out of the seven detected haplotypes was confirmed by their basal positions in relation to other Ellobius haplotypes in the phylogenetic tree. Additionally, 'haplotype B' was basal in relation to other E. talpinus haplotypes and found present in very distant sampling sites. BLASTN search revealed 195 NUMTs in the E. talpinus nuclear genome, including fragments of all four PCR amplified pseudogenes. Although the majority of the NUMTs studied were short, the entire mtDNA had copies in the nuclear genome. The most numerous NUMTs were found for rrnL, COXI, and D-loop. CONCLUSIONS Numerous NUMTs are present in E. talpinus and can be difficult to discriminate against mtDNA sequences. Thus, in future population or phylogenetic studies in E. talpinus, the possibility of cryptic NUMTs amplification should always be taken into account.
Collapse
Affiliation(s)
- Kristina Kuprina
- Institute of Botany and Landscape Ecology, University of Greifswald, Soldmannstr. 15, Greifswald, 17489, Germany.
- Department of Vertebrate Zoology, Saint Petersburg State University, Universitetskaya nab. 7/9, Saint Petersburg, 199034, Russia.
| | - Antonina Smorkatcheva
- Department of Vertebrate Zoology, Saint Petersburg State University, Universitetskaya nab. 7/9, Saint Petersburg, 199034, Russia
| | - Anna Rudyk
- Department of Vertebrate Zoology, Saint Petersburg State University, Universitetskaya nab. 7/9, Saint Petersburg, 199034, Russia
| | - Svetlana Galkina
- Department of Genetics and Biotechnology, Saint Petersburg State University, Universitetskaya nab. 7/9, Saint Petersburg, 199034, Russia
| |
Collapse
|
9
|
Tao Y, He C, Lin D, Gu Z, Pu W. Comprehensive Identification of Mitochondrial Pseudogenes (NUMTs) in the Human Telomere-to-Telomere Reference Genome. Genes (Basel) 2023; 14:2092. [PMID: 38003036 PMCID: PMC10671835 DOI: 10.3390/genes14112092] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2023] [Revised: 11/09/2023] [Accepted: 11/12/2023] [Indexed: 11/26/2023] Open
Abstract
Practices related to mitochondrial research have long been hindered by the presence of mitochondrial pseudogenes within the nuclear genome (NUMTs). Even though partially assembled human reference genomes like hg38 have included NUMTs compilation, the exhaustive NUMTs within the only complete reference genome (T2T-CHR13) remain unknown. Here, we comprehensively identified the fixed NUMTs within the reference genome using human pan-mitogenome (HPMT) from GeneBank. The inclusion of HPMT serves the purpose of establishing an authentic mitochondrial DNA (mtDNA) mutational spectrum for the identification of NUMTs, distinguishing it from the polymorphic variations found in NUMTs. Using HPMT, we identified approximately 10% of additional NUMTs in three human reference genomes under stricter thresholds. And we also observed an approximate 6% increase in NUMTs in T2T-CHR13 compared to hg38, including NUMTs on the short arms of chromosomes 13, 14, and 15 that were not assembled previously. Furthermore, alignments based on 20-mer from mtDNA suggested the presence of more mtDNA-like short segments within the nuclear genome, which should be avoided for short amplicon or cell free mtDNA detection. Finally, through the assay of transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) on cell lines before and after mtDNA elimination, we concluded that NUMTs have a minimal impact on bulk ATAC-seq, even though 16% of sequencing data originated from mtDNA.
Collapse
Affiliation(s)
- Yichen Tao
- MOE Key Laboratory of Contemporary Anthropology, School of Life Sciences, Fudan University, Shanghai 200438, China; (Y.T.); (D.L.)
| | - Chengpeng He
- Greater Bay Area Institute of Precision Medicine (Guangzhou), Fudan University, Nansha District, Guangzhou 511458, China;
| | - Deng Lin
- MOE Key Laboratory of Contemporary Anthropology, School of Life Sciences, Fudan University, Shanghai 200438, China; (Y.T.); (D.L.)
| | - Zhenglong Gu
- MOE Key Laboratory of Contemporary Anthropology, School of Life Sciences, Fudan University, Shanghai 200438, China; (Y.T.); (D.L.)
- Greater Bay Area Institute of Precision Medicine (Guangzhou), Fudan University, Nansha District, Guangzhou 511458, China;
| | - Weilin Pu
- Greater Bay Area Institute of Precision Medicine (Guangzhou), Fudan University, Nansha District, Guangzhou 511458, China;
| |
Collapse
|
10
|
Jiang L, Liu J, Li S, Wen Y, Zheng X, Qin L, Hou Y, Wang Z. CmVCall: An automated and adjustable nanopore analysis pipeline for heteroplasmy detection of the control region in human mitochondrial genome. Forensic Sci Int Genet 2023; 67:102930. [PMID: 37595417 DOI: 10.1016/j.fsigen.2023.102930] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2023] [Revised: 08/09/2023] [Accepted: 08/10/2023] [Indexed: 08/20/2023]
Abstract
Genetic associations between human mitochondrial DNA (mtDNA) heteroplasmy and mitochondrial diseases, aging, and cancer have been elaborated, contributing a lot to the further understanding of mtDNA polymorphic spectrum in anthropology, population, and forensic genetics. In the past decade, heteroplasmy detection using Sanger sequencing and next generation sequencing (NGS) was hampered by the former's inefficiency and the latter's inherent bias due to amplification and mapping of short reads, respectively. Nanopore sequencing stands out for its ability to yield long contiguous segments of DNA, providing a new insight into heterogeneity authentication. In addition to MinION from Oxford Nanopore Technologies, an alternative nanopore sequencer QNome (Qitan Technology) has also been applied to various biological research and the forensic applicability of this platform has been proved recently. In this study, we evaluated the performance of four commonly used variant callers in the heterogeneity authentication of the control region of human mtDNA based on simulations of different ratios generated by mixing QNome nanopore sequencing reads of two synthetic sequences. Then, an open-source and python-based nanopore analytics pipeline, CmVCall was developed and incorporated multiple programs including reads filtering, removal of nuclear mitochondrial sequences (NUMTs), alignment, optional 'Correction' mode, and heterogeneity identification. CmVCall can achieve high precision, accuracy, and recall of 100%, 99.9%, and 92.3% with a 5% heteroplasmy level in 'Correction' mode. Moreover, blood, saliva, and hair shaft samples from monozygotic (MZ) twins were used for heterogeneity evaluation and comparison with the NGS data. Results of MZ twin samples showed that CmVCall could identify more point heteroplasmy sites, revealing significant levels of inter- and intra-individual mtDNA polymorphism. In conclusion, we believe that this analysis pipeline will lay a solid foundation for the development of a comprehensive nanopore analysis pipeline targeting the whole mitochondrial genome.
Collapse
Affiliation(s)
- Lirong Jiang
- Institute of Forensic Medicine, West China School of Basic Medical Sciences & Forensic Medicine, Sichuan University, Chengdu 610041, China
| | - Jing Liu
- Institute of Forensic Medicine, West China School of Basic Medical Sciences & Forensic Medicine, Sichuan University, Chengdu 610041, China
| | - Suyu Li
- Institute of Forensic Medicine, West China School of Basic Medical Sciences & Forensic Medicine, Sichuan University, Chengdu 610041, China
| | - Yufeng Wen
- School of Life Sciences, Jilin University, Changchun 130012, China
| | - Xinyue Zheng
- Institute of Forensic Medicine, West China School of Basic Medical Sciences & Forensic Medicine, Sichuan University, Chengdu 610041, China
| | - Liu Qin
- Qitan Technology Ltd., Chengdu, Chengdu 610044, China.
| | - Yiping Hou
- Institute of Forensic Medicine, West China School of Basic Medical Sciences & Forensic Medicine, Sichuan University, Chengdu 610041, China.
| | - Zheng Wang
- Institute of Forensic Medicine, West China School of Basic Medical Sciences & Forensic Medicine, Sichuan University, Chengdu 610041, China.
| |
Collapse
|
11
|
Gupta R, Kanai M, Durham TJ, Tsuo K, McCoy JG, Kotrys AV, Zhou W, Chinnery PF, Karczewski KJ, Calvo SE, Neale BM, Mootha VK. Nuclear genetic control of mtDNA copy number and heteroplasmy in humans. Nature 2023; 620:839-848. [PMID: 37587338 PMCID: PMC10447254 DOI: 10.1038/s41586-023-06426-5] [Citation(s) in RCA: 22] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2022] [Accepted: 07/11/2023] [Indexed: 08/18/2023]
Abstract
Mitochondrial DNA (mtDNA) is a maternally inherited, high-copy-number genome required for oxidative phosphorylation1. Heteroplasmy refers to the presence of a mixture of mtDNA alleles in an individual and has been associated with disease and ageing. Mechanisms underlying common variation in human heteroplasmy, and the influence of the nuclear genome on this variation, remain insufficiently explored. Here we quantify mtDNA copy number (mtCN) and heteroplasmy using blood-derived whole-genome sequences from 274,832 individuals and perform genome-wide association studies to identify associated nuclear loci. Following blood cell composition correction, we find that mtCN declines linearly with age and is associated with variants at 92 nuclear loci. We observe that nearly everyone harbours heteroplasmic mtDNA variants obeying two principles: (1) heteroplasmic single nucleotide variants tend to arise somatically and accumulate sharply after the age of 70 years, whereas (2) heteroplasmic indels are maternally inherited as mixtures with relative levels associated with 42 nuclear loci involved in mtDNA replication, maintenance and novel pathways. These loci may act by conferring a replicative advantage to certain mtDNA alleles. As an illustrative example, we identify a length variant carried by more than 50% of humans at position chrM:302 within a G-quadruplex previously proposed to mediate mtDNA transcription/replication switching2,3. We find that this variant exerts cis-acting genetic control over mtDNA abundance and is itself associated in-trans with nuclear loci encoding machinery for this regulatory switch. Our study suggests that common variation in the nuclear genome can shape variation in mtCN and heteroplasmy dynamics across the human population.
Collapse
Affiliation(s)
- Rahul Gupta
- Howard Hughes Medical Institute and Department of Molecular Biology, Massachusetts General Hospital, Boston, MA, USA.
- Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Analytic and Translational Genetics Unit, Center for Genomic Medicine, Massachusetts General Hospital, Boston, MA, USA.
| | - Masahiro Kanai
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Analytic and Translational Genetics Unit, Center for Genomic Medicine, Massachusetts General Hospital, Boston, MA, USA
| | - Timothy J Durham
- Howard Hughes Medical Institute and Department of Molecular Biology, Massachusetts General Hospital, Boston, MA, USA
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Kristin Tsuo
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Analytic and Translational Genetics Unit, Center for Genomic Medicine, Massachusetts General Hospital, Boston, MA, USA
| | - Jason G McCoy
- Howard Hughes Medical Institute and Department of Molecular Biology, Massachusetts General Hospital, Boston, MA, USA
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Anna V Kotrys
- Howard Hughes Medical Institute and Department of Molecular Biology, Massachusetts General Hospital, Boston, MA, USA
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Wei Zhou
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Analytic and Translational Genetics Unit, Center for Genomic Medicine, Massachusetts General Hospital, Boston, MA, USA
| | - Patrick F Chinnery
- Department of Clinical Neurosciences, University of Cambridge, Cambridge, UK
- MRC Mitochondrial Biology Unit, University of Cambridge, Cambridge, UK
| | - Konrad J Karczewski
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Analytic and Translational Genetics Unit, Center for Genomic Medicine, Massachusetts General Hospital, Boston, MA, USA
| | - Sarah E Calvo
- Howard Hughes Medical Institute and Department of Molecular Biology, Massachusetts General Hospital, Boston, MA, USA
- Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Benjamin M Neale
- Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Analytic and Translational Genetics Unit, Center for Genomic Medicine, Massachusetts General Hospital, Boston, MA, USA.
| | - Vamsi K Mootha
- Howard Hughes Medical Institute and Department of Molecular Biology, Massachusetts General Hospital, Boston, MA, USA.
- Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Department of Systems Biology, Harvard Medical School, Boston, MA, USA.
| |
Collapse
|
12
|
Giosa D, Lombardo D, Musolino C, Chines V, Raffa G, Casuscelli di Tocco F, D'Aliberti D, Caminiti G, Saitta C, Alibrandi A, Aiese Cigliano R, Romeo O, Navarra G, Raimondo G, Pollicino T. Mitochondrial DNA is a target of HBV integration. Commun Biol 2023; 6:684. [PMID: 37400627 DOI: 10.1038/s42003-023-05017-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Accepted: 06/05/2023] [Indexed: 07/05/2023] Open
Abstract
Hepatitis B virus (HBV) may integrate into the genome of infected cells and contribute to hepatocarcinogenesis. However, the role of HBV integration in hepatocellular carcinoma (HCC) development remains unclear. In this study, we apply a high-throughput HBV integration sequencing approach that allows sensitive identification of HBV integration sites and enumeration of integration clones. We identify 3339 HBV integration sites in paired tumour and non-tumour tissue samples from 7 patients with HCC. We detect 2107 clonally expanded integrations (1817 in tumour and 290 in non-tumour tissues), and a significant enrichment of clonal HBV integrations in mitochondrial DNA (mtDNA) preferentially occurring in the oxidative phosphorylation genes (OXPHOS) and D-loop region. We also find that HBV RNA sequences are imported into the mitochondria of hepatoma cells with the involvement of polynucleotide phosphorylase (PNPASE), and that HBV RNA might have a role in the process of HBV integration into mtDNA. Our results suggest a potential mechanism by which HBV integration may contribute to HCC development.
Collapse
Affiliation(s)
- Domenico Giosa
- Department of Clinical and Experimental Medicine, University Hospital of Messina, Messina, Italy
- Laboratory of Molecular Hepatology, University Hospital of Messina, Messina, Italy
| | - Daniele Lombardo
- Department of Clinical and Experimental Medicine, University Hospital of Messina, Messina, Italy
- Laboratory of Molecular Hepatology, University Hospital of Messina, Messina, Italy
| | - Cristina Musolino
- Laboratory of Molecular Hepatology, University Hospital of Messina, Messina, Italy
- Department of Human Pathology, University Hospital of Messina, Messina, Italy
| | - Valeria Chines
- Department of Clinical and Experimental Medicine, University Hospital of Messina, Messina, Italy
- Laboratory of Molecular Hepatology, University Hospital of Messina, Messina, Italy
| | - Giuseppina Raffa
- Department of Clinical and Experimental Medicine, University Hospital of Messina, Messina, Italy
- Laboratory of Molecular Hepatology, University Hospital of Messina, Messina, Italy
| | - Francesca Casuscelli di Tocco
- Department of Clinical and Experimental Medicine, University Hospital of Messina, Messina, Italy
- Laboratory of Molecular Hepatology, University Hospital of Messina, Messina, Italy
| | - Deborah D'Aliberti
- Department of Clinical and Experimental Medicine, University Hospital of Messina, Messina, Italy
- Laboratory of Molecular Hepatology, University Hospital of Messina, Messina, Italy
| | - Giuseppe Caminiti
- Laboratory of Molecular Hepatology, University Hospital of Messina, Messina, Italy
| | - Carlo Saitta
- Department of Clinical and Experimental Medicine, University Hospital of Messina, Messina, Italy
| | | | | | - Orazio Romeo
- Department of ChiBioFarAm, University of Messina, Messina, Italy
| | - Giuseppe Navarra
- Department of Human Pathology, University Hospital of Messina, Messina, Italy
| | - Giovanni Raimondo
- Department of Clinical and Experimental Medicine, University Hospital of Messina, Messina, Italy
| | - Teresa Pollicino
- Department of Clinical and Experimental Medicine, University Hospital of Messina, Messina, Italy.
- Laboratory of Molecular Hepatology, University Hospital of Messina, Messina, Italy.
| |
Collapse
|
13
|
Frascarelli C, Zanetti N, Nasca A, Izzo R, Lamperti C, Lamantea E, Legati A, Ghezzi D. Nanopore long-read next-generation sequencing for detection of mitochondrial DNA large-scale deletions. Front Genet 2023; 14:1089956. [PMID: 37456669 PMCID: PMC10344361 DOI: 10.3389/fgene.2023.1089956] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2022] [Accepted: 06/13/2023] [Indexed: 07/18/2023] Open
Abstract
Primary mitochondrial diseases are progressive genetic disorders affecting multiple organs and characterized by mitochondrial dysfunction. These disorders can be caused by mutations in nuclear genes coding proteins with mitochondrial localization or by genetic defects in the mitochondrial genome (mtDNA). The latter include point pathogenic variants and large-scale deletions/rearrangements. MtDNA molecules with the wild type or a variant sequence can exist together in a single cell, a condition known as mtDNA heteroplasmy. MtDNA single point mutations are typically detected by means of Next-Generation Sequencing (NGS) based on short reads which, however, are limited for the identification of structural mtDNA alterations. Recently, new NGS technologies based on long reads have been released, allowing to obtain sequences of several kilobases in length; this approach is suitable for detection of structural alterations affecting the mitochondrial genome. In the present work we illustrate the optimization of two sequencing protocols based on long-read Oxford Nanopore Technology to detect mtDNA structural alterations. This approach presents strong advantages in the analysis of mtDNA compared to both short-read NGS and traditional techniques, potentially becoming the method of choice for genetic studies on mtDNA.
Collapse
Affiliation(s)
- Chiara Frascarelli
- Unit of Medical Genetics and Neurogenetics, Fondazione IRCCS Istituto Neurologico Carlo Besta, Milan, Italy
| | - Nadia Zanetti
- Unit of Medical Genetics and Neurogenetics, Fondazione IRCCS Istituto Neurologico Carlo Besta, Milan, Italy
| | - Alessia Nasca
- Unit of Medical Genetics and Neurogenetics, Fondazione IRCCS Istituto Neurologico Carlo Besta, Milan, Italy
| | - Rossella Izzo
- Unit of Medical Genetics and Neurogenetics, Fondazione IRCCS Istituto Neurologico Carlo Besta, Milan, Italy
| | - Costanza Lamperti
- Unit of Medical Genetics and Neurogenetics, Fondazione IRCCS Istituto Neurologico Carlo Besta, Milan, Italy
| | - Eleonora Lamantea
- Unit of Medical Genetics and Neurogenetics, Fondazione IRCCS Istituto Neurologico Carlo Besta, Milan, Italy
| | - Andrea Legati
- Unit of Medical Genetics and Neurogenetics, Fondazione IRCCS Istituto Neurologico Carlo Besta, Milan, Italy
| | - Daniele Ghezzi
- Unit of Medical Genetics and Neurogenetics, Fondazione IRCCS Istituto Neurologico Carlo Besta, Milan, Italy
- Department of Pathophysiology and Transplantation (DEPT), University of Milan, Milan, Italy
| |
Collapse
|
14
|
Hebert PDN, Bock DG, Prosser SWJ. Interrogating 1000 insect genomes for NUMTs: A risk assessment for estimates of species richness. PLoS One 2023; 18:e0286620. [PMID: 37289794 PMCID: PMC10249859 DOI: 10.1371/journal.pone.0286620] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2023] [Accepted: 05/22/2023] [Indexed: 06/10/2023] Open
Abstract
The nuclear genomes of most animal species include NUMTs, segments of the mitogenome incorporated into their chromosomes. Although NUMT counts are known to vary greatly among species, there has been no comprehensive study of their frequency/attributes in the most diverse group of terrestrial organisms, insects. This study examines NUMTs derived from a 658 bp 5' segment of the cytochrome c oxidase I (COI) gene, the barcode region for the animal kingdom. This assessment is important because unrecognized NUMTs can elevate estimates of species richness obtained through DNA barcoding and derived approaches (eDNA, metabarcoding). This investigation detected nearly 10,000 COI NUMTs ≥ 100 bp in the genomes of 1,002 insect species (range = 0-443). Variation in nuclear genome size explained 56% of the mitogenome-wide variation in NUMT counts. Although insect orders with the largest genome sizes possessed the highest NUMT counts, there was considerable variation among their component lineages. Two thirds of COI NUMTs possessed an IPSC (indel and/or premature stop codon) allowing their recognition and exclusion from downstream analyses. The remainder can elevate species richness as they showed 10.1% mean divergence from their mitochondrial homologue. The extent of exposure to "ghost species" is strongly impacted by the target amplicon's length. NUMTs can raise apparent species richness by up to 22% when a 658 bp COI amplicon is examined versus a doubling of apparent richness when 150 bp amplicons are targeted. Given these impacts, metabarcoding and eDNA studies should target the longest possible amplicons while also avoiding use of 12S/16S rDNA as they triple NUMT exposure because IPSC screens cannot be employed.
Collapse
Affiliation(s)
- Paul D. N. Hebert
- Centre for Biodiversity Genomics, University of Guelph, Guelph, ON, Canada
| | - Dan G. Bock
- Centre for Biodiversity Genomics, University of Guelph, Guelph, ON, Canada
| | - Sean W. J. Prosser
- Centre for Biodiversity Genomics, University of Guelph, Guelph, ON, Canada
| |
Collapse
|
15
|
Vandiver AR, Hoang AN, Herbst A, Lee CC, Aiken JM, McKenzie D, Teitell MA, Timp W, Wanagat J. Nanopore sequencing identifies a higher frequency and expanded spectrum of mitochondrial DNA deletion mutations in human aging. Aging Cell 2023; 22:e13842. [PMID: 37132288 PMCID: PMC10265159 DOI: 10.1111/acel.13842] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2023] [Revised: 03/23/2023] [Accepted: 03/27/2023] [Indexed: 05/04/2023] Open
Abstract
Mitochondrial DNA (mtDNA) deletion mutations cause many human diseases and are linked to age-induced mitochondrial dysfunction. Mapping the mutation spectrum and quantifying mtDNA deletion mutation frequency is challenging with next-generation sequencing methods. We hypothesized that long-read sequencing of human mtDNA across the lifespan would detect a broader spectrum of mtDNA rearrangements and provide a more accurate measurement of their frequency. We employed nanopore Cas9-targeted sequencing (nCATS) to map and quantitate mtDNA deletion mutations and develop analyses that are fit-for-purpose. We analyzed total DNA from vastus lateralis muscle in 15 males ranging from 20 to 81 years of age and substantia nigra from three 20-year-old and three 79-year-old men. We found that mtDNA deletion mutations detected by nCATS increased exponentially with age and mapped to a wider region of the mitochondrial genome than previously reported. Using simulated data, we observed that large deletions are often reported as chimeric alignments. To address this, we developed two algorithms for deletion identification which yield consistent deletion mapping and identify both previously reported and novel mtDNA deletion breakpoints. The identified mtDNA deletion frequency measured by nCATS correlates strongly with chronological age and predicts the deletion frequency as measured by digital PCR approaches. In substantia nigra, we observed a similar frequency of age-related mtDNA deletions to those observed in muscle samples, but noted a distinct spectrum of deletion breakpoints. NCATS-mtDNA sequencing allows the identification of mtDNA deletions on a single-molecule level, characterizing the strong relationship between mtDNA deletion frequency and chronological aging.
Collapse
Affiliation(s)
- Amy R. Vandiver
- Division of Dermatology, Department of MedicineUCLALos AngelesCaliforniaUSA
- Veterans Administration Greater Los Angeles Healthcare SystemLos AngelesCaliforniaUSA
| | - Austin N. Hoang
- Division of Geriatrics, Department of MedicineUCLALos AngelesCaliforniaUSA
| | - Allen Herbst
- US Geological Survey National Wildlife Health CenterMadisonWisconsinUSA
| | - Cathy C. Lee
- Veterans Administration Greater Los Angeles Healthcare SystemLos AngelesCaliforniaUSA
- Division of Geriatrics, Department of MedicineUCLALos AngelesCaliforniaUSA
| | - Judd M. Aiken
- Department of Agricultural, Food and Nutritional SciencesUniversity of AlbertaEdmontonAlbertaCanada
| | - Debbie McKenzie
- Department of Biological SciencesUniversity of AlbertaEdmontonAlbertaCanada
| | - Michael A. Teitell
- Molecular Biology InstituteUniversity of California at Los AngelesLos AngelesCaliforniaUSA
- Department of Pathology and Laboratory Medicine, David Geffen School of MedicineUniversity of California at Los AngelesLos AngelesCaliforniaUSA
- Department of Bioengineering, California NanoSystems Institute, Broad Center for Regenerative Medicine and Stem Cell ResearchUniversity of California at Los AngelesLos AngelesCaliforniaUSA
- Department of Pediatrics, David Geffen School of MedicineUniversity of California at Los AngelesLos AngelesCaliforniaUSA
- Jonsson Comprehensive Cancer Center, David Geffen School of MedicineUniversity of California at Los AngelesLos AngelesCaliforniaUSA
| | - Winston Timp
- Department of Molecular Biology and GeneticsJohns Hopkins University School of MedicineBaltimoreMarylandUSA
- Department of Biomedical EngineeringJohns Hopkins UniversityBaltimoreMarylandUSA
- Department of Genetic MedicineJohns Hopkins University School of MedicineBaltimoreMarylandUSA
| | - Jonathan Wanagat
- Veterans Administration Greater Los Angeles Healthcare SystemLos AngelesCaliforniaUSA
- Division of Geriatrics, Department of MedicineUCLALos AngelesCaliforniaUSA
| |
Collapse
|
16
|
Xue L, Moreira JD, Smith KK, Fetterman JL. The Mighty NUMT: Mitochondrial DNA Flexing Its Code in the Nuclear Genome. Biomolecules 2023; 13:753. [PMID: 37238623 PMCID: PMC10216076 DOI: 10.3390/biom13050753] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Revised: 04/07/2023] [Accepted: 04/24/2023] [Indexed: 05/28/2023] Open
Abstract
Nuclear-mitochondrial DNA segments (NUMTs) are mitochondrial DNA (mtDNA) fragments that have been inserted into the nuclear genome. Some NUMTs are common within the human population but most NUMTs are rare and specific to individuals. NUMTs range in size from 24 base pairs to encompassing nearly the entire mtDNA and are found throughout the nuclear genome. Emerging evidence suggests that the formation of NUMTs is an ongoing process in humans. NUMTs contaminate sequencing results of the mtDNA by introducing false positive variants, particularly heteroplasmic variants present at a low variant allele frequency (VAF). In our review, we discuss the prevalence of NUMTs in the human population, the potential mechanisms of de novo NUMT insertion via DNA repair mechanisms, and provide an overview of the existing approaches for minimizing NUMT contamination. Apart from filtering known NUMTs, both wet lab-based and computational methods can be used to minimize the contamination of NUMTs in analyses of human mtDNA. Current approaches include: (1) isolating mitochondria to enrich for mtDNA; (2) applying basic local alignment to identify NUMTs for subsequent filtering; (3) bioinformatic pipelines for NUMT detection; (4) k-mer-based NUMT detection; and (5) filtering candidate false positive variants by mtDNA copy number, VAF, or sequence quality score. Multiple approaches must be applied in order to effectively identify NUMTs in samples. Although next-generation sequencing is revolutionizing our understanding of heteroplasmic mtDNA, it also raises new challenges with the high prevalence and individual-specific NUMTs that need to be handled with care in studies of mitochondrial genetics.
Collapse
Affiliation(s)
- Liying Xue
- Evans Department of Medicine, Boston University Chobanian & Avedisian School of Medicine, Boston, MA 02118, USA
| | - Jesse D. Moreira
- Department of Health Sciences, Programs in Human Physiology, Boston University Sargent College, Boston, MA 02215, USA
| | - Karan K. Smith
- Evans Department of Medicine, Boston University Chobanian & Avedisian School of Medicine, Boston, MA 02118, USA
| | - Jessica L. Fetterman
- Evans Department of Medicine, Boston University Chobanian & Avedisian School of Medicine, Boston, MA 02118, USA
| |
Collapse
|
17
|
Zhou W, Karan KR, Gu W, Klein HU, Sturm G, De Jager PL, Bennett DA, Hirano M, Picard M, Mills RE. Somatic nuclear mitochondrial DNA insertions are prevalent in the human brain and accumulate over time in fibroblasts. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.02.03.527065. [PMID: 36778249 PMCID: PMC9915708 DOI: 10.1101/2023.02.03.527065] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
The transfer of mitochondrial DNA into the nuclear genomes of eukaryotes (Numts) has been linked to lifespan in non-human species 1-3 and recently demonstrated to occur in rare instances from one human generation to the next 4. Here we investigated numtogenesis dynamics in humans in two ways. First, we quantified Numts in 1,187 post-mortem brain and blood samples from different individuals. Compared to circulating immune cells (n=389), post-mitotic brain tissue (n=798) contained more Numts, consistent with their potential somatic accumulation. Within brain samples we observed a 5.5-fold enrichment of somatic Numt insertions in the dorsolateral prefrontal cortex compared to cerebellum samples, suggesting that brain Numts arose spontaneously during development or across the lifespan. Moreover, more brain Numts was linked to earlier mortality. The brains of individuals with no cognitive impairment who died at younger ages carried approximately 2 more Numts per decade of life lost than those who lived longer. Second, we tested the dynamic transfer of Numts using a repeated-measures WGS design in a human fibroblast model that recapitulates several molecular hallmarks of aging 5. These longitudinal experiments revealed a gradual accumulation of one Numt every ~13 days. Numtogenesis was independent of large-scale genomic instability and unlikely driven cell clonality. Targeted pharmacological perturbations including chronic glucocorticoid signaling or impairing mitochondrial oxidative phosphorylation (OxPhos) only modestly increased the rate of numtogenesis, whereas patient-derived SURF1-mutant cells exhibiting mtDNA instability accumulated Numts 4.7-fold faster than healthy donors. Combined, our data document spontaneous numtogenesis in human cells and demonstrate an association between brain cortical somatic Numts and human lifespan. These findings open the possibility that mito-nuclear horizontal gene transfer among human post-mitotic tissues produce functionally-relevant human Numts over timescales shorter than previously assumed.
Collapse
Affiliation(s)
- Weichen Zhou
- Department of Computational Medicine and Bioinformatics, University of Michigan Medical School, Ann Arbor, MI 48109, USA
| | - Kalpita R. Karan
- Department of Psychiatry, Division of Behavioral Medicine, Columbia University Irving Medical Center, New York, USA
| | - Wenjin Gu
- Department of Computational Medicine and Bioinformatics, University of Michigan Medical School, Ann Arbor, MI 48109, USA
| | - Hans-Ulrich Klein
- Center for Translational & Computational Neuroimmunology, Department of Neurology, Columbia University Irving Medical Center, New York, NY 10032 USA
- Taub Institute for Research on Alzheimer’s Disease and the Aging Brain, Columbia University Irving Medical Center, New York, NY 10032 USA
| | - Gabriel Sturm
- Department of Psychiatry, Division of Behavioral Medicine, Columbia University Irving Medical Center, New York, USA
- Department of Biochemistry and Biophysics, University of California, San Francisco, San Francisco, CA 94158, USA
| | - Philip L. De Jager
- Center for Translational & Computational Neuroimmunology, Department of Neurology, Columbia University Irving Medical Center, New York, NY 10032 USA
- Taub Institute for Research on Alzheimer’s Disease and the Aging Brain, Columbia University Irving Medical Center, New York, NY 10032 USA
| | - David A. Bennett
- Rush Alzheimer’s Disease Center, Rush University Medical Center, Chicago, IL 60612 USA
| | - Michio Hirano
- Center for Translational & Computational Neuroimmunology, Department of Neurology, Columbia University Irving Medical Center, New York, NY 10032 USA
| | - Martin Picard
- Department of Psychiatry, Division of Behavioral Medicine, Columbia University Irving Medical Center, New York, USA
- Department of Neurology, H. Houston Merritt Center, Columbia University Translational Neuroscience Initiative, Columbia University Irving Medical Center, New York, USA
- New York State Psychiatric Institute, New York, USA
| | - Ryan E Mills
- Department of Computational Medicine and Bioinformatics, University of Michigan Medical School, Ann Arbor, MI 48109, USA
- Department of Human Genetics, University of Michigan Medical School, Ann Arbor, MI 48109, USA
| |
Collapse
|
18
|
Legati A, Ghezzi D, Viscomi C. Mitochondrial DNA Sequencing and Heteroplasmy Quantification by Next Generation Sequencing. Methods Mol Biol 2023; 2615:381-395. [PMID: 36807805 DOI: 10.1007/978-1-0716-2922-2_26] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/23/2023]
Abstract
Over the last 10 years, next generation sequencing (NGS) became the gold standard for both diagnosis and discovery of new disease genes responsible for heterogeneous disorders, such as mitochondrial encephalomyopathies. The application of this technology to mtDNA mutations poses extra challenges compared to other genetic conditions because of the peculiarities of mitochondrial genetics and the requirement for proper NGS data management and analysis. Here, we describe a detailed, clinically relevant protocol to sequence the whole mtDNA and quantify heteroplasmy levels of mtDNA variants, starting from total DNA through the generation of a single PCR amplicon.
Collapse
Affiliation(s)
- Andrea Legati
- Unit of Medical Genetics and Neurogenetics, Fondazione IRCCS Istituto Neurologico Carlo Besta, Milan, Italy
| | - Daniele Ghezzi
- Unit of Medical Genetics and Neurogenetics, Fondazione IRCCS Istituto Neurologico Carlo Besta, Milan, Italy.,Department of Pathophysiology and Transplantation, University of Milan, Milan, Italy.,Lab of Neurogenetics and Mitochondrial Disorders, Fondazione IRCCS Istituto Neurologico Carlo Besta/Università degli Studi di Milano, Milan, Italy
| | - Carlo Viscomi
- Department of Biomedical Sciences, University of Padova, Padova, Italy.
| |
Collapse
|
19
|
Gupta R, Kanai M, Durham TJ, Tsuo K, McCoy JG, Chinnery PF, Karczewski KJ, Calvo SE, Neale BM, Mootha VK. Nuclear genetic control of mtDNA copy number and heteroplasmy in humans. MEDRXIV : THE PREPRINT SERVER FOR HEALTH SCIENCES 2023:2023.01.19.23284696. [PMID: 36711677 PMCID: PMC9882621 DOI: 10.1101/2023.01.19.23284696] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
Human mitochondria contain a high copy number, maternally transmitted genome (mtDNA) that encodes 13 proteins required for oxidative phosphorylation. Heteroplasmy arises when multiple mtDNA variants co-exist in an individual and can exhibit complex dynamics in disease and in aging. As all proteins involved in mtDNA replication and maintenance are nuclear-encoded, heteroplasmy levels can, in principle, be under nuclear genetic control, however this has never been shown in humans. Here, we develop algorithms to quantify mtDNA copy number (mtCN) and heteroplasmy levels using blood-derived whole genome sequences from 274,832 individuals of diverse ancestry and perform GWAS to identify nuclear loci controlling these traits. After careful correction for blood cell composition, we observe that mtCN declines linearly with age and is associated with 92 independent nuclear genetic loci. We find that nearly every individual carries heteroplasmic variants that obey two key patterns: (1) heteroplasmic single nucleotide variants are somatic mutations that accumulate sharply after age 70, while (2) heteroplasmic indels are maternally transmitted as mtDNA mixtures with resulting levels influenced by 42 independent nuclear loci involved in mtDNA replication, maintenance, and novel pathways. These nuclear loci do not appear to act by mtDNA mutagenesis, but rather, likely act by conferring a replicative advantage to specific mtDNA molecules. As an illustrative example, the most common heteroplasmy we identify is a length variant carried by >50% of humans at position m.302 within a G-quadruplex known to serve as a replication switch. We find that this heteroplasmic variant exerts cis -acting genetic control over mtDNA abundance and is itself under trans -acting genetic control of nuclear loci encoding protein components of this regulatory switch. Our study showcases how nuclear haplotype can privilege the replication of specific mtDNA molecules to shape mtCN and heteroplasmy dynamics in the human population.
Collapse
Affiliation(s)
- Rahul Gupta
- Howard Hughes Medical Institute and Department of Molecular Biology, Massachusetts General Hospital, United States
- Broad Institute of MIT and Harvard, United States
- Analytic and Translational Genetics Unit, Center for Genomic Medicine, Massachusetts General Hospital, United States
| | - Masahiro Kanai
- Broad Institute of MIT and Harvard, United States
- Analytic and Translational Genetics Unit, Center for Genomic Medicine, Massachusetts General Hospital, United States
| | - Timothy J Durham
- Howard Hughes Medical Institute and Department of Molecular Biology, Massachusetts General Hospital, United States
- Broad Institute of MIT and Harvard, United States
| | - Kristin Tsuo
- Broad Institute of MIT and Harvard, United States
- Analytic and Translational Genetics Unit, Center for Genomic Medicine, Massachusetts General Hospital, United States
| | - Jason G McCoy
- Howard Hughes Medical Institute and Department of Molecular Biology, Massachusetts General Hospital, United States
- Broad Institute of MIT and Harvard, United States
| | - Patrick F Chinnery
- Department of Clinical Neurosciences & MRC Mitochondrial Biology Unit, University of Cambridge, United Kingdom
| | - Konrad J Karczewski
- Broad Institute of MIT and Harvard, United States
- Analytic and Translational Genetics Unit, Center for Genomic Medicine, Massachusetts General Hospital, United States
| | - Sarah E Calvo
- Howard Hughes Medical Institute and Department of Molecular Biology, Massachusetts General Hospital, United States
- Broad Institute of MIT and Harvard, United States
| | - Benjamin M Neale
- Broad Institute of MIT and Harvard, United States
- Analytic and Translational Genetics Unit, Center for Genomic Medicine, Massachusetts General Hospital, United States
| | - Vamsi K Mootha
- Howard Hughes Medical Institute and Department of Molecular Biology, Massachusetts General Hospital, United States
- Broad Institute of MIT and Harvard, United States
- Department of Systems Biology, Harvard Medical School, United States
| |
Collapse
|
20
|
Triant DA, Pearson WR. Comparison of detection methods and genome quality when quantifying nuclear mitochondrial insertions in vertebrate genomes. Front Genet 2022; 13:984513. [PMID: 36482890 PMCID: PMC9723244 DOI: 10.3389/fgene.2022.984513] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2022] [Accepted: 11/03/2022] [Indexed: 01/27/2024] Open
Abstract
The integration of mitochondrial genome fragments into the nuclear genome is well documented, and the transfer of these mitochondrial nuclear pseudogenes (numts) is thought to be an ongoing evolutionary process. With the increasing number of eukaryotic genomes available, genome-wide distributions of numts are often surveyed. However, inconsistencies in genome quality can reduce the accuracy of numt estimates, and methods used for identification can be complicated by the diverse sizes and ages of numts. Numts have been previously characterized in rodent genomes and it was postulated that they might be more prevalent in a group of voles with rapidly evolving karyotypes. Here, we examine 37 rodent genomes, and an additional 26 vertebrate genomes, while also considering numt detection methods. We identify numts using DNA:DNA and protein:translated-DNA similarity searches and compare numt distributions among rodent and vertebrate taxa to assess whether some groups are more susceptible to transfer. A combination of protein sequence comparisons (protein:translated-DNA) and BLASTN genomic DNA searches detect 50% more numts than genomic DNA:DNA searches alone. In addition, higher-quality RefSeq genomes produce lower estimates of numts than GenBank genomes, suggesting that lower quality genome assemblies can overestimate numts abundance. Phylogenetic analysis shows that mitochondrial transfers are not associated with karyotypic diversity among rodents. Surprisingly, we did not find a strong correlation between numt counts and genome size. Estimates using DNA: DNA analyses can underestimate the amount of mitochondrial DNA that is transferred to the nucleus.
Collapse
Affiliation(s)
- Deborah A. Triant
- Department of Biochemistry and Molecular Genetics, University of Virginia, Charlottesville, VA, United States
| | | |
Collapse
|
21
|
Mito-SiPE is a sequence-independent and PCR-free mtDNA enrichment method for accurate ultra-deep mitochondrial sequencing. Commun Biol 2022; 5:1269. [PMID: 36402890 PMCID: PMC9675811 DOI: 10.1038/s42003-022-04182-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2021] [Accepted: 10/27/2022] [Indexed: 11/21/2022] Open
Abstract
The analysis of somatic variation in the mitochondrial genome requires deep sequencing of mitochondrial DNA. This is ordinarily achieved by selective enrichment methods, such as PCR amplification or probe hybridization. These methods can introduce bias and are prone to contamination by nuclear-mitochondrial sequences (NUMTs), elements that can introduce artefacts into heteroplasmy analysis. We isolated intact mitochondria using differential centrifugation and alkaline lysis and subjected purified mitochondrial DNA to a sequence-independent and PCR-free method to obtain ultra-deep (>80,000X) sequencing coverage of the mitochondrial genome. This methodology avoids false-heteroplasmy calls that occur when long-range PCR amplification is used for mitochondrial DNA enrichment. Previously published methods employing mitochondrial DNA purification did not measure mitochondrial DNA enrichment or utilise high coverage short-read sequencing. Here, we describe a protocol that yields mitochondrial DNA and have quantified the increased level of mitochondrial DNA post-enrichment in 7 different mouse tissues. This method will enable researchers to identify changes in low frequency heteroplasmy without introducing PCR biases or NUMT contamination that are incorrectly identified as heteroplasmy when long-range PCR is used.
Collapse
|
22
|
Wei W, Schon KR, Elgar G, Orioli A, Tanguy M, Giess A, Tischkowitz M, Caulfield MJ, Chinnery PF. Nuclear-embedded mitochondrial DNA sequences in 66,083 human genomes. Nature 2022; 611:105-114. [PMID: 36198798 PMCID: PMC9630118 DOI: 10.1038/s41586-022-05288-7] [Citation(s) in RCA: 58] [Impact Index Per Article: 29.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2022] [Accepted: 08/29/2022] [Indexed: 02/02/2023]
Abstract
DNA transfer from cytoplasmic organelles to the cell nucleus is a legacy of the endosymbiotic event-the majority of nuclear-mitochondrial segments (NUMTs) are thought to be ancient, preceding human speciation1-3. Here we analyse whole-genome sequences from 66,083 people-including 12,509 people with cancer-and demonstrate the ongoing transfer of mitochondrial DNA into the nucleus, contributing to a complex NUMT landscape. More than 99% of individuals had at least one of 1,637 different NUMTs, with 1 in 8 individuals having an ultra-rare NUMT that is present in less than 0.1% of the population. More than 90% of the extant NUMTs that we evaluated inserted into the nuclear genome after humans diverged from apes. Once embedded, the sequences were no longer under the evolutionary constraint seen within the mitochondrion, and NUMT-specific mutations had a different mutational signature to mitochondrial DNA. De novo NUMTs were observed in the germline once in every 104 births and once in every 103 cancers. NUMTs preferentially involved non-coding mitochondrial DNA, linking transcription and replication to their origin, with nuclear insertion involving multiple mechanisms including double-strand break repair associated with PR domain zinc-finger protein 9 (PRDM9) binding. The frequency of tumour-specific NUMTs differed between cancers, including a probably causal insertion in a myxoid liposarcoma. We found evidence of selection against NUMTs on the basis of size and genomic location, shaping a highly heterogenous and dynamic human NUMT landscape.
Collapse
Affiliation(s)
- Wei Wei
- Department of Clinical Neuroscience, School of Clinical Medicine, University of Cambridge, Cambridge, UK
- Medical Research Council Mitochondrial Biology Unit, University of Cambridge, Cambridge, UK
| | - Katherine R Schon
- Department of Clinical Neuroscience, School of Clinical Medicine, University of Cambridge, Cambridge, UK
- Medical Research Council Mitochondrial Biology Unit, University of Cambridge, Cambridge, UK
- Academic Department of Medical Genetics, School of Clinical Medicine, University of Cambridge, Cambridge, UK
| | | | | | | | | | - Marc Tischkowitz
- Academic Department of Medical Genetics, School of Clinical Medicine, University of Cambridge, Cambridge, UK
| | - Mark J Caulfield
- William Harvey Research Institute, Queen Mary University of London, London, UK
| | - Patrick F Chinnery
- Department of Clinical Neuroscience, School of Clinical Medicine, University of Cambridge, Cambridge, UK.
- Medical Research Council Mitochondrial Biology Unit, University of Cambridge, Cambridge, UK.
| |
Collapse
|
23
|
Chiaratti MR, Chinnery PF. Modulating mitochondrial DNA mutations: factors shaping heteroplasmy in the germ line and somatic cells. Pharmacol Res 2022; 185:106466. [PMID: 36174964 DOI: 10.1016/j.phrs.2022.106466] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/31/2022] [Revised: 09/21/2022] [Accepted: 09/23/2022] [Indexed: 11/30/2022]
Abstract
Until recently it was thought that most humans only harbor one type of mitochondrial DNA (mtDNA), however, deep sequencing and single-cell analysis has shown the converse - that mixed populations of mtDNA (heteroplasmy) are the norm. This is important because heteroplasmy levels can change dramatically during transmission in the female germ line, leading to high levels causing severe mitochondrial diseases. There is also emerging evidence that low level mtDNA mutations contribute to common late onset diseases such as neurodegenerative disorders and cardiometabolic diseases because the inherited mutation levels can change within developing organs and non-dividing cells over time. Initial predictions suggested that the segregation of mtDNA heteroplasmy was largely stochastic, with an equal tendency for levels to increase or decrease. However, transgenic animal work and single-cell analysis have shown this not to be the case during germ-line transmission and in somatic tissues during life. Mutation levels in specific mtDNA regions can increase or decrease in different contexts and the underlying molecular mechanisms are starting to be unraveled. In this review we provide a synthesis of recent literature on the mechanisms of selection for and against mtDNA variants. We identify the most pertinent gaps in our understanding and suggest ways these could be addressed using state of the art techniques.
Collapse
Affiliation(s)
- Marcos R Chiaratti
- Departamento de Genética e Evolução, Centro de Ciências Biológicas e da Saúde, Universidade Federal de São Carlos, São Carlos, Brazil.
| | - Patrick F Chinnery
- Department of Clinical Neurosciences, School of Clinical Medicine, University of Cambridge, Cambridge Biomedical Campus, Cambridge, UK; Medical Research Council Mitochondrial Biology Unit, University of Cambridge, Cambridge Biomedical Campus, Cambridge, UK.
| |
Collapse
|
24
|
Keraite I, Becker P, Canevazzi D, Frias-López C, Dabad M, Tonda-Hernandez R, Paramonov I, Ingham MJ, Brun-Heath I, Leno J, Abulí A, Garcia-Arumí E, Heath SC, Gut M, Gut IG. A method for multiplexed full-length single-molecule sequencing of the human mitochondrial genome. Nat Commun 2022; 13:5902. [PMID: 36202811 PMCID: PMC9537161 DOI: 10.1038/s41467-022-33530-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2022] [Accepted: 09/21/2022] [Indexed: 11/09/2022] Open
Abstract
Methods to reconstruct the mitochondrial DNA (mtDNA) sequence using short-read sequencing come with an inherent bias due to amplification and mapping. They can fail to determine the phase of variants, to capture multiple deletions and to cover the mitochondrial genome evenly. Here we describe a method to target, multiplex and sequence at high coverage full-length human mitochondrial genomes as native single-molecules, utilizing the RNA-guided DNA endonuclease Cas9. Combining Cas9 induced breaks, that define the mtDNA beginning and end of the sequencing reads, as barcodes, we achieve high demultiplexing specificity and delineation of the full-length of the mtDNA, regardless of the structural variant pattern. The long-read sequencing data is analysed with a pipeline where our custom-developed software, baldur, efficiently detects single nucleotide heteroplasmy to below 1%, physically determines phase and can accurately disentangle complex deletions. Our workflow is a tool for studying mtDNA variation and will accelerate mitochondrial research. Accurate analysis of mitochondrial DNA is important for mitochondrial disease clinical research and diagnostics. Here, authors present a method using Cas9 cleavage, nanopore sequencing and a custom pipeline to identify pathogenic variants, deletions and accurately quantify heteroplasmy to below 1%.
Collapse
Affiliation(s)
- Ieva Keraite
- CNAG-CRG, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona, Spain
| | - Philipp Becker
- CNAG-CRG, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona, Spain.,Qiagen, Hilden, Germany
| | - Davide Canevazzi
- CNAG-CRG, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona, Spain
| | - Cristina Frias-López
- CNAG-CRG, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona, Spain
| | - Marc Dabad
- CNAG-CRG, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona, Spain
| | - Raúl Tonda-Hernandez
- CNAG-CRG, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona, Spain
| | - Ida Paramonov
- CNAG-CRG, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona, Spain
| | - Matthew John Ingham
- CNAG-CRG, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona, Spain
| | - Isabelle Brun-Heath
- Institute for Research in Biomedicine (IRB Barcelona) - The Barcelona Institute of Science and Technology (BIST), Barcelona, Spain.,Joint IRB-BSC Program in Computational Biology, Barcelona, Spain
| | - Jordi Leno
- Department of Clinical and Molecular Genetics and Rare Disease, Hospital Universitari Vall d'Hebron, Barcelona, Spain.,Medicine Genetics Group, VHIR, Hospital Universitari Vall d'Hebron, Barcelona, Spain
| | - Anna Abulí
- Department of Clinical and Molecular Genetics and Rare Disease, Hospital Universitari Vall d'Hebron, Barcelona, Spain.,Medicine Genetics Group, VHIR, Hospital Universitari Vall d'Hebron, Barcelona, Spain
| | - Elena Garcia-Arumí
- Department of Clinical and Molecular Genetics and Rare Disease, Hospital Universitari Vall d'Hebron, Barcelona, Spain.,Research Group on Neuromuscular and Mitochondrial Disorders, VHIR, Hospital Universitari Vall d'Hebron, Barcelona, Spain.,Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER), Instituto de Salud Carlos III, Barcelona, Spain
| | - Simon Charles Heath
- CNAG-CRG, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona, Spain.,Universitat Pompeu Fabra, Barcelona, Spain
| | - Marta Gut
- CNAG-CRG, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona, Spain. .,Universitat Pompeu Fabra, Barcelona, Spain.
| | - Ivo Glynne Gut
- CNAG-CRG, Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Barcelona, Spain. .,Universitat Pompeu Fabra, Barcelona, Spain.
| |
Collapse
|
25
|
Daly GT, Pastukh VM, Tan YB, Francis CM, Aggen CZ, Groark SC, Edwards C, Mulekar MS, Hamo M, Simmons JD, Kutcher ME, Hartsell EM, Dinwiddie DL, Turpin ZM, Bass HW, Roberts JT, Gillespie MN, Langley RJ. Novel attributes of cell-free plasma mitochondrial DNA in traumatic injury. Clin Transl Med 2022; 12:e1055. [PMID: 36245326 PMCID: PMC9574491 DOI: 10.1002/ctm2.1055] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2022] [Revised: 08/22/2022] [Accepted: 09/29/2022] [Indexed: 01/28/2023] Open
Affiliation(s)
- Grant T. Daly
- Department of PharmacologyUniversity of South Alabama College of MedicineMobileAlabamaUSA
| | - Viktor M. Pastukh
- Department of PharmacologyUniversity of South Alabama College of MedicineMobileAlabamaUSA
| | - Yong B. Tan
- Department of SurgeryUniversity of South Alabama Colleges of MedicineMobileAlabamaUSA
| | - C. Michael Francis
- Department of Physiology and Cell BiologyUniversity of South Alabama College of MedicineMobileAlabamaUSA
| | - C. Zack Aggen
- Department of PharmacologyUniversity of South Alabama College of MedicineMobileAlabamaUSA
| | - S. Chris Groark
- Department of PharmacologyUniversity of South Alabama College of MedicineMobileAlabamaUSA
| | - Carson Edwards
- Department of PharmacologyUniversity of South Alabama College of MedicineMobileAlabamaUSA
| | - Madhuri S. Mulekar
- Department of Mathematics and StatisticsUniversity of South Alabama Colleges of Medicine and Arts and SciencesMobileAlabamaUSA
| | - Mohammad Hamo
- Department of PharmacologyUniversity of South Alabama College of MedicineMobileAlabamaUSA
| | - Jon D. Simmons
- Department of SurgeryUniversity of South Alabama Colleges of MedicineMobileAlabamaUSA
| | - Matthew E. Kutcher
- Department of SurgeryUniversity of Mississippi School of MedicineJacksonMississippiUSA
| | - Emily M. Hartsell
- Department of PharmacologyUniversity of South Alabama College of MedicineMobileAlabamaUSA
| | - Darrell L. Dinwiddie
- Department of PediatricsUniversity of New Mexico School of MedicineAlbuquerqueNew MexicoUSA
| | - Zachary M. Turpin
- Department of Biological Science, College of Arts and SciencesFlorida State UniversityTallahasseeFloridaUSA
| | - Hank W. Bass
- Department of Biological Science, College of Arts and SciencesFlorida State UniversityTallahasseeFloridaUSA
| | - Justin T. Roberts
- Department of PharmacologyUniversity of South Alabama College of MedicineMobileAlabamaUSA
| | - Mark N. Gillespie
- Department of PharmacologyUniversity of South Alabama College of MedicineMobileAlabamaUSA
| | - Raymond J. Langley
- Department of PharmacologyUniversity of South Alabama College of MedicineMobileAlabamaUSA
| |
Collapse
|
26
|
Development and validation of a SYBR green-based mitochondrial DNA quantification method by following the MIQE and other guidelines. Leg Med (Tokyo) 2022; 58:102096. [PMID: 35689884 DOI: 10.1016/j.legalmed.2022.102096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2022] [Revised: 04/11/2022] [Accepted: 05/27/2022] [Indexed: 01/28/2023]
Abstract
In forensic mitochondrial DNA (mtDNA) analysis, quantitative PCR (qPCR) is usually performed to obtain high-quality sequence data for subsequent Sanger or massively parallel sequencing. Unlike methods for nuclear DNA quantification using qPCR, a calibrator is necessary to obtain mtDNA concentrations (i.e., copies/µL). Herein, we developed and validated a mtDNA quantification method based on a SYBR Green assay by following MIQE [Bustin et al., Clin. Chem. 55 (2009) 611-22] and other guidelines. Primers were designed to amplify nucleotide positions 16,190-16,420 in hypervariable region 1 for qPCR using PowerUp SYBR Green and QuantStudio 5. The optimized conditions were 0.3 µM each primer and an annealing temperature of 60 °C under a 2-step cycling protocol. K562 DNA at 100 pg/µL was converted into a mtDNA concentration of 16,400 copies/µL using linearized plasmid DNA. This mtDNA calibrator was obtained by cloning the synthesized DNA fragments of mtDNA (positions 16,140-16,470) containing a 100-bp inversion. The linear dynamic range of the K562 standard curve was 10,000-0.1 pg/µL (r2 ≥ 0.999). The accuracy was examined using NIST SRM 2372a, and its components A, B, and C were quantified with differences of -29.4%, -35.0%, and -22.0%, respectively, against the mtDNA concentrations calculated from published NIST data. We also examined the specificity of the primers, stability of the reaction mix, precision, tolerance against PCR inhibitors, and cross-reactivity against DNA from various animal taxa. Our newly developed mtDNA quantification method is expected to be useful for forensic mtDNA analysis.
Collapse
|
27
|
Post hoc deconvolution of human mitochondrial DNA mixtures by EMMA 2 using fine-tuned Phylotree nomenclature. Comput Struct Biotechnol J 2022; 20:3630-3638. [PMID: 35860401 PMCID: PMC9283771 DOI: 10.1016/j.csbj.2022.06.053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2022] [Revised: 06/24/2022] [Accepted: 06/25/2022] [Indexed: 11/23/2022] Open
Abstract
MtDNA mixtures are observed frequently and difficult to deconvolute. Most previous methods require raw data or quantitative information. EMMA 2 produces valid splittings from consensus sequences of any sequencing technology. EMMA 2 can deconvolute 2 and 3 person mixtures in a fast and traceable way.
In this paper we present a new algorithm for splitting (partial) human mitogenomes into components with high similarity to haplogroup motifs of Phylotree. The algorithm reads a (partial) mitogenome coded by the differences to the reference (rCRS) and outputs the estimated haplogroups of the putative components. The algorithm requires no special information on the raw data of the sequencing process and is therefore suited for the post hoc analysis of mixtures of any sequencing technology. The software EMMA 2 implementing the algorithm will be made available via the EMPOP (https://empop.online) database and extends the nine years old software EMMA for haplogrouping single mitogenomes to mixtures with at most three components.
Collapse
|
28
|
Parakatselaki ME, Zhu CT, Rand D, Ladoukakis ED. NUMTs Can Imitate Biparental Transmission of mtDNA-A Case in Drosophila melanogaster. Genes (Basel) 2022; 13:genes13061023. [PMID: 35741785 PMCID: PMC9222939 DOI: 10.3390/genes13061023] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2022] [Revised: 05/27/2022] [Accepted: 06/02/2022] [Indexed: 11/16/2022] Open
Abstract
mtDNA sequences can be incorporated into the nuclear genome and produce nuclear mitochondrial fragments (NUMTs), which resemble mtDNA in their sequence but are transmitted biparentally, like the nuclear genome. NUMTs can be mistaken as real mtDNA and may lead to the erroneous impression that mtDNA is biparentally transmitted. Here, we report a case of mtDNA heteroplasmy in a Drosophila melanogaster DGRP line, in which the one haplotype was biparentally transmitted in an autosomal manner. Given the sequence identity of this haplotype with the mtDNA, the crossing experiments led to uncertainty about whether heteroplasmy was real or an artifact due to a NUMT. More specific experiments revealed that there is a large NUMT insertion in the X chromosome of a specific DGRP line, imitating biparental inheritance of mtDNA. Our result suggests that studies on mtDNA heteroplasmy and on mtDNA inheritance should first exclude the possibility of NUMT interference in their data.
Collapse
Affiliation(s)
| | - Chen-Tseh Zhu
- Department of Ecology, Evolutionary and Organismal Biology, Brown University, Providence, RI 02912, USA; (C.-T.Z.); (D.R.)
| | - David Rand
- Department of Ecology, Evolutionary and Organismal Biology, Brown University, Providence, RI 02912, USA; (C.-T.Z.); (D.R.)
| | - Emmanuel D. Ladoukakis
- Department of Biology, University of Crete, Voutes University Campus, 70013 Iraklio, Greece;
- Correspondence: ; Tel.: +30-281-039-4067
| |
Collapse
|
29
|
Battle SL, Puiu D, Verlouw J, Broer L, Boerwinkle E, Taylor KD, Rotter JI, Rich SS, Grove ML, Pankratz N, Fetterman JL, Liu C, Arking D. A bioinformatics pipeline for estimating mitochondrial DNA copy number and heteroplasmy levels from whole genome sequencing data. NAR Genom Bioinform 2022; 4:lqac034. [PMID: 35591888 PMCID: PMC9112767 DOI: 10.1093/nargab/lqac034] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2021] [Revised: 03/16/2022] [Accepted: 04/25/2022] [Indexed: 12/19/2022] Open
Abstract
Mitochondrial diseases are a heterogeneous group of disorders that can be caused by mutations in the nuclear or mitochondrial genome. Mitochondrial DNA (mtDNA) variants may exist in a state of heteroplasmy, where a percentage of DNA molecules harbor a variant, or homoplasmy, where all DNA molecules have the same variant. The relative quantity of mtDNA in a cell, or copy number (mtDNA-CN), is associated with mitochondrial function, human disease, and mortality. To facilitate accurate identification of heteroplasmy and quantify mtDNA-CN, we built a bioinformatics pipeline that takes whole genome sequencing data and outputs mitochondrial variants, and mtDNA-CN. We incorporate variant annotations to facilitate determination of variant significance. Our pipeline yields uniform coverage by remapping to a circularized chrM and by recovering reads falsely mapped to nuclear-encoded mitochondrial sequences. Notably, we construct a consensus chrM sequence for each sample and recall heteroplasmy against the sample's unique mitochondrial genome. We observe an approximately 3-fold increased association with age for heteroplasmic variants in non-homopolymer regions and, are better able to capture genetic variation in the D-loop of chrM compared to existing software. Our bioinformatics pipeline more accurately captures features of mitochondrial genetics than existing pipelines that are important in understanding how mitochondrial dysfunction contributes to disease.
Collapse
Affiliation(s)
- Stephanie L Battle
- McKusick-Nathans Institute, Department of Genetic Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Daniela Puiu
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | | | - Joost Verlouw
- Department of Internal Medicine, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Linda Broer
- Department of Internal Medicine, Erasmus Medical Center, Rotterdam, The Netherlands
| | - Eric Boerwinkle
- Human Genetics Center, Department of Epidemiology, Human Genetics, and Environmental Sciences, School of Public Health, The University of Texas Health Science Center at Houston, Houston, TX, USA
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX, USA
| | - Kent D Taylor
- The Institute for Translational Genomics and Population Sciences, Department of Pediatrics, The Lundquist Institute for Biomedical Innovation at Harbor-UCLA Medical Center, Torrance, CA, USA
| | - Jerome I Rotter
- The Institute for Translational Genomics and Population Sciences, Department of Pediatrics, The Lundquist Institute for Biomedical Innovation at Harbor-UCLA Medical Center, Torrance, CA, USA
| | - Stephan S Rich
- Center for Public Health Genomics, University of Virginia, Charlottesville, VA, USA
| | - Megan L Grove
- Human Genetics Center, Department of Epidemiology, Human Genetics, and Environmental Sciences, School of Public Health, The University of Texas Health Science Center at Houston, Houston, TX, USA
| | - Nathan Pankratz
- Department of Laboratory Medicine and Pathology, University of Minnesota, Minneapolis, MN, USA
| | - Jessica L Fetterman
- Evans Department of Medicine and the Whitaker Cardiovascular Institute, Boston University School of Medicine, Boston, MA, USA
| | - Chunyu Liu
- Framingham Heart Study, Boston University School of Medicine, Boston, MA, USA
- Department of Biostatistics, Boston University School of Public Health, Boston, MA, USA
| | - Dan E Arking
- McKusick-Nathans Institute, Department of Genetic Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| |
Collapse
|
30
|
Kołomański M, Szyda J, Frąszczak M, Mielczarek M. DNA sequence features underlying large-scale duplications and deletions in human. J Appl Genet 2022; 63:527-533. [PMID: 35590085 PMCID: PMC9365719 DOI: 10.1007/s13353-022-00704-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2021] [Revised: 03/22/2022] [Accepted: 05/05/2022] [Indexed: 11/25/2022]
Abstract
Copy number variants (CNVs) may cover up to 12% of the whole genome and have substantial impact on phenotypes. We used 5867 duplications and 33,181 deletions available from the 1000 Genomes Project to characterise genomic regions vulnerable to CNV formation and to identify sequence features characteristic for those regions. The GC content for deletions was lower and for duplications was higher than for randomly selected regions. In regions flanking deletions and downstream of duplications, content was higher than in the random sequences, but upstream of duplication content was lower. In duplications and downstream of deletion regions, the percentage of low-complexity sequences was not different from the randomised data. In deletions and upstream of CNVs, it was higher, while for downstream of duplications, it was lower as compared to random sequences. The majority of CNVs intersected with genic regions — mainly with introns. GC content may be associated with CNV formation and CNVs, especially duplications are initiated in low-complexity regions. Moreover, CNVs located or overlapped with introns indicate their role in shaping intron variability. Genic CNV regions were enriched in many essential biological processes such as cell adhesion, synaptic transmission, transport, cytoskeleton organization, immune response and metabolic mechanisms, which indicates that these large-scaled variants play important biological roles.
Collapse
Affiliation(s)
- Mateusz Kołomański
- Biostatistics Group, Department of Genetics, Wroclaw University of Environmental and Life Sciences, Wroclaw, Poland
| | - Joanna Szyda
- Biostatistics Group, Department of Genetics, Wroclaw University of Environmental and Life Sciences, Wroclaw, Poland
| | - Magdalena Frąszczak
- Biostatistics Group, Department of Genetics, Wroclaw University of Environmental and Life Sciences, Wroclaw, Poland
| | - Magda Mielczarek
- Biostatistics Group, Department of Genetics, Wroclaw University of Environmental and Life Sciences, Wroclaw, Poland.
| |
Collapse
|
31
|
Lüth T, Schaake S, Grünewald A, May P, Trinh J, Weissensteiner H. Benchmarking Low-Frequency Variant Calling With Long-Read Data on Mitochondrial DNA. Front Genet 2022; 13:887644. [PMID: 35664331 PMCID: PMC9161029 DOI: 10.3389/fgene.2022.887644] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2022] [Accepted: 04/18/2022] [Indexed: 11/13/2022] Open
Abstract
Background: Sequencing quality has improved over the last decade for long-reads, allowing for more accurate detection of somatic low-frequency variants. In this study, we used mixtures of mitochondrial samples with different haplogroups (i.e., a specific set of mitochondrial variants) to investigate the applicability of nanopore sequencing for low-frequency single nucleotide variant detection. Methods: We investigated the impact of base-calling, alignment/mapping, quality control steps, and variant calling by comparing the results to a previously derived short-read gold standard generated on the Illumina NextSeq. For nanopore sequencing, six mixtures of four different haplotypes were prepared, allowing us to reliably check for expected variants at the predefined 5%, 2%, and 1% mixture levels. We used two different versions of Guppy for base-calling, two aligners (i.e., Minimap2 and Ngmlr), and three variant callers (i.e., Mutserve2, Freebayes, and Nanopanel2) to compare low-frequency variants. We used F1 score measurements to assess the performance of variant calling. Results: We observed a mean read length of 11 kb and a mean overall read quality of 15. Ngmlr showed not only higher F1 scores but also higher allele frequencies (AF) of false-positive calls across the mixtures (mean F1 score = 0.83; false-positive allele frequencies < 0.17) compared to Minimap2 (mean F1 score = 0.82; false-positive AF < 0.06). Mutserve2 had the highest F1 scores (5% level: F1 score >0.99, 2% level: F1 score >0.54, and 1% level: F1 score >0.70) across all callers and mixture levels. Conclusion: We here present the benchmarking for low-frequency variant calling with nanopore sequencing by identifying current limitations.
Collapse
Affiliation(s)
- Theresa Lüth
- Institute of Neurogenetics, University of Lübeck and University Hospital Schleswig-Holstein, Lübeck, Germany
| | - Susen Schaake
- Institute of Neurogenetics, University of Lübeck and University Hospital Schleswig-Holstein, Lübeck, Germany
| | - Anne Grünewald
- Institute of Neurogenetics, University of Lübeck and University Hospital Schleswig-Holstein, Lübeck, Germany
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Patrick May
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Joanne Trinh
- Institute of Neurogenetics, University of Lübeck and University Hospital Schleswig-Holstein, Lübeck, Germany
- *Correspondence: Joanne Trinh, ; Hansi Weissensteiner,
| | - Hansi Weissensteiner
- Institute of Genetic Epidemiology, Medical University of Innsbruck, Innsbruck, Austria
- *Correspondence: Joanne Trinh, ; Hansi Weissensteiner,
| |
Collapse
|
32
|
Fujii K, Mita Y, Watahiki H, Fukagawa T, Kitayama T, Mizuno N, Nakahara H, Sekiguchi K. Poly_NumtS_430 or HSA_NumtS_587 observed in massively parallel sequencing of the mitochondrial HV1 and HV2 regions. Forensic Sci Int Genet 2022; 59:102717. [DOI: 10.1016/j.fsigen.2022.102717] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2021] [Revised: 03/12/2022] [Accepted: 05/01/2022] [Indexed: 11/04/2022]
|
33
|
Maniego J, Pesko B, Habershon-Butcher J, Hincks P, Taylor P, Tozaki T, Ohnuma A, Stewart G, Proudman C, Ryder E. Use of mitochondrial sequencing to detect gene doping in horses via gene editing and somatic cell nuclear transfer. Drug Test Anal 2022; 14:1429-1437. [PMID: 35362263 DOI: 10.1002/dta.3267] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2022] [Revised: 03/24/2022] [Accepted: 03/29/2022] [Indexed: 11/09/2022]
Abstract
Gene editing and subsequent cloning techniques offer great potential not only in genetic disease correction in domestic animals, but also in livestock production by enhancement of desirable traits. The existence of the technology, however, leaves it open to potential misuse in performance-led sports such as horseracing and other equestrian events. Recent advances in equine gene editing, regarding the generation of gene-edited embryos using CRISPR/Cas9 technology and somatic cell nuclear transfer, has highlighted the need to develop tools to detect potential prohibited use of the technology. One possible method involves the characterisation of the mitochondrial genome (which is not routinely preserved during cloning) and comparing it to the sequence of the registered dam. We present here our approach to whole-mitochondrial sequencing using tiled long-range PCR and next-generation sequencing. To determine whether the background mutation rate in the mitochondrial genome could potentially confound results, we sequenced ten sets of dam and foal duos. We found variation between duos but none within duos, indicating that this method is feasible for future screening systems. Analysis of WGS data from over one hundred Thoroughbred horses revealed wide variation in the mitochondria sequence within the breed, further displaying the utility of this approach.
Collapse
Affiliation(s)
- Jillian Maniego
- Sport and Specialised Analytical Services, LGC, Newmarket Road, Fordham, Cambridgeshire, UK
| | - Bogusia Pesko
- Sport and Specialised Analytical Services, LGC, Newmarket Road, Fordham, Cambridgeshire, UK
| | | | - Pamela Hincks
- Sport and Specialised Analytical Services, LGC, Newmarket Road, Fordham, Cambridgeshire, UK
| | - Polly Taylor
- Sport and Specialised Analytical Services, LGC, Newmarket Road, Fordham, Cambridgeshire, UK
| | - Teruaki Tozaki
- Genetic Analysis Department, Laboratory of Racing Chemistry, Utsunomiya, Japan
| | - Aoi Ohnuma
- Genetic Analysis Department, Laboratory of Racing Chemistry, Utsunomiya, Japan
| | - Graham Stewart
- School of Biosciences and Medicine, University of Surrey, Guildford, UK
| | - Christopher Proudman
- School of Veterinary Medicine, Daphne Jackson Road, University of Surrey, Guildford, UK
| | - Edward Ryder
- Sport and Specialised Analytical Services, LGC, Newmarket Road, Fordham, Cambridgeshire, UK
| |
Collapse
|
34
|
Sturk-Andreaggi K, Ring JD, Ameur A, Gyllensten U, Bodner M, Parson W, Marshall C, Allen M. The Value of Whole-Genome Sequencing for Mitochondrial DNA Population Studies: Strategies and Criteria for Extracting High-Quality Mitogenome Haplotypes. Int J Mol Sci 2022; 23:ijms23042244. [PMID: 35216360 PMCID: PMC8876724 DOI: 10.3390/ijms23042244] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2022] [Revised: 02/08/2022] [Accepted: 02/09/2022] [Indexed: 02/04/2023] Open
Abstract
Whole-genome sequencing (WGS) data present a readily available resource for mitochondrial genome (mitogenome) haplotypes that can be utilized for genetics research including population studies. However, the reconstruction of the mitogenome is complicated by nuclear mitochondrial DNA (mtDNA) segments (NUMTs) that co-align with the mtDNA sequences and mimic authentic heteroplasmy. Two minimum variant detection thresholds, 5% and 10%, were assessed for the ability to produce authentic mitogenome haplotypes from a previously generated WGS dataset. Variants associated with NUMTs were detected in the mtDNA alignments for 91 of 917 (~8%) Swedish samples when the 5% frequency threshold was applied. The 413 observed NUMT variants were predominantly detected in two regions (nps 12,612–13,105 and 16,390–16,527), which were consistent with previously documented NUMTs. The number of NUMT variants was reduced by ~97% (400) using a 10% frequency threshold. Furthermore, the 5% frequency data were inconsistent with a platinum-quality mitogenome dataset with respect to observed heteroplasmy. These analyses illustrate that a 10% variant detection threshold may be necessary to ensure the generation of reliable mitogenome haplotypes from WGS data resources.
Collapse
Affiliation(s)
- Kimberly Sturk-Andreaggi
- Department of Immunology Genetics and Pathology, Uppsala University, Uppsala 751 08, Sweden; (A.A.); (U.G.)
- Armed Forces Medical Examiner System’s Armed Forces DNA Identification Laboratory (AFMES-AFDIL), Dover Air Force Base, Dover, DE 19902, USA; (J.D.R.); (C.M.)
- SNA International, LLC, Alexandria, VA 22314, USA
- Correspondence: (K.S.-A.); (M.A.)
| | - Joseph D. Ring
- Armed Forces Medical Examiner System’s Armed Forces DNA Identification Laboratory (AFMES-AFDIL), Dover Air Force Base, Dover, DE 19902, USA; (J.D.R.); (C.M.)
- SNA International, LLC, Alexandria, VA 22314, USA
| | - Adam Ameur
- Department of Immunology Genetics and Pathology, Uppsala University, Uppsala 751 08, Sweden; (A.A.); (U.G.)
| | - Ulf Gyllensten
- Department of Immunology Genetics and Pathology, Uppsala University, Uppsala 751 08, Sweden; (A.A.); (U.G.)
| | - Martin Bodner
- Institute of Legal Medicine, Medical University of Innsbruck, Innsbruck 6020, Austria; (M.B.); (W.P.)
| | - Walther Parson
- Institute of Legal Medicine, Medical University of Innsbruck, Innsbruck 6020, Austria; (M.B.); (W.P.)
- Forensic Science Program, The Pennsylvania State University, University Park, PA 16801, USA
| | - Charla Marshall
- Armed Forces Medical Examiner System’s Armed Forces DNA Identification Laboratory (AFMES-AFDIL), Dover Air Force Base, Dover, DE 19902, USA; (J.D.R.); (C.M.)
- SNA International, LLC, Alexandria, VA 22314, USA
- Forensic Science Program, The Pennsylvania State University, University Park, PA 16801, USA
| | - Marie Allen
- Department of Immunology Genetics and Pathology, Uppsala University, Uppsala 751 08, Sweden; (A.A.); (U.G.)
- Correspondence: (K.S.-A.); (M.A.)
| |
Collapse
|
35
|
Hazkani-Covo E. A Burst of Numt Insertion in the Dasyuridae Family During Marsupial Evolution. Front Ecol Evol 2022. [DOI: 10.3389/fevo.2022.844443] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Nuclear pseudogenes of mitochondrial origin (numts) are common in all eukaryotes. Our previous scan of numts in sequenced nuclear genomes suggested that the highest numt content currently known in animals is that in the gray short-tailed opossum. The present work sought to determine numt content in marsupials and to compare it to those in placental and monothematic mammals as well as in non-mammalian vertebrates. To achieve this, 70 vertebrate species with available nuclear and mitochondrial genomes were scanned for numt content. An extreme numt content was found in the Dasyuridae, with 3,450 in Sarcophilus harrisii (1,955 kb) and 2,813 in Antechinus flavipes (847 kb). The evolutionarily closest species analyzed, the extinct Thylacinus cynocephalus belonging to the Thylacindae family, had only 435 numts (238 kb). These two Dasyuridae genomes featured the highest numt content identified in animals to date. A phylogenetic analysis of numts longer than 300 bp, using a Diprotodonita mitochondrial tree, indicated a burst of numt insertion that began before the divergence of the Dasyurini and Phascogalini, reaching a peak in the early evolution of the two tribes. No comparable increase was found in the early divergent species T. cynocephalus. Divergence of the Dasyuridae tribes has been previously dated to shortly after the Miocene climate transition, characterized by a rapid temperature decline. Interestingly, deviation from optimal growth temperature is one of the environmental factors reported to increase numt insertions in a laboratory setting.
Collapse
|
36
|
Laricchia KM, Lake NJ, Watts NA, Shand M, Haessly A, Gauthier L, Benjamin D, Banks E, Soto J, Garimella K, Emery J, Rehm HL, MacArthur DG, Tiao G, Lek MV, Mootha VK, Calvo SE. Mitochondrial DNA variation across 56,434 individuals in gnomAD. Genome Res 2022; 32:569-582. [PMID: 35074858 PMCID: PMC8896463 DOI: 10.1101/gr.276013.121] [Citation(s) in RCA: 27] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2021] [Accepted: 01/19/2022] [Indexed: 11/25/2022]
Abstract
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: (i) detecting homoplasmic and heteroplasmic variants, present respectively in all or a fraction of mtDNA molecules, (ii) circular mtDNA genome, and (iii) 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 greater than 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. 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.
Collapse
|
37
|
Bicci I, Calabrese C, Golder ZJ, Gomez-Duran A, Chinnery PF. Single-molecule mitochondrial DNA sequencing shows no evidence of CpG methylation in human cells and tissues. Nucleic Acids Res 2021; 49:12757-12768. [PMID: 34850165 PMCID: PMC8682748 DOI: 10.1093/nar/gkab1179] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2021] [Revised: 10/29/2021] [Accepted: 11/18/2021] [Indexed: 01/11/2023] Open
Abstract
Methylation on CpG residues is one of the most important epigenetic modifications of nuclear DNA, regulating gene expression. Methylation of mitochondrial DNA (mtDNA) has been studied using whole genome bisulfite sequencing (WGBS), but recent evidence has uncovered technical issues which introduce a potential bias during methylation quantification. Here, we validate the technical concerns of WGBS, and develop and assess the accuracy of a new protocol for mtDNA nucleotide variant-specific methylation using single-molecule Oxford Nanopore Sequencing (ONS). Our approach circumvents confounders by enriching for full-length molecules over nuclear DNA. Variant calling analysis against showed that 99.5% of homoplasmic mtDNA variants can be reliably identified providing there is adequate sequencing depth. We show that some of the mtDNA methylation signal detected by ONS is due to sequence-specific false positives introduced by the technique. The residual signal was observed across several human primary and cancer cell lines and multiple human tissues, but was always below the error threshold modelled using negative controls. We conclude that there is no evidence for CpG methylation in human mtDNA, thus resolving previous controversies. Additionally, we developed a reliable protocol to study epigenetic modifications of mtDNA at single-molecule and single-base resolution, with potential applications beyond CpG methylation.
Collapse
Affiliation(s)
- Iacopo Bicci
- MRC-Mitochondrial Biology Unit, The Keith Peters Building, Cambridge CB2 0XY, UK.,Department of Clinical Neurosciences, University of Cambridge, Cambridge Biomedical Campus, Hills Road, Cambridge CB2 0XY, UK
| | - Claudia Calabrese
- MRC-Mitochondrial Biology Unit, The Keith Peters Building, Cambridge CB2 0XY, UK.,Department of Clinical Neurosciences, University of Cambridge, Cambridge Biomedical Campus, Hills Road, Cambridge CB2 0XY, UK
| | - Zoe J Golder
- MRC-Mitochondrial Biology Unit, The Keith Peters Building, Cambridge CB2 0XY, UK.,Department of Clinical Neurosciences, University of Cambridge, Cambridge Biomedical Campus, Hills Road, Cambridge CB2 0XY, UK
| | - Aurora Gomez-Duran
- MRC-Mitochondrial Biology Unit, The Keith Peters Building, Cambridge CB2 0XY, UK.,Department of Clinical Neurosciences, University of Cambridge, Cambridge Biomedical Campus, Hills Road, Cambridge CB2 0XY, UK.,Centro de Investigaciones Biológicas Margarita Salas. Spanish National Research Council, Madrid, Spain
| | - Patrick F Chinnery
- MRC-Mitochondrial Biology Unit, The Keith Peters Building, Cambridge CB2 0XY, UK.,Department of Clinical Neurosciences, University of Cambridge, Cambridge Biomedical Campus, Hills Road, Cambridge CB2 0XY, UK
| |
Collapse
|
38
|
Emser SV, Schaschl H, Millesi E, Steinborn R. Extension of Mitogenome Enrichment Based on Single Long-Range PCR: mtDNAs and Putative Mitochondrial-Derived Peptides of Five Rodent Hibernators. Front Genet 2021; 12:685806. [PMID: 35027919 PMCID: PMC8749263 DOI: 10.3389/fgene.2021.685806] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2021] [Accepted: 11/10/2021] [Indexed: 12/14/2022] Open
Abstract
Enriching mitochondrial DNA (mtDNA) for sequencing entire mitochondrial genomes (mitogenomes) can be achieved by single long-range PCR. This avoids interference from the omnipresent nuclear mtDNA sequences (NUMTs). The approach is currently restricted to the use of samples collected from humans and ray-finned fishes. Here, we extended the use of single long-range PCR by introducing back-to-back oligonucleotides that target a sequence of extraordinary homology across vertebrates. The assay was applied to five hibernating rodents, namely alpine marmot, Arctic and European ground squirrels, and common and garden dormice, four of which have not been fully sequenced before. Analysis of the novel mitogenomes focussed on the prediction of mitochondrial-derived peptides (MDPs) providing another level of information encoded by mtDNA. The comparison of MOTS-c, SHLP4 and SHLP6 sequences across vertebrate species identified segments of high homology that argue for future experimentation. In addition, we evaluated four candidate polymorphisms replacing an amino acid in mitochondrially encoded subunits of the oxidative phosphorylation (OXPHOS) system that were reported in relation to cold-adaptation. No obvious pattern was found for the diverse sets of mammalian species that either apply daily or multiday torpor or otherwise cope with cold. In summary, our single long-range PCR assay applying a pair of back-to-back primers that target a consensus sequence motif of Vertebrata has potential to amplify (intact) mitochondrial rings present in templates from a taxonomically diverse range of vertebrates. It could be promising for studying novel mitogenomes, mitotypes of a population and mitochondrial heteroplasmy in a sensitive, straightforward and flexible manner.
Collapse
Affiliation(s)
- Sarah V. Emser
- Genomics Core Facility, VetCore, University of Veterinary Medicine, Vienna, Austria
- Department of Behavioral and Cognitive Biology, University of Vienna, Vienna, Austria
| | - Helmut Schaschl
- Department of Evolutionary Anthropology, University of Vienna, Vienna, Austria
| | - Eva Millesi
- Department of Behavioral and Cognitive Biology, University of Vienna, Vienna, Austria
| | - Ralf Steinborn
- Genomics Core Facility, VetCore, University of Veterinary Medicine, Vienna, Austria
| |
Collapse
|
39
|
Singh LN, Kao SH, Wallace DC. Unlocking the Complexity of Mitochondrial DNA: A Key to Understanding Neurodegenerative Disease Caused by Injury. Cells 2021; 10:cells10123460. [PMID: 34943968 PMCID: PMC8715673 DOI: 10.3390/cells10123460] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2021] [Revised: 12/03/2021] [Accepted: 12/06/2021] [Indexed: 12/12/2022] Open
Abstract
Neurodegenerative disorders that are triggered by injury typically have variable and unpredictable outcomes due to the complex and multifactorial cascade of events following the injury and during recovery. Hence, several factors beyond the initial injury likely contribute to the disease progression and pathology, and among these are genetic factors. Genetics is a recognized factor in determining the outcome of common neurodegenerative diseases. The role of mitochondrial genetics and function in traditional neurodegenerative diseases, such as Alzheimer’s and Parkinson’s diseases, is well-established. Much less is known about mitochondrial genetics, however, regarding neurodegenerative diseases that result from injuries such as traumatic brain injury and ischaemic stroke. We discuss the potential role of mitochondrial DNA genetics in the progression and outcome of injury-related neurodegenerative diseases. We present a guide for understanding mitochondrial genetic variation, along with the nuances of quantifying mitochondrial DNA variation. Evidence supporting a role for mitochondrial DNA as a risk factor for neurodegenerative disease is also reviewed and examined. Further research into the impact of mitochondrial DNA on neurodegenerative disease resulting from injury will likely offer key insights into the genetic factors that determine the outcome of these diseases together with potential targets for treatment.
Collapse
Affiliation(s)
- Larry N. Singh
- Center for Mitochondrial and Epigenomic Medicine, Division of Human Genetics, The Children’s Hospital of Philadelphia, Philadelphia, PA 19104, USA;
- Correspondence:
| | - Shih-Han Kao
- Resuscitation Science Center, The Children’s Hospital of Philadelphia, Philadelphia, PA 19104, USA;
| | - Douglas C. Wallace
- Center for Mitochondrial and Epigenomic Medicine, Division of Human Genetics, The Children’s Hospital of Philadelphia, Philadelphia, PA 19104, USA;
- Department of Pediatrics, Division of Human Genetics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104, USA
| |
Collapse
|
40
|
Laaksonen J, Mishra PP, Seppälä I, Raitoharju E, Marttila S, Mononen N, Lyytikäinen LP, Kleber ME, Delgado GE, Lepistö M, Almusa H, Ellonen P, Lorkowski S, März W, Hutri-Kähönen N, Raitakari O, Kähönen M, Salonen JT, Lehtimäki T. Mitochondrial genome-wide analysis of nuclear DNA methylation quantitative trait loci. Hum Mol Genet 2021; 31:1720-1732. [PMID: 35077545 PMCID: PMC9122653 DOI: 10.1093/hmg/ddab339] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2021] [Revised: 11/12/2021] [Accepted: 11/16/2021] [Indexed: 11/13/2022] Open
Abstract
Abstract
Mitochondria have a complex communication network with the surrounding cell and can alter nuclear DNA methylation (DNAm). Variation in the mitochondrial DNA (mtDNA) has also been linked to differential DNAm. Genome-wide association studies have identified numerous DNAm quantitative trait loci, but these studies have not examined the mitochondrial genome. Herein, we quantified nuclear DNAm from blood and conducted a mitochondrial genome-wide association study of DNAm, with an additional emphasis on sex- and prediabetes-specific heterogeneity. We used the Young Finns Study (n = 926) with sequenced mtDNA genotypes as a discovery sample and sought replication in the Ludwigshafen Risk and Cardiovascular Health study (n = 2317). We identified numerous significant associations in the discovery phase (P < 10−9), but they were not replicated when accounting for multiple testing. In total, 27 associations were nominally replicated with a P < 0.05. The replication analysis presented no evidence of sex- or prediabetes-specific heterogeneity. The 27 associations were included in a joint meta-analysis of the two cohorts, and 19 DNAm sites associated with mtDNA variants, while four other sites showed haplogroup associations. An expression quantitative trait methylation analysis was performed for the identified DNAm sites, pinpointing two statistically significant associations. This study provides evidence of a mitochondrial genetic control of nuclear DNAm with little evidence found for sex- and prediabetes-specific effects. The lack of a comparable mtDNA data set for replication is a limitation in our study and further studies are needed to validate our results.
Collapse
Affiliation(s)
- Jaakko Laaksonen
- To whom correspondence should be addressed at: Department of Clinical Chemistry, Fimlab Laboratories and Finnish Cardiovascular Research Center Tampere, Faculty of Medicine and Health Technology, Tampere University, Arvo Ylpön katu 34, PO Box 100, Tampere FI-33014, Finland. Tel: +358 504080774; E-mail:
| | - Pashupati P Mishra
- Department of Clinical Chemistry, Fimlab Laboratories and Finnish Cardiovascular Research Center Tampere, Faculty of Medicine and Health Technology, Tampere University, Tampere 33520, Finland
| | - Ilkka Seppälä
- Department of Clinical Chemistry, Fimlab Laboratories and Finnish Cardiovascular Research Center Tampere, Faculty of Medicine and Health Technology, Tampere University, Tampere 33520, Finland
| | - Emma Raitoharju
- Department of Clinical Chemistry, Fimlab Laboratories and Finnish Cardiovascular Research Center Tampere, Faculty of Medicine and Health Technology, Tampere University, Tampere 33520, Finland
- Molecular Epidemiology, Faculty of Medicine and Health Technology, Tampere University, Tampere 33520, Finland
| | - Saara Marttila
- Molecular Epidemiology, Faculty of Medicine and Health Technology, Tampere University, Tampere 33520, Finland
- Gerontology Research Center, Tampere University, Tampere 33520, Finland
| | - Nina Mononen
- Department of Clinical Chemistry, Fimlab Laboratories and Finnish Cardiovascular Research Center Tampere, Faculty of Medicine and Health Technology, Tampere University, Tampere 33520, Finland
| | - Leo-Pekka Lyytikäinen
- Department of Clinical Chemistry, Fimlab Laboratories and Finnish Cardiovascular Research Center Tampere, Faculty of Medicine and Health Technology, Tampere University, Tampere 33520, Finland
| | - Marcus E Kleber
- Vth Department of Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim 68167, Germany
| | - Graciela E Delgado
- Vth Department of Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim 68167, Germany
| | - Maija Lepistö
- Institute for Molecular Medicine (FIMM), University of Helsinki, Helsinki 00290, Finland
| | - Henrikki Almusa
- Institute for Molecular Medicine (FIMM), University of Helsinki, Helsinki 00290, Finland
| | - Pekka Ellonen
- Institute for Molecular Medicine (FIMM), University of Helsinki, Helsinki 00290, Finland
| | - Stefan Lorkowski
- Institute of Nutritional Sciences, Friedrich Schiller University Jena, Jena 07743, Germany
- Competence Cluster for Nutrition and Cardiovascular Health (nutriCARD) Halle-Jena-Leipzig, Jena 07743, Germany
| | - Winfried März
- Vth Department of Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim 68167, Germany
- Competence Cluster for Nutrition and Cardiovascular Health (nutriCARD) Halle-Jena-Leipzig, Jena 07743, Germany
- SYNLAB Academy, SYNLAB Holding Deutschland GmbH, Augsburg 86156, Germany
- Clinical Institute of Medical and Chemical Laboratory Diagnostics, Medical University of Graz, Graz 8010, Austria
| | - Nina Hutri-Kähönen
- Tampere Centre for Skills Training and Simulation, Tampere University, Tampere 33520, Finland
| | - Olli Raitakari
- Centre for Population Health Research, University of Turku and Turku University Hospital, Turku 20520, Finland
- Research Centre for Applied and Preventive Cardiovascular Medicine, University of Turku, Turku 20520, Finland
- Department of Clinical Physiology and Nuclear Medicine, Turku University Hospital, Turku 20520, Finland
| | - Mika Kähönen
- Department of Clinical Physiology, Tampere University Hospital, Tampere 33520, Finland
- Finnish Cardiovascular Research Center Tampere, Faculty of Medicine and Health Technology, Tampere University, Tampere 33520, Finland
| | - Jukka T Salonen
- Department of Public Health, Faculty of Medicine, University of Helsinki, Helsinki 00014, Finland
- MAS-Metabolic Analytical Services Oy, Helsinki 00990, Finland
| | - Terho Lehtimäki
- Department of Clinical Chemistry, Fimlab Laboratories and Finnish Cardiovascular Research Center Tampere, Faculty of Medicine and Health Technology, Tampere University, Tampere 33520, Finland
| |
Collapse
|
41
|
Singh LN, Ennis B, Loneragan B, Tsao NL, Lopez Sanchez MIG, Li J, Acheampong P, Tran O, Trounce IA, Zhu Y, Potluri P, Emanuel BS, Rader DJ, Arany Z, Damrauer SM, Resnick AC, Anderson SA, Wallace DC. MitoScape: A big-data, machine-learning platform for obtaining mitochondrial DNA from next-generation sequencing data. PLoS Comput Biol 2021; 17:e1009594. [PMID: 34762648 PMCID: PMC8610268 DOI: 10.1371/journal.pcbi.1009594] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2021] [Revised: 11/23/2021] [Accepted: 10/27/2021] [Indexed: 11/18/2022] Open
Abstract
The growing number of next-generation sequencing (NGS) data presents a unique opportunity to study the combined impact of mitochondrial and nuclear-encoded genetic variation in complex disease. Mitochondrial DNA variants and in particular, heteroplasmic variants, are critical for determining human disease severity. While there are approaches for obtaining mitochondrial DNA variants from NGS data, these software do not account for the unique characteristics of mitochondrial genetics and can be inaccurate even for homoplasmic variants. We introduce MitoScape, a novel, big-data, software for extracting mitochondrial DNA sequences from NGS. MitoScape adopts a novel departure from other algorithms by using machine learning to model the unique characteristics of mitochondrial genetics. We also employ a novel approach of using rho-zero (mitochondrial DNA-depleted) data to model nuclear-encoded mitochondrial sequences. We showed that MitoScape produces accurate heteroplasmy estimates using gold-standard mitochondrial DNA data. We provide a comprehensive comparison of the most common tools for obtaining mtDNA variants from NGS and showed that MitoScape had superior performance to compared tools in every statistically category we compared, including false positives and false negatives. By applying MitoScape to common disease examples, we illustrate how MitoScape facilitates important heteroplasmy-disease association discoveries by expanding upon a reported association between hypertrophic cardiomyopathy and mitochondrial haplogroup T in men (adjusted p-value = 0.003). The improved accuracy of mitochondrial DNA variants produced by MitoScape will be instrumental in diagnosing disease in the context of personalized medicine and clinical diagnostics. Recent studies have highlighted the importance of mitochondrial DNA variation in both primary mitochondrial disease and complex, human pathology including COVID-19, and space-flight stress. The vast amount of existing, next-generation sequencing (NGS) data can be leveraged to interrogate both nuclear and mitochondrial DNA (mtDNA) sequence simultaneously, allowing for analysis of the interplay between mitochondrial and nuclear encoded genes in mitochondrial function. Identifying mtDNA sequence accurately is complicated by the presence of nuclear encoded mitochondrial sequences (NUMTs), which are homologous to mtDNA. Current software for analyzing mtDNA from NGS do not accurately model the unique characteristics of mitochondrial genetics. We introduce MitoScape, a novel, big-data, software which models mitochondrial genetics through machine learning to accurately identify mtDNA sequence from NGS data. MitoScape takes advantage of rho-zero cell data to model the characteristics of NUMTs. We show that MitoScape produces more accurate heteroplasmy estimates compared to published software. We provide an example of applying MitoScape in replicating an association between hypertrophic cardiomyopathy and mitochondrial haplogroup T in men. MitoScape is an important contribution to mitochondrial genomics allowing for accurate mtDNA variants, and the ability to tailor mtDNA analysis in different population and disease contexts, which is not available in other software.
Collapse
Affiliation(s)
- Larry N. Singh
- Center for Mitochondrial and Epigenomic Medicine, Division of Human Genetics, The Children’s Hospital of Philadelphia, Philadelphia, Pennsylvania, United States of America
- * E-mail:
| | - Brian Ennis
- Center for Data-Driven Discovery in Biomedicine (D3b), The Children’s Hospital of Philadelphia, Philadelphia, Pennsylvania, United States of America
| | - Bryn Loneragan
- Center for Eye Research Australia, Ophthalmology, Department of Surgery, University of Melbourne, Melbourne, Australia
| | - Noah L. Tsao
- Department of Surgery, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| | - M. Isabel G. Lopez Sanchez
- Center for Eye Research Australia, Ophthalmology, Department of Surgery, University of Melbourne, Melbourne, Australia
| | - Jianping Li
- Department of Psychiatry, The Children’s Hospital of Philadelphia and the University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| | - Patrick Acheampong
- Center for Mitochondrial and Epigenomic Medicine, Division of Human Genetics, The Children’s Hospital of Philadelphia, Philadelphia, Pennsylvania, United States of America
| | - Oanh Tran
- 22q and You Center, Division of Human Genetics, The Children’s Hospital of Philadelphia and the University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| | - Ian A. Trounce
- Center for Eye Research Australia, Ophthalmology, Department of Surgery, University of Melbourne, Melbourne, Australia
| | - Yuankun Zhu
- Center for Data-Driven Discovery in Biomedicine (D3b), The Children’s Hospital of Philadelphia, Philadelphia, Pennsylvania, United States of America
| | - Prasanth Potluri
- Center for Mitochondrial and Epigenomic Medicine, Division of Human Genetics, The Children’s Hospital of Philadelphia, Philadelphia, Pennsylvania, United States of America
| | | | - Beverly S. Emanuel
- 22q and You Center, Division of Human Genetics, The Children’s Hospital of Philadelphia and the University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| | - Daniel J. Rader
- Department of Genetics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| | - Zoltan Arany
- Cardiovascular Institute, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| | - Scott M. Damrauer
- Department of Surgery, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| | - Adam C. Resnick
- Center for Data-Driven Discovery in Biomedicine (D3b), The Children’s Hospital of Philadelphia, Philadelphia, Pennsylvania, United States of America
| | - Stewart A. Anderson
- Department of Psychiatry, The Children’s Hospital of Philadelphia and the University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| | - Douglas C. Wallace
- Center for Mitochondrial and Epigenomic Medicine, Division of Human Genetics, The Children’s Hospital of Philadelphia, Philadelphia, Pennsylvania, United States of America
| |
Collapse
|
42
|
Lüth T, Wasner K, Klein C, Schaake S, Tse R, Pereira SL, Laß J, Sinkkonen L, Grünewald A, Trinh J. Nanopore Single-Molecule Sequencing for Mitochondrial DNA Methylation Analysis: Investigating Parkin-Associated Parkinsonism as a Proof of Concept. Front Aging Neurosci 2021; 13:713084. [PMID: 34650424 PMCID: PMC8506010 DOI: 10.3389/fnagi.2021.713084] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2021] [Accepted: 09/01/2021] [Indexed: 01/08/2023] Open
Abstract
Objective: To establish a workflow for mitochondrial DNA (mtDNA) CpG methylation using Nanopore whole-genome sequencing and perform first pilot experiments on affected Parkin biallelic mutation carriers (Parkin-PD) and healthy controls. Background: Mitochondria, including mtDNA, are established key players in Parkinson's disease (PD) pathogenesis. Mutations in Parkin, essential for degradation of damaged mitochondria, cause early-onset PD. However, mtDNA methylation and its implication in PD is understudied. Herein, we establish a workflow using Nanopore sequencing to directly detect mtDNA CpG methylation and compare mtDNA methylation between Parkin-related PD and healthy individuals. Methods: To obtain mtDNA, whole-genome Nanopore sequencing was performed on blood-derived from five Parkin-PD and three control subjects. In addition, induced pluripotent stem cell (iPSC)-derived midbrain neurons from four of these patients with PD and the three control subjects were investigated. The workflow was validated, using methylated and unmethylated 897 bp synthetic DNA samples at different dilution ratios (0, 50, 100% methylation) and mtDNA without methylation. MtDNA CpG methylation frequency (MF) was detected using Nanopolish and Megalodon. Results: Across all blood-derived samples, we obtained a mean coverage of 250.3X (SD ± 80.5X) and across all neuron-derived samples 830X (SD ± 465X) of the mitochondrial genome. We detected overall low-level CpG methylation from the blood-derived DNA (mean MF ± SD = 0.029 ± 0.041) and neuron-derived DNA (mean MF ± SD = 0.019 ± 0.035). Validation of the workflow, using synthetic DNA samples showed that highly methylated DNA molecules were prone to lower Guppy Phred quality scores and thereby more likely to fail Guppy base-calling. CpG methylation in blood- and neuron-derived DNA was significantly lower in Parkin-PD compared to controls (Mann-Whitney U-test p < 0.05). Conclusion: Nanopore sequencing is a useful method to investigate mtDNA methylation architecture, including Guppy-failed reads is of importance when investigating highly methylated sites. We present a mtDNA methylation workflow and suggest methylation variability across different tissues and between Parkin-PD patients and controls as an initial model to investigate.
Collapse
Affiliation(s)
- Theresa Lüth
- Institute of Neurogenetics BMF, University of Lübeck and University Hospital Schleswig-Holstein Campus Lübeck, Lübeck, Germany
| | - Kobi Wasner
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Christine Klein
- Institute of Neurogenetics BMF, University of Lübeck and University Hospital Schleswig-Holstein Campus Lübeck, Lübeck, Germany
| | - Susen Schaake
- Institute of Neurogenetics BMF, University of Lübeck and University Hospital Schleswig-Holstein Campus Lübeck, Lübeck, Germany
| | - Ronnie Tse
- Institute of Neurogenetics BMF, University of Lübeck and University Hospital Schleswig-Holstein Campus Lübeck, Lübeck, Germany
| | - Sandro L. Pereira
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Joshua Laß
- Institute of Neurogenetics BMF, University of Lübeck and University Hospital Schleswig-Holstein Campus Lübeck, Lübeck, Germany
| | - Lasse Sinkkonen
- Department of Life Sciences and Medicine, University of Luxembourg, Belvaux, Luxembourg
| | - Anne Grünewald
- Institute of Neurogenetics BMF, University of Lübeck and University Hospital Schleswig-Holstein Campus Lübeck, Lübeck, Germany
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Joanne Trinh
- Institute of Neurogenetics BMF, University of Lübeck and University Hospital Schleswig-Holstein Campus Lübeck, Lübeck, Germany
| |
Collapse
|
43
|
Zhang G, Geng D, Guo Q, Liu W, Li S, Gao W, Wang Y, Zhang M, Wang Y, Bu Y, Niu H. Genomic landscape of mitochondrial DNA insertions in 23 bat genomes: characteristics, loci, phylogeny, and polymorphism. Integr Zool 2021; 17:890-903. [PMID: 34496458 DOI: 10.1111/1749-4877.12582] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The transfer of mitochondrial DNA to the nuclear genome gives rise to the nuclear DNA sequences of mitochondrial origin (NUMTs), considered as a driving force in genome evolution. In this study, NUMTs in 23 bat genomes were investigated and compared systematically. The results showed that NUMTs existed in 22 genomes except for Noctilio leporinus, suggesting that mitochondrial fragment insertion in the nuclear genome was a common event in bat genomes. However, remarkable variations in NUMTs number, cumulative length, and proportion in the nuclear genome were discovered across bat species. Further orthologous NUMT loci analysis of the Phyllostomidae family indicated that the NUMTs insertion events in bat genomes were homoplasy-free. The NUMTs were mainly inserted into the intergenic regions, particularly, co-localized with repetitive sequences (especially transposable elements). However, several NUMTs were inserted into genes, some of which were in the exon region of functional genes. One NUMT in the genome of Pteropus alecto surprisingly matched with cDNA of ATP8B3 that provided evidence of NUMTs with coding function. Phylogenic analysis on NUMTs originating from COXI and COXII loci highlighted 2 clusters of Yinpterochiroptera and Yangochiroptera for Chiroptera. Seven NUMTs from Rhinolophus ferrumequinum were amplified, and the sequencing results confirmed the reliability of the NUMT analysis. One of them was polymorphic for the presence or absence of the NUMT insertion, and each genotype of NUMT loci showed a distinct regional distribution pattern. The information obtained in this study provides novel insights into the NUMT organization and features in bat genomes and establishes a basis for further studying of the evolution of bat species.
Collapse
Affiliation(s)
- Guojun Zhang
- College of Life Sciences, Henan Normal University, Xinxiang, China.,School of Basic Medical Sciences, Xinxiang Medical University, Xinxiang, China
| | - Deqi Geng
- College of Life Sciences, Henan Normal University, Xinxiang, China
| | - Qiulin Guo
- College of Life Sciences, Henan Normal University, Xinxiang, China
| | - Wei Liu
- College of Life Sciences, Henan Normal University, Xinxiang, China
| | - Shufen Li
- College of Life Sciences, Henan Normal University, Xinxiang, China
| | - Wujun Gao
- College of Life Sciences, Henan Normal University, Xinxiang, China
| | - Yongfei Wang
- College of Life Sciences, Henan Normal University, Xinxiang, China
| | - Min Zhang
- College of Life Sciences, Henan Normal University, Xinxiang, China
| | - Yilin Wang
- College of Life Sciences, Henan Normal University, Xinxiang, China
| | - Yanzhen Bu
- College of Life Sciences, Henan Normal University, Xinxiang, China
| | - Hongxing Niu
- College of Life Sciences, Henan Normal University, Xinxiang, China
| |
Collapse
|
44
|
The Isolation and Deep Sequencing of Mitochondrial DNA. Methods Mol Biol 2021. [PMID: 34080167 DOI: 10.1007/978-1-0716-1270-5_27] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
Abstract
In recent years, next-generation sequencing (NGS) has become a powerful tool for studying both inherited and somatic heteroplasmic mitochondrial DNA (mtDNA) variation. NGS has proved particularly powerful when combined with single-cell isolation techniques, allowing the investigation of low-level heteroplasmic variants both between cells and within tissues. Nevertheless, there remain significant challenges, especially around the selective enrichment of mtDNA from total cellular DNA and the avoidance of nuclear pseudogenes. This chapter summarizes the techniques needed to enrich, amplify, sequence, and analyse mtDNA using NGS .
Collapse
|
45
|
Abstract
Variation in the mitochondrial DNA (mtDNA) sequence is common in certain tumours. Two classes of cancer mtDNA variants can be identified: de novo mutations that act as 'inducers' of carcinogenesis and functional variants that act as 'adaptors', permitting cancer cells to thrive in different environments. These mtDNA variants have three origins: inherited variants, which run in families, somatic mutations arising within each cell or individual, and variants that are also associated with ancient mtDNA lineages (haplogroups) and are thought to permit adaptation to changing tissue or geographic environments. In addition to mtDNA sequence variation, mtDNA copy number and perhaps transfer of mtDNA sequences into the nucleus can contribute to certain cancers. Strong functional relevance of mtDNA variation has been demonstrated in oncocytoma and prostate cancer, while mtDNA variation has been reported in multiple other cancer types. Alterations in nuclear DNA-encoded mitochondrial genes have confirmed the importance of mitochondrial metabolism in cancer, affecting mitochondrial reactive oxygen species production, redox state and mitochondrial intermediates that act as substrates for chromatin-modifying enzymes. Hence, subtle changes in the mitochondrial genotype can have profound effects on the nucleus, as well as carcinogenesis and cancer progression.
Collapse
Affiliation(s)
- Piotr K Kopinski
- Howard Hughes Medical Institute, University of Pennsylvania, Philadelphia, PA, USA
- Center for Mitochondrial and Epigenomic Medicine, The Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - Larry N Singh
- Center for Mitochondrial and Epigenomic Medicine, The Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - Shiping Zhang
- Center for Mitochondrial and Epigenomic Medicine, The Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - Marie T Lott
- Center for Mitochondrial and Epigenomic Medicine, The Children's Hospital of Philadelphia, Philadelphia, PA, USA
| | - Douglas C Wallace
- Center for Mitochondrial and Epigenomic Medicine, The Children's Hospital of Philadelphia, Philadelphia, PA, USA.
- Department of Pediatrics, Division of Human Genetics, The Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA.
| |
Collapse
|
46
|
Rossmann MP, Dubois SM, Agarwal S, Zon LI. Mitochondrial function in development and disease. Dis Model Mech 2021; 14:269120. [PMID: 34114603 PMCID: PMC8214736 DOI: 10.1242/dmm.048912] [Citation(s) in RCA: 43] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Mitochondria are organelles with vital functions in almost all eukaryotic cells. Often described as the cellular ‘powerhouses’ due to their essential role in aerobic oxidative phosphorylation, mitochondria perform many other essential functions beyond energy production. As signaling organelles, mitochondria communicate with the nucleus and other organelles to help maintain cellular homeostasis, allow cellular adaptation to diverse stresses, and help steer cell fate decisions during development. Mitochondria have taken center stage in the research of normal and pathological processes, including normal tissue homeostasis and metabolism, neurodegeneration, immunity and infectious diseases. The central role that mitochondria assume within cells is evidenced by the broad impact of mitochondrial diseases, caused by defects in either mitochondrial or nuclear genes encoding for mitochondrial proteins, on different organ systems. In this Review, we will provide the reader with a foundation of the mitochondrial ‘hardware’, the mitochondrion itself, with its specific dynamics, quality control mechanisms and cross-organelle communication, including its roles as a driver of an innate immune response, all with a focus on development, disease and aging. We will further discuss how mitochondrial DNA is inherited, how its mutation affects cell and organismal fitness, and current therapeutic approaches for mitochondrial diseases in both model organisms and humans. Summary: Mitochondria have a plethora of functions beyond metabolism. This Review discusses the emerging and multifaceted roles of mitochondria in different model organisms and human disease biology.
Collapse
Affiliation(s)
- Marlies P Rossmann
- Department of Stem Cell and Regenerative Biology, Harvard University, Cambridge, MA 01238, USA.,Stem Cell Program and Division of Hematology/Oncology, Boston Children's Hospital, Harvard Medical School, Howard Hughes Medical Institute, Boston, MA 02115, USA
| | - Sonia M Dubois
- Division of Hematology/Oncology, Boston Children's Hospital, Harvard Medical School, Boston, MA 02115, USA
| | - Suneet Agarwal
- Division of Hematology/Oncology, Boston Children's Hospital, Harvard Medical School, Boston, MA 02115, USA
| | - Leonard I Zon
- Department of Stem Cell and Regenerative Biology, Harvard University, Cambridge, MA 01238, USA.,Stem Cell Program and Division of Hematology/Oncology, Boston Children's Hospital, Harvard Medical School, Howard Hughes Medical Institute, Boston, MA 02115, USA
| |
Collapse
|
47
|
Hofreiter M, Sneberger J, Pospisek M, Vanek D. Progress in forensic bone DNA analysis: Lessons learned from ancient DNA. Forensic Sci Int Genet 2021; 54:102538. [PMID: 34265517 DOI: 10.1016/j.fsigen.2021.102538] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2020] [Revised: 03/07/2021] [Accepted: 05/25/2021] [Indexed: 01/18/2023]
Abstract
Research on ancient and forensic DNA is related in many ways, and the two fields must deal with similar obstacles. Therefore, communication between these two communities has the potential to improve results in both research fields. Here, we present the insights gained in the ancient DNA community with regard to analyzing DNA from aged skeletal material and the potential use of the developed protocols in forensic work. We discuss the various steps, from choosing samples for DNA extraction to deciding between classical PCR amplification and massively parallel sequencing approaches. Based on the progress made in ancient DNA analyses combined with the requirements of forensic work, we suggest that there is substantial potential for incorporating ancient DNA approaches into forensic protocols, a process that has already begun to a considerable extent. However, taking full advantage of the experiences gained from ancient DNA work will require comparative studies by the forensic DNA community to tailor the methods developed for ancient samples to the specific needs of forensic studies and case work. If successful, in our view, the benefits for both communities would be considerable.
Collapse
Affiliation(s)
- Michael Hofreiter
- Institute for Biochemistry and Biology, University of Potsdam, Karl-Liebknecht-Str. 24-25, 14476 Potsdam, Germany.
| | - Jiri Sneberger
- Department of Genetics and Microbiology, Faculty of Science, Charles University in Prague, Vinicna 5, Prague 2 12843, Czech Republic; Department of the History of the Middle Ages of Museum of West Bohemia, Kopeckeho sady 2, Pilsen 30100, Czech Republic; Nuclear Physics Institute of the CAS, Na Truhlarce 39/64, Prague 18086, Czech Republic
| | - Martin Pospisek
- Department of Genetics and Microbiology, Faculty of Science, Charles University in Prague, Vinicna 5, Prague 2 12843, Czech Republic; Biologicals s.r.o., Sramkova 315, Ricany 25101, Czech Republic
| | - Daniel Vanek
- Forensic DNA Service, Janovskeho 18, Prague 7 17000, Czech Republic; Institute of Legal Medicine, Bulovka Hospital, Prague, Czech Republic; Charles University in Prague, 2nd Faculty of Medicine, Prague, Czech Republic.
| |
Collapse
|
48
|
Formenti G, Rhie A, Balacco J, Haase B, Mountcastle J, Fedrigo O, Brown S, Capodiferro MR, Al-Ajli FO, Ambrosini R, Houde P, Koren S, Oliver K, Smith M, Skelton J, Betteridge E, Dolucan J, Corton C, Bista I, Torrance J, Tracey A, Wood J, Uliano-Silva M, Howe K, McCarthy S, Winkler S, Kwak W, Korlach J, Fungtammasan A, Fordham D, Costa V, Mayes S, Chiara M, Horner DS, Myers E, Durbin R, Achilli A, Braun EL, Phillippy AM, Jarvis ED. Complete vertebrate mitogenomes reveal widespread repeats and gene duplications. Genome Biol 2021; 22:120. [PMID: 33910595 PMCID: PMC8082918 DOI: 10.1186/s13059-021-02336-9] [Citation(s) in RCA: 50] [Impact Index Per Article: 16.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2020] [Accepted: 03/31/2021] [Indexed: 01/22/2023] Open
Abstract
BACKGROUND Modern sequencing technologies should make the assembly of the relatively small mitochondrial genomes an easy undertaking. However, few tools exist that address mitochondrial assembly directly. RESULTS As part of the Vertebrate Genomes Project (VGP) we develop mitoVGP, a fully automated pipeline for similarity-based identification of mitochondrial reads and de novo assembly of mitochondrial genomes that incorporates both long (> 10 kbp, PacBio or Nanopore) and short (100-300 bp, Illumina) reads. Our pipeline leads to successful complete mitogenome assemblies of 100 vertebrate species of the VGP. We observe that tissue type and library size selection have considerable impact on mitogenome sequencing and assembly. Comparing our assemblies to purportedly complete reference mitogenomes based on short-read sequencing, we identify errors, missing sequences, and incomplete genes in those references, particularly in repetitive regions. Our assemblies also identify novel gene region duplications. The presence of repeats and duplications in over half of the species herein assembled indicates that their occurrence is a principle of mitochondrial structure rather than an exception, shedding new light on mitochondrial genome evolution and organization. CONCLUSIONS Our results indicate that even in the "simple" case of vertebrate mitogenomes the completeness of many currently available reference sequences can be further improved, and caution should be exercised before claiming the complete assembly of a mitogenome, particularly from short reads alone.
Collapse
Affiliation(s)
- Giulio Formenti
- The Vertebrate Genome Lab, Rockefeller University, New York, NY, USA.
- Laboratory of Neurogenetics of Language, Rockefeller University, New York, NY, USA.
- The Howards Hughes Medical Institute, Chevy Chase, MD, USA.
| | - Arang Rhie
- Genome Informatics Section, Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Jennifer Balacco
- The Vertebrate Genome Lab, Rockefeller University, New York, NY, USA
| | - Bettina Haase
- The Vertebrate Genome Lab, Rockefeller University, New York, NY, USA
| | | | - Olivier Fedrigo
- The Vertebrate Genome Lab, Rockefeller University, New York, NY, USA
| | - Samara Brown
- Laboratory of Neurogenetics of Language, Rockefeller University, New York, NY, USA
- The Howards Hughes Medical Institute, Chevy Chase, MD, USA
| | | | - Farooq O Al-Ajli
- Monash University Malaysia Genomics Facility, School of Science, Bandar Sunway, Selangor Darul Ehsan, Malaysia
- Tropical Medicine and Biology Multidisciplinary Platform, Monash University Malaysia, Bandar Sunway, Selangor Darul Ehsan, Malaysia
- Qatar Falcon Genome Project, Doha, State of Qatar
| | - Roberto Ambrosini
- Department of Environmental Science and Policy, University of Milan, Milan, Italy
| | - Peter Houde
- Department of Biology, New Mexico State University, Las Cruces, NM, USA
| | - Sergey Koren
- Genome Informatics Section, Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | | | | | | | | | | | | | - Iliana Bista
- Wellcome Sanger Institute, Cambridge, UK
- Department of Genetics, University of Cambridge, Cambridge, UK
| | | | | | | | | | | | - Shane McCarthy
- Wellcome Sanger Institute, Cambridge, UK
- Department of Genetics, University of Cambridge, Cambridge, UK
| | - Sylke Winkler
- Max Planck Institute of Molecular Cell Biology & Genetics, Dresden, Germany
| | | | | | | | - Daniel Fordham
- Oxford Nanopore Technologies Ltd, Oxford Science Park, Oxford, UK
| | - Vania Costa
- Oxford Nanopore Technologies Ltd, Oxford Science Park, Oxford, UK
| | - Simon Mayes
- Oxford Nanopore Technologies Ltd, Oxford Science Park, Oxford, UK
| | - Matteo Chiara
- Department of Biosciences, University of Milan, Milan, Italy
| | - David S Horner
- Department of Biosciences, University of Milan, Milan, Italy
| | - Eugene Myers
- Max Planck Institute of Molecular Cell Biology & Genetics, Dresden, Germany
| | - Richard Durbin
- Wellcome Sanger Institute, Cambridge, UK
- Department of Genetics, University of Cambridge, Cambridge, UK
| | - Alessandro Achilli
- Department of Biology and Biotechnology "L. Spallanzani", University of Pavia, Pavia, Italy
| | - Edward L Braun
- Department of Biology, University of Florida, Gainesville, FL, USA
| | - Adam M Phillippy
- Genome Informatics Section, Computational and Statistical Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA
| | - Erich D Jarvis
- The Vertebrate Genome Lab, Rockefeller University, New York, NY, USA.
- Laboratory of Neurogenetics of Language, Rockefeller University, New York, NY, USA.
- The Howards Hughes Medical Institute, Chevy Chase, MD, USA.
| |
Collapse
|
49
|
Akimova E, Gassner FJ, Schubert M, Rebhandl S, Arzt C, Rauscher S, Tober V, Zaborsky N, Greil R, Geisberger R. SAMHD1 restrains aberrant nucleotide insertions at repair junctions generated by DNA end joining. Nucleic Acids Res 2021; 49:2598-2608. [PMID: 33591315 PMCID: PMC7969033 DOI: 10.1093/nar/gkab051] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2020] [Revised: 01/15/2021] [Accepted: 01/19/2021] [Indexed: 12/15/2022] Open
Abstract
Aberrant end joining of DNA double strand breaks leads to chromosomal rearrangements and to insertion of nuclear or mitochondrial DNA into breakpoints, which is commonly observed in cancer cells and constitutes a major threat to genome integrity. However, the mechanisms that are causative for these insertions are largely unknown. By monitoring end joining of different linear DNA substrates introduced into HEK293 cells, as well as by examining end joining of CRISPR/Cas9 induced DNA breaks in HEK293 and HeLa cells, we provide evidence that the dNTPase activity of SAMHD1 impedes aberrant DNA resynthesis at DNA breaks during DNA end joining. Hence, SAMHD1 expression or low intracellular dNTP levels lead to shorter repair joints and impede insertion of distant DNA regions prior end repair. Our results reveal a novel role for SAMHD1 in DNA end joining and provide new insights into how loss of SAMHD1 may contribute to genome instability and cancer development.
Collapse
Affiliation(s)
- Ekaterina Akimova
- Department of Internal Medicine III with Haematology, Medical Oncology, Haemostaseology, Infectiology and Rheumatology, Oncologic Center, Paracelsus Medical University, Salzburg, Austria.,Salzburg Cancer Research Institute - Laboratory for Immunological and Molecular Cancer Research (SCRI-LIMCR); Cancer Cluster Salzburg, 5020 Salzburg, Austria.,Department of Biosciences, Paris Lodron University of Salzburg, 5020 Salzburg, Austria
| | - Franz Josef Gassner
- Department of Internal Medicine III with Haematology, Medical Oncology, Haemostaseology, Infectiology and Rheumatology, Oncologic Center, Paracelsus Medical University, Salzburg, Austria.,Salzburg Cancer Research Institute - Laboratory for Immunological and Molecular Cancer Research (SCRI-LIMCR); Cancer Cluster Salzburg, 5020 Salzburg, Austria
| | - Maria Schubert
- Department of Internal Medicine III with Haematology, Medical Oncology, Haemostaseology, Infectiology and Rheumatology, Oncologic Center, Paracelsus Medical University, Salzburg, Austria.,Salzburg Cancer Research Institute - Laboratory for Immunological and Molecular Cancer Research (SCRI-LIMCR); Cancer Cluster Salzburg, 5020 Salzburg, Austria
| | - Stefan Rebhandl
- Department of Internal Medicine III with Haematology, Medical Oncology, Haemostaseology, Infectiology and Rheumatology, Oncologic Center, Paracelsus Medical University, Salzburg, Austria.,Salzburg Cancer Research Institute - Laboratory for Immunological and Molecular Cancer Research (SCRI-LIMCR); Cancer Cluster Salzburg, 5020 Salzburg, Austria
| | - Claudia Arzt
- Department of Internal Medicine III with Haematology, Medical Oncology, Haemostaseology, Infectiology and Rheumatology, Oncologic Center, Paracelsus Medical University, Salzburg, Austria.,Salzburg Cancer Research Institute - Laboratory for Immunological and Molecular Cancer Research (SCRI-LIMCR); Cancer Cluster Salzburg, 5020 Salzburg, Austria
| | - Stefanie Rauscher
- Department of Internal Medicine III with Haematology, Medical Oncology, Haemostaseology, Infectiology and Rheumatology, Oncologic Center, Paracelsus Medical University, Salzburg, Austria.,Salzburg Cancer Research Institute - Laboratory for Immunological and Molecular Cancer Research (SCRI-LIMCR); Cancer Cluster Salzburg, 5020 Salzburg, Austria.,Department of Biosciences, Paris Lodron University of Salzburg, 5020 Salzburg, Austria
| | - Vanessa Tober
- Department of Internal Medicine III with Haematology, Medical Oncology, Haemostaseology, Infectiology and Rheumatology, Oncologic Center, Paracelsus Medical University, Salzburg, Austria.,Salzburg Cancer Research Institute - Laboratory for Immunological and Molecular Cancer Research (SCRI-LIMCR); Cancer Cluster Salzburg, 5020 Salzburg, Austria.,Department of Biosciences, Paris Lodron University of Salzburg, 5020 Salzburg, Austria
| | - Nadja Zaborsky
- Department of Internal Medicine III with Haematology, Medical Oncology, Haemostaseology, Infectiology and Rheumatology, Oncologic Center, Paracelsus Medical University, Salzburg, Austria.,Salzburg Cancer Research Institute - Laboratory for Immunological and Molecular Cancer Research (SCRI-LIMCR); Cancer Cluster Salzburg, 5020 Salzburg, Austria
| | - Richard Greil
- Department of Internal Medicine III with Haematology, Medical Oncology, Haemostaseology, Infectiology and Rheumatology, Oncologic Center, Paracelsus Medical University, Salzburg, Austria.,Salzburg Cancer Research Institute - Laboratory for Immunological and Molecular Cancer Research (SCRI-LIMCR); Cancer Cluster Salzburg, 5020 Salzburg, Austria
| | - Roland Geisberger
- Department of Internal Medicine III with Haematology, Medical Oncology, Haemostaseology, Infectiology and Rheumatology, Oncologic Center, Paracelsus Medical University, Salzburg, Austria.,Salzburg Cancer Research Institute - Laboratory for Immunological and Molecular Cancer Research (SCRI-LIMCR); Cancer Cluster Salzburg, 5020 Salzburg, Austria
| |
Collapse
|
50
|
Leroux É, Brosseau C, Angers B, Angers A, Breton S. [Mitochondrial DNA methylation: Controversies, issues and perspectives]. Med Sci (Paris) 2021; 37:258-264. [PMID: 33739273 DOI: 10.1051/medsci/2021011] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Abstract
DNA methylation is an epigenetic mechanism that has been largely probed regarding eukaryotic nuclear genome and bacteria, and its role is especially crucial in the regulation of gene expression. In mammals, it is almost exclusively acting on a cytosine preceding a guanine (CpG), whereas it presents itself mainly in a non-CpG context in bacteria's DNA. Conversely to nuclear and bacterial genomes, the existence of methylation in the mitochondrial genome is still widely debated. This controversy has been attributed to structural differences between the nuclear and mitochondrial genomes, and to the techniques used to study methylation of cytosines, which were rather optimized for the study of nuclear DNA. However, novel studies suggest that cytosine methylation is truly existing in mitochondria, and that it is mostly found in a non-CpG context, just like in their evolutionary relative, the bacteria.
Collapse
Affiliation(s)
- Émélie Leroux
- Département de sciences biologiques, Université de Montréal, Campus MIL, Faculté des Arts et des Sciences, CP 6128, Succursale Centre-Ville, Montréal QC, H3C 3J7, Canada
| | - Cindy Brosseau
- Département de sciences biologiques, Université de Montréal, Campus MIL, Faculté des Arts et des Sciences, CP 6128, Succursale Centre-Ville, Montréal QC, H3C 3J7, Canada
| | - Bernard Angers
- Département de sciences biologiques, Université de Montréal, Campus MIL, Faculté des Arts et des Sciences, CP 6128, Succursale Centre-Ville, Montréal QC, H3C 3J7, Canada
| | - Annie Angers
- Département de sciences biologiques, Université de Montréal, Campus MIL, Faculté des Arts et des Sciences, CP 6128, Succursale Centre-Ville, Montréal QC, H3C 3J7, Canada
| | - Sophie Breton
- Département de sciences biologiques, Université de Montréal, Campus MIL, Faculté des Arts et des Sciences, CP 6128, Succursale Centre-Ville, Montréal QC, H3C 3J7, Canada
| |
Collapse
|