Quantifying mitochondrial DNA copy number using robust regression to interpret real time PCR results

Biological and translational attributes of mitochondrial DNA copy number: Laboratory perspective to clinical relevance

The mitochondrial DNA copy number (mtDNAcn) plays a vital role in cellular energy metabolism and mitochondrial health. As mitochondria are responsible for adenosine triphosphate production through oxidative phosphorylation, maintaining an appropriate mtDNAcn level is vital for the overall cellular function. Alterations in mtDNAcn have been linked to various diseases, including neurodegenerative disorders, metabolic conditions, and cancers, making it an important biomarker for understanding the disease pathogenesis. The accurate estimation of mtDNAcn is essential for clinical applications. Quantitative polymerase chain reaction and next-generation sequencing are commonly employed techniques with distinct advantages and limitations. Clinically, mtDNAcn serves as a valuable indicator for early diagnosis, disease progression, and treatment response. For instance, in oncology, elevated mtDNAcn levels in blood samples are associated with tumor aggressiveness and can aid in monitoring treatment efficacy. In neurodegenerative diseases such as Alzheimer’s and Parkinson’s, altered mtDNAcn patterns provide insights into disease mechanisms and progression. Understanding and estimating mtDNAcn are critical for advancing diagnostic and therapeutic strategies in various medical fields. As research continues to uncover the implications of mtDNAcn alterations, its potential as a clinical biomarker is likely to expand, thereby enhancing our ability to diagnose and manage complex diseases.

Valproic Acid and Lamotrigine Differentially Modulate the Telomere Length in Epilepsy Patients

Background/Objectives: Antiseizure drugs (ASDs) are the primary therapy for epilepsy, and the choice varies according to seizure type. Epilepsy patients experience chronic mitochondrial oxidative stress and increased levels of pro-inflammatory mediators, recognizable hallmarks of biological aging; however, few studies have explored aging markers in epilepsy. Herein, we addressed for the first time the impact of ASDs on molecular aging by measuring the telomere length (TL) and mtDNA copy number (mtDNA-CN). Methods: We used real-time quantitative PCR (QPCR) in epilepsy patients compared to matched healthy controls (CTs) and assessed the association with plasma levels of ASDs and other clinical variables. The sample comprised 64 epilepsy patients and 64 CTs. Patients were grouped based on monotherapy with lamotrigine (LTG) or valproic acid (VPA), and those treated with a combination therapy (LTG + VPA). Multivariable logistic regression was applied to analyze the obtained data. Results: mtDNA-CN was similar between patients and controls, and none of the comparisons were significant for this marker. TL was shorter in not seizure-free patients than in CTs (1.50 ± 0.35 vs. 1.68 ± 0.34; p < 0.05), regardless of the ASD therapy. These patients exhibited the highest proportion of adverse drug reactions. TL was longer in patients on VPA monotherapy, followed by patients on LTG monotherapy and patients on an LTG + VPA combined scheme (1.77 ± 0.24; 1.50 ± 0.32; 1.36 ± 0.37, respectively; p < 0.05), suggesting that ASD treatment differentially modulates TL. Conclusions: Our findings suggest that clinicians could consider TL measurements to decide the best ASD treatment option (VPA and/or LTG) to help predict ASD responses in epilepsy patients.

Mdivi-1: Effective but complex mitochondrial fission inhibitor

Mdivi-1, Mitochondrial DIVIsion inhibitor 1, has been widely employed in research under the assumption that it exclusively influences mitochondrial fusion, but effects other than mitochondrial dynamics have been underinvestigated. This paper provides transcriptome and DNA methylome-wide analysis for Mdivi-1 treated SH-SY5Y human neuroblastoma cells using RNA sequencing (RNA-seq) and methyl capture sequencing (MC-seq) methods. Gene ontology analysis of RNA sequences revealed that p53 transcriptional gene network and DNA replication initiation-related genes were significantly up and down-regulated, respectively, showing the correlation with the arrest cell cycle in the G1 phase. MC-seq, a powerful sequencing method for capturing DNA methylation status in CpG sites, revealed that although Mdivi-1 does not induce dramatic DNA methylation change, the subtle alterations were concentrated within the CpG island. Integrative analysis of both sequencing data disclosed that the p53 transcriptional network was activated while the Parkinson's disease pathway was halted. Next, we investigated several changes in mitochondria in response to Mdivi-1. Copy number and transcription of mitochondrial DNA were suppressed. ROS levels increased, and elevated ROS triggered mitochondrial retrograde signaling rather than inducing direct DNA damage. In this study, we could better understand the molecular network of Mdivi-1 by analyzing DNA methylation and mRNA transcription in the nucleus and further investigating various changes in mitochondria, providing inspiration for studying nuclear-mitochondrial communications.

