Multidimensional vector space representation for convergent evolution and molecular phylogeny

Mitochondrial determinants of mammalian longevity

Current ageing theories are far from satisfactory because of the many determinants involved in ageing. The well-known rate-of-living theory assumes that the product (lifetime energy expenditure, LEE) of maximum lifespan (MLS) and mass-specific basal metabolic rate (msBMR) is approximately constant. Although this theory provides a significant inverse correlation between msBMR and MLS as a whole for mammals, it remains problematic for two reasons. First, several interspecies studies within respective orders (typically within rodents) have shown no inverse relationships between msBMR and MLS. Second, LEE values widely vary in mammals and birds. Here, to solve these two problems, we introduced a new quantity designated as mitochondrial (mt) lifetime energy output, mtLEO = MLS × mtMR, in place of LEE, by using the mt metabolic rate (mtMR) per mitochondrion. Thereby, we found that mtLEO values were distributed more narrowly than LEE ones, and strongly correlated with the four amino-acid variables (AAVs) of Ser, Thr and Cys contents and hydrophobicity of mtDNA-encoded membrane proteins (these variables were related to the stability of these proteins). Consequently, only these two mt items, mtMR and the AAVs, solved the above-mentioned problems in the rate-of-living theory, and thus extensively improved the correlation with MLS compared with that given by LEE.

Evolution of mitochondrial power in vertebrate metazoans

BACKGROUND

Basal metabolic rate (BMR) has a very strong body-mass (M) dependence in an individual animal group, and BMR per unit mass (msBMR) converges on a markedly narrow range even across major taxonomic groups. However, it is here a basic question in metazoan biology how much BMR per unit mitochondrion (mtBMR) changes, and then whether mtBMR can be related to the original molecular mechanism of action of mt-encoded membrane proteins (MMPs) playing a central role in cellular energy production.

METHODOLOGY/PRINCIPAL FINDINGS

Analyzing variations of amino-acid compositions of MMPs across 13 metazoan animal groups, incorporating 2022 sequences, we found a strong inverse correlation between Ser/Thr composition (STC) and hydrophobicity (HYD). A majority of animal groups showed an evolutionary pathway of a gradual increase in HYD and decrease in STC, whereas only the deuterostome lineage revealed a rapid decrease in HYD and increase in STC. The strongest correlations appeared in 5 large subunits (ND4, ND5, ND2, CO1, and CO3) undergoing dynamic conformational changes for the proton-pumping function. The pathway of the majority groups is well understood as reflecting natural selection to reduce mtBMR, since simply raising HYD in MMPs (surrounded by the lipid bilayer) weakens their mobility and strengthens their stability. On the other hand, the marked decrease in HYD of the deuterostome elevates mtBMR, but is accompanied with their instability heightening a turnover rate of mitochondria and then cells. Interestingly, cooperative networks of interhelical hydrogen-bonds between motifs involving Ser and Thr residues can enhance MMP stability.

CONCLUSION/SIGNIFICANCE

This stability enhancement lowers turnover rates of mitochondria/cells and may prolong even longevity, and was indeed founded by strong positive correlations of STC with both mtBMR and longevity. The lowest HYD and highest STC in Aves and Mammals are congruent with their very high mtBMR and long longevity.

Evidence for an ancient adaptive episode of convergent molecular evolution

Documented cases of convergent molecular evolution due to selection are fairly unusual, and examples to date have involved only a few amino acid positions. However, because convergence mimics shared ancestry and is not accommodated by current phylogenetic methods, it can strongly mislead phylogenetic inference when it does occur. Here, we present a case of extensive convergent molecular evolution between snake and agamid lizard mitochondrial genomes that overcomes an otherwise strong phylogenetic signal. Evidence from morphology, nuclear genes, and most sites in the mitochondrial genome support one phylogenetic tree, but a subset of mostly amino acid-altering substitutions (primarily at the first and second codon positions) across multiple mitochondrial genes strongly supports a radically different phylogeny. The relevant sites generally evolved slowly but converged between ancient lineages of snakes and agamids. We estimate that approximately 44 of 113 predicted convergent changes distributed across all 13 mitochondrial protein-coding genes are expected to have arisen from nonneutral causes-a remarkably large number. Combined with strong previous evidence for adaptive evolution in snake mitochondrial proteins, it is likely that much of this convergent evolution was driven by adaptation. These results indicate that nonneutral convergent molecular evolution in mitochondria can occur at a scale and intensity far beyond what has been documented previously, and they highlight the vulnerability of standard phylogenetic methods to the presence of nonneutral convergent sequence evolution.

