The case for the continuing use of the revised Cambridge Reference Sequence (rCRS) and the standardization of notation in human mitochondrial DNA studies

The discovery of a ten-generation m.C1494T pedigree in the east of England with probable links to King Richard III

This paper reports the discovery of a m.C1494T pedigree in the east of England made during a search for matrilineal relations of King Richard III. The mitochondrial DNA variant m.C1494T has been associated with aminoglycoside-induced deafness. This variant is very uncommon. although pedigrees with this variant have previously been found in China and Spain. The members of the newly identified pedigree all belong to the mitochondrial haplogroup J1c2c3, which is also the haplogroup of King Richard III. The presence of a few people in the USA from the same haplogroup has previously been noted, and it is now known that one of the people can show his descent from a couple who lived in Nottinghamshire, England, in the late 1700's. The mitochondrial DNA sequence of this man, at present living in the USA, and of his 4th cousin, twice removed, living in Lincoln, England, has shown they belong to haplogroup J1c2c3 and both have the variant m.C1494T; thereby, allowing the production of a multi-generational pedigree originating in the east of England. Fortunately, deafness has not been found in any living member of this large pedigree. It was also noted that the link to the family of King Richard III has not been firmly defined; however the circumstantial evidence is strong as many of his family members lived in this part of England.

Mitochondrial DNA: Inherent Complexities Relevant to Genetic Analyses

Mitochondrial DNA (mtDNA) exhibits distinct characteristics distinguishing it from the nuclear genome, necessitating specific analytical methods in genetic studies. This comprehensive review explores the complex role of mtDNA in a variety of genetic studies, including genome-wide, epigenome-wide, and phenome-wide association studies, with a focus on its implications for human traits and diseases. Here, we discuss the structure and gene-encoding properties of mtDNA, along with the influence of environmental factors and epigenetic modifications on its function and variability. Particularly significant are the challenges posed by mtDNA's high mutation rate, heteroplasmy, and copy number variations, and their impact on disease susceptibility and population genetic analyses. The review also highlights recent advances in methodological approaches that enhance our understanding of mtDNA associations, advocating for refined genetic research techniques that accommodate its complexities. By providing a comprehensive overview of the intricacies of mtDNA, this paper underscores the need for an integrated approach to genetic studies that considers the unique properties of mitochondrial genetics. Our findings aim to inform future research and encourage the development of innovative methodologies to better interpret the broad implications of mtDNA in human health and disease.

Clinical and genetic analysis of essential hypertension with mitochondrial tRNAMet 4435A>G and YARS2 mutation

OBJECTIVES

To investigate the role of m.4435A>G and YARS2 c.572G>T (p.G191V) mutations in the development of essential hypertension.

METHODS

A hypertensive patient with m.4435A>G and YARS2 p.G191V mutations was identified from previously collected mitochondrial genome and exon sequencing data. Clinical data were collected, and a molecular genetic study was conducted in the proband and his family members. Peripheral venous blood was collected, and immortalized lymphocyte lines constructed. The mitochondrial transfer RNA (tRNA), mitochondrial protein, adenosine triphosphate (ATP), mitochondrial membrane potential (MMP), and reactive oxygen species (ROS) in the constructed lymphocyte cell lines were measured.

RESULTS

Mitochondrial genome sequencing showed that all maternal members carried a highly conserved m.4435A>G mutation. The m.4435A>G mutation might affect the secondary structure and folding free energy of mitochondrial tRNA and change its stability, which may influence the anticodon ring structure. Compared with the control group, the cell lines carrying m.4435A>G and YARS2 p.G191V mutations had decreased mitochondrial tRNA homeostasis, mitochondrial protein expression, ATP production and MMP levels, as well as increased ROS levels (all P<0.05).

CONCLUSIONS

The YARS2 p.G191V mutation aggravates the changes in mitochondrial translation and mitochondrial function caused by m.4435A>G through affecting the steady-state level of mitochondrial tRNA and further leads to cell dysfunction, indicating that YARS2 p.G191V and m.4435A>G mutations have a synergistic effect in this family and jointly participate in the occurrence and development of essential hypertension.

CD4+ T cell mitochondrial genotype in Multiple Sclerosis: a cross-sectional and longitudinal analysis

Multiple Sclerosis (MS) is a chronic autoimmune demyelinating disease of the central nervous system (CNS), with a largely unknown etiology, where mitochondrial dysfunction likely contributes to neuroaxonal loss and brain atrophy. Mirroring the CNS, peripheral immune cells from patients with MS, particularly CD4+ T cells, show inappropriate mitochondrial phenotypes and/or oxidative phosphorylation (OxPhos) insufficiency, with a still unknown contribution of mitochondrial DNA (mtDNA). We hypothesized that mitochondrial genotype in CD4+ T cells might influence MS disease activity and progression. Thus, we performed a retrospective cross-sectional and longitudinal study on patients with a recent diagnosis of either Clinically Isolated Syndrome (CIS) or Relapsing-Remitting MS (RRMS) at two timepoints: 6 months (VIS1) and 36 months (VIS2) after disease onset. Our primary outcomes were the differences in mtDNA extracted from CD4+ T cells between: (I) patients with CIS/RRMS (PwMS) at VIS1 and age- and sex-matched healthy controls (HC), in the cross-sectional analysis, and (II) different diagnostic evolutions in PwMS from VIS1 to VIS2, in the longitudinal analysis. We successfully performed mtDNA whole genome sequencing (mean coverage: 2055.77 reads/base pair) in 183 samples (61 triplets). Nonetheless, mitochondrial genotype was not associated with a diagnosis of CIS/RRMS, nor with longitudinal diagnostic evolution.

PCR in Forensic Science: A Critical Review

The polymerase chain reaction (PCR) has played a fundamental role in our understanding of the world, and has applications across a broad range of disciplines. The introduction of PCR into forensic science marked the beginning of a new era of DNA profiling. This era has pushed PCR to its limits and allowed genetic data to be generated from trace DNA. Trace samples contain very small amounts of degraded DNA associated with inhibitory compounds and ions. Despite significant development in the PCR process since it was first introduced, the challenges of profiling inhibited and degraded samples remain. This review examines the evolution of the PCR from its inception in the 1980s, through to its current application in forensic science. The driving factors behind PCR evolution for DNA profiling are discussed along with a critical comparison of cycling conditions used in commercial PCR kits. Newer PCR methods that are currently used in forensic practice and beyond are examined, and possible future directions of PCR for DNA profiling are evaluated.

Comprehensive analysis of mitochondrial DNA variants, mitochondrial DNA copy number and oxidative damage in psoriatic arthritis