Longitudinal Changes in Mitochondrial DNA Copy Number and Telomere Length in Patients with Parkinson's Disease

Parkinson's disease (PD) pathophysiology includes mitochondrial dysfunction, neuroinflammation, and aging as its biggest risk factors. Mitochondrial DNA copy number (mtDNA-CN) and telomere length (TL) are biological aging markers with inconclusive results regarding their association with PD. A case-control study was used to measure TL and mtDNA-CN using qPCR in PBMCs. PD patients were naive at baseline (T0) and followed-up at one (T1) and two (T2) years after the dopaminergic treatment (DRT). Plasmatic cytokines were determined by ELISA in all participants, along with clinical parameters of patients at T0. While TL was shorter in patients vs. controls at all time points evaluated (p < 0.01), mtDNA-CN showed no differences. An increase in mtDNA-CN and TL was observed in treated patients vs. naive ones (p < 0.001). Our statistical model analyzed both aging markers with covariates, showing a strong correlation between them (r = 0.57, p < 0.01), and IL-17A levels positively correlating with mtDNA-CN only in untreated patients (r = 0.45, p < 0.05). TL and mtDNA-CN could be useful markers for monitoring inflammation progression or treatment response in PD. DRT might modulate TL and mtDNA-CN, reflecting a compensatory mechanism to counteract mitochondrial dysfunction in PD, but this needs further investigation.

Mitochondria as the Target of Hepatotoxicity and Drug-Induced Liver Injury: Molecular Mechanisms and Detection Methods

One of the major mechanisms of drug-induced liver injury includes mitochondrial perturbation and dysfunction. This is not a surprise, given that mitochondria are essential organelles in most cells, which are responsible for energy homeostasis and the regulation of cellular metabolism. Drug-induced mitochondrial dysfunction can be influenced by various factors and conditions, such as genetic predisposition, the presence of metabolic disorders and obesity, viral infections, as well as drugs. Despite the fact that many methods have been developed for studying mitochondrial function, there is still a need for advanced and integrative models and approaches more closely resembling liver physiology, which would take into account predisposing factors. This could reduce the costs of drug development by the early prediction of potential mitochondrial toxicity during pre-clinical tests and, especially, prevent serious complications observed in clinical settings.

Emerging methods for and novel insights gained by absolute quantification of mitochondrial DNA copy number and its clinical applications

The past thirty years have seen a surge in interest in pathophysiological roles of mitochondria, and the accurate quantification of mitochondrial DNA copy number (mCN) in cells and tissue samples is a fundamental aspect of assessing changes in mitochondrial health and biogenesis. Quantification of mCN between studies is surprisingly variable due to a combination of physiological variability and diverse protocols being used to measure this endpoint. The advent of novel methods to quantify nucleic acids like digital polymerase chain reaction (dPCR) and high throughput sequencing offer the ability to measure absolute values of mCN. We conducted an in-depth survey of articles published between 1969 -- 2020 to create an overview of mCN values, to assess consensus values of tissue-specific mCN, and to evaluate consistency between methods of assessing mCN. We identify best practices for methods used to assess mCN, and we address the impact of using specific loci on the mitochondrial genome to determine mCN. Current data suggest that clinical measurement of mCN can provide diagnostic and prognostic value in a range of diseases and health conditions, with emphasis on cancer and cardiovascular disease, and the advent of means to measure absolute mCN should improve future clinical applications of mCN measurements.

Characterization of mitogenomes from four Mucorales species and insights into pathogenicity

BACKGROUND

Mucorales, as one major order of Zygomycetes fungi, can infect human beings and cause serious consequence. We have noticed the pathogenicity of Mucorales is closely related to energy metabolism, while mitochondria play the role of energy factories in almost all biological activities.

METHODS

Virulence of M irregularis, M hiemalis, L corymbifera and R arrhizus were verified in Galleria mellonella larvae, as well as mitochondrial gene copies analysed with RT-qPCR. Mitogenomes of the four Mucorales species were sequenced based on illumina NovaSeq technology to study their characteristic features and functional regions.