Adaptive threonine increase in transmembrane regions of mitochondrial proteins in higher primates

BACKGROUND

The mitochondrial (mt) gene tree of placental mammals reveals a very strong acceleration of the amino acid (AA) replacement rate and a change in AA compositional bias in the lineage leading to the higher primates (simians), in contrast to the nuclear gene tree. Whether this acceleration and compositional bias were caused by adaptive evolution at the AA level or directional mutation pressure at the DNA level has been vigorously debated.

METHODOLOGY/PRINCIPAL FINDINGS

Our phylogenetic analysis indicates that the rate acceleration in the simian lineage is accompanied by a marked increase in threonine (Thr) residues in the transmembrane helix regions of mt DNA-encoded proteins. This Thr increase involved the replacement of hydrophobic AAs in the membrane interior. Even after accounting for lack of independence due to phylogeny, a regression analysis reveals a statistical significant positive correlation between Thr composition and longevity in primates.

CONCLUSION/SIGNIFICANCE

Because crucial roles of Thr and Ser in membrane proteins have been proposed to be the formation of hydrogen bonds enhancing helix-helix interactions, the Thr increase detected in the higher primates might be adaptive by serving to reinforce stability of mt proteins in the inner membrane. The correlation between Thr composition in the membrane interior and the longevity of animals is striking, especially because some mt functions are thought to be involved in aging.

Complex population history of two Anopheles dirus mosquito species in Southeast Asia suggests the influence of Pleistocene climate change rather than human-mediated effects

Anopheles dirus and Anopheles baimaii are closely related species which feed on primates, particularly humans, and transmit malaria in the tropical forests of mainland Southeast Asia. Here, we report an in-depth phylogeographic picture based on 269 individuals from 21 populations from mainland Southeast Asia. Analysis of 1537 bp of mtDNA sequence revealed that the population history of A. baimaii is far more complex than previously thought. An old expansion (pre-300 kyr BP) was inferred in northern India/Bangladesh with a wave of south-eastwards expansion arriving at the Thai border (ca 135-173 kyr BP) followed by leptokurtic dispersal very recently (ca 16 kyr BP) into peninsular Thailand. The long and complex population history of these anthropophilic species suggests their expansions are not in response to the relatively recent (ca 40 kyr BP) human expansions in mainland Southeast Asia but, rather, fit well with our understanding of Pleistocene climatic change there.

Mammalian evolution and biomedicine: new views from phylogeny

Recent progress resolving the phylogenetic relationships of the major lineages of mammals has had a broad impact in evolutionary biology, comparative genomics and the biomedical sciences. Novel insights into the timing and historical biogeography of early mammalian diversification have resulted from a new molecular tree for placental mammals coupled with dating approaches that relax the assumption of the molecular clock. We highlight the numerous applications to come from a well-resolved phylogeny and genomic prospecting in multiple lineages of mammals, from identifying regulatory elements in mammalian genomes to assessing the functional consequences of mutations in human disease loci and those driving adaptive evolution.

Robust time estimation reconciles views of the antiquity of placental mammals

BACKGROUND

Molecular studies have reported divergence times of modern placental orders long before the Cretaceous-Tertiary boundary and far older than paleontological data. However, this discrepancy may not be real, but rather appear because of the violation of implicit assumptions in the estimation procedures, such as non-gradual change of evolutionary rate and failure to correct for convergent evolution.

METHODOLOGY/PRINCIPAL FINDINGS