Growing evidence suggests that abnormalities in mitochondrial DNA (mtDNA) are involved in the pathogenesis of various inflammatory and immuno-mediated diseases. The present study analysed the entire mitochondrial genome by next-generation sequencing (NGS) in 23 patients with psoriatic arthritis (PsA) and 20 healthy controls to identify PsA-related variants. Changes in mtDNA copy number (mtDNAcn) were also evaluated by quantitative polymerase chain reaction (qPCR) and mtDNA oxidative damage was measured using an 8-hydroxy-2'-deoxyguanosine assay. NGS analysis revealed a total of 435 variants including 187 in patients with PsA only and 122 in controls only. Additionally, 126 common variants were found, of which 2 variants differed significantly in their frequencies among patients and controls (P<0.05), and may be associated with susceptibility to PsA. A total of 33 missense variants in mtDNA-encoded genes for complexes I, III, IV and V were identified only in patients with PsA. Of them, 25 variants were predicted to be deleterious by affecting the functions and structures of encoded proteins, and 13 variants were predicted to affect protein's stability. mtDNAcn analysis revealed decreased mtDNA content in patients with PsA compared with controls (P=0.0001) but the decrease in mtDNAcn was not correlated with patients' age or inflammatory biomarkers (P>0.05). Moreover, a higher level of oxidative damage was observed in patients with PsA compared with controls (P=0.03). The results of the present comprehensive analysis of mtDNA in PsA revealed that certain mtDNA variants may be implicated in the predisposition/pathogenesis of PsA, highlighting the importance of NGS in the identification of mtDNA variants in PsA. The current results also demonstrated that decreased mtDNAcn in PsA may be a consequence of increased oxidative stress. These data provide valuable insights into the contribution of mtDNA defects to the pathogenesis of PsA. Additional studies in larger cohorts are needed to elucidate the role of mtDNA defects in PsA.

A systematic comparison of human mitochondrial genome assembly tools

BACKGROUND

Mitochondria are the cell organelles that produce most of the chemical energy required to power the cell's biochemical reactions. Despite being a part of a eukaryotic host cell, the mitochondria contain a separate genome whose origin is linked with the endosymbiosis of a prokaryotic cell by the host cell and encode independent genomic information throughout their genomes. Mitochondrial genomes accommodate essential genes and are regularly utilized in biotechnology and phylogenetics. Various assemblers capable of generating complete mitochondrial genomes are being continuously developed. These tools often use whole-genome sequencing data as an input containing reads from the mitochondrial genome. Till now, no published work has explored the systematic comparison of all the available tools for assembling human mitochondrial genomes using short-read sequencing data. This evaluation is required to identify the best tool that can be well-optimized for small-scale projects or even national-level research.

RESULTS

In this study, we have tested the mitochondrial genome assemblers for both simulated datasets and whole genome sequencing (WGS) datasets of humans. For the highest computational setting of 16 computational threads with the simulated dataset having 1000X read depth, MitoFlex took the least execution time of 69 s, and IOGA took the longest execution time of 1278 s. NOVOPlasty utilized the least computational memory of approximately 0.098 GB for the same setting, whereas IOGA utilized the highest computational memory of 11.858 GB. In the case of WGS datasets for humans, GetOrganelle and MitoFlex performed the best in capturing the SNPs information with a mean F1-score of 0.919 at the sequencing depth of 10X. MToolBox and NOVOPlasty performed consistently across all sequencing depths with a mean F1 score of 0.897 and 0.890, respectively.

CONCLUSIONS

Based on the overall performance metrics and consistency in assembly quality for all sequencing data, MToolBox performed the best. However, NOVOPlasty was the second fastest tool in execution time despite being single-threaded, and it utilized the least computational resources among all the assemblers when tested on simulated datasets. Therefore, NOVOPlasty may be more practical when there is a significant sample size and a lack of computational resources. Besides, as long-read sequencing gains popularity, mitochondrial genome assemblers must be developed to use long-read sequencing data.

Possible Role of Mitochondrial Transfer RNA Gene 5816 A > G Genetic Polymorphism (m.5816A > G) in a 3-Year-Old Child with Dystonia: Report of a Case

Background  Mutations in the mitochondrial transfer RNA (mt-tRNA) gene are a hotspot for mitochondrial DNA (mtDNA) mutations and are most common in mitochondrial diseases. Methods  We identified the mt-tRNA gene 5816 A > G (m.5816 A > G) mutation in a 3-year-old child with dystonia who died. We performed clinical evaluation, genetic analysis, and biochemical investigation with mitochondrial function testing. Results  Our patient was found to have dystonia with hyperlactatemia. Electroencephalogram findings were abnormal in children with numerous multifocal spikes, multispike, spikes and slow waves, slow waves and low amplitude fast waves, more pronounced in the occipital region bilaterally, and occurring continuously during sleep. One year later, the preexisting patient had seizures lasting 1 to 2 hours and subsequently died. mtDNA sequencing revealed that the proband, her mother, and her grandmother all carried the m.5816A > G mutation. Oxygen consumption rate (OCR) assays revealed that the proband's basal resting OCR, adenosine triphosphate production, proton leak, maximal respiration, and spare capacity OCR were all significantly lower compared with healthy children of the same age. Conclusion  The present case demonstrates a childhood dystonia caused by a mt-tRNA gene 5816 A > G mutation, which has never been reported before. Our findings provide valuable new insights into the pathogenic mechanism and function of the m.5816A > G mutation.

A Customized Human Mitochondrial DNA Database (hMITO DB v1.0) for Rapid Sequence Analysis, Haplotyping and Geo-Mapping

The field of mitochondrial genomics has advanced rapidly and has revolutionized disciplines such as molecular anthropology, population genetics, and medical genetics/oncogenetics. However, mtDNA next-generation sequencing (NGS) analysis for matrilineal haplotyping and phylogeographic inference remains hindered by the lack of a consolidated mitogenome database and an efficient bioinformatics pipeline. To address this, we developed a customized human mitogenome database (hMITO DB) embedded in a CLC Genomics workflow for read mapping, variant analysis, haplotyping, and geo-mapping. The database was constructed from 4286 mitogenomes. The macro-haplogroup (A to Z) distribution and representative phylogenetic tree were found to be consistent with published literature. The hMITO DB automated workflow was tested using mtDNA-NGS sequences derived from Pap smears and cervical cancer cell lines. The auto-generated read mapping, variants track, and table of haplotypes and geo-origins were completed in 15 min for 47 samples. The mtDNA workflow proved to be a rapid, efficient, and accurate means of sequence analysis for translational mitogenomics.

Clinical and genetic analysis of essential hypertension with CYB gene m.15024G>A mutation

OBJECTIVES

To explore the role of mitochondrial CYB 15024G>A mutation in the development of essential hypertension.

METHODS

Mitochondrial genome sequences of hypertensive patients were obtained from previous studies. Clinical and genetic data of a hypertensive patient with mitochondrial CYB 15024G>A mutation and its pedigree were analyzed. Lymphocytes derived from patient and family members were transformed into immortalized lymphoblastoid cell lines, and the levels of adenosine triphosphate (ATP), mitochondrial membrane potential and intracellular reactive oxygen species (ROS) were detected.

RESULTS

The penetrance of this essential hypertension family was 42.9%, and the age of onset was 46-68 years old. Mitochondrial genome sequencing results showed that all maternal members carried a highly conserved mitochondrial CYB 15024G>A mutation. This mutation could affect the free energy of mitochondrial CYB for secondary and tertiary structure and protein folding, thereby changing its structural stability and the structure of the electron transfer function area around the mutation site. Compared with the control, the cell line carrying the mitochondrial CYB 15024G>A mutation showed significantly decreased levels of mitochondrial CYB, ATP and mitochondrial membrane potential, and increased levels of ROS (P<0.01).

CONCLUSIONS