RESULTS

Variant virulence of M irregularis, M hiemalis, L corymbifera and R arrhizu were verified by clinical retrospective data and our G mellonella infection models, also copies of mitochondrial genes indicated the significant associations with pathogenicity. A total of 274.18 clean reads were generated to be assembled; the complete mitogenomes of the four Mucorales species were obtained with totally different length. After the genomes annotated and compared, M irregularis was found more similar with M hiemalis than those of L corymbifera and R arrhizus, especially the small (rrns) and large (rrnl) subunits of mitochondrial ribosomal RNA (rRNA) genes. The GC content, ncRNAs and the distribution of the SNPs and InDels were also compared, and the GC content rate of fungi seems to be related to the fungal thermal adaptability. In addition, linear mitogenomes of these four Mucorales showed diverse arrangements of orf genes and directionality of some conserved gene elements.

CONCLUSION

This study uncovered the pathogenicity variances among the four Mucorales species and the relationship between their mitogenomic features and clinical pathogenicity. Further studies like spatial structure of mitochondrial genomes and the comprehensive analysis of transcription regulation are needed.

Targeted epigenetic induction of mitochondrial biogenesis enhances antitumor immunity in mouse model

Considering the potential of combinatorial therapies in overcoming existing limitations of cancer immunotherapy, there is an increasing need to identify small-molecule modulators of immune cells capable of augmenting the effect of programmed cell death protein 1 (PD-1) blockade, leading to better cancer treatment. Although epigenetic drugs showed potential in combination therapy, the lack of sequence specificity is a major concern. Here, we identify and develop a DNA-based epigenetic activator with tri-arginine vector called EnPGC-1 that can trigger the targeted induction of the peroxisome proliferator-activated receptor-gamma coactivator 1 alpha/beta (PGC-1α/β), a regulator of mitochondrial biogenesis. EnPGC-1 enhances mitochondrial activation, energy metabolism, proliferation of CD8⁺ T cells in vitro, and, in particular, enhances oxidative phosphorylation, a feature of long-lived memory T cells. Genome-wide gene analysis suggests that EnPGC-1 and not the control compounds can regulate T cell activation as a major biological process. EnPGC-1 also synergizes with PD-1 blockade to enhance antitumor immunity and improved host survival.

Somatic Mitochondrial DNA Point Mutations Used as Biomarkers to Demonstrate Genomic Heterogeneity in Primary Prostate Cancer

Primary prostate tumor heterogeneity is poorly understood, leaving research efforts with challenges regarding the initiation and advancement of the disease. The growth of tumor cells is accompanied by mutations in nuclear and in mitochondrial genomes. Thus, mitochondrial DNA mutations may be used as tumor cell markers. By the use of laser capture microdissection coupled with assays for mitochondrial point mutation detection, mtDNA mutations were used to trace mutated cells at a histological level. Point mutations in mtDNA were determined in 12 primary prostate cancers. The tumors represent different pathology-prognostic grade groups. Known mutational hotspots of the mtDNA were scanned for heteroplasmy. All specimens with mtDNA heteroplasmy were subsequently subsampled by laser capture microdissection. From a total number of 1728 microsamples, mitochondrial DNA target sequences were amplified and base substitutions detected by cycling temperature capillary electrophoresis. Real-time PCR was used as a quantitative assay to determine the relative mtDNA copy number of 12 tumors studied, represented by two samples from each (N = 24); a high degree (75%) demonstrated tumor specimen heterogeneity. A grid of 96 spots isolated by laser capture microdissection demonstrated interfocal sample heterogeneity and increased the limit of detection. The spots demonstrated a wide range of mutant fractions from 0 to 100% mutant copies. The mitochondrial DNA copy number in the samples was determined by real-time PCR. No correlation between copy number and pathology-prognostic grade groups was observed. Somatic mitochondrial DNA point mutations represent traceable biomarkers demonstrating heterogeneity in primary prostate cancer. Mutations can be detected in areas before changes in tissue histopathology are evident to the pathologist.

Role of mitochondrial fission-related genes in mitochondrial morphology and energy metabolism in mouse embryonic stem cells