New procedures for divergence-time estimation robust to abrupt changes in the rate of molecular evolution are described. We used a variant of the multidimensional vector space (MVS) procedure to take account of possible convergent evolution. Numerical simulations of abrupt rate change and convergent evolution showed good performance of the new procedures in contrast to current methods. Application to complete mitochondrial genomes identified marked rate accelerations and decelerations, which are not obtained with current methods. The root of placental mammals is estimated to be approximately 18 million years more recent than when assuming a log Brownian motion model. Correcting the pairwise distances for convergent evolution using MVS lowers the age of the root about another 20 million years compared to using standard maximum likelihood tree branch lengths. These two procedures combined revise the root time of placental mammals from around 122 million years ago to close to 84 million years ago. As a result, the estimated distribution of molecular divergence times is broadly consistent with quantitative analysis of the North American fossil record and traditional morphological views.

CONCLUSIONS/SIGNIFICANCE

By including the dual effects of abrupt rate change and directly accounting for convergent evolution at the molecular level, these estimates provide congruence between the molecular results, paleontological analyses and morphological expectations. The programs developed here are provided along with sample data that reproduce the results of this study and are especially applicable studies using genome-scale sequence lengths.

Core set approach to reduce uncertainty of gene trees

Background

A genealogy based on gene sequences within a species plays an essential role in the estimation of the character, structure, and evolutionary history of that species. Because intraspecific sequences are more closely related than interspecific ones, detailed information on the evolutionary process may be available by determining all the node sequences of trees and provide insight into functional constraints and adaptations. However, strong evolutionary correlations on a few lineages make this determination difficult as a whole, and the maximum parsimony (MP) method frequently allows a number of topologies with a same total branching length.

Results

Kitazoe et al. developed multidimensional vector-space representation of phylogeny. It converts additivity of evolutionary distances to orthogonality among the vectors expressing branches, and provides a unified index to measure deviations from the orthogoality. In this paper, this index is used to detect and exclude sequences with large deviations from orthogonality, and then selects a maximum subset ("core set") of sequences for which MP generates a single solution. Once the core set tree is formed whose all the node sequences are given, the excluded sequences are found to have basically two phylogenetic positions on this tree, respectively. Fortunately, since multiple substitutions are rare in intra-species sequences, the variance of nucleotide transitions is confined to a small range. By applying the core set approach to 38 partial env sequences of HIV-1 in a single patient and also 198 mitochondrial COI and COII DNA sequences of Anopheles dirus, we demonstrate how consistently this approach constructs the tree.

Conclusion

In the HIV dataset, we confirmed that the obtained core set tree is the unique maximum set for which MP proposes a single tree. In the mosquito data set, the fluctuation of nucleotide transitions caused by the sequences excluded from the core set was very small. We reproduced this core-set tree by simulation based on random process, and applied our approach to many sets of the obtained endpoint sequences. Consequently, the ninety percent of the endpoint sequences was identified as the core sets and the obtained node sequences were perfectly identical to the true ones.

Evolution of Structural Shape in Bacterial Globin-Related Proteins

The globin family of proteins has a characteristic structural pattern of helix interactions that nonetheless exhibits some variation. A simplified model for globin structural evolution was developed in which protein shape evolved by random change of contacts between helices. A conserved globin domain of 15 bacterial proteins representing four structural families was studied. Using a parsimony approach ancestral structural states could be reconstructed. The distribution of number of contact changes per site for a fixed topology tree fit a gamma distribution. Homoplasy was high, with multiple changes per site and no support for an invariant class of residue-residue contacts. Contacts changed more slowly than sequence. A phylogenetic reconstruction using a distance measure based on the proportion of shared contacts was generally consistent with a sequence-based phylogeny but not highly resolved. Contact pattern convergence between members of different globin family proteins could not be detected. Simulation studies indicated the convergence test was sensitive enough to have detected convergence involving only 10% of the contacts, suggesting a limit on the extent of selection for a specific contact pattern. Contact site methods may provide additional approaches to study the relationship between protein structure and sequence evolution.

Retroposed elements as archives for the evolutionary history of placental mammals

