Comparison of whole genome sequencing and targeted sequencing for mitochondrial DNA

Genetic landscape of primary mitochondrial diseases in children and adults using molecular genetics and genomic investigations of mitochondrial and nuclear genome

BACKGROUND

Primary mitochondrial diseases (PMD) are one of the most common metabolic genetic disorders. They are due to pathogenic variants in the mitochondrial genome (mtDNA) or nuclear genome (nDNA) that impair mitochondrial function and/or structure. We hypothesize that there is overlap between PMD and other genetic diseases that are mimicking PMD. For this reason, we performed a retrospective cohort study.

METHODS

All individuals with suspected PMD that underwent molecular genetic and genomic investigations were included. Individuals were grouped for comparison: (1) individuals with mtDNA-PMD; (2) individuals with nDNA-PMD; (3) individuals with other genetic diseases mimicking PMD (non-PMD); (4) individuals without a confirmed genetic diagnosis.

RESULTS

297 individuals fulfilled inclusion criteria. The diagnostic yield of molecular genetics and genomic investigations was 31.3%, including 37% for clinical exome sequencing and 15.8% for mitochondrial genome sequencing. We identified 71 individuals with PMD (mtDNA n = 41, nDNA n = 30) and 22 individuals with non-PMD. Adults had higher percentage of mtDNA-PMD compared to children (p-value = 0.00123). There is a statistically significant phenotypic difference between children and adults with PMD.

CONCLUSION

We report a large cohort of individuals with PMD and the diagnostic yield of urine mitochondrial genome sequencing (16.1%). We think liver phenotype might be progressive and should be studied further in PMD. We showed a relationship between non-PMD genes and their indirect effects on mitochondrial machinery. Differentiation of PMD from non-PMD can be achieved using specific phenotypes as there was a statistically significant difference for muscular, cardiac, and ophthalmologic phenotypes, seizures, hearing loss, peripheral neuropathy in PMD group compared to non-PMD group.

Genetic diversity and origin of Kazakh Tobet Dogs

The Kazakh Tobet is an indigenous Kazakh dog breed that has been used to guard livestock since ancient times. To understand the genetic structure and phylogenetic relationship of the Kazakh Tobet breed with other herding and livestock guarding dog breeds, we analysed short tandem repeat data of 107 Kazakh Tobet dogs from different regions of Kazakhstan and Mongolia, as well as whole genome sequencing data from two Kazakh Tobet dogs and 43 dogs from 24 working breeds. Our results indicate a high genetic diversity of the Kazakh Tobet, with the average number of alleles per locus ranging from 6.00 to 10.22 and observed heterozygosity ranging from 76 to 78%. The breed has a complex genetic structure characterised by seven different clusters. The neighbour-joining tree constructed based on 14,668,406 autosomal and the maximum likelihood tree based on mitochondrial D-loop sequences indicate a common genetic heritage between the Kazakh Tobet, the Central Asian Shepherd Dog and the Turkish Akbash. The presence of haplotype A18 in the Kazakh Tobets supports the hypothesis of the ancient origin of the breed, which was previously suggested by archaeological finds and written sources. These results provide an important genetic basis for the ongoing efforts to improve the Kazakh Tobet breed, to ensure its preservation as an independent genetic lineage and to recognise a breed on an international level.

Two complete mitochondrial genomes of Boulenophrys (Anura: Megophryidae: Megophryinae): characteristics and phylogenetic implications

The Chinese horned toads, Boulenophrys boettgeri (Boulenger, 1899) and Boulenophrys kuatunensis (Pope, 1929), are two captivating species within the family Megophryidae, which inhabit the mountainous streams in the Eastern of China. In this study, two new complete mitochondrial genomes of B. boettgeri and B. kuatunensis were sequenced, assembled, and annotated using next-generation sequencing. The length of mitochondrial genomes of B. boettgeri and B. kuatunensis was 16,597 and 17,921 bp, respectively, with both containing 13 protein coding genes, 22 tRNA genes, two rRNA genes, and one putative control region. Phylogenetic relationships based on protein-coding mitochondrial genes showed that the two Boulenophrys species formed a cluster with other Boulenophrys species. The two new sequences provide valuable insights into the mitochondrial genomes of these two species, offering important data for understanding the phylogenetic relationships of the genus Boulenophrys.