Mitochondrial CYB 15024G>A mutation may affect the structure of respiratory chain subunits and mitochondrial function, leading to cell dysfunction, which suggests that the mutation may play a synergistic role in essential hypertension.

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.

Parallel sequencing of extrachromosomal circular DNAs and transcriptomes in single cancer cells

Extrachromosomal DNAs (ecDNAs) are common in cancer, but many questions about their origin, structural dynamics and impact on intratumor heterogeneity are still unresolved. Here we describe single-cell extrachromosomal circular DNA and transcriptome sequencing (scEC&T-seq), a method for parallel sequencing of circular DNAs and full-length mRNA from single cells. By applying scEC&T-seq to cancer cells, we describe intercellular differences in ecDNA content while investigating their structural heterogeneity and transcriptional impact. Oncogene-containing ecDNAs were clonally present in cancer cells and drove intercellular oncogene expression differences. In contrast, other small circular DNAs were exclusive to individual cells, indicating differences in their selection and propagation. Intercellular differences in ecDNA structure pointed to circular recombination as a mechanism of ecDNA evolution. These results demonstrate scEC&T-seq as an approach to systematically characterize both small and large circular DNA in cancer cells, which will facilitate the analysis of these DNA elements in cancer and beyond.

Forensic evaluation of mitochondrial DNA heteroplasmy in Gujarat population, India

BACKGROUND

Owing to its high copy number and its small size, mtDNA analysis is the most reliable choice when biological materials from crime scenes are degraded or have mixed STR profiles.

AIM

To examine the occurrence of heteroplasmy along with its frequency and pattern in both HV1 and HV2 regions of the mtDNA among unrelated individuals from India.

SUBJECTS AND METHODS

Mitochondrial DNA control region [hypervariable region one (HV1) and hypervariable region two (HV2)] were analysed in blood and buccal tissues of 104 unrelated individuals from the Indian state of Gujarat.

RESULTS

A high frequency of point heteroplasmy (PH) and length heteroplasmy (LH) was revealed. PH was detected in 7.69% of the population, with a higher frequency observed in blood than in buccal samples. However, there were no statistically significant differences in PH between the two tissues (Chi-square = 0.552, p ≥ 0.05). A total of six PH positions were detected: three at HV1, and another three at HV2. The studied population showed 46.15% LH in the HV1 and HV2 regions of both tissues. The LH positions observed in the Gujarat population were the same as those previously reported at HV1 np16184-16193 and HV2 np303-315.

CONCLUSIONS

Our findings suggest that differences in the pattern of heteroplasmy found in different tissues can complicate the forensic analysis, on the other hand, the probability of a match between the questioned and reference samples increases when the heteroplasmy is identical in both tissues. Variability of PH among persons and even within tissues recommends analysing multiple tissue samples before drawing a conclusion in forensic mtDNA analyses.

Mitochondrial DNA variation across 56,434 individuals in gnomAD

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.

Mitochondrial Disorders in Alzheimer's Disease

Alzheimer's disease is the most common age-related neurodegenerative disease. Understanding of its etiology and pathogenesis is constantly expanding. Thus, the increasing attention of researchers is directed to the study of the role of mitochondrial disorders. In addition, in recent years, the concept of Alzheimer's disease as a stress-induced disease has begun to form more and more actively. The stress-induced damage to the neuronal system can trigger a vicious circle of pathological processes, among which mitochondrial dysfunctions have a significant place, since mitochondria represent a substantial component in the anti-stress activity of the cell. The study of mitochondrial disorders in Alzheimer's disease is relevant for at least two reasons: first, as important pathogenetic component in this disease; second, due to vital role of mitochondria in formation of the body resistance to various conditions, including stressful ones, throughout the life. This literature review analyzes the results of a number of recent studies assessing potential significance of the mitochondrial disorders in Alzheimer's disease. The probable mechanisms of mitochondrial disorders associated with the development of this disease are considered: bioenergetic dysfunctions, changes in mitochondrial DNA (including assessment of the significance of its haplogroup features), disorders in the dynamics of these organelles, oxidative damage to calcium channels, damage to MAM complexes (membranes associated with mitochondria; mitochondria-associated membranes), disruptions of the mitochondrial quality control system, mitochondrial permeability, etc. The issues of the "primary" or "secondary" mitochondrial damage in Alzheimer's disease are discussed. Potentials for the development of new methods for diagnosis and therapy of mitochondrial disorders in Alzheimer's disease are considered.

A homozygous nonsense mutation in DCBLD2 is a candidate cause of developmental delay, dysmorphic features and restrictive cardiomyopathy

DCBLD2 encodes discodin, CUB and LCCL domain-containing protein 2, a type-I transmembrane receptor that is involved in intracellular receptor signalling pathways and the regulation of cell growth. In this report, we describe a 5-year-old female who presented severe clinical features, including restrictive cardiomyopathy, developmental delay, spasticity and dysmorphic features. Trio-whole-exome sequencing and segregation analysis were performed to identify the genetic cause of the disease within the family. A novel homozygous nonsense variant in the DCBLD2 gene (c.80G > A, p.W27*) was identified as the most likely cause of the patient's phenotype. This nonsense variant falls in the extracellular N-terminus of DCBLD2 and thus might affect proper protein function of the transmembrane receptor. A number of in vitro investigations were performed on the proband's skin fibroblasts compared to normal fibroblasts, which allowed a comprehensive assessment resulting in the functional characterization of the identified DCBLD2 nonsense variant in different cellular processes. Our data propose a significant association between the identified variant and the observed reduction in cell proliferation, cell cycle progression, intracellular ROS, and Ca2 + levels, which would likely explain the phenotypic presentation of the patient as associated with lethal restrictive cardiomyopathy.

A history of you, me, and humanity: mitochondrial DNA in anthropological research

Within genetic anthropology, mitochondrial DNA (mtDNA) has garnered a prominent if not enduring place within the anthropological toolkit. MtDNA has provided new and innovative perspectives on the emergence and dispersal of our species, interactions with extinct human species, and illuminated relationships between human groups. In this paper, I provide a brief overview of the major findings ascertained from mtDNA about human origins, human dispersal across the globe, interactions with other hominin species, and the more recent uses of mtDNA in direct to consumer ancestry tests. Relative to nuclear DNA, mtDNA is a small section of the genome and due to its inheritance pattern provides a limited resolution of population history and an individual's genetic ancestry. Consequently, some scholars dismiss mtDNA as insignificant due to the limited inferences that may be made using the locus. Regardless, mtDNA provides some useful insights to understanding how social, cultural, and environmental factors have shaped patterns of genetic variability. Furthermore, with regard to the experiences of historically marginalized groups, in particular those of African descent throughout the Americas, mtDNA has the potential to fill gaps in knowledge that would otherwise remain unknown. Within anthropological sciences, the value of this locus for understanding human experience is maximized when contextualized with complementary lines of evidence.

Effect of cell microenvironment on the drug sensitivity of hepatocellular cancer cells