Reconstruction of the placental mammalian (eutherian) evolutionary tree has undergone diverse revisions, and numerous aspects remain hotly debated. Initial hierarchical divisions based on morphology contained many misgroupings due to features that evolved independently by similar selection processes. Molecular analyses corrected many of these misgroupings and the superordinal hierarchy of placental mammals was recently assembled into four clades. However, long or rapid evolutionary periods, as well as directional mutation pressure, can produce molecular homoplasies, similar characteristics lacking common ancestors. Retroposed elements, by contrast, integrate randomly into genomes with negligible probabilities of the same element integrating independently into orthologous positions in different species. Thus, presence/absence analyses of these elements are a superior strategy for molecular systematics. By computationally scanning more than 160,000 chromosomal loci and judiciously selecting from only phylogenetically informative retroposons for experimental high-throughput PCR applications, we recovered 28 clear, independent monophyly markers that conclusively verify the earliest divergences in placental mammalian evolution. Using tests that take into account ancestral polymorphisms, multiple long interspersed elements and long terminal repeat element insertions provide highly significant evidence for the monophyletic clades Boreotheria (synonymous with Boreoeutheria), Supraprimates (synonymous with Euarchontoglires), and Laurasiatheria. More importantly, two retropositions provide new support for a prior scenario of early mammalian evolution that places the basal placental divergence between Xenarthra and Epitheria, the latter comprising all remaining placentals. Due to its virtually homoplasy-free nature, the analysis of retroposon presence/absence patterns avoids the pitfalls of other molecular methodologies and provides a rapid, unequivocal means for revealing the evolutionary history of organisms.

The authors identified and sequenced retroposons in mammalian genomes. The presence and absence of these retroposons provided evolutionary markers from which the authors reconstructed the phylogenetic history of placental mammals.

The Platypus Is in Its Place: Nuclear Genes and Indels Confirm the Sister Group Relation of Monotremes and Therians

Morphological data supports monotremes as the sister group of Theria (extant marsupials + eutherians), but phylogenetic analyses of 12 mitochondrial protein-coding genes have strongly supported the grouping of monotremes with marsupials: the Marsupionta hypothesis. Various nuclear genes tend to support Theria, but a comprehensive study of long concatenated sequences and broad taxon sampling is lacking. We therefore determined sequences from six nuclear genes and obtained additional sequences from the databases to create two large and independent nuclear data sets. One (data set I) emphasized taxon sampling and comprised five genes, with a concatenated length of 2,793 bp, from 21 species (two monotremes, six marsupials, nine placentals, and four outgroups). The other (data set II) emphasized gene sampling and comprised eight genes and three proteins, with a concatenated length of 10,773 bp or 3,669 amino acids, from five taxa (a monotreme, a marsupial, a rodent, human, and chicken). Both data sets were analyzed by parsimony, minimum evolution, maximum likelihood, and Bayesian methods using various models and data partitions. Data set I gave bootstrap support values for Theria between 55% and 100%, while support for Marsupionta was at most 12.3%. Taking base compositional bias into account generally increased the support for Theria. Data set II exclusively supported Theria, with the highest possible values and significantly rejected Marsupionta. Independent phylogenetic evidence in support of Theria was obtained from two single amino acid deletions and one insertion, while no supporting insertions and deletions were found for Marsupionta. On the basis of our data sets, the time of divergence between Monotremata and Theria was estimated at 231-217 MYA and between Marsupialia and Eutheria at 193-186 MYA. The morphological evidence for a basal position of Monotremata, well separated from Theria, is thus fully supported by the available molecular data from nuclear genes.

Orthologous repeats and mammalian phylogenetic inference

Determining phylogenetic relationships between species is a difficult problem, and many phylogenetic relationships remain unresolved, even among eutherian mammals. Repetitive elements provide excellent markers for phylogenetic analysis, because their mode of evolution is predominantly homoplasy-free and unidirectional. Historically, phylogenetic studies using repetitive elements have relied on biological methods such as PCR analysis, and computational inference is limited to a few isolated repeats. Here, we present a novel computational method for inferring phylogenetic relationships from partial sequence data using orthologous repeats. We apply our method to reconstructing the phylogeny of 28 mammals, using more than 1000 orthologous repeats obtained from sequence data available from the NISC Comparative Sequencing Program. The resulting phylogeny has robust bootstrap numbers, and broadly matches results from previous studies which were obtained using entirely different data and methods. In addition, we shed light on some of the debatable aspects of the phylogeny. With rapid expansion of available partial sequence data, computational analysis of repetitive elements holds great promise for the future of phylogenetic inference.