Marcella Attimonelli

The genetic and metabolic signature of oncocytic transformation implicates HIF1α destabilization

We previously showed that disruptive complex I mutations in mitochondrial DNA are the main genetic hallmark of oncocytic tumors of the thyroid and kidney. We here report a high frequency of homoplasmic disruptive mutations in a large panel of oncocytic pituitary and head-and-neck tumors. The presence of such mutations implicates disassembly of respiratory complex I in vivo which in turn contributes to the inability of oncocytic tumors to stabilize HIF1alpha and to display pseudo-hypoxia. By utilizing transmitochondrial cytoplasmic hybrids (cybrids), we induced the shift to homoplasmy of a truncating mutation in the mitochondria-coded MTND1 gene. Such shift is associated with a profound metabolic impairment leading to the imbalance of alpha-ketoglutarate and succinate, the Krebs cycle metabolites which are the main responsible for HIF1alpha stabilization. We conclude that the main hallmarks of oncocytic transformation, namely the occurrence of homoplasmic disruptive mutations and complex I disassembly, may explain the benign nature of oncocytic neoplasms through lack of HIF1alpha stabilization.

HmtDB, a genomic resource for mitochondrion-based human variability studies

HmtDB (http://www.hmtdb.uniba.it:8080/hmdb) is a open resource created to support population genetics and mitochondrial disease studies. The database hosts human mitochondrial genome sequences annotated with population and variability data, the latter being estimated through the application of the SiteVar software based on site-specific nucleotide and amino acid variability calculations. The annotations are manually curated thus adding value to the quality of the information provided to the end-user. Classifier tools implemented in HmtDB allow the prediction of the haplogroup for any human mitochondrial genome currently stored in HmtDB or externally submitted de novo by an end-user. Haplogroup definition is based on the Phylotree system. End-users accessing HmtDB are hence allowed to (i) browse the database through the use of a multi-criterion ‘query’ system; (ii) analyze their own human mitochondrial sequences via the ‘classify’ tool (for complete genomes) or by downloading the ‘fragment-classifier’ tool (for partial sequences); (iii) download multi-alignments with reference genomes as well as variability data.

Prospects of barcoding the Italian wild dendroflora: oaks reveal severe limitations to tracking species identity.

DNA barcoding may be particularly important in influencing ecology, economic issues, and the fundamental crisis facing biodiversity as a standardized, species‐level identification tool for taxonomy assessment. Trees play important roles in the conservation of many land ecosystems, the wood trade, and the definition of biogeographical processes; nevertheless, peculiar biological, evolutionary and taxonomical features will probably constitute an intriguing challenge to barcoders. We examined whether four marker regions (trnh‐psba, rbcL, rpoc1, matK) proposed by the Consortium for the Barcode of Life (CBOL) matched species taxonomy in a preliminary tree biodiversity survey of Italian forested land. Our objective was to provide a test of future in situ applications of DNA barcodes by evaluating the efficacy of species discrimination under the criteria of uniformity of methods and natural co‐occurrence of the species in the main forest ecosystems. Fifty‐two species were included in a floristic study. We obtained 73% total discrimination success, with trnH‐psbA as the best performing marker and oaks as the least responsive plants to the markers used. A further taxon‐based study of Quercus (thirty specimens, 12 species) revealed that this genus is refractory to barcoding (0% discrimination success), a probable consequence of low variation rate at the plastid genome level, hybridization, and the incidence of biogeography. We conclude that some species‐rich tree genera in small geographical regions may prove exceptionally difficult to barcode. Until more efficient markers are developed, we recommend that improved and diversified sampling (multiple locations of sympatric and co‐occurring congenerics) be embraced as a timely and important goal for the precise assessment of haplotype specificity to facilitate the productive application of barcoding in practice.

HmtDB, a Human Mitochondrial Genomic Resource Based on Variability Studies Supporting Population Genetics and Biomedical Research

BackgroundPopulation genetics studies based on the analysis of mtDNA and mitochondrial disease studies have produced a huge quantity of sequence data and related information. These data are at present worldwide distributed in differently organised databases and web sites not well integrated among them. Moreover it is not generally possible for the user to submit and contemporarily analyse its own data comparing them with the content of a given database, both for population genetics and mitochondrial disease data.ResultsHmtDB is a well-integrated web-based human mitochondrial bioinformatic resource aimed at supporting population genetics and mitochondrial disease studies, thanks to a new approach based on site-specific nucleotide and aminoacid variability estimation. HmtDB consists of a database of Human Mitochondrial Genomes, annotated with population data, and a set of bioinformatic tools, able to produce site-specific variability data and to automatically characterize newly sequenced human mitochondrial genomes. A query system for the retrieval of genomes and a web submission tool for the annotation of new genomes have been designed and will soon be implemented. The first release contains 1255 fully annotated human mitochondrial genomes. Nucleotide site-specific variability data and multialigned genomes can be downloaded. Intra-human and inter-species aminoacid variability data estimated on the 13 coding for proteins genes of the 1255 human genomes and 60 mammalian species are also available. HmtDB is freely available, upon registration, at http://www.hmdb.uniba.it.ConclusionThe HmtDB project will contribute towards completing and/or refining haplogroup classification and revealing the real pathogenic potential of mitochondrial mutations, on the basis of variability estimation.

MToolBox: a highly automated pipeline for heteroplasmy annotation and prioritization analysis of human mitochondrial variants in high-throughput sequencing

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.

MitoDrome: a database of Drosophila melanogaster nuclear genes encoding proteins targeted to the mitochondrion.

Mitochondria are organelles present in the cytoplasm of most eukaryotic cells; although they have their own DNA, the majority of the proteins necessary for a functional mitochondrion are coded by the nuclear DNA and only after transcription and translation they are imported in the mitochondrion as proteins. The primary role of the mitochondrion is electron transport and oxidative phosphorylation. Although it has been studied for a long time, the interest of researchers in mitochondria is still alive thanks to the discovery of mitochondrial role in apoptosis, aging and cancer. Aim of the MitoDrome database is to annotate the Drosophila melanogaster nuclear genes coding for mitochondrial proteins in order to contribute to the functional characterization of nuclear genes coding for mitochondrial proteins and to knowledge of gene diseases related to mitochondrial dysfunctions. Indeed D. melanogaster is one of the most studied organisms and a model for the Human genome. Data are derived from the comparison of Human mitochondrial proteins versus the Drosophila genome, ESTs and cDNA sequence data available in the FlyBase database. Links from the MitoDrome entries to the related homologous entries available in MitoNuC will be soon imple-mented. The MitoDrome database is available at http://bighost.area.ba.cnr.it/BIG/MitoDrome. Data are organised in a flat-file format and can be retrieved using the SRS system.

Extraction and annotation of human mitochondrial genomes from 1000 Genomes Whole Exome Sequencing data

BackgroundWhole 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.ResultsA previously published pipeline aimed at assembling mitochondrial genomes from off-target WES reads and further improved to detect insertions and deletions (indels) and heteroplasmy in a dataset of 1242 samples from the 1000 Genomes project, enabled to obtain a nearly complete mitochondrial genome from 943 samples (76% analyzed exomes). The robustness of our computational strategy was highlighted by the reduction of reads amount recognized as mitochondrial in the original annotation produced by the Consortium, due to NumtS filtering.An accurate survey was carried out on 1242 individuals. 215 indels, mostly heteroplasmic, and 3407 single base variants were mapped. A homogeneous mismatches distribution was observed along the whole mitochondrial genome, while a lower frequency of indels was found within protein-coding regions, where frameshift mutations may be deleterious. The majority of indels and mismatches found were not previously annotated in mitochondrial databases since conventional sequencing methods were limited to homoplasmy or quasi-homoplasmy detection. Intriguingly, upon filtering out non haplogroup-defining variants, we detected a widespread population occurrence of rare events predicted to be damaging. Eventually, samples were stratified into blood- and lymphoblastoid-derived to detect possibly different trends of mutability in the two datasets, an analysis which did not yield significant discordances.ConclusionsTo the best of our knowledge, this is likely the most extended population-scale mitochondrial genotyping in humans enriched with the estimation of heteroplasmies.

Respiratory complex I is essential to induce a Warburg profile in mitochondria-defective tumor cells.

BackgroundAerobic 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.MethodsComplex 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.ResultsComplementation 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.ConclusionRespiratory 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.

Placing mitochondrial DNA mutations within the progression model of type I endometrial carcinoma.

Mitochondrial DNA (mtDNA) mutations have been described in almost all types of cancer. However, their exact role and timing of occurrence during tumor development and progression are still a matter of debate. A Vogelstein-like model of progression is well established for endometrial carcinoma (EC), however, mtDNA has been scarcely investigated in these tumors despite the fact that mitochondrial biogenesis increase has been shown to be a hallmark of type I EC. Here, we screened a panel of 23 type I EC tissues and matched typical hyperplasia for mutations in mtDNA and in four oncosupressors/oncogenes, namely PTEN, KRAS, CTNNB1 and TP53. Overall, mtDNA mutations were identified in 69% of cases, while mutational events in nuclear genes occurred in 56% of the cases, indicating that mtDNA mutations may precede the genetic instability of these genes canonically involved in progression from hyperplasia to tumor. Protein expression analysis revealed an increase in mitochondrial biogenesis and activation of oxidative stress response mechanisms in tumor tissues, but not in hyperplasia, in correlation with the occurrence of pathogenic mtDNA mutations. Our results point out an involvement of mtDNA mutations in EC progression and explain the increase in mitochondrial biogenesis of type I EC. Last, since mtDNA mutations occur after hyperplasia, their potential role in contributing to genetic instability may be envisioned.

The reference human nuclear mitochondrial sequences compilation validated and implemented on the UCSC genome browser

BackgroundEukaryotic nuclear genomes contain fragments of mitochondrial DNA called NumtS (Nuclear mitochondrial Sequences), whose mode and time of insertion, as well as their functional/structural role within the genome are debated issues. Insertion sites match with chromosomal breaks, revealing that micro-deletions usually occurring at non-homologous end joining loci become reduced in presence of NumtS. Some NumtS are involved in recombination events leading to fragment duplication. Moreover, NumtS are polymorphic, a feature that renders them candidates as population markers. Finally, they are a cause of contamination during human mtDNA sequencing, leading to the generation of false heteroplasmies.ResultsHere we present RHNumtS.2, the most exhaustive human NumtSome catalogue annotating 585 NumtS, 97% of which were here validated in a European individual and in HapMap samples. The NumtS complete dataset and related features have been made available at the UCSC Genome Browser. The produced sequences have been submitted to INSDC databases. The implementation of the RHNumtS.2 tracks within the UCSC Genome Browser has been carried out with the aim to facilitate browsing of the NumtS tracks to be exploited in a wide range of research applications.ConclusionsWe aimed at providing the scientific community with the most exhaustive overview on the human NumtSome, a resource whose aim is to support several research applications, such as studies concerning human structural variation, diversity, and disease, as well as the detection of false heteroplasmic mtDNA variants. Upon implementation of the NumtS tracks, the application of the BLAT program on the UCSC Genome Browser has now become an additional tool to check for heteroplasmic artefacts, supported by data available through the NumtS tracks.