A PCR-independent approach for mtDNA enrichment and next-generation sequencing: comprehensive evaluation and clinical application

BACKGROUND

Sequencing the mitochondrial genome has been increasingly important for the investigation of primary mitochondrial diseases (PMD) and mitochondrial genetics. To overcome the limitations originating from PCR-based mtDNA enrichment, we set out to develop and evaluate a PCR-independent approach in this study, named Pime-Seq (PCR-independent mtDNA enrichment and next generation Sequencing).

RESULTS

By using the optimized mtDNA enrichment procedure, the mtDNA reads ratio reached 88.0 ± 7.9% in the sequencing library when applied on human PBMC samples. We found the variants called by Pime-Seq were highly consistent among technical repeats. To evaluate the accuracy and reliability of this method, we compared Pime-Seq with lrPCR based NGS by performing both methods simultaneously on 45 samples, yielding 1677 concordant variants, as well as 146 discordant variants with low-level heteroplasmic fraction, in which Pime-Seq showed higher reliability. Furthermore, we applied Pime-Seq on 4 samples of PMD patients retrospectively, and successfully detected all the pathogenic mtDNA variants. In addition, we performed a prospective study on 192 apparently healthy pregnant women during prenatal screening, in which Pime-Seq identified pathogenic mtDNA variants in 4 samples, providing extra information for better health monitoring in these cases.

CONCLUSIONS

Pime-Seq can obtain highly enriched mtDNA in a PCR-independent manner for high quality and reliable mtDNA deep-sequencing, which provides us an effective and promising tool for detecting mtDNA variants for both clinical and research purposes.

Applications of genome sequencing as a single platform for clinical constitutional genetic testing

The number of human disease genes has dramatically increased over the past decade, largely fueled by ongoing advances in sequencing technologies. In parallel, the number of available clinical genetic tests has also increased, including the utilization of exome sequencing for undiagnosed diseases. Although most clinical sequencing tests have been centered on enrichment-based multigene panels and exome sequencing, the continued improvements in performance and throughput of genome sequencing suggest that this technology is emerging as a potential platform for routine clinical genetic testing. A notable advantage is a single workflow with the opportunity to reflexively interrogate content as clinically indicated; however, challenges with implementing routine clinical genome sequencing still remain. This review is centered on evaluating the applications of genome sequencing as a single platform for clinical constitutional genetic testing, including its potential utility for diagnostic testing, carrier screening, cytogenomic molecular karyotyping, prenatal testing, mitochondrial genome interrogation, and pharmacogenomic and polygenic risk score testing.

A Practical Guide to Whole Genome Sequencing in the NICU

The neonatal period is a peak time for the presentation of genetic disorders that can be diagnosed using whole genome sequencing (WGS). While any one genetic disorder is individually rare, they collectively contribute to significant morbidity, mortality, and health-care costs. As the cost of WGS continues to decline and becomes increasingly available, the ordering of rapid WGS for NICU patients with signs or symptoms of an underlying genetic condition is now feasible. However, many neonatal clinicians are not comfortable with the testing, and unfortunately, there is a dearth of geneticists to facilitate testing for every patient that needs it. Here, we will review the science behind WGS, diagnostic capabilities, limitations of testing, time to consider testing, test initiation, interpretation of results, developing a plan of care that incorporates genomic information, and returning WGS results to families.

Revisiting mitogenome evolution in Medusozoa with eight new mitochondrial genomes

