Gene families are growing rapidly, but standard methods for inferring phylogenies do not scale to alignments with over 10,000 sequences. We present FastTree, a method for constructing large phylogenies and for estimating their reliability. Instead of storing a distance matrix, FastTree stores sequence profiles of internal nodes in the tree. FastTree uses these profiles to implement neighbor-joining and uses heuristics to quickly identify candidate joins. FastTree then uses nearest-neighbor interchanges to reduce the length of the tree. For an alignment with N sequences, L sites, and a different characters, a distance matrix requires O(N^2) space and O(N^2 L) time, but FastTree requires just O( NLa + N sqrt(N) ) memory and O( N sqrt(N) log(N) L a ) time. To estimate the tree's reliability, FastTree uses local bootstrapping, which gives another 100-fold speedup over a distance matrix. For example, FastTree computed a tree and support values for 158,022 distinct 16S ribosomal RNAs in 17 hours and 2.4 gigabytes of memory. Just computing pairwise Jukes-Cantor distances and storing them, without inferring a tree or bootstrapping, would require 17 hours and 50 gigabytes of memory. In simulations, FastTree was slightly more accurate than neighbor joining, BIONJ, or FastME; on genuine alignments, FastTree's topologies had higher likelihoods. FastTree is available at http://microbesonline.org/fasttree.

/pdf/fast-tree-computing-large-minimum-evolution-trees-with-4grq0gwfpy.pdf

Fast Tree: Computing Large Minimum-Evolution Trees with Profiles instead of a Distance Matrix

Mitochondria cooperate with their host cells by contributing to bioenergetics, metabolism, biosynthesis, and cell death or survival functions. Reactive oxygen species (ROS) generated by mitochondria participate in stress signalling in normal cells but also contribute to the initiation of nuclear or mitochondrial DNA mutations that promote neoplastic transformation. In cancer cells, mitochondrial ROS amplify the tumorigenic phenotype and accelerate the accumulation of additional mutations that lead to metastatic behaviour. As mitochondria carry out important functions in normal cells, disabling their function is not a feasible therapy for cancer. However, ROS signalling contributes to proliferation and survival in many cancers, so the targeted disruption of mitochondria-to-cell redox communication represents a promising avenue for future therapy.

/pdf/mitochondrial-ros-in-cancer-initiators-amplifiers-or-an-kvgp64ttyp.pdf

Mitochondrial ROS in cancer: initiators, amplifiers or an Achilles' heel?

Study of monogenic mitochondrial cardiomyopathies may yield insights into mitochondrial roles in cardiac development and disease. Here, we combined patient-derived and genetically engineered induced pluripotent stem cells (iPSCs) with tissue engineering to elucidate the pathophysiology underlying the cardiomyopathy of Barth syndrome (BTHS), a mitochondrial disorder caused by mutation of the gene encoding tafazzin (TAZ ). Using BTHS iPSC-derived cardiomyocytes (iPSC-CMs), we defined metabolic, structural and functional abnormalities associated with TAZ mutation. BTHS iPSC-CMs assembled sparse and irregular sarcomeres, and engineered BTHS 'heart-on-chip' tissues contracted weakly. Gene replacement and genome editing demonstrated that TAZ mutation is necessary and sufficient for these phenotypes. Sarcomere assembly and myocardial contraction abnormalities occurred in the context of normal whole-cell ATP levels. Excess levels of reactive oxygen species mechanistically linked TAZ mutation to impaired cardiomyocyte function. Our study provides new insights into the pathogenesis of Barth syndrome, suggests new treatment strategies and advances iPSC-based in vitro modeling of cardiomyopathy.

/pdf/modeling-the-mitochondrial-cardiomyopathy-of-barth-syndrome-2w94o1kzz6.pdf

Modeling the mitochondrial cardiomyopathy of Barth syndrome with induced pluripotent stem cell and heart-on-chip technologies

The Great Hungarian Plain was a crossroads of cultural transformations that have shaped European prehistory. Here we analyse a 5,000-year transect of human genomes, sampled from petrous bones giving consistently excellent endogenous DNA yields, from 13 Hungarian Neolithic, Copper, Bronze and Iron Age burials including two to high (~22 × ) and seven to ~1 × coverage, to investigate the impact of these on Europe’s genetic landscape. These data suggest genomic shifts with the advent of the Neolithic, Bronze and Iron Ages, with interleaved periods of genome stability. The earliest Neolithic context genome shows a European hunter-gatherer genetic signature and a restricted ancestral population size, suggesting direct contact between cultures after the arrival of the first farmers into Europe. The latest, Iron Age, sample reveals an eastern genomic influence concordant with introduced Steppe burial rites. We observe transition towards lighter pigmentation and surprisingly, no Neolithic presence of lactase persistence.

/pdf/genome-flux-and-stasis-in-a-five-millennium-transect-of-5628lmd1wn.pdf

Genome flux and stasis in a five millennium transect of European prehistory