The native hepatocellular cancer (HCC) microenvironment is characterized by more hypoxic, hypoglycemic, and acidic conditions than those used in standard cell culture. This study aimed to investigate whether HCC cells cultured in more native conditions have an altered phenotype and drug sensitivity compared to those cultured in standard conditions. Six HCC cell lines were cultured in “standard” (21% O₂, 25 mM glucose) or more “native” (1% O₂, 5 mM glucose, 10 mM lactate) conditions. Cells were assessed for growth rates, cell cycle distribution, relevant metabolite and protein levels, genome-wide gene expression, mitochondrial DNA sequence and sensitivity to relevant drugs. Many differences in cellular and molecular phenotypes and drug sensitivity were observed between the cells. HCC cells cultured in native conditions had slower doubling times, increased HK2 and GLUT, lower PHDA and ATP levels, and mutations in mitochondrial DNA. Thirty-one genes, including the hypoxia-associated NDRG1, were differentially expressed between the cells. HCC patients in The Cancer Genome Atlas (TCGA) with tumors with a high score based on these 31 genes had a poorer prognosis than those with a low score (p = 0.002). From 90 comparisons of drug sensitivity, increased resistance and sensitivity for cells cultured in native conditions was observed in 14 (16%) and 8 (9%) comparisons respectively. In conclusion, cells cultured in more native conditions can have a more glycolytic and aggressive phenotype and varied drug sensitivity to those cultured in standard conditions, and may provide new insights to understanding tumor biology and drug development.

Next-generation sequencing of the whole mitochondrial genome identifies novel and common variants in patients with psoriasis, type 2 diabetes mellitus and psoriasis with comorbid type 2 diabetes mellitus

Recent studies have shown the role of mitochondrial DNA (mtDNA) variants in the pathogenesis of both psoriasis (Ps) and type 2 diabetes (T2D) amongst different ethnicities. However, no studies have investigated the mtDNA variants present in patients with Ps, T2D, and both Ps and T2D (Ps-T2D) in the Arab population. The entire mitochondrial genomes of Kuwaiti subjects with Ps, T2D, Ps-T2D and healthy controls were sequenced using Ion Torrent next-generation sequencing. A total of 36 novel mutations and 51 previously reported mutations were identified in the patient groups that were absent in the controls. Amongst the novel mutations, eight were non-synonymous and exhibited amino acid changes. Of these, two missense mutations (G5262A and A12397G) in the ND genes were detected in the Ps group and a C15735T missense mutation in the CYB gene was detected in Ps-T2D. Other known sequence variations were seen more frequently in all or certain patient groups compared with the controls (P<0.05). Additionally, the A8701G missense mutation in the ATPase 6 gene missense mutation was also observed in a higher frequency in the Ps group compared with the control. The present study is the first to perform a complete mitochondrial genome sequence analysis of Kuwaiti subjects with Ps, T2D and Ps-T2D, and both novel and known mtDNA variants were discovered. The patient-specific novel non-synonymous mutations may be co-responsible in the determination of these diseases. The higher frequency of certain mtDNA variants in the patients compared with the controls may suggest a role in predisposing patients to these diseases. Further functional analyses are required to reveal the role of the identified mutations in these disease conditions.

Mitochondrial Haplogroup Reveals the Genetic Basis of Diabetes Mellitus Type 2 Comorbidity in Psoriasis

OBJECTIVE

Published data show a clear link between psoriasis (Ps) and the increasing prevalence of comorbid conditions, such as diabetes mellitus type 2 (DM2). The role of the mitochondrial genomic haplogroup in the potential coexistence of Ps and DM2 comorbidity is the subject of this study.

MATERIAL AND METHODS

Ninety-eight Kuwaiti individuals were recruited in 4 cohorts (20 healthy controls, 15 with DM2, 34 with Ps, and 29 with Ps and diabetes mellitus). An Ion Torrent S5XL was used to sequence mitochondrial DNA (mtDNA). χ2 test was used to assess differences in the distribution of each haplogroup between cases and controls (p < 0.05). The Bonferroni correction was applied (p < 0.004). The mtDNA haplogroups were analyzed by HaploGrep.

RESULTS

Haplogroups R0, U, J, T, N, L3, M, H, X, HV, R, and K were detected in the studied population. Haplogroup M had a high risk for Ps (odds ratio (OR) 4.0, p = 0.003). Haplogroup R0 and J had decreased the risk of DM2 (OR 0.28, p = 0.007).

CONCLUSION

Our results indicated that mtDNA haplogroups have a potential contribution to the pathogenesis of Ps and DM2 comorbidity. We show for the first time that the comorbidity of diabetes in Ps may be related to mitochondrial dysfunction.

Potentially adaptive SARS-CoV-2 mutations discovered with novel spatiotemporal and explainable AI models

BACKGROUND

A mechanistic understanding of the spread of SARS-CoV-2 and diligent tracking of ongoing mutagenesis are of key importance to plan robust strategies for confining its transmission. Large numbers of available sequences and their dates of transmission provide an unprecedented opportunity to analyze evolutionary adaptation in novel ways. Addition of high-resolution structural information can reveal the functional basis of these processes at the molecular level. Integrated systems biology-directed analyses of these data layers afford valuable insights to build a global understanding of the COVID-19 pandemic.

RESULTS

Here we identify globally distributed haplotypes from 15,789 SARS-CoV-2 genomes and model their success based on their duration, dispersal, and frequency in the host population. Our models identify mutations that are likely compensatory adaptive changes that allowed for rapid expansion of the virus. Functional predictions from structural analyses indicate that, contrary to previous reports, the Asp614Gly mutation in the spike glycoprotein (S) likely reduced transmission and the subsequent Pro323Leu mutation in the RNA-dependent RNA polymerase led to the precipitous spread of the virus. Our model also suggests that two mutations in the nsp13 helicase allowed for the adaptation of the virus to the Pacific Northwest of the USA. Finally, our explainable artificial intelligence algorithm identified a mutational hotspot in the sequence of S that also displays a signature of positive selection and may have implications for tissue or cell-specific expression of the virus.

CONCLUSIONS

These results provide valuable insights for the development of drugs and surveillance strategies to combat the current and future pandemics.

Phylogenomic study and classification of mitochondrial DNA through virtual genomic fingerprints

In the present study, we evaluated the ability of the Virtual Analysis Method for Phylogenomic fingerprint Estimation (VAMPhyRE) toolkit to classify human mitochondrial DNA (mtDNA) haplogroups. In total, 357 random mtDNA sequences were obtained from different haplogroups, based on the classification of PhyloTree. Additionally, we included a control group of five sequences (Pan paniscus, Pan troglodytes, Homo sapiens neanderthalensis, Yoruba15, and the revised Cambridge reference sequence). VAMPhyRE employs a virtual hybridization technique, using probes that specifically bind to their complementary sequences in the genome. We used 65,536 probes of 8 nucleotides to identify potential sites where hybridization occurs between the mtDNA and the specific probe, forming different heteroduplexes and thus, creating a unique and specific genomic fingerprint for each sequence. Genomic fingerprints were compared, and a table of distances was calculated to obtain a mitochondrial phylogenomic tree with the macrohaplogroups, L, N, M, and R, and their corresponding haplogroups, according to universal nomenclature. The results obtained suggest an accuracy of 97.25% for the distribution of the 357 mtDNA sequences in the four macrohaplogroups and their corresponding haplogroups when compared with other mtDNA classification tools that require reference sequences and do not offer an analysis based on an evolutionary approach. These data are available online at http://biomedbiotec.encb.ipn.mx/VAMPhyRE/.