Mitogenomics has improved our understanding of medusozoan phylogeny. However, sequenced medusozoan mitogenomes remain scarce, and Medusozoa phylogeny studies often analyze mitogenomic sequences without incorporating mitogenome rearrangements. To better understand medusozoan evolution, we analyzed Medusozoa mitogenome phylogeny by sequencing and assembling eight mitogenomes from three classes (Cubozoa, Hydrozoa, and Scyphozoa). We reconstructed the mitogenome phylogeny using these mitogenomes and 84 other existing cnidarian mitogenomes to study mitochondrial gene rearrangements. All reconstructed mitogenomes had 13 mitochondrial protein-coding genes and two ribosomal genes typical for Medusozoa. Non-cubozoan mitogenomes were all linear and had typical gene orders, while arrangement of genes in the fragmented Cubozoa (Morbakka sp.) mitogenome differed from other Cubozoa mitogenomes. Gene order comparisons and ancestral state reconstruction suggest minimal rearrangements within medusozoan classes except for Hydrozoa. Our findings support a staurozoan ancestral medusozoan gene order, expand the pool of available medusozoan mitogenomes, and enhance our understanding of medusozoan phylogenetic relationships.

Mitochondrial DNA in Human Diversity and Health: From the Golden Age to the Omics Era

Mitochondrial DNA (mtDNA) is a small fraction of our hereditary material. However, this molecule has had an overwhelming presence in scientific research for decades until the arrival of high-throughput studies. Several appealing properties justify the application of mtDNA to understand how human populations are-from a genetic perspective-and how individuals exhibit phenotypes of biomedical importance. Here, I review the basics of mitochondrial studies with a focus on the dawn of the field, analysis methods and the connection between two sides of mitochondrial genetics: anthropological and biomedical. The particularities of mtDNA, with respect to inheritance pattern, evolutionary rate and dependence on the nuclear genome, explain the challenges of associating mtDNA composition and diseases. Finally, I consider the relevance of this single locus in the context of omics research. The present work may serve as a tribute to a tool that has provided important insights into the past and present of humankind.

At-Risk Genomic Findings for Pediatric-Onset Disorders From Genome Sequencing vs Medically Actionable Gene Panel in Proactive Screening of Newborns and Children

Importance

Although the clinical utility of genome sequencing for critically ill children is well recognized, its utility for proactive pediatric screening is not well explored.

Objective

To evaluate molecular findings from screening ostensibly healthy children with genome sequencing compared with a gene panel for medically actionable pediatric conditions.

Design, Setting, and Participants

This case series study was conducted among consecutive, apparently healthy children undergoing proactive genetic screening for pediatric disorders by genome sequencing (n = 562) or an exome-based panel of 268 genes (n = 606) from March 1, 2018, through July 31, 2022.

Exposures

Genetic screening for pediatric-onset disorders using genome sequencing or an exome-based panel of 268 genes.

Main Outcomes and Measures

Molecular findings indicative of genetic disease risk.

Results

Of 562 apparently healthy children (286 girls [50.9%]; median age, 29 days [IQR, 9-117 days]) undergoing screening by genome sequencing, 46 (8.2%; 95% CI, 5.9%-10.5%) were found to be at risk for pediatric-onset disease, including 22 children (3.9%) at risk for high-penetrance disorders. Sequence analysis uncovered molecular diagnoses among 32 individuals (5.7%), while copy number variant analysis uncovered molecular diagnoses among 14 individuals (2.5%), including 4 individuals (0.7%) with chromosome scale abnormalities. Overall, there were 47 molecular diagnoses, with 1 individual receiving 2 diagnoses; of the 47 potential diagnoses, 22 (46.8%) were associated with high-penetrance conditions. Pathogenic variants in medically actionable pediatric genes were found in 6 individuals (1.1%), constituting 12.8% (6 of 47) of all diagnoses. At least 1 pharmacogenomic variant was reported for 89.0% (500 of 562) of the cohort. In contrast, of 606 children (293 girls [48.3%]; median age, 26 days [IQR, 10-67 days]) undergoing gene panel screening, only 13 (2.1%; 95% CI, 1.0%-3.3%) resulted in potential childhood-onset diagnoses, a significantly lower rate than those screened by genome sequencing (P < .001).