Recent sequencing studies have extensively explored the somatic alterations present in the nuclear genomes of cancers. Although mitochondria control energy metabolism and apoptosis, the origins and impact of cancer-associated mutations in mtDNA are unclear. In this study, we analyzed somatic alterations in mtDNA from 1675 tumors. We identified 1907 somatic substitutions, which exhibited dramatic replicative strand bias, predominantly C > T and A > G on the mitochondrial heavy strand. This strand-asymmetric signature differs from those found in nuclear cancer genomes but matches the inferred germline process shaping primate mtDNA sequence content. A number of mtDNA mutations showed considerable heterogeneity across tumor types. Missense mutations were selectively neutral and often gradually drifted towards homoplasmy over time. In contrast, mutations resulting in protein truncation undergo negative selection and were almost exclusively heteroplasmic. Our findings indicate that the endogenous mutational mechanism has far greater impact than any other external mutagens in mitochondria and is fundamentally linked to mtDNA replication.

/pdf/origins-and-functional-consequences-of-somatic-mitochondrial-3coyl4kyob.pdf

Origins and functional consequences of somatic mitochondrial DNA mutations in human cancer

Heteroplasmy, the mixture of mitochondrial genomes (mtDNA), varies among individuals and cells. Heteroplasmy levels alter the penetrance of pathological mtDNA mutations, and the susceptibility to age-related diseases such as Parkinson's disease. Although mitochondrial dysfunction occurs in age-related type 2 diabetes mellitus (T2DM), the involvement of heteroplasmy in diabetes is unclear. We hypothesized that the heteroplasmic mutational (HM) pattern may change in T2DM. To test this, we used next-generation sequencing, i.e. massive parallel sequencing (MPS), along with PCR–cloning–Sanger sequencing to analyze HM in blood and skeletal muscle DNA samples from monozygotic (MZ) twins either concordant or discordant for T2DM. Great variability was identified in the repertoires and amounts of HMs among individuals, with a tendency towards more mutations in skeletal muscle than in blood. Whereas many HMs were unique, many were either shared among twin pairs or among tissues of the same individual, regardless of their prevalence. This suggested a heritable influence on even low abundance HMs. We found no clear differences between T2DM and controls. However, we found ∼5-fold increase of HMs in non-coding sequences implying the influence of negative selection (P 5% of the molecules per sample) and in low abundance heteroplasmy (<5% of the molecules). Although our study found no evidence supporting the involvement of HMs in the etiology of T2DM, the twin study found clear evidence of a heritable influence on the accumulation of HMs as well as the signatures of selection in heteroplasmic mutations.

/pdf/mitochondrial-dna-heteroplasmy-in-diabetes-and-normal-adults-3p1xwv38y3.pdf

Mitochondrial DNA heteroplasmy in diabetes and normal adults: role of acquired and inherited mutational patterns in twins

Mutations frequently reoccur in the human mitochondrial DNA (mtDNA). However, it is unclear whether recurrent mtDNA nodal mutations (RNMs), that is, recurrent mutations in stems of unrelated phylogenetic nodes, are functional and hence selectively constrained. To answer this question, we performed comprehensive parsimony and maximum likelihood analyses of 9,868 publicly available whole human mtDNAs revealing 1,606 single nodal mutations (SNMs) and 679 RNMs. We then evaluated the potential functionality of synonymous, nonsynonymous and RNA SNMs and RNMs. For synonymous mutations, we have implemented the Codon Adaptation Index. For nonsynonymous mutations, we assessed evolutionary conservation, and employed previously described pathogenicity score assessment tools. For RNA genes’ mutations, we designed a bioinformatic tool which compiled evolutionary conservation and potential effect on RNA structure. While comparing the functionality scores of nonsynonymous and RNA SNMs and RNMs with those of disease-causing mtDNA mutations, we found significant difference (P < 0.001). However, 24 RNMs and 67 SNMs had comparable values with disease-causing mutations reflecting their potential function thus being the best candidates to participate in adaptive events of unrelated lineages. Strikingly, some functional RNMs occurred in unrelated mtDNA lineages that independently altered susceptibility to the same diseases, thus suggesting common functionality. To our knowledge, this is the most comprehensive analysis of selective signatures in the mtDNA not only within proteins but also within RNA genes. For the first time, we discover virtually all positively selected RNMs in our phylogeny while emphasizing their dual role in past evolutionary events and in disease today.

Functional Recurrent Mutations in the Human Mitochondrial Phylogeny: Dual Roles in Evolution and Disease

Multiple human mutational landscapes of normal and cancer conditions are currently available. However, while the unique mutational patterns of tumors have been extensively studied, little attention has been paid to similarities between malignant and normal conditions. Here we compared the pattern of mutations in the mitochondrial genomes (mtDNAs) of cancer (98 sequences) and natural populations (2400 sequences). De novo mtDNA mutations in cancer preferentially colocalized with ancient variants in human phylogeny. A significant portion of the cancer mutations was organized in recurrent combinations (COMs), reaching a length of seven mutations, which also colocalized with ancient variants. Thus, by analyzing similarities rather than differences in patterns of mtDNA mutations in tumor and human evolution, we discovered evidence for similar selective constraints, suggesting a functional potential for these mutations.

/pdf/mtdna-mutation-pattern-in-tumors-and-human-evolution-are-2mj0kj90qd.pdf

Ilia Zhidkov

Papers

Mitochondrial DNA heteroplasmy in diabetes and normal adults: role of acquired and inherited mutational patterns in twins

Functional Recurrent Mutations in the Human Mitochondrial Phylogeny: Dual Roles in Evolution and Disease

mtDNA mutation pattern in tumors and human evolution are shaped by similar selective constraints

MitoBamAnnotator: A web-based tool for detecting and annotating heteroplasmy in human mitochondrial DNA sequences.

Evolution and disease converge in the mitochondrion.