Mitochondrial Genetic Research Articles

A Successful Approach to Extracted Human mtDNA and Amplify some Genes from FFPE Samples

Tissues that have been paraffin-embedded and fixed in formalin (FFPE) are a significant diagnostic resource for molecular and pathophysiological methods of cancer and numerous other illnesses. The most difficult obstacle is obtaining high-quality of DNA from FFPE tissues, especially mitochondrial DNA (mtDNA), as formalin fixation significantly compromises the integrity of DNA. To extract mtDNA, the Geneaid gSYNCTM DNA kit (cat#GS100) was used with fifty breast cancer tissues (block). To check the mtDNA quality, specific primers for the following genes, ATP6 gene (8527-9024) and D-Loop region (256-16277), were used to run a PCR reaction. Our data showed that the mtDNA was successfully extracted with high quantity and quality. The two genes were amplified via PCR with very high rate (97% average). In conclusion, the method described here has many possibilities for applications in various molecular research projects involving FFPET substances. It can be used to obtain a high quantity and quality of mtDNA from FFPE samples for different molecular diagnostic and research methods such as PCR, qPCR, Sanger sequencing, and whole genome sequencing.

Chimeric mitochondrial RNA transcripts predict mitochondrial genome deletion mutations in mitochondrial genetic diseases and aging.

While it is well understood that mitochondrial DNA (mtDNA) deletion mutations cause incurable diseases and contribute to aging, little is known about the transcriptional products that arise from these DNA structural variants. We hypothesized that mitochondrial genomes containing deletion mutations express chimeric mitochondrial RNAs. To test this, we analyzed human and rat RNA sequencing data to identify, quantitate, and characterize chimeric mitochondrial RNAs. We observed increased chimeric mitochondrial RNA frequency in samples from patients with mitochondrial genetic diseases and in samples from aged humans. The spectrum of chimeric mitochondrial transcripts reflected the known pattern of mtDNA deletion mutations. To test the hypothesis that mtDNA deletions induce chimeric RNA transcripts, we treated 18 mo and 34 mo rats with guanidinopropionic acid to induce high levels of skeletal muscle mtDNA deletion mutations. With mtDNA deletion induction, we demonstrate that the chimeric mitochondrial transcript frequency also increased and correlated strongly with an orthogonal DNA-based mutation assay performed on identical samples. Further, we show that the frequency of chimeric mitochondrial transcripts predicts expression of both nuclear and mitochondrial genes central to mitochondrial function, demonstrating the utility of these events as metrics of age-induced metabolic change. Mapping and quantitation of chimeric mitochondrial RNAs provides an accessible, orthogonal approach to DNA-based mutation assays, offers a potential method for identifying mitochondrial pathology in widely accessible datasets, and opens a new area of study in mitochondrial genetics and transcriptomics.

TreeWave: command line tool for alignment-free phylogeny reconstruction based on graphical representation of DNA sequences and genomic signal processing

BackgroundGenomic sequence similarity comparison is a crucial research area in bioinformatics. Multiple Sequence Alignment (MSA) is the basic technique used to identify regions of similarity between sequences, although MSA tools are widely used and highly accurate, they are often limited by computational complexity, and inaccuracies when handling highly divergent sequences, which leads to the development of alignment-free (AF) algorithms.ResultsThis paper presents TreeWave, a novel AF approach based on frequency chaos game representation and discrete wavelet transform of sequences for phylogeny inference. We validate our method on various genomic datasets such as complete virus genome sequences, bacteria genome sequences, human mitochondrial genome sequences, and rRNA gene sequences. Compared to classical methods, our tool demonstrates a significant reduction in running time, especially when analyzing large datasets. The resulting phylogenetic trees show that TreeWave has similar classification accuracy to the classical MSA methods based on the normalized Robinson-Foulds distances and Baker’s Gamma coefficients.ConclusionsTreeWave is an open source and user-friendly command line tool for phylogeny reconstruction. It is a faster and more scalable tool that prioritizes computational efficiency while maintaining accuracy. TreeWave is freely available at https://github.com/nasmaB/TreeWave.

Similarity of Human Mitochondrial DNA Nucleotide Substitution Spectra Reconstructed Over One and Many Generations

Using phylogenetic analysis of mitochondrial whole genome nucleotide sequences (mtDNA), allowing the study of genetic changes over many generations, a spectrum of nucleotide substitutions (along the L-strand of mtDNA) was reconstructed in European populations. The spectra of mtDNA nucleotide substitutions observed in a heteroplasmic state (at ≥1% and ≥5% levels) in first generation children were also analyzed. It was found that the spectra of nucleotide substitutions reconstructed over one and many generations do not differ practically in the main parameters: the distribution of pyrimidine and purine substitutions (with predominance of transitions TC), the ratio of the number of transitions and transversions. Analysis of the phylogenetic tree of mtDNA haplotypes in Europeans clearly revealed the influence of negative (purifying) selection on mitochondrial gene pools. It is suggested that the selective processes guiding the mtDNA evolution in one and many generations are of a similar nature, i.e., caused by negative selection. The problem of how mutations occur and spread in mitochondria of germ line cells is discussed.

My Early Years of Yeast Mitochondrial Genetics.

There have been massive technological advances in molecular biology and genetics over the past five decades. I have personally experienced these advances and here I reflect on those origins, from my perspective, studying yeast mitochondrial genetics leading up to deciphering the functions of the mitochondrial genome. The yeast contributions commenced in the middle of the last century with pure genetics, correlating mutants with phenotypes, in order to discover genes, just like the early explorations to discover new lands. The quest was to explore the mitochondrial genome and find its genes and their products. It was most fortunate that DNA sequencing technologies became available in the late 1970s, and laboratories were restructured enormously to keep pace with the emerging technologies. There were considerable costs in equipping laboratories, purchasing ultracentrifuges and restriction endonucleases, and undertaking DNA sequencing; additionally, workers required special safety gear.

Quantifying constraint in the human mitochondrial genome.

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

Copy Number Analysis of Whole Exome Sequencing Data

Abstract Background Genetic alterations implicated in constitutional disorders range in size from single nucleotide substitutions to losses and gains of millions of bases. Whole exome sequencing (WES) captures many clinically relevant variants, but an analytic pipeline optimized for detection of small substitutions, insertions, and deletions may miss pathogenic copy number variants (CNVs). Historically, our laboratory performed separate cytogenetic tests to evaluate for larger genomic imbalances. Addition of copy number analysis to our WES pipeline could improve the assay’s diagnostic yield. This study describes part of our strategy to validate copy number analysis of our constitutional whole exome sequencing data. Methods We identified 18 CNV specimens that underwent both chromosomal microarray analysis (CMA) and WES between 2017 and 2023. Variants included 13 deletions (size range 14Kb-2Mb) and 5 duplications (size range 325Kb-10Mb). The standard method of constitutional copy number analysis in our clinical laboratory is CMA using the Affymetrix Cytoscan HD Array and analysis software. CNVs identified on CMA were reviewed with reference to the Database of Genomic Variants and multiple academic hospital databases. WES was performed on an Illumina sequencing system following capture-based library enrichment with the Integrated DNA Technologies xGen Exome Research, CNV Backbone, Human mtDNA Research, and Human ID Research Panel probes. CNVs identified on WES were called using the Broad Institute’s Genomic Analysis Toolkit (GATK) with a quality score cutoff of 120. Whole exome sequencing data was visualized on the Integrative Genomics Viewer. Results 14 out of 18 CNVs reported on CMA were also detected on WES with a quality score at or above 120 (78%). This includes 10 out of 13 deletions (77%) and 3 out of 5 duplications (60%). For both deletions and duplications, the size ranges of detected variants at or above the cutoff score and undetected variants overlapped. One duplication reported on CMA was detected on WES with a quality score below 120. Visualization of whole exome sequencing data showed variation in probe density that may have affected CNV calling. Conclusions The majority of CNVs reported after CMA (14/18, 78%) were also detected on WES data using a GATK quality score cutoff of 120. The data does not suggest an association between CNV size and detection by WES, although interpretation is limited by the small number of cases. Data visualization may be required to identify some CNVs. Future steps to characterize the limitations of copy number analysis of WES data include assessments of variance in probe density and read depth across the exome. Copy number analysis of WES data from cases with normal CMA results will also be performed for assay validation.

The 9 bp Deletion between the Mitochondrial COII and Lysine tRNA Genes in a Caucasian Population with Cognitive Disorders: An Observational Study.

The loss of one of the two copies of the 9 bp tandem repeat sequence (CCCCCTCTA) located in the small non-coding region between the cytochrome oxidase II (COII) and the lysine tRNA genes in human mtDNA has been reported to be polymorphic in Asian, Oceanian and Sub-Saharan African populations, but it has rarely been observed in Europe. In this study, we will evaluate the possible association between the MIC9D polymorphism and cognitive disorders. A genetic analysis of unrelated Sicilian patients with cognitive deficits was performed to identify the 9 bp deletion MIC9D polymorphism. The MIC9D polymorphism was found in six patients, whereas this variant was absent in control individuals without cognitive deficits. The patients with the MIC9D polymorphism exhibited more complex clinical presentations; in particular, all had neuromuscular disorders and five also presented with behavioral disorders. The present study suggests a potential association between the MIC9D polymorphism and cognitive impairment with concurrent neuromuscular and behavioral involvement.

A preliminary study on detecting human DNA in aquatic environments: Potential of eDNA in forensics

Human environmental DNA (eDNA) application have not been fully applied or adequately considered in the fields of eDNA and forensics. Nonetheless, this technique holds great potential as a complementary tool for detecting human DNA in aquatic environments, particularly in cases involving crimes connect to such environments. However, the detectability or stability of eDNA can vary depending on several factors. Therefore, this preliminary study investigates the detection and degradation rates of human eDNA, as well as the recovery of nuclear short tandem repeat (STR) profiles and mitochondrial DNA (mtDNA) sequencing, using water samples from both saltwater and freshwater sources. To conduct the experiment, whole human blood was spiked into the water samples. Water samples were then filtered using a 5 µm pore size filter, and samples were collected at various time intervals up to 23 days. A human specific qPCR assay targeting HV1 region of human mtDNA was used to detect human eDNA. Results demonstrated that human eDNA remains detectable for up to 36 hours in freshwater samples and up 84 hours in saltwater samples. The limit of detection (LOD) of human eDNA, (205 copies/µl), was achieved after 60 hours in freshwater and 180 hours in saltwater samples. Partial STR profiles could be recovered up to 24 hours for freshwater and saltwater. Results from mtDNA sequencing indicate that full mtDNA profiles could be recovered from freshwater samples up to 48 hours and remained detectable up to 72 hours in saltwater. Overall, the findings of this study underscore the importance of considering and incorporating human eDNA analysis as a valuable tool in forensic practice. By harnessing the power of eDNA, law enforcement agencies can enhance their investigation capabilities, improve the accuracy of forensic reconstructions, and ultimately contribute to the resolution of cases involving aquatic environments. Further research and validation are needed to optimize and expand the utilization of eDNA techniques in forensic investigations.

Unconstrained Precision Mitochondrial Genome Editing with αDdCBEs.

DddA-derived cytosine base editors (DdCBEs) enable the targeted introduction of C•G-to-T•A conversions in mitochondrial DNA (mtDNA). DdCBEs work in pairs, with each arm composed of a transcription activator-like effector (TALE), a split double-stranded DNA deaminase half, and a uracil glycosylase inhibitor. This pioneering technology has helped improve our understanding of cellular processes involving mtDNA and has paved the way for the development of models and therapies for genetic disorders caused by pathogenic mtDNA variants. Nonetheless, given the intrinsic properties of TALE proteins, several target sites in human mtDNA are predicted to remain out of reach to DdCBEs and other TALE-based technologies. Specifically, due to the conventional requirement for a thymine immediately upstream of the TALE target sequences (i.e., the 5'-T constraint), over 150 loci in the human mitochondrial genome are presumed to be inaccessible to DdCBEs. Previous attempts at circumventing this requirement, either by developing monomeric DdCBEs or utilizing DNA-binding domains alternative to TALEs, have resulted in suboptimal specificity profiles with reduced therapeutic potential. Here, aiming to challenge and elucidate the relevance of the 5'-T constraint in the context of DdCBE-mediated mtDNA editing, and to expand the range of motifs that are editable by this technology, we generated DdCBEs containing TALE proteins engineered to recognize all 5' bases. These modified DdCBEs are herein referred to as αDdCBEs. Notably, 5'-T-noncompliant canonical DdCBEs efficiently edited mtDNA at diverse loci. However, they were frequently outperformed by αDdCBEs, which exhibited significant improvements in activity and specificity, regardless of the most 5' bases of their TALE binding sites. Furthermore, we showed that αDdCBEs are compatible with the enhanced DddAtox variants DddA6 and DddA11, and we validated TALE shifting with αDdCBEs as an effective approach to optimize base editing outcomes. Overall, αDdCBEs enable efficient, specific, and unconstrained mitochondrial base editing.

Human sperm mitochondrial DNA copy numbers and deletion rates: Comparing persons living in two urban industrial agglomerations differing in sources of air pollution

Persons living in industrial environments are exposed to levels of air pollution that can affect their health and fertility. The Czech capital city, Prague, and the Ostrava industrial agglomeration differ in their major sources of air pollution. In Prague, heavy traffic produces high levels of nitrogen oxides throughout the year. In the Ostrava region, an iron industry and local heating are sources of particulate matter (PM) and benzo[a]pyrene (B[a]P), especially in the winter. We evaluated the effects of air pollution on human sperm mitochondrial DNA (mtDNA). Using real-time PCR, we analysed sperm mtDNA copy number and deletion rate in Prague city policemen in two seasons (spring and autumn) and compared the results with those from Ostrava. In Prague, the sperm mtDNA deletion rate was significantly higher in autumn than in spring, which is the opposite of the results from Ostrava. The sperm mtDNA copy number did not show any seasonal differences in either of the cities; it was correlated negatively with sperm concentration, motility, and viability, and with sperm chromatin integrity (assessed with the Sperm Chromatin Structure Assay). The comparison between the two cities showed that the sperm mtDNA deletion rate in spring and the sperm mtDNA copy number in autumn were significantly lower in Prague vs. Ostrava. Our study supports the hypothesis that sperm mtDNA deletion rate is affected by the composition of air pollution. Sperm mtDNA abundance is closely associated with chromatin damage and standard semen characteristics.

Rapamycin's lifespan effect is modulated by mito-nuclear epistasis in Drosophila.

The macrolide drug rapamycin is a benchmark anti-ageing drug, which robustly extends lifespan of diverse organisms. For any health intervention, it is paramount to establish whether benefits are distributed equitably among individuals and populations, and ideally to match intervention to recipients' needs. However, how responses to rapamycin vary is surprisingly understudied. Here we investigate how among-population variation in both mitochondrial and nuclear genetics shapes rapamycin's effects on lifespan. We show that epistatic "mito-nuclear" interactions, between mitochondria and nuclei, modulate the response to rapamycin treatment. Differences manifest as differential demographic effects of rapamycin, with altered age-specific mortality rate. However, a fitness cost of rapamycin early in life does not show a correlated response, suggesting that mito-nuclear epistasis can decouple costs and benefits of treatment. These findings suggest that a deeper understanding of how variation in mitochondrial and nuclear genomes shapes physiology may facilitate tailoring of anti-ageing therapy to individual need.

Gulf war illness: a tale of two genomes

IntroductionGulf War illness (GWI) is an environmentally-triggered chronic multisymptom illness typified by protean symptoms, in which mitochondrial impairment is evident. It has been likened to accelerated aging. Nuclear genetics of detoxification have been linked to GWI.ObjectiveTo see whether mitochondrial (mt) haplogroup U – a heritable profile of mitochondrial DNA that has been tied to aging-related conditions – significantly predicts greater GWI severity; and to assess whether GWI severity is influenced by mitochondrial as well as nuclear genetics. 54 consenting Gulf War veterans gave information on GWI severity, of whom 52 had nuclear DNA assessment; and 45 had both nuclear and mitochondrial DNA assessments. Regression with robust standard errors assessed prediction of GWI severity as a function of nuclear genetics (butyrylcholinesterase variants), mitochondrial genetics (haplogroup U, previously tied to aging-related conditions); or both.ResultsBChE “adverse” variants significantly predicted GWI severity (β(SE) = 23.4(11.4), p = 0.046), as did mt haplogroup U (β(SE) = 36.4(13.6), p = 0.010). In a model including both, BChE was no longer significant, but mt haplogroup U retained significance (β(SE) = 36.7(13.0), p = 0.007). This is the first study to show that mitochondrial genetics are tied to GWI severity in Gulf-deployed veterans. Other data affirm a tie to nuclear genetics, making GWI indeed a “tale of two genomes.”

Exploring Mitochondrial Heterogeneity and Evolutionary Dynamics in Thelephora ganbajun through Population Genomics.

Limited exploration in fungal mitochondrial genetics has uncovered diverse inheritance modes. The mitochondrial genomes are inherited uniparentally in the majority of sexual eukaryotes, our discovery of persistent mitochondrial heterogeneity within the natural population of the basidiomycete fungus Thelephora ganbajun represents a significant advance in understanding mitochondrial inheritance and evolution in eukaryotes. Here, we present a comprehensive analysis by sequencing and assembling the complete mitogenomes of 40 samples exhibiting diverse cox1 heterogeneity patterns from various geographical origins. Additionally, we identified heterogeneous variants in the nad5 gene, which, similar to cox1, displayed variability across multiple copies. Notably, our study reveals a distinct prevalence of introns and homing endonucleases in these heterogeneous genes. Furthermore, we detected potential instances of horizontal gene transfer involving homing endonucleases. Population genomic analyses underscore regional variations in mitochondrial genome composition among natural samples exhibiting heterogeneity. Thus, polymorphisms in heterogeneous genes, introns, and homing endonucleases significantly influence mitochondrial structure, structural variation, and evolutionary dynamics in this species. This study contributes valuable insights into mitochondrial genome architecture, population dynamics, and the evolutionary implications of mitochondrial heterogeneity in sexual eukaryotes.

Highly accurate single-color fluorogenic DNA decoding sequencing for mutational genotyping

We proposed a single-color fluorogenic DNA decoding sequencing method designed to improve sequencing accuracy, increase read length and throughput, as well as decrease scanning time. This method involves the incorporation of a mixture of four types of 3’-O-modified nucleotide reversible terminators into each reaction. Among them, two nucleotides are labeled with the same fluorophore, while the remaining two are unlabeled. Only one nucleotide can be extended in each reaction, and an encoding that partially defines base composition can be obtained. Through cyclic interrogation of a template twice with different nucleotide combinations, two sets of encodings are sequentially obtained, enabling the determination of the sequence. We demonstrate the feasibility of this method using established sequencing chemistry, achieving a cycle efficiency of approximately 99.5 %. Notably, this strategy exhibits remarkable efficacy in the detection and correction of sequencing errors, achieving a theoretical error rate of 0.00016 % at a sequencing depth of ×2, which is lower than Sanger sequencing. This method is theoretically compatible with the existing sequencing-by-synthesis (SBS) platforms, and the instrument is simpler, which may facilitate further reductions in sequencing costs, thereby broadening its applications in biology and medicine. Moreover, we demonstrate the capability to detect known mutation sites using information from only a single sequencing run. We validate this approach by accurately identifying a mutation site in the human mitochondrial DNA.

Replication and Transcription of Human Mitochondrial DNA.

Mammalian mitochondrial DNA (mtDNA) is replicated and transcribed by phage-like DNA and RNA polymerases, and our understanding of these processes has progressed substantially over the last several decades. Molecular mechanisms have been elucidated by biochemistry and structural biology and essential in vivo roles established by cell biology and mouse genetics. Single molecules of mtDNA are packaged by mitochondrial transcription factor A into mitochondrial nucleoids, and their level of compaction influences the initiation of both replication and transcription. Mutations affecting the molecular machineries replicating and transcribing mtDNA are important causes of human mitochondrial disease, reflecting the critical role of the genome in oxidative phosphorylation system biogenesis. Mechanisms controlling mtDNA replication and transcription still need to be clarified, and future research in this area is likely to open novel therapeutic possibilities for treating mitochondrial dysfunction.

Similarity of Human Mitochondrial DNA Nucleotide Substitution Spectra Reconstructed over One and Many Generations

Similarity of Human Mitochondrial DNA Nucleotide Substitution Spectra Reconstructed over One and Many Generations

The human mitochondrial mRNA structurome reveals mechanisms of gene expression.

The human mitochondrial genome encodes crucial oxidative phosphorylation system proteins, pivotal for aerobic energy transduction. They are translated from nine monocistronic and two bicistronic transcripts whose native structures remain unexplored, posing a gap in understanding mitochondrial gene expression. In this work, we devised the mitochondrial dimethyl sulfate mutational profiling with sequencing (mitoDMS-MaPseq) method and applied detection of RNA folding ensembles using expectation-maximization (DREEM) clustering to unravel the native mitochondrial messenger RNA (mt-mRNA) structurome in wild-type (WT) and leucine-rich pentatricopeptide repeat-containing protein (LRPPRC)-deficient cells. Our findings elucidate LRPPRC's role as a holdase contributing to maintaining mt-mRNA folding and efficient translation. mt-mRNA structural insights in WT mitochondria, coupled with metabolic labeling, unveil potential mRNA-programmed translational pausing and a distinct programmed ribosomal frameshifting mechanism. Our data define a critical layer of mitochondrial gene expression regulation. These mt-mRNA folding maps provide a reference for studying mt-mRNA structures in diverse physiological and pathological contexts.

Mitoclone2: an R package for elucidating clonal structure in single-cell RNA-sequencing data using mitochondrial variants.

Clonal cell population dynamics play a critical role in both disease and development. Due to high mitochondrial mutation rates under both healthy and diseased conditions, mitochondrial genomic variability is a particularly useful resource in facilitating the identification of clonal population structure. Here we present mitoClone2, an all-inclusive R package allowing for the identification of clonal populations through integration of mitochondrial heteroplasmic variants discovered from single-cell sequencing experiments. Our package streamlines the investigation of this phenomenon by providing: built-in compatibility with commonly used tools for the delineation of clonal structure, the ability to directly use multiplexed BAM files as input, annotations for both human and mouse mitochondrial genomes, and helper functions for calling, filtering, clustering, and visualizing variants.

Mitochondrial genetics through the lens of single-cell multi-omics.

Mitochondria carry their own genetic information encoding for a subset of protein-coding genes and translational machinery essential for cellular respiration and metabolism. Despite its small size, the mitochondrial genome, its natural genetic variation and molecular phenotypes have been challenging to study using bulk sequencing approaches, due to its variation in cellular copy number, non-Mendelian modes of inheritance and propensity for mutations. Here we highlight emerging strategies designed to capture mitochondrial genetic variation across individual cells for lineage tracing and studying mitochondrial genetics in primary human cells and clinical specimens. We review recent advances surrounding single-cell mitochondrial genome sequencing and its integration with functional genomic readouts, including leveraging somatic mitochondrial DNA mutations as clonal markers that can resolve cellular population dynamics in complex human tissues. Finally, we discuss how single-cell whole mitochondrial genome sequencing approaches can be utilized to investigate mitochondrial genetics and its contribution to cellular heterogeneity and disease.