Conclusions and Relevance

In this case series study, genome sequencing as a proactive screening approach for children, due to its unrestrictive gene content and technical advantages in comparison with an exome-based gene panel for medically actionable childhood conditions, uncovered a wide range of heterogeneous high-penetrance pediatric conditions that could guide early interventions and medical management.

Targeting mitochondrial impairment for the treatment of cardiovascular diseases: From hypertension to ischemia-reperfusion injury, searching for new pharmacological targets

Mitochondria and mitochondrial proteins represent a group of promising pharmacological target candidates in the search of new molecular targets and drugs to counteract the onset of hypertension and more in general cardiovascular diseases (CVDs). Indeed, several mitochondrial pathways result impaired in CVDs, showing ATP depletion and ROS production as common traits of cardiac tissue degeneration. Thus, targeting mitochondrial dysfunction in cardiomyocytes can represent a successful strategy to prevent heart failure. In this context, the identification of new pharmacological targets among mitochondrial proteins paves the way for the design of new selective drugs. Thanks to the advances in omics approaches, to a greater availability of mitochondrial crystallized protein structures and to the development of new computational approaches for protein 3D-modelling and drug design, it is now possible to investigate in detail impaired mitochondrial pathways in CVDs. Furthermore, it is possible to design new powerful drugs able to hit the selected pharmacological targets in a highly selective way to rescue mitochondrial dysfunction and prevent cardiac tissue degeneration. The role of mitochondrial dysfunction in the onset of CVDs appears increasingly evident, as reflected by the impairment of proteins involved in lipid peroxidation, mitochondrial dynamics, respiratory chain complexes, and membrane polarization maintenance in CVD patients. Conversely, little is known about proteins responsible for the cross-talk between mitochondria and cytoplasm in cardiomyocytes. Mitochondrial transporters of the SLC25A family, in particular, are responsible for the translocation of nucleotides (e.g., ATP), amino acids (e.g., aspartate, glutamate, ornithine), organic acids (e.g. malate and 2-oxoglutarate), and other cofactors (e.g., inorganic phosphate, NAD⁺, FAD, carnitine, CoA derivatives) between the mitochondrial and cytosolic compartments. Thus, mitochondrial transporters play a key role in the mitochondria-cytosol cross-talk by leading metabolic pathways such as the malate/aspartate shuttle, the carnitine shuttle, the ATP export from mitochondria, and the regulation of permeability transition pore opening. Since all these pathways are crucial for maintaining healthy cardiomyocytes, mitochondrial carriers emerge as an interesting class of new possible pharmacological targets for CVD treatments.

Adaptive sequencing using nanopores and deep learning of mitochondrial DNA

Nanopore sequencing is an emerging technology that reads DNA by utilizing a unique method of detecting nucleic acid sequences and identifies the various chemical modifications they carry. Deep learning has increased in popularity as a useful technique to solve many complex computational tasks. 'Adaptive sequencing' is an implementation of selective sequencing, intended for use on the nanopore sequencing platform. In this study, we demonstrated an alternative method of software-based selective sequencing that is performed in real time by combining nanopore sequencing and deep learning. Our results showed the feasibility of using deep learning for classifying signals from only the first 200 nucleotides in a raw nanopore sequencing signal format. This was further demonstrated by comparing the accuracy of our deep learning classification model across data from several human cell lines and other eukaryotic organisms. We used custom deep learning models and a script that utilizes a 'Read Until' framework to target mitochondrial molecules in real time from a human cell line sample. This achieved a significant separation and enrichment ability of 2.3-fold. In a series of very short sequencing experiments (10, 30 and 120 min), we identified genomic and mitochondrial reads with accuracy above 90%, although mitochondrial DNA comprised only 0.1% of the total input material. The uniqueness of our method is the ability to distinguish two groups of DNA even without a labeled reference. This contrasts with studies that required a well-defined reference, whether of a DNA sequence or of another type of representation. Additionally, our method showed higher correlation to the theoretically possible enrichment factor, compared with other published methods. We believe that our results will lay the foundation for rapid and selective sequencing using nanopore technology and will pave the approach for clinical applications that use nanopore sequencing data.

