Claudia Calabrese

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.

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.

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.

Searching for a needle in the haystack: Comparing six methods to evaluate heteroplasmy in difficult sequence context

Mitochondrial DNA (mtDNA) mutations have been involved in disease, aging and cancer and furthermore exploited for evolutionary and forensic investigation. When investigating mtDNA mutations the peculiar aspects of mitochondrial genetics, such as heteroplasmy and threshold effect, require suitable approaches which must be sensitive enough to detect low-level heteroplasmy and, precise enough to quantify the exact mutational load. In order to establish the optimal approach for the evaluation of heteroplasmy, six methods were experimentally compared for their capacity to reveal and quantify mtDNA variants. Drawbacks and advantages of cloning, Fluorescent PCR (F-PCR), denaturing High Performance Liquid Chromatography (dHPLC), quantitative Real-Time PCR (qRTPCR), High Resolution Melting (HRM) and 454 pyrosequencing were determined. In particular, detection and quantification of a mutation in a difficult sequence context were investigated, through analysis of an insertion in a homopolymeric stretch (m.3571insC).

A multi-parametric workflow for the prioritization of mitochondrial DNA variants of clinical interest.

Assigning a pathogenic role to mitochondrial DNA (mtDNA) variants and unveiling the potential involvement of the mitochondrial genome in diseases are challenging tasks in human medicine. Assuming that rare variants are more likely to be damaging, we designed a phylogeny-based prioritization workflow to obtain a reliable pool of candidate variants for further investigations. The prioritization workflow relies on an exhaustive functional annotation through the mtDNA extraction pipeline MToolBox and includes Macro Haplogroup Consensus Sequences to filter out fixed evolutionary variants and report rare or private variants, the nucleotide variability as reported in HmtDB and the disease score based on several predictors of pathogenicity for non-synonymous variants. Cutoffs for both the disease score as well as for the nucleotide variability index were established with the aim to discriminate sequence variants contributing to defective phenotypes. The workflow was validated on mitochondrial sequences from Leber’s Hereditary Optic Neuropathy affected individuals, successfully identifying 23 variants including the majority of the known causative ones. The application of the prioritization workflow to cancer datasets allowed to trim down the number of candidate for subsequent functional analyses, unveiling among these a high percentage of somatic variants. Prioritization criteria were implemented in both standalone (http://sourceforge.net/projects/mtoolbox/) and web version (https://mseqdr.org/mtoolbox.php) of MToolBox.

A comprehensive characterization of mitochondrial DNA mutations in glioblastoma multiforme.

Glioblastoma multiforme (GBM) is the most malignant brain cancer in adults, with a poor prognosis, whose molecular stratification still represents a challenge in pathology and clinics. On the other hand, mitochondrial DNA (mtDNA) mutations have been found in most tumors as modifiers of the bioenergetics state, albeit in GBM a characterization of the mtDNA status is lacking to date. Here, a characterization of the burden of mtDNA mutations in GBM samples was performed. First, investigation of tumor-specific vs. non tumor-specific mutations was carried out with the MToolBox bioinformatics pipeline by analyzing 45 matched tumor/blood samples, from whole genome or whole exome sequencing datasets obtained from The Cancer Genome Atlas (TCGA) consortium. Additionally, the entire mtDNA sequence was obtained in a dataset of 104 fresh-frozen GBM samples. Mitochondrial mutations with potential pathogenic interest were prioritized based on heteroplasmic fraction, nucleotide variability, and in silico prediction of pathogenicity. A preliminary biochemical analysis of the activity of mitochondrial respiratory complexes was also performed on fresh-frozen GBM samples. Although a high number of mutations was detected, we report that the large majority of them does not pass the prioritization filters. Therefore, a relatively limited burden of pathogenic mutations is indeed carried by GBM, which did not appear to determine a general impairment of the respiratory chain. This article is part of a Directed Issue entitled: Energy Metabolism Disorders and Therapies.

Evidence of association of human papillomavirus with prognosis worsening in glioblastoma multiforme

BACKGROUND Glioblastoma multiforme (GBM) is the most malignant brain tumor in adults, but its etiology still remains unknown. Recently, a role of viruses such as cytomegalovirus and JC virus in gliomagenesis has been suggested. Since human papillomavirus (HPV) is considered the most common oncogenic virus in humans, we evaluated its occurrence in GBM samples. MATERIAL AND METHODS Fifty-two formalin-fixed paraffin-embedded primary glioblastoma specimens were retrospectively analyzed. The presence of HPV genome on tumor DNA was assessed by MY/GP nested PCR. Confirmation of HPV detection was obtained by chromogenic in situ hybridization (CISH) and immunohistochemistry (IHC) with an antibody directed against the L1 capsidic protein. Finally, univariate and multivariate proportional-hazards models were used to compare the risk of death among HPV-positive and HPV-negative patients. RESULTS Strikingly, viral DNA was detected after PCR in 12 cases (23%). HPV16 genome was present in 25% infected samples, whereas the remaining samples tested positive for HPV6. CISH confirmed positivity in all infected samples for which enough material was available. Moreover, IHC positivity suggested that production of viral proteins from HPV genome is an ongoing process in GBM cancer cells. Finally an association between HPV infection and a worse prognosis was found in patients upon age stratification with a univariate analysis (HR, 2.10; 95% CI, 1.00-4.44; log-rank P = .045). CONCLUSIONS HPV infection status may be considered an independent prognostic factor in GBM patients and suggests that prevention may be considered, should HPV be recognized as a causative agent in gliomagenesis.

Deep sequencing unearths Nuclear mitochondrial Sequences under Leber's hereditary optic neuropathy-associated false heteroplasmic mitochondrial DNA variants

Lebers hereditary optic neuropathy (LHON) is associated with mitochondrial DNA (mtDNA) ND mutations that are mostly homoplasmic. However, these mutations are not sufficient to explain the peculiar features of penetrance and the tissue-specific expression of the disease and are believed to be causative in association with unknown environmental or other genetic factors. Discerning between clear-cut pathogenetic variants, such as those that appear to be heteroplasmic, and less penetrant variants, such as the homoplasmic, remains a challenging issue that we have addressed here using next-generation sequencing approach. We set up a protocol to quantify MTND5 heteroplasmy levels in a family in which the proband manifests a LHON phenotype. Furthermore, to study this mtDNA haplotype, we applied the cybridization protocol. The results demonstrate that the mutations are mostly homoplasmic, whereas the suspected heteroplasmic feature of the observed mutations is due to the co-amplification of Nuclear mitochondrial Sequences.

A platform independent RNA-Seq protocol for the detection of transcriptome complexity.

BackgroundRecent studies have demonstrated an unexpected complexity of transcription in eukaryotes. The majority of the genome is transcribed and only a little fraction of these transcripts is annotated as protein coding genes and their splice variants. Indeed, most transcripts are the result of antisense, overlapping and non-coding RNA expression. In this frame, one of the key aims of high throughput transcriptome sequencing is the detection of all RNA species present in the cell and the first crucial step for RNA-seq users is represented by the choice of the strategy for cDNA library construction. The protocols developed so far provide the utilization of the entire library for a single sequencing run with a specific platform.ResultsWe set up a unique protocol to generate and amplify a strand-specific cDNA library representative of all RNA species that may be implemented with all major platforms currently available on the market (Roche 454, Illumina, ABI/SOLiD). Our method is reproducible, fast, easy-to-perform and even allows to start from low input total RNA. Furthermore, we provide a suitable bioinformatics tool for the analysis of the sequences produced following this protocol.ConclusionWe tested the efficiency of our strategy, showing that our method is platform-independent, thus allowing the simultaneous analysis of the same sample with different NGS technologies, and providing an accurate quantitative and qualitative portrait of complex whole transcriptomes.