Specifications of the ACMG/AMP standards and guidelines for mitochondrial DNA variant interpretation

Mitochondrial DNA (mtDNA) variant pathogenicity interpretation has special considerations given unique features of the mtDNA genome, including maternal inheritance, variant heteroplasmy, threshold effect, absence of splicing, and contextual effects of haplogroups. Currently, there are insufficient standardized criteria for mtDNA variant assessment, which leads to inconsistencies in clinical variant pathogenicity reporting. An international working group of mtDNA experts was assembled within the Mitochondrial Disease Sequence Data Resource Consortium and obtained Expert Panel status from ClinGen. This group reviewed the 2015 American College of Medical Genetics and Association of Molecular Pathology standards and guidelines that are widely used for clinical interpretation of DNA sequence variants and provided further specifications for additional and specific guidance related to mtDNA variant classification. These Expert Panel consensus specifications allow for consistent consideration of the unique aspects of the mtDNA genome that directly influence variant assessment, including addressing mtDNA genome composition and structure, haplogroups and phylogeny, maternal inheritance, heteroplasmy, and functional analyses unique to mtDNA, as well as specifications for utilization of mtDNA genomic databases and computational algorithms.

Clinical use of NGS data from the targeted gene panel for mitochondrial diseases screening

BACKGROUND AND OBJECTIVE

Mitochondrial diseases are a frequent cause of inherited genetic disorders caused by mutations in both the mitochondrial and nuclear human genome. The new technique of high-throughput sequencing, which is used more and more frequently around the world, is most often focused on nuclear DNA. In some cases, such data after proper IT processing could also allow to determine alterations in mtDNA genome. In our work, we want to verify that off-target reads from targeted gene panels are sufficient data to determine pathogenic variants in the mitochondrial genome.

METHODS

We analyzed 50 patients who underwent routine diagnostics with the Illumina's TruSight One Sequencing Panel. In the entire bioinformatic analysis process, we have used only free, user-friendly and generally available online tools that do not require specialized IT knowledge.

RESULTS

Most of the data obtained were suitable for determining the presence of homoplasmic variants in mtDNA; 84% of the data met the minimum 20-fold coverage requirement as defined in the scientific literature for clinical data. We managed to identify 16 pathogenic variants in the examined genetic material (mtDNA) according to the ClinVar database.

CONCLUSIONS

In conclusion, we have outlined that off-target reads from targeted gene panel (TruSight One Sequencing Panel) may also be suitable for determining potentially pathogenic homoplasmic variants in mtDNA. We also described a simple pipeline based only on free tools available online. Introducing such a pipeline into a standard procedure of clinical units which carry out such research undoubtedly can extend the diagnostic potential by information about mtDNA, especially when it is based on purely free tools that do not require specialized bioinformatic knowledge.

A novel droplet digital PCR human mtDNA assay for fecal source tracking

Human mitochondrial DNA provides a promising target for fecal source tracking because it is unique and intrinsic to humans. We developed a TaqMan chemistry assay, hCYTB484, targeting the cytochrome b gene of the human mitochondrial genome on a droplet digital PCR (ddPCR) platform and compared the performance of hCYTB484 with the HF183/BacR287 assay, a widely used assay targeting human-associated Bacteroides. For both assays, we defined the analytical limit of detection and analytical lower limit of quantification using frequency of detection and imprecision goals, respectively. We then established these analytical limits using empirical ddPCR data, presenting a novel approach to determining the analytical lower limit of quantification. We evaluated assay sensitivity using individual human feces from US, Bangladesh, and Mozambique and evaluated assay specificity using cow, pig, chicken, and goat samples collected from the US. To compare assay performance across a range of thresholds, we utilized receiver operating characteristic curves. The hCYTB484 marker was detected and quantifiable in 100% of the human feces from the 3 geographical distant regions whereas the HF183/BacR287 marker was detectable and quantifiable in 51% and 31% (respectively) of human feces samples. The hCYTB484 marker also was more specific (97%), having fewer detections in pig, chicken, and goat samples than the HF183/BacR287 marker (80%). The higher performance of the hCYTB484 marker in individual feces from geographically distant regions is desirable in the detection of fecal pollution from sources to which fewer individuals contribute, such as the non-sewered forms of sanitation (e.g. pit latrines and septic tanks) that serve most of Earth's population and carry the highest risk of exposure to fecal-oral pathogens.

What can a comparative genomics approach tell us about the pathogenicity of mtDNA mutations in human populations?

Mitochondrial disorders are heterogeneous, showing variable presentation and penetrance. Over the last three decades, our ability to recognize mitochondrial patients and diagnose these mutations, linking genotype to phenotype, has greatly improved. However, it has become increasingly clear that these strides in diagnostics have not benefited all population groups. Recent studies have demonstrated that patients from genetically understudied populations, in particular those of black African heritage, are less likely to receive a diagnosis of mtDNA disease. It has been suggested that haplogroup context might influence the presentation and penetrance of mtDNA disease; thus, the spectrum of mutations that are associated with disease in different populations. However, to date there is only one well-established example of such an effect: the increased penetrance of two Leber's hereditary optic neuropathy mutations on a haplogroup J background. This paper conducted the most extensive investigation to date into the importance of haplogroup context on the pathogenicity of mtDNA mutations. We searched for proven human point mutations across 726 multiple sequence alignments derived from 33 non-human species absent of disease. A total of 58 pathogenic point mutations arise in the sequences of these species. We assessed the sequence context and found evidence of population variants that could modulate the phenotypic expression of these point mutations masking the pathogenic effects seen in humans. This supports the theory that sequence context is influential in the presentation of mtDNA disease and has implications for diagnostic practices. We have shown that our current understanding of the pathogenicity of mtDNA point mutations, primarily built on studies of individuals with haplogroups HVUKTJ, will not present a complete picture. This will have the effect of creating a diagnostic inequality, whereby individuals who do not belong to these lineages are less likely to receive a genetic diagnosis.

NUMT Confounding Biases Mitochondrial Heteroplasmy Calls in Favor of the Reference Allele

Homology between mitochondrial DNA (mtDNA) and nuclear DNA of mitochondrial origin (nuMTs) causes confounding when aligning short sequence reads to the reference human genome, as the true sequence origin cannot be determined. Using a systematic in silico approach, we here report the impact of all potential mitochondrial variants on alignment accuracy and variant calling. A total of 49,707 possible mutations were introduced across the 16,569 bp reference mitochondrial genome (16,569 × 3 alternative alleles), one variant at-at-time. The resulting in silico fragmentation and alignment to the entire reference genome (GRCh38) revealed preferential mapping of mutated mitochondrial fragments to nuclear loci, as variants increased loci similarity to nuMTs, for a total of 807, 362, and 41 variants at 333, 144, and 27 positions when using 100, 150, and 300 bp single-end fragments. We subsequently modeled these affected variants at 50% heteroplasmy and carried out variant calling, observing bias in the reported allele frequencies in favor of the reference allele. Four variants (chrM:6023A, chrM:4456T, chrM:5147A, and chrM:7521A) including a possible hypertension factor, chrM:4456T, caused 100% loss of coverage at the mutated position (with all 100 bp single-end fragments aligning to homologous, nuclear positions instead of chrM), rendering these variants undetectable when aligning to the entire reference genome. Furthermore, four mitochondrial variants reported to be pathogenic were found to cause significant loss of coverage and select haplogroup-defining SNPs were shown to exacerbate the loss of coverage caused by surrounding variants. Increased fragment length and use of paired-end reads both improved alignment accuracy.