Examining the Association between Mitochondrial Genome Variation and Coronary Artery Disease

Large-scale genome-wide association studies have identified hundreds of single-nucleotide variants (SNVs) significantly associated with coronary artery disease (CAD). However, collectively, these explain <20% of the heritability. Hypothesis: Here, we hypothesize that mitochondrial (MT)-SNVs might present one potential source of this “missing heritability”. Methods: We analyzed 265 MT-SNVs in ~500,000 UK Biobank individuals, exploring two different CAD definitions: a more stringent (myocardial infarction and/or revascularization; HARD = 20,405), and a more inclusive (angina and chronic ischemic heart disease; SOFT = 34,782). Results: In HARD cases, the most significant (p < 0.05) associations were for m.295C>T (control region) and m.12612A>G (ND5), found more frequently in cases (OR = 1.05), potentially related to reduced cardiorespiratory fitness in response to exercise, as well as for m.12372G>A (ND5) and m.11467A>G (ND4), present more frequently in controls (OR = 0.97), previously associated with lower ROS production rate. In SOFT cases, four MT-SNVs survived multiple testing corrections (at FDR < 5%), all potentially conferring increased CAD risk. Of those, m.11251A>G (ND4) and m.15452C>A (CYB) have previously shown significant associations with body height. In line with this, we observed that CAD cases were slightly less physically active, and their average body height was ~2.00 cm lower compared to controls; both traits are known to be related to increased CAD risk. Gene-based tests identified CO2 associated with HARD/SOFT CAD, whereas ND3 and CYB associated with SOFT cases (p < 0.05), dysfunction of which has been related to MT oxidative stress, obesity/T2D (CO2), BMI (ND3), and angina/exercise intolerance (CYB). Finally, we observed that macro-haplogroup I was significantly (p < 0.05) more frequent in HARD cases vs. controls (3.35% vs. 3.08%), potentially associated with response to exercise. Conclusions: We found only spurious associations between MT genome variation and HARD/SOFT CAD and conclude that more MT-SNV data in even larger study cohorts may be needed to conclusively determine the role of MT DNA in CAD.

Advanced approach for comprehensive mtDNA genome testing in mitochondrial disease

Mitochondrial disease diagnosis requires interrogation of both nuclear and mitochondrial (mtDNA) genomes for single-nucleotide variants (SNVs) and copy number alterations, both in the proband and often maternal relatives, together with careful phenotype correlation. We developed a comprehensive mtDNA sequencing test ('MitoGenome') using long-range PCR (LR-PCR) to amplify the full length of the mtDNA genome followed by next generation sequencing (NGS) to accurately detect SNVs and large-scale mtDNA deletions (LSMD), combined with droplet digital PCR (ddPCR) for LSMD heteroplasmy quantification. Overall, MitoGenome tests were performed on 428 samples from 394 patients with suspected or confirmed mitochondrial disease. The positive yield was 11% (43/394), including 34 patients with pathogenic or likely pathogenic SNVs (the most common being m.3243A > G in 8/34 (24%) patients), 8 patients with single LSMD, and 3 patients with multiple LSMD exceeding 10% heteroplasmy levels. Two patients with both LSMD and pathogenic SNV were detected. Overall, this LR-PCR/NGS assay provides a highly accurate and comprehensive diagnostic method for simultaneous mtDNA SNV detection at heteroplasmy levels as low as 1% and LSMD detection at heteroplasmy levels below 10%. Inclusion of maternal samples for variant classification and ddPCR to quantify LSMD heteroplasmy levels further enables accurate pathogenicity assessment and clinical correlation interpretation of mtDNA genome sequence variants and copy number alterations.