Motivation: The increasing availability of mitochondria-targeted and off-target sequencing data in whole-exome and whole-genome sequencing studies (WXS and WGS) has risen the demand of effective pipelines to accurately measure heteroplasmy and to easily recognize the most functionally important mitochondrial variants among a huge number of candidates. To this purpose, we developed MToolBox, a highly automated pipeline to reconstruct and analyze human mitochondrial DNA from high-throughput sequencing data. Results: MToolBox implements an effective computational strategy for mitochondrial genomes assembling and haplogroup assignment also including a prioritization analysis of detected variants. MToolBox provides a Variant Call Format file featuring, for the first time, allele-specific heteroplasmy and annotation files with prioritized variants. MToolBox was tested on simulated samples and applied on 1000 Genomes WXS datasets. Availability and implementation: MToolBox package is available at https://sourceforge.net/projects/mtoolbox/ . Contact:marcella.attimonelli@uniba.it Supplementary information: Supplementary data are available at Bioinformatics online.

Background: Aerobic glycolysis, namely the Warburg effect, is the main hallmark of cancer cells. Mitochondrial respiratory dysfunction has been proposed to be one of the major causes for such glycolytic shift. This hypothesis has been revisited as tumors appear to undergo waves of gene regulation during progression, some of which rely on functional mitochondria. In this framework, the role of mitochondrial complex I is still debated, in particular with respect to the effect of mitochondrial DNA mutations in cancer metabolism. The aim of this work is to provide the proof of concept that functional complex I is necessary to sustain tumor progression. Methods: Complex I-null osteosarcoma cells were complemented with allotopically expressed complex I subunit 1 (MT-ND1). Complex I re-assembly and function recovery, also in terms of NADH consumption, were assessed. Clones were tested for their ability to grow in soft agar and to generate tumor masses in nude mice. Hypoxia levels were evaluated via pimonidazole staining and hypoxia-inducible factor-1α (HIF-1α) immunoblotting and histochemical staining. 454-pyrosequencing was implemented to obtain global transcriptomic profiling of allotopic and non-allotopic xenografts. Results: Complementation of a truncative mutation in the gene encoding MT-ND1, showed that a functional enzyme was required to perform the glycolytic shift during the hypoxia response and to induce a Warburg profile in vitro and in vivo, fostering cancer progression. Such trigger was mediated by HIF-1α, whose stabilization was regulated after recovery of the balance between α-ketoglutarate and succinate due to a recuperation of NADH consumption that followed complex I rescue. Conclusion: Respiratory complex I is essential for the induction of Warburg effect and adaptation to hypoxia of cancer cells, allowing them to sustain tumor growth. Differently from other mitochondrial tumor suppressor genes, therefore, a complex I severe mutation such as the one here reported may confer anti-tumorigenic properties, highlighting the prognostic values of such genetic markers in cancer.

Interest in human mitochondrial genetic data is constantly increasing among both clinicians and researchers, due to the involvement of mitochondrial DNA (mtDNA) in a number of physiological and pathological processes. Thanks to new sequencing technologies and modern databases, the large amount of information on mtDNA variability may be exploited to gain insights into the relationship between mtDNA variants, phenotypes and diseases. To facilitate this process, we have developed the HmtVar resource, a variant-focused database that allows the exploration of a dataset of over 40 000 human mitochondrial variants. Mitochondrial variation data, initially gathered from the HmtDB platform, are integrated with in-house pathogenicity assessments based on various evaluation criteria and with a set of additional annotations from third-party resources. The result is a comprehensive collection of information of crucial importance for human mitochondrial variation studies and investigation of common and rare diseases in which the mitochondrion may be involved. HmtVar is accessible at https://www.hmtvar.uniba.it and data may be retrieved using either a web interface through the Query page or a state-of-the-art API for programmatic access.

Whole Exome Sequencing (WES) is one of the most used and cost-effective next generation technologies that allows sequencing of all nuclear exons. Off-target regions may be captured if they present high sequence similarity with baits. Bioinformatics tools have been optimized to retrieve a large amount of WES off-target mitochondrial DNA (mtDNA), by exploiting the aspecificity of probes, partially overlapping to Nuclear mitochondrial Sequences (NumtS). The 1000 Genomes project represents one of the widest resources to extract mtDNA sequences from WES data, considering the large effort the scientific community is undertaking to reconstruct human population history using mtDNA as marker, and the involvement of mtDNA in pathology.

The colonization of the nuclear genome by mitochondrial DNA is an ongoing process in eukaryotes and plays an important role in genomic variability. Notwithstanding the DNA sequence availability of about 100 complete eukaryotic genomes, up to now NumtS distribution has been fully reported for a small number of sequenced eukaryotic species. With the aim to clarify the time and way of NumtS evolution, we explored the genomic distribution of NumtS in 23 eukaryotic species using an intra/interspecies in silico approach based on a cross-species similarity search and deeply investigate the evolution of NumtS in mammals. The intra- and interspecies analysis underlined how some mitochondrial regions that populated nuclear genomes can be considered as hotspots. Considering the large amount of NumtS we found in platypus and opossum genomes, we hypothesized the occurrence of an earlier colonization that happened prior to the Prototherian/Therian mammal divergence, approximately 160–210 million years ago. These events are still detectable due to the species-specific dynamics that have affected these genomes. Phylogenetic analyses of NumtS derived from two different mitochondrial DNA loci allowed us to recognize the unusual NumtS evolution that acted differently on primate and non-primate species’ genomes.