Frequency and spectrum of MT-TT variants associated with Leber's hereditary optic neuropathy in a Chinese cohort of subjects

Leber’s hereditary optic neuropathy (LHON) is a maternally inherited eye disease. In our previous investigations, we have reported the spectrum and frequency of mitochondrial MT-ND1, MT-ND4 and MT-ND6 gene in Chinese LHON population. This study aimed to assess the molecular epidemiology of MT-TT mutations in Chinese families with LHON. A cohort of 352 Chinese Han probands lacking the known LHON-associated mtDNA mutations and 376 control subjects underwent molecular analysis of mtDNA. All variants were evaluated for evolutionary conservation, structural and functional consequences. Fifteen variants were identified in the MT-TT gene by mitochondrial genome analysis of LHON pedigrees, which was substantially higher than that of individuals from general Chinese populations. The incidences of the two known LHON-associated mutations, m.15927G > A and m.15951A > G, were 2.27% and 1.14%, respectively. Nine putative LHON-associated variants were identified in 20 probands, translated into 2.1% cases of this cohort. Moreover, mtDNAs in 41 probands carrying the MT-TT mutation(s) were widely dispersed among nine Eastern Asian haplogroups. Our results suggest that the MT-TT gene is a mutational hotspot for these 352 Chinese families lacking the known LHON-associated mutations. These data further showed the molecular epidemiology of MT-TT mutations in Chinese Han LHON pedigrees.

A descriptive pilot study of mitochondrial mutations & clinical phenotype in fibromyalgia syndrome

Background & objectives:

Fibromyalgia syndrome (FMS) is one of the most common chronic pain conditions of unknown aetiology. Mitochondrial dysfunction has been reported in FMS with some studies reporting the presence of mitochondrial mutation namely A3243G, which also causes mitochondrial encephalomyopathy, lactic acidosis and stroke-like episodes. This pilot study was conducted to assess this mutation and also detect large deletions in mitochondrial DNA (mtDNA) in patients with FMS.

Methods:

Thirty female patients with FMS participated and 30 matched controls were included. Genomic DNA was subjected to polymerase chain reaction (PCR) amplification using specific primers followed by restriction digestion with ApaI enzyme to detect the specific A3243G mtDNA mutation. Long-range PCR was done in two sets to detect the large deletions in the mtDNA. Biochemical parameters including thyroid-stimulating hormone and vitamin D levels were also looked at.

Results:

None of the patients were found to carry the common mutation or large deletions. Low vitamin D level was a common finding. Hypothyroidism was found in a few patients.

Interpretation & conclusions:

Although the common mutation or large mtDNA deletions were not detected in blood mtDNA in the FMS patients, mutations in the muscle and sequence variation in mtDNA remained a possibility. Future studies in both blood and muscle tissue including mtDNA sequencing are warranted in such patients to determine if a subset of FMS patients have mitochondrial myopathy.

Identification of World War II bone remains found in Ukraine using classical anthropological and mitochondrial DNA results

Gyula Ágner was a Royal Hungarian First Lieutenant (1st Lt.) during the World War II and died at 30 years old due to a mine shrapnel injury on 27 April 1944 in Luczky, Ukraine. In October 2014, the Hungarian Ministry of Defence exhumated the remains then transported them to Budapest in Hungary. Classical anthropological methods were used to determine morphological gender, height and age at death; furthermore, metrical and pathological characters were also analysed. Determination of maternal lineage was the only solution to examine the possible relationship of the bone fragments. Gyula Ágner did not have direct descendants, thus the living niece of the deceased (his sister’s daughter) served as the reference person during the investigations. Hypervariable regions of the mtDNA control region (HV1, HV2 and HV3) were amplified by Qiagen® Multiplex PCR Kit in different monoplex reactions. The results of the anthropological and genetical analysis supported the hypothesis that the bone remains belong to Gyula Ágner.

Molecular markers used in forensic genetics

Forensic genetics is a field that has become subject to increasing interest in recent years. Both the technology and the markers used for forensic purposes have changed since the 1980s. The minisatellite sequences used in the famous Pitchfork case introduced genetics to the forensic sciences. Minisatellite sequences have now been replaced by more sensitive microsatellite markers, which have become the basis for the creation of genetic profile databases. Modern molecular methods also exploit single nucleotide polymorphisms, which are often the only way to identify degraded DNA samples. The same type of variation is taken into consideration in attempting to establish the ethnicity of a perpetrator and to determine phenotypic traits such as the eye or hair colour of the individual who is the source of the genetic material. This paper contains a review of the techniques and molecular markers used in human and animal forensic genetics, and also presents the potential trends in forensic genetics such as phenotyping.

DNA Methylation-a Potential Source of Mitochondria DNA Base Mismatch in the Development of Diabetic Retinopathy

In the development of diabetic retinopathy, retinal mitochondria are dysfunctional, and mitochondrial DNA (mtDNA) is damaged with increased base mismatches and hypermethylated cytosines. DNA methylation is also a potential source of mutation, and in diabetes, the noncoding region, the displacement loop (D-loop), experiences more methylation and base mismatches than other regions of the mtDNA. Our aim was to investigate a possible crosstalk between mtDNA methylation and base mismatches in the development of diabetic retinopathy. The effect of inhibition of Dnmts (by 5-aza-2'-deoxycytidine or Dnmt1-siRNA) on glucose-induced mtDNA base mismatches was investigated in human retinal endothelial cells by surveyor endonuclease digestion and validated by Sanger sequencing. The role of deamination factors on increased base mismatches was determined in the cells genetically modulated for mitochondrial superoxide dismutase (Sod2) or cytidine-deaminase (APOBEC3A). The results were confirmed in an in vivo model using retinal microvasculature from diabetic mice overexpressing Sod2. Inhibition of DNA methylation, or regulation of cytosine deamination, significantly inhibited an increase in base mismatches at the D-loop and prevented mitochondrial dysfunction. Overexpression of Sod2 in mice also prevented diabetes-induced D-loop hypermethylation and increase in base mismatches. The crosstalk between DNA methylation and base mismatches continued even after termination of hyperglycemia, suggesting its role in the metabolic memory phenomenon associated with the progression of diabetic retinopathy. Inhibition of DNA methylation limits the availability of methylated cytosine for deamination, suggesting a crosstalk between DNA methylation and base mismatches. Thus, regulation of DNA methylation, or its deamination, should impede the development of diabetic retinopathy by preventing formation of base mismatches and mitochondrial dysfunction.

Pseudogenization of the Humanin gene is common in the mitochondrial DNA of many vertebrates

In the human the peptide Humanin is produced from the small Humanin gene which is embedded as a gene-within-a-gene in the 16S ribosomal molecule of the mitochondrial DNA (mtDNA). The peptide itself appears to be significant in the prevention of cell death in many tissues and improve cognition in animal models. By using simple data mining techniques, it is possible to show that 99.4% of the human Humanin sequences in the GenBank database are unaffected by mutations. However, in other vertebrates, pseudogenization of the Humanin gene is a common feature; occurring apparently randomly in some species and not others. The persistence, or loss, of a functional Humanin gene may be an important factor in laboratory animals, especially if they are being used as animal models in studies of Alzheimer's disease (AD). The exact reason why Humanin underwent pseudogenization in some vertebrate species during their evolution remains to be determined. This study was originally planned to review the available information about Humanin and it was a surprise to be able to show that pseudogenization has occurred in a gene in the mtDNA and is not restricted solely to chromosomal genes.

Acquired resistance to PI3K/mTOR inhibition is associated with mitochondrial DNA mutation and glycolysis

Acquired resistance (AQR) to drug treatment occurs frequently in cancer patients and remains an impediment to successful therapy. The aim of this study was to gain insight into how AQR arises following the application of PI3K/mTOR inhibitors. H1975 lung cancer cells with EGFR T790M mutations that confer resistance to EGFR inhibitors underwent prolonged treatment with the PI3K/mTOR inhibitor, BEZ235. Monoclonal cells with stable and increased resistance to BEZ235 were obtained after 8 months treatment. These AQR clones showed class-specific resistance to PI3K/mTOR inhibitors, reduced G1 cell cycle arrest and impedance of migration following PI3K/mTOR inhibition, reduced PTEN expression and increased Akt and S6RP phosphorylation. Transcriptome analysis revealed the AQR clones had increased expression of the metabolite transporters SLC16A9 and SLC16A7, suggestive of altered cell metabolism. Subsequent experiments revealed that AQR clones possess features consistent with elevated glycolysis, including increased levels of glucose, lactate, glutamine, glucose dependence, GLUT1 expression, and rates of post-glucose extracellular acidification, and decreased levels of reactive oxygen species and rates of oxygen consumption. Combination treatment of BEZ235 with the glycolysis inhibitor 3-bromopyruvate was synergistic in AQR clones, but only additive in parental cells. DNA sequencing revealed the presence of a mitochondrial DNA (mtDNA) MT-C01 variant in AQR but not parental cells. Depletion of mitochondrial DNA in parental cells induced resistance to BEZ235 and other PI3K/mTOR inhibitors, and was accompanied by increased glycolysis. The results of this study provide the first evidence that a metabolic switch associated with mtDNA mutation can be an underlying mechanism for AQR.

The Landscape of mtDNA Modifications in Cancer: A Tale of Two Cities

Mitochondria from normal and cancerous cells represent a tale of two cities, wherein both execute similar processes but with different cellular and molecular effects. Given the number of reviews currently available which describe the functional implications of mitochondrial mutations in cancer, this article focuses on documenting current knowledge in the abundance and distribution of somatic mitochondrial mutations, followed by elucidation of processes which affect the fate of mutations in cancer cells. The conclusion includes an overview of translational implications for mtDNA mutations, as well as recommendations for future research uniting mitochondrial variants and tumorigenesis.

Sequence variations of mitochondrial DNA D-loop region in patients with acute myeloid leukemia

The aim of the present study was to explore variations of the displacement (D)-loop region in patients with acute myeloid leukemia (AML) and their possible associations with AML pathogenesis. Blood or bone marrow samples from 216 patients with AML (158 AML patients in the first stage, and 58 more patients with AML-M3 for further verification), and 146 healthy controls were collected. Sanger sequencing was performed for the D-loop region ranging between nucleotide (nt)15811 and nt 775. With the exception of mitochondrial microsatellite instability (mtMSI) variations, a total of 2,630 variations in 232 loci were identified with similar variation rates/person in patients with AML and controls when compared with the revised Cambridge reference sequence (8.54±2.14 vs. 8.77±2.15; P=0.366). A positive association between AML and variation-T152C was identified, which occurred more frequently in patients with AML compared with in controls [26.6 vs. 17.1%; P=0.048; odds ratio (OR), 1.752; 95% confidence interval (CI), 1.004-3.058]. Furthermore, T152C was identified to be associated with promyelocytic leukemia-retinoic acid receptor α(PML-RARα) and French-American-British AML subtypes, with a tendency to occur in patients with AML-M3. The AML-M3 sample size was extended by 58 cases, and it was identified that the T152C variation rate was significantly higher in patients with AML-M3 compared with that of controls (41.0 vs. 17.1%; P<0.001; OR, 3.228; 95% CI, 1.714-6.079). However, no association was identified between the T152C variation and clinical characteristics, or chemotherapy response in patients with AML-M3. In addition, the mtMSIs, including D310, mt514-523 (CA)_n and T16189C, demonstrated no association with AML risk. Together, the results of the present study suggest that the mitochondrial DNA D-loop region is high variable, and that T152C is associated with AML risk, particularly regarding the M3 subtype. T152C mayparticipate in AML pathogenesis and may be a diagnostic biomarker; however further studies with larger sample sizes are required in order to verify its value.

Mitochondrial DNA sequencing and large-scale genotyping identifies MT-ND4 gene mutation m.11696G>A associated with idiopathic oligoasthenospermia

Genetic variants of mitochondrial DNA (mtDNA) were implicated to be associated with male infertility. Our previous whole mitochondrial genome sequencing and association study has identified two susceptibility mtDNA variants for oligoasthenospermia in Han Chinese men. In this study, we tested promising associations in an extended validation using 670 idiopathic oligoasthenospermia cases and 793 healthy controls to identify additional risk variants. We found that the genetic variant of m.11696G>A showed significantly higher frequency in the case group than that in the control group (odds ratio (OR) 2.21, 95% CI 1.21-4.04) (P=7.90×10⁻³). To elucidate the exact role of the genetic variants in spermatogenesis, two main sperm parameters (sperm count and motility) were taken into account. We found that m.11696G>A was associated with low sperm motility, with the OR of 2.38 (95 % CI 1.27-4.46) (P =5.22×10⁻³). These results advance our understanding of the genetic susceptibility to oligoasthenospermia and more functional studies are needed to provide insights into its pathogenic mechanism.

Scanning the mitochondrial genome for mutations by cycling temperature capillary electrophoresis

To bypass possible nuclear contamination and to exclusively amplify DNA from the mitochondrion, a set of 23 primers was selected. On the mitochondrial DNA selection fragments, a second set of fragments was used to amplify and identify mutant fractions with a detection limit of 1% . This mutation scanning method analyzed 76% of the mitochondrial genome and was used to examine 94 tumours from different tissues of origin. In all, 87 tumours had one or more mutations, leaving seven samples without observed mutations. Sanger sequencing verified samples carrying mutations with a mutant fraction exceeding 30%. The generated data validate that several regions of the mitochondrial DNA have more mutations than others.

HaploGrep 2: mitochondrial haplogroup classification in the era of high-throughput sequencing

Mitochondrial DNA (mtDNA) profiles can be classified into phylogenetic clusters (haplogroups), which is of great relevance for evolutionary, forensic and medical genetics. With the extensive growth of the underlying phylogenetic tree summarizing the published mtDNA sequences, the manual process of haplogroup classification would be too time-consuming. The previously published classification tool HaploGrep provided an automatic way to address this issue. Here, we present the completely updated version HaploGrep 2 offering several advanced features, including a generic rule-based system for immediate quality control (QC). This allows detecting artificial recombinants and missing variants as well as annotating rare and phantom mutations. Furthermore, the handling of high-throughput data in form of VCF files is now directly supported. For data output, several graphical reports are generated in real time, such as a multiple sequence alignment format, a VCF format and extended haplogroup QC reports, all viewable directly within the application. In addition, HaploGrep 2 generates a publication-ready phylogenetic tree of all input samples encoded relative to the revised Cambridge Reference Sequence. Finally, new distance measures and optimizations of the algorithm increase accuracy and speed-up the application. HaploGrep 2 can be accessed freely and without any registration at http://haplogrep.uibk.ac.at.

From NGS assembly challenges to instability of fungal mitochondrial genomes: A case study in genome complexity

The presence of repetitive or non-unique DNA persisting over sizable regions of a eukaryotic genome can hinder the genome's successful de novo assembly from short reads: ambiguities in assigning genome locations to the non-unique subsequences can result in premature termination of contigs and thus overfragmented assemblies. Fungal mitochondrial (mtDNA) genomes are compact (typically less than 100 kb), yet often contain short non-unique sequences that can be shown to impede their successful de novo assembly in silico. Such repeats can also confuse processes in the cell in vivo. A well-studied example is ectopic (out-of-register, illegitimate) recombination associated with repeat pairs, which can lead to deletion of functionally important genes that are located between the repeats. Repeats that remain conserved over micro- or macroevolutionary timescales despite such risks may indicate functionally or structurally (e.g., for replication) important regions. This principle could form the basis of a mining strategy for accelerating discovery of function in genome sequences. We present here our screening of a sample of 11 fully sequenced fungal mitochondrial genomes by observing where exact k-mer repeats occurred several times; initial analyses motivated us to focus on 17-mers occurring more than three times. Based on the diverse repeats we observe, we propose that such screening may serve as an efficient expedient for gaining a rapid but representative first insight into the repeat landscapes of sparsely characterized mitochondrial chromosomes. Our matching of the flagged repeats to previously reported regions of interest supports the idea that systems of persisting, non-trivial repeats in genomes can often highlight features meriting further attention.

Mitochondrial DNA variants correlate with symptoms in myalgic encephalomyelitis/chronic fatigue syndrome

Background

Mitochondrial dysfunction has been hypothesized to occur in Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS), a disease characterized by fatigue, cognitive difficulties, pain, malaise, and exercise intolerance. We investigated whether haplogroup, single nucleotide polymorphisms (SNPs), or heteroplasmy of mitochondrial DNA (mtDNA) were associated with health status and/or symptoms.