Mitochondria, the major organelles that produce energy for cell survival and function, dynamically change their morphology via fusion and fission, a process called mitochondrial dynamics. The details of the underlying mechanism of mitochondrial dynamics have not yet been elucidated. Here, we aimed to investigate the function of mitochondrial fission genes in embryonic stem cells (ESCs). To this end, we generated homozygous knockout ESC lines, namely, Fis1^-/-, Mff^-/-, and Dnm1l^-/- ESCs, using the CRISPR-Cas9 system. Interestingly, the Fis1^-/-, Mff^-/-, and Dnm1l^-/- ESCs showed normal morphology, self-renewal, and the ability to differentiate into all three germ layers in vitro. However, transmission electron microscopy showed a significant increase in the cytoplasm to nucleus ratio and mitochondrial elongation in Dnm1l^-/- ESCs, which was due to incomplete fission. To assess the change in metabolic energy, we analyzed oxidative phosphorylation (OXPHOS), glycolysis, and the intracellular ATP concentration. The ESC knockout lines showed an increase in OXPHOS, decrease in glycolysis, and an increase in intracellular ATP concentration, which was related to mitochondrial elongation. In particular, the Dnm1l knockout most significantly affected mitochondrial morphology, energy metabolism, and ATP production in ESCs. Furthermore, RNA sequencing and gene ontology analysis showed that the differentially expressed genes in Mff^-/- ESCs were distinct from those in Dnm1l^-/- or Fis1^-/- ESCs. In total, five metabolism-related genes, namely, Aass, Cdo1, Cyp2b23, Nt5e, and Pck2, were expressed in all three knockout ESC lines, and three of them were associated with regulation of ATP generation.

GATA-1 isoforms differently contribute to the production and compartmentation of reactive oxygen species in the myeloid leukemia cell line K562

Maintenance of a balanced expression of the two isoforms of the transcription factor GATA‐1, the full‐length protein (GATA‐1_FL) and a shorter isoform (GATA‐1 _S), contributes to control hematopoiesis, whereas their dysregulation can alter the differentiation/proliferation potential of hematopoietic precursors thereby eventually leading to a variety of hematopoietic disorders. Although it is well established that these isoforms play opposite roles in these remarkable processes, most of the molecular pathways involved remain unknown.

Here, we demonstrate that GATA‐1_FL and GATA‐1_S are able to differently influence intracellular redox states and reactive oxygen species (ROS) compartmentation in the erythroleukemic K562 cell line, thus shedding novel mechanistic insights into the processes of cell proliferation and apoptosis resistance in myeloid precursors. Furthermore, given the role played by ROS signaling as a strategy to escape apoptosis and evade cell‐mediated immunity in myeloid cells, this study highlights a mechanism through which aberrant expression of GATA‐1 isoforms could play a role in the leukemogenic process.

Developmental validation of a Nextera XT mitogenome Illumina MiSeq sequencing method for high-quality samples

Generating mitochondrial genome (mitogenome) data from reference samples in a rapid and efficient manner is critical to harnessing the greater power of discrimination of the entire mitochondrial DNA (mtDNA) marker. The method of long-range target enrichment, Nextera XT library preparation, and Illumina sequencing on the MiSeq is a well-established technique for generating mitogenome data from high-quality samples. To this end, a validation was conducted for this mitogenome method processing up to 24 samples simultaneously along with analysis in the CLC Genomics Workbench and utilizing the AQME (AFDIL-QIAGEN mtDNA Expert) tool to generate forensic profiles. This validation followed the Federal Bureau of Investigation's Quality Assurance Standards (QAS) for forensic DNA testing laboratories and the Scientific Working Group on DNA Analysis Methods (SWGDAM) validation guidelines. The evaluation of control DNA, non-probative samples, blank controls, mixtures, and nonhuman samples demonstrated the validity of this method. Specifically, the sensitivity was established at ≥25 pg of nuclear DNA input for accurate mitogenome profile generation. Unreproducible low-level variants were observed in samples with low amplicon yields. Further, variant quality was shown to be a useful metric for identifying sequencing error and crosstalk. Success of this method was demonstrated with a variety of reference sample substrates and extract types. These studies further demonstrate the advantages of using NGS techniques by highlighting the quantitative nature of heteroplasmy detection. The results presented herein from more than 175 samples processed in ten sequencing runs, show this mitogenome sequencing method and analysis strategy to be valid for the generation of reference data.