Methods

Illumina sequencing of PCR-amplified mtDNA was performed to analyze sequence and extent of heteroplasmy of mtDNAs of 193 cases and 196 age- and gender-matched controls from DNA samples collected by the Chronic Fatigue Initiative. Association testing was carried out to examine possible correlations of mitochondrial sequences with case/control status and symptom constellation and severity as reported by subjects on Short Form-36 and DePaul Symptom Questionnaires.

Results

No ME/CFS subject exhibited known disease-causing mtDNA mutations. Extent of heteroplasmy was low in all subjects. Although no association between mtDNA SNPs and ME/CFS vs. healthy status was observed, haplogroups J, U and H as well as eight SNPs in ME/CFS cases were significantly associated with individual symptoms, symptom clusters, or symptom severity.

Conclusions

Analysis of mitochondrial genomes in ME/CFS cases indicates that individuals of a certain haplogroup or carrying specific SNPs are more likely to exhibit certain neurological, inflammatory, and/or gastrointestinal symptoms. No increase in susceptibility to ME/CFS of individuals carrying particular mitochondrial genomes or SNPs was observed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12967-016-0771-6) contains supplementary material, which is available to authorized users.

An original phylogenetic approach identified mitochondrial haplogroup T1a1 as inversely associated with breast cancer risk in BRCA2 mutation carriers

INTRODUCTION

Individuals carrying pathogenic mutations in the BRCA1 and BRCA2 genes have a high lifetime risk of breast cancer. BRCA1 and BRCA2 are involved in DNA double-strand break repair, DNA alterations that can be caused by exposure to reactive oxygen species, a main source of which are mitochondria. Mitochondrial genome variations affect electron transport chain efficiency and reactive oxygen species production. Individuals with different mitochondrial haplogroups differ in their metabolism and sensitivity to oxidative stress. Variability in mitochondrial genetic background can alter reactive oxygen species production, leading to cancer risk. In the present study, we tested the hypothesis that mitochondrial haplogroups modify breast cancer risk in BRCA1/2 mutation carriers.

METHODS

We genotyped 22,214 (11,421 affected, 10,793 unaffected) mutation carriers belonging to the Consortium of Investigators of Modifiers of BRCA1/2 for 129 mitochondrial polymorphisms using the iCOGS array. Haplogroup inference and association detection were performed using a phylogenetic approach. ALTree was applied to explore the reference mitochondrial evolutionary tree and detect subclades enriched in affected or unaffected individuals.

RESULTS

We discovered that subclade T1a1 was depleted in affected BRCA2 mutation carriers compared with the rest of clade T (hazard ratio (HR) = 0.55; 95% confidence interval (CI), 0.34 to 0.88; P = 0.01). Compared with the most frequent haplogroup in the general population (that is, H and T clades), the T1a1 haplogroup has a HR of 0.62 (95% CI, 0.40 to 0.95; P = 0.03). We also identified three potential susceptibility loci, including G13708A/rs28359178, which has demonstrated an inverse association with familial breast cancer risk.

CONCLUSIONS

This study illustrates how original approaches such as the phylogeny-based method we used can empower classical molecular epidemiological studies aimed at identifying association or risk modification effects.

DomeTree: a canonical toolkit for mitochondrial DNA analyses in domesticated animals

Mitochondrial DNA (mtDNA) is widely used in various genetic studies of domesticated animals. Many applications require comprehensive knowledge about the phylogeny of mtDNA variants. Herein, we provide the most up-to-date mtDNA phylogeny (i.e. haplogroup tree or matrilineal genealogy) and a standardized hierarchical haplogroup nomenclature system for domesticated cattle, dogs, goats, horses, pigs, sheep, yaks and chickens. These high-resolution mtDNA haplogroup trees based on 1240 complete or near-complete mtDNA genome sequences are available in open resource DomeTree (http://www.dometree.org). In addition, we offer the software MitoToolPy (http://www.mitotool.org/mp.html) to facilitate the mtDNA data analyses. We will continuously and regularly update DomeTree and MitoToolPy.

Mitochondrial genomes of domestic animals need scrutiny

More than 1000 complete or near-complete mitochondrial DNA (mtDNA) sequences have been deposited in GenBank for eight common domestic animals (cattle, dog, goat, horse, pig, sheep, yak and chicken) and their close wild ancestors or relatives, as well. Nevertheless, few efforts have been performed to evaluate the sequence data quality. Herein, we conducted a phylogenetic survey of these complete or near-complete mtDNA sequences based on mtDNA haplogroup trees for the eight animals. We show that errors due to artificial recombination, surplus of mutations and phantom mutations do exist in 14.5% (194/1342) of mtDNA sequences and all of them should be treated with wide caution. We propose some caveats for future mtDNA studies of domestic animals.