Reconstructing the phylogenetic history of long-term effective population size and life-history traits using patterns of amino acid replacement in mitochondrial genomes of mammals and birds

Half a Century of Controversy: The Neutralist/Selectionist Debate in Molecular Evolution

The neutral and nearly neutral theories, introduced more than 50 yr ago, have raised and still raise passionate discussion regarding the forces governing molecular evolution and their relative importance. The debate, initially focused on the amount of within-species polymorphism and constancy of the substitution rate, has spread, matured, and now underlies a wide range of topics and questions. The neutralist/selectionist controversy has structured the field and influences the way molecular evolutionary scientists conceive their research.

OUP accepted manuscript

The nearly neutral theory is a common framework to describe natural selection at the molecular level. This theory emphasizes the importance of slightly deleterious mutations by recognizing their ability to segregate and eventually get fixed due to genetic drift in spite of the presence of purifying selection. As genetic drift is stronger in smaller than in larger populations, a correlation between population size and molecular measures of natural selection is expected within the nearly neutral theory. However, this hypothesis was originally formulated under equilibrium conditions. As most natural populations are not in equilibrium, testing the relationship empirically may lead to confounded outcomes. Demographic nonequilibria, for instance following a change in population size, are common scenarios that are expected to push the selection–drift relationship off equilibrium. By explicitly modeling the effects of a change in population size on allele frequency trajectories in the Poisson random field framework, we obtain analytical solutions of the nonstationary allele frequency spectrum. This enables us to derive exact results of measures of natural selection and effective population size in a demographic nonequilibrium. The study of their time-dependent relationship reveals a substantial deviation from the equilibrium selection–drift balance after a change in population size. Moreover, we show that the deviation is sensitive to the combination of different measures. These results therefore constitute relevant tools for empirical studies to choose suitable measures for investigating the selection–drift relationship in natural populations. Additionally, our new modeling approach extends existing population genetics theory and can serve as foundation for methodological developments.

Protein Structure, Models of Sequence Evolution, and Data Type Effects in Phylogenetic Analyses of Mitochondrial Data: A Case Study in Birds

Phylogenomic analyses have revolutionized the study of biodiversity, but they have revealed that estimated tree topologies can depend, at least in part, on the subset of the genome that is analyzed. For example, estimates of trees for avian orders differ if protein-coding or non-coding data are analyzed. The bird tree is a good study system because the historical signal for relationships among orders is very weak, which should permit subtle non-historical signals to be identified, while monophyly of orders is strongly corroborated, allowing identification of strong non-historical signals. Hydrophobic amino acids in mitochondrially-encoded proteins, which are expected to be found in transmembrane helices, have been hypothesized to be associated with non-historical signals. We tested this hypothesis by comparing the evolution of transmembrane helices and extramembrane segments of mitochondrial proteins from 420 bird species, sampled from most avian orders. We estimated amino acid exchangeabilities for both structural environments and assessed the performance of phylogenetic analysis using each data type. We compared those relative exchangeabilities with values calculated using a substitution matrix for transmembrane helices estimated using a variety of nuclear- and mitochondrially-encoded proteins, allowing us to compare the bird-specific mitochondrial models with a general model of transmembrane protein evolution. To complement our amino acid analyses, we examined the impact of protein structure on patterns of nucleotide evolution. Models of transmembrane and extramembrane sequence evolution for amino acids and nucleotides exhibited striking differences, but there was no evidence for strong topological data type effects. However, incorporating protein structure into analyses of mitochondrially-encoded proteins improved model fit. Thus, we believe that considering protein structure will improve analyses of mitogenomic data, both in birds and in other taxa.

Quantifying the impact of changes in effective population size and expression level on the rate of coding sequence evolution

Molecular sequences are shaped by selection, where the strength of selection relative to drift is determined by effective population size (N_e). Populations with high N_e are expected to undergo stronger purifying selection, and consequently to show a lower substitution rate for selected mutations relative to the substitution rate for neutral mutations (ω). However, computational models based on biophysics of protein stability have suggested that ω can also be independent of N_e. Together, the response of ω to changes in N_e depends on the specific mapping from sequence to fitness. Importantly, an increase in protein expression level has been found empirically to result in decrease of ω, an observation predicted by theoretical models assuming selection for protein stability. Here, we derive a theoretical approximation for the response of ω to changes in N_e and expression level, under an explicit genotype-phenotype-fitness map. The method is generally valid for additive traits and log-concave fitness functions. We applied these results to protein undergoing selection for their conformational stability and corroborate out findings with simulations under more complex models. We predict a weak response of ω to changes in either N_e or expression level, which are interchangeable. Based on empirical data, we propose that fitness based on the conformational stability may not be a sufficient mechanism to explain the empirically observed variation in ω across species. Other aspects of protein biophysics might be explored, such as protein-protein interactions, which can lead to a stronger response of ω to changes in N_e.

Slow crabs - fast genomes: Locomotory capacity predicts skew magnitude in crustacean mitogenomes

Base composition skews (G-C/G+C) of mitochondrial genomes are believed to be primarily driven by mutational pressure, which is positively correlated with metabolic rate. In marine animals, metabolic rate is also positively correlated with locomotory capacity. Given the central role of mitochondria in energy metabolism, we hypothesised that selection for locomotory capacity should be positively correlated with the strength of purifying selection (dN/dS), and thus be negatively correlated with the skew magnitude. Therefore, these two models assume diametrically opposite associations between the metabolic rate and skew magnitude: positive correlation in the prevailing paradigm, and negative in our working hypothesis. We examined correlations between the skew magnitude, metabolic rate, locomotory capacity, and several other variables previously associated with mitochondrial evolution on 287 crustacean mitogenomes. Weakly locomotory taxa had higher skew magnitude and ω (dN/dS) values, but not the gene order rearrangement rate. Skew and ω magnitudes were correlated. Multilevel regression analyses indicated that three competing variables, body size, gene order rearrangement rate, and effective population size, had negligible impacts on the skew magnitude. In most crustacean lineages selection for locomotory capacity appears to be the primary factor determining the skew magnitude. Contrary to the prevailing paradigm, this implies that adaptive selection outweighs nonadaptive selection (mutation pressure) in crustaceans. However, we found indications that effective population size (nonadaptive factor) may outweigh the impact of locomotory capacity in sessile crustaceans (Thecostraca). In conclusion, skew magnitude is a product of the interplay between adaptive and nonadaptive factors, the balance of which varies among lineages.

Reconstructing the History of Variation in Effective Population Size along Phylogenies

The nearly neutral theory predicts specific relations between effective population size (Ne) and patterns of divergence and polymorphism, which depend on the shape of the distribution of fitness effects (DFE) of new mutations. However, testing these relations is not straightforward, owing to the difficulty in estimating Ne. Here, we introduce an integrative framework allowing for an explicit reconstruction of the phylogenetic history of Ne, thus leading to a quantitative test of the nearly neutral theory and an estimation of the allometric scaling of the ratios of nonsynonymous over synonymous polymorphism (πN/πS) and divergence (dN/dS) with respect to Ne. As an illustration, we applied our method to primates, for which the nearly neutral predictions were mostly verified. Under a purely nearly neutral model with a constant DFE across species, we find that the variation in πN/πS and dN/dS as a function of Ne is too large to be compatible with current estimates of the DFE based on site frequency spectra. The reconstructed history of Ne shows a 10-fold variation across primates. The mutation rate per generation u, also reconstructed over the tree by the method, varies over a 3-fold range and is negatively correlated with Ne. As a result of these opposing trends for Ne and u, variation in πS is intermediate, primarily driven by Ne but substantially influenced by u. Altogether, our integrative framework provides a quantitative assessment of the role of Ne and u in modulating patterns of genetic variation, while giving a synthetic picture of their history over the clade.

Inferring long-term effective population size with Mutation-Selection Models

Mutation–selection phylogenetic codon models are grounded on population genetics first principles and represent a principled approach for investigating the intricate interplay between mutation, selection, and drift. In their current form, mutation–selection codon models are entirely characterized by the collection of site-specific amino-acid fitness profiles. However, thus far, they have relied on the assumption of a constant genetic drift, translating into a unique effective population size (Ne) across the phylogeny, clearly an unrealistic assumption. This assumption can be alleviated by introducing variation in Ne between lineages. In addition to Ne, the mutation rate (μ) is susceptible to vary between lineages, and both should covary with life-history traits (LHTs). This suggests that the model should more globally account for the joint evolutionary process followed by all of these lineage-specific variables (Ne, μ, and LHTs). In this direction, we introduce an extended mutation–selection model jointly reconstructing in a Bayesian Monte Carlo framework the fitness landscape across sites and long-term trends in Ne, μ, and LHTs along the phylogeny, from an alignment of DNA coding sequences and a matrix of observed LHTs in extant species. The model was tested against simulated data and applied to empirical data in mammals, isopods, and primates. The reconstructed history of Ne in these groups appears to correlate with LHTs or ecological variables in a way that suggests that the reconstruction is reasonable, at least in its global trends. On the other hand, the range of variation in Ne inferred across species is surprisingly narrow. This last point suggests that some of the assumptions of the model, in particular concerning the assumed absence of epistatic interactions between sites, are potentially problematic.

Long-Lived Species of Bivalves Exhibit Low MT-DNA Substitution Rates

Bivalves represent valuable taxonomic group for aging studies given their wide variation in longevity (from 1–2 to >500 years). It is well known that aging is associated to the maintenance of Reactive Oxygen Species homeostasis and that mitochondria phenotype and genotype dysfunctions accumulation is a hallmark of these processes. Previous studies have shown that mitochondrial DNA mutation rates are linked to lifespan in vertebrate species, but no study has explored this in invertebrates. To this end, we performed a Bayesian Phylogenetic Covariance model of evolution analysis using 12 mitochondrial protein-coding genes of 76 bivalve species. Three life history traits (maximum longevity, generation time and mean temperature tolerance) were tested against 1) synonymous substitution rates (dS), 2) conservative amino acid replacement rates (Kc) and 3) ratios of radical over conservative amino acid replacement rates (Kr/Kc). Our results confirm the already known correlation between longevity and generation time and show, for the first time in an invertebrate class, a significant negative correlation between dS and longevity. This correlation was not as strong when generation time and mean temperature tolerance variations were also considered in our model (marginal correlation), suggesting a confounding effect of these traits on the relationship between longevity and mtDNA substitution rate. By confirming the negative correlation between dS and longevity previously documented in birds and mammals, our results provide support for a general pattern in substitution rates.

Noncoding regions underpin avian bill shape diversification at macroevolutionary scales

Recent progress has been made in identifying genomic regions implicated in trait evolution on a microevolutionary scale in many species, but whether these are relevant over macroevolutionary time remains unclear. Here, we directly address this fundamental question using bird beak shape, a key evolutionary innovation linked to patterns of resource use, divergence, and speciation, as a model trait. We integrate class-wide geometric-morphometric analyses with evolutionary sequence analyses of 10,322 protein-coding genes as well as 229,001 genomic regions spanning 72 species. We identify 1434 protein-coding genes and 39,806 noncoding regions for which molecular rates were significantly related to rates of bill shape evolution. We show that homologs of the identified protein-coding genes as well as genes in close proximity to the identified noncoding regions are involved in craniofacial embryo development in mammals. They are associated with embryonic stem cell pathways, including BMP and Wnt signaling, both of which have repeatedly been implicated in the morphological development of avian beaks. This suggests that identifying genotype-phenotype association on a genome-wide scale over macroevolutionary time is feasible. Although the coding and noncoding gene sets are associated with similar pathways, the actual genes are highly distinct, with significantly reduced overlap between them and bill-related phenotype associations specific to noncoding loci. Evidence for signatures of recent diversifying selection on our identified noncoding loci in Darwin finch populations further suggests that regulatory rather than coding changes are major drivers of morphological diversification over macroevolutionary times.

Molecular Evolution in Small Steps under Prevailing Negative Selection: A Nearly Universal Rule of Codon Substitution

The widely accepted view that evolution proceeds in small steps is based on two premises: 1) negative selection acts strongly against large differences and 2) positive selection favors small-step changes. The two premises are not biologically connected and should be evaluated separately. We now extend a previous approach to studying codon evolution in the entire genome. Codon substitution rate is a function of the physicochemical distance between amino acids (AAs), equated with the step size of evolution. Between nine pairs of closely related species of plants, invertebrates, and vertebrates, the evolutionary rate is strongly and negatively correlated with a set of AA distances (Δ_U, scaled to [0, 1]). Δ_U, a composite measure of evolutionary rates across diverse taxa, is influenced by almost all of the 48 physicochemical properties used here. The new analyses reveal a crucial trend hidden from previous studies: Δ_U is strongly correlated with the evolutionary rate (R² > 0.8) only when the genes are predominantly under negative selection. Because most genes in most taxa are strongly constrained by negative selection, Δ_U has indeed appeared to be a nearly universal measure of codon evolution. In conclusion, molecular evolution at the codon level generally takes small steps due to the prevailing negative selection. Whether positive selection may, or may not, follow the small-step rule is addressed in a companion study.

Polymorphism Data Assist Estimation of the Nonsynonymous over Synonymous Fixation Rate Ratio ω for Closely Related Species

The ratio of nonsynonymous over synonymous sequence divergence, dN/dS, is a widely used estimate of the nonsynonymous over synonymous fixation rate ratio ω, which measures the extent to which natural selection modulates protein sequence evolution. Its computation is based on a phylogenetic approach and computes sequence divergence of protein-coding DNA between species, traditionally using a single representative DNA sequence per species. This approach ignores the presence of polymorphisms and relies on the indirect assumption that new mutations fix instantaneously, an assumption which is generally violated and reasonable only for distantly related species. The violation of the underlying assumption leads to a time-dependence of sequence divergence, and biased estimates of ω in particular for closely related species, where the contribution of ancestral and lineage-specific polymorphisms to sequence divergence is substantial. We here use a time-dependent Poisson random field model to derive an analytical expression of dN/dS as a function of divergence time and sample size. We then extend our framework to the estimation of the proportion of adaptive protein evolution α. This mathematical treatment enables us to show that the joint usage of polymorphism and divergence data can assist the inference of selection for closely related species. Moreover, our analytical results provide the basis for a protocol for the estimation of ω and α for closely related species. We illustrate the performance of this protocol by studying a population data set of four corvid species, which involves the estimation of ω and α at different time-scales and for several choices of sample sizes.

LMAP_S: Lightweight Multigene Alignment and Phylogeny eStimation

Background

Recent advances in genome sequencing technologies and the cost drop in high-throughput sequencing continue to give rise to a deluge of data available for downstream analyses. Among others, evolutionary biologists often make use of genomic data to uncover phenotypic diversity and adaptive evolution in protein-coding genes. Therefore, multiple sequence alignments (MSA) and phylogenetic trees (PT) need to be estimated with optimal results. However, the preparation of an initial dataset of multiple sequence file(s) (MSF) and the steps involved can be challenging when considering extensive amount of data. Thus, it becomes necessary the development of a tool that removes the potential source of error and automates the time-consuming steps of a typical workflow with high-throughput and optimal MSA and PT estimations.

Results

We introduce LMAP_S (Lightweight Multigene Alignment and Phylogeny eStimation), a user-friendly command-line and interactive package, designed to handle an improved alignment and phylogeny estimation workflow: MSF preparation, MSA estimation, outlier detection, refinement, consensus, phylogeny estimation, comparison and editing, among which file and directory organization, execution, manipulation of information are automated, with minimal manual user intervention. LMAP_S was developed for the workstation multi-core environment and provides a unique advantage for processing multiple datasets. Our software, proved to be efficient throughout the workflow, including, the (unlimited) handling of more than 20 datasets.

Conclusions

We have developed a simple and versatile LMAP_S package enabling researchers to effectively estimate multiple datasets MSAs and PTs in a high-throughput fashion. LMAP_S integrates more than 25 software providing overall more than 65 algorithm choices distributed in five stages. At minimum, one FASTA file is required within a single input directory. To our knowledge, no other software combines MSA and phylogeny estimation with as many alternatives and provides means to find optimal MSAs and phylogenies. Moreover, we used a case study comparing methodologies that highlighted the usefulness of our software. LMAP_S has been developed as an open-source package, allowing its integration into more complex open-source bioinformatics pipelines. LMAP_S package is released under GPLv3 license and is freely available at https://lmap-s.sourceforge.io/.

Simulation data for the estimation of numerical constants for approximating pairwise evolutionary distances between amino acid sequences

Estimating the number of substitution events per site that have occurred during the evolution of a pair of amino acid sequences is a common task in phylogenetics and comparative genomics that often requires quite slow maximum-likelihood procedures when taking into account explicit evolutionary models. Data presented in this article are large sets of numbers of substitution events and associated numbers of observed differences between pairs of aligned amino acid sequences that have been generated through a simulation procedure of sequence evolution under a broad range of evolutionary models. These data are available at https://zenodo.org/record/2653704 (doi:10.5281/zenodo.2653704). They are accompanied in this paper by figures showing the strong relationship between the corresponding evolutionary and uncorrected distances, as well as estimated numerical constants that determine non-linear functions that fit the simulated data. These numerical constants can be useful to quickly estimate pairwise evolutionary distances directly from uncorrected distances between aligned amino acid sequences.

Physicochemical Amino Acid Properties Better Describe Substitution Rates in Large Populations

Substitutions between chemically distant amino acids are known to occur less frequently than those between more similar amino acids. This knowledge, however, is not reflected in most codon substitution models, which treat all nonsynonymous changes as if they were equivalent in terms of impact on the protein. A variety of methods for integrating chemical distances into models have been proposed, with a common approach being to divide substitutions into radical or conservative categories. Nevertheless, it remains unclear whether the resulting models describe sequence evolution better than their simpler counterparts. We propose a parametric codon model that distinguishes between radical and conservative substitutions, allowing us to assess if radical substitutions are preferentially removed by selection. Applying our new model to a range of phylogenomic data, we find differentiating between radical and conservative substitutions provides significantly better fit for large populations, but see no equivalent improvement for smaller populations. Comparing codon and amino acid models using these same data shows that alignments from large populations tend to select phylogenetic models containing information about amino acid exchangeabilities, whereas the structure of the genetic code is more important for smaller populations. Our results suggest selection against radical substitutions is, on average, more pronounced in large populations than smaller ones. The reduced observable effect of selection in smaller populations may be due to stronger genetic drift making it more challenging to detect preferences. Our results imply an important connection between the life history of a phylogenetic group and the model that best describes its evolution.

GC-biased gene conversion conceals the prediction of the nearly neutral theory in avian genomes

BACKGROUND

The nearly neutral theory of molecular evolution predicts that the efficacy of natural selection increases with the effective population size. This prediction has been verified by independent observations in diverse taxa, which show that life-history traits are strongly correlated with measures of the efficacy of selection, such as the dN/dS ratio. Surprisingly, avian taxa are an exception to this theory because correlations between life-history traits and dN/dS are apparently absent. Here we explore the role of GC-biased gene conversion on estimates of substitution rates as a potential driver of these unexpected observations.

RESULTS

We analyze the relationship between dN/dS estimated from alignments of 47 avian genomes and several proxies for effective population size. To distinguish the impact of GC-biased gene conversion from selection, we use an approach that accounts for non-stationary base composition and estimate dN/dS separately for changes affected or unaffected by GC-biased gene conversion. This analysis shows that the impact of GC-biased gene conversion on substitution rates can explain the lack of correlations between life-history traits and dN/dS. Strong correlations between life-history traits and dN/dS are recovered after accounting for GC-biased gene conversion. The correlations are robust to variation in base composition and genomic location.

CONCLUSIONS

Our study shows that gene sequence evolution across a wide range of avian lineages meets the prediction of the nearly neutral theory, the efficacy of selection increases with effective population size. Moreover, our study illustrates that accounting for GC-biased gene conversion is important to correctly estimate the strength of selection.

Determinants of the Efficacy of Natural Selection on Coding and Noncoding Variability in Two Passerine Species

Population genetic theory predicts that selection should be more effective when the effective population size (N_e) is larger, and that the efficacy of selection should correlate positively with recombination rate. Here, we analyzed the genomes of ten great tits and ten zebra finches. Nucleotide diversity at 4-fold degenerate sites indicates that zebra finches have a 2.83-fold larger N_e. We obtained clear evidence that purifying selection is more effective in zebra finches. The proportion of substitutions at 0-fold degenerate sites fixed by positive selection (α) is high in both species (great tit 48%; zebra finch 64%) and is significantly higher in zebra finches. When α was estimated on GC-conservative changes (i.e., between A and T and between G and C), the estimates reduced in both species (great tit 22%; zebra finch 53%). A theoretical model presented herein suggests that failing to control for the effects of GC-biased gene conversion (gBGC) is potentially a contributor to the overestimation of α, and that this effect cannot be alleviated by first fitting a demographic model to neutral variants. We present the first estimates in birds for α in the untranslated regions, and found evidence for substantial adaptive changes. Finally, although purifying selection is stronger in high-recombination regions, we obtained mixed evidence for α increasing with recombination rate, especially after accounting for gBGC. These results highlight that it is important to consider the potential confounding effects of gBGC when quantifying selection and that our understanding of what determines the efficacy of selection is incomplete.

Extreme Deviations from Expected Evolutionary Rates in Archaeal Protein Families

Origin of new biological functions is a complex phenomenon ranging from single-nucleotide substitutions to the gain of new genes via horizontal gene transfer or duplication. Neofunctionalization and subfunctionalization of proteins is often attributed to the emergence of paralogs that are subject to relaxed purifying selection or positive selection and thus evolve at accelerated rates. Such phenomena potentially could be detected as anomalies in the phylogenies of the respective gene families. We developed a computational pipeline to search for such anomalies in 1,834 orthologous clusters of archaeal genes, focusing on lineage-specific subfamilies that significantly deviate from the expected rate of evolution. Multiple potential cases of neofunctionalization and subfunctionalization were identified, including some ancient, house-keeping gene families, such as ribosomal protein S10, general transcription factor TFIIB and chaperone Hsp20. As expected, many cases of apparent acceleration of evolution are associated with lineage-specific gene duplication. On other occasions, long branches in phylogenetic trees correspond to horizontal gene transfer across long evolutionary distances. Significant deceleration of evolution is less common than acceleration, and the underlying causes are not well understood; functional shifts accompanied by increased constraints could be involved. Many gene families appear to be “highly evolvable,” that is, include both long and short branches. Even in the absence of precise functional predictions, this approach allows one to select targets for experimentation in search of new biology.

In Cold Blood: Compositional Bias and Positive Selection Drive the High Evolutionary Rate of Vampire Bats Mitochondrial Genomes

Mitochondrial genomes of animals have long been considered to evolve under the action of purifying selection. Nevertheless, there is increasing evidence that they can also undergo episodes of positive selection in response to shifts in physiological or environmental demands. Vampire bats experienced such a shift, as they are the only mammals feeding exclusively on blood and possessing anatomical adaptations to deal with the associated physiological requirements (e.g., ingestion of high amounts of liquid water and iron). We sequenced eight new chiropteran mitogenomes including two species of vampire bats, five representatives of other lineages of phyllostomids and one close outgroup. Conducting detailed comparative mitogenomic analyses, we found evidence for accelerated evolutionary rates at the nucleotide and amino acid levels in vampires. Moreover, the mitogenomes of vampire bats are characterized by an increased cytosine (C) content mirrored by a decrease in thymine (T) compared with other chiropterans. Proteins encoded by the vampire bat mitogenomes also exhibit a significant increase in threonine (Thr) and slight reductions in frequency of the hydrophobic residues isoleucine (Ile), valine (Val), methionine (Met), and phenylalanine (Phe). We show that these peculiar substitution patterns can be explained by the co-occurrence of both neutral (mutational bias) and adaptive (positive selection) processes. We propose that vampire bat mitogenomes may have been impacted by selection on mitochondrial proteins to accommodate the metabolism and nutritional qualities of blood meals.

Analysis of mitochondrial DNA sequence and copy number variation across five high-altitude species and their low-altitude relatives

High-altitude inhospitable environments impose a formidable life challenge for the local animals. Training and exposure to high-altitude environments produce both distinct physiological and phenotypic characteristics. The mitochondrion, an organelle crucial for the energy production, plays an important role in hypoxia adaptation. In this study, we investigated the mitochondrial DNA (mtDNA) polymorphism and copy number variation between the population pairs from distinct altitudes across the multi-species. Higher mitochondrial DNA control region's genetic diversity is conspicuous in high-altitude animals versus low-altitude relatives. We also found an accordant decrease of mtDNA copy number in most of the tissues from high-altitude animals. Compared to mammals, chickens have significantly distinct mitogenomic characteristics, and more significant changes in the skeletal muscle mtDNA copy number between high- and low-altitude individuals. Our study catches a snapshot of the biological similarities and differences in the mitochondrial high-altitude acclimation across the species.

Human long intrinsically disordered protein regions are frequent targets of positive selection

Intrinsically disordered regions occur frequently in proteins and are characterized by a lack of a well-defined three-dimensional structure. Although these regions do not show a higher order of structural organization, they are known to be functionally important. Disordered regions are rapidly evolving, largely attributed to relaxed purifying selection and an increased role of genetic drift. It has also been suggested that positive selection might contribute to their rapid diversification. However, for our own species, it is currently unknown whether positive selection has played a role during the evolution of these protein regions. Here, we address this question by investigating the evolutionary pattern of more than 6600 human proteins with intrinsically disordered regions and their ordered counterparts. Our comparative approach with data from more than 90 mammalian genomes uses a priori knowledge of disordered protein regions, and we show that this increases the power to detect positive selection by an order of magnitude. We can confirm that human intrinsically disordered regions evolve more rapidly, not only within humans but also across the entire mammalian phylogeny. They have, however, experienced substantial evolutionary constraint, hinting at their fundamental functional importance. We find compelling evidence that disordered protein regions are frequent targets of positive selection and estimate that the relative rate of adaptive substitutions differs fourfold between disordered and ordered protein regions in humans. Our results suggest that disordered protein regions are important targets of genetic innovation and that the contribution of positive selection in these regions is more pronounced than in other protein parts.

Genomic Signature of an Avian Lilliput Effect across the K-Pg Extinction

Survivorship following major mass extinctions may be associated with a decrease in body size-a phenomenon called the Lilliput Effect. Body size is a strong predictor of many life history traits (LHTs), and is known to influence demography and intrinsic biological processes. Pronounced changes in organismal size throughout Earth history are therefore likely to be associated with concomitant genome-wide changes in evolutionary rates. Here, we report pronounced heterogeneity in rates of molecular evolution (varying up to $\sim$20-fold) across a large-scale avian phylogenomic data set and show that nucleotide substitution rates are strongly correlated with body size and metabolic rate. We also identify potential body size reductions associated with the Cretaceous-Paleogene (K-Pg) transition, consistent with a Lilliput Effect in the wake of that mass extinction event. We posit that selection for reduced body size across the K-Pg extinction horizon may have resulted in transient increases in substitution rate along the deepest branches of the extant avian tree of life. This "hidden" rate acceleration may result in both strict and relaxed molecular clocks over-estimating the age of the avian crown group through the relationship between life history and demographic parameters that scale with molecular substitution rate. If reductions in body size (and/or selection for related demographic parameters like short generation times) are a common property of lineages surviving mass extinctions, this phenomenon may help resolve persistent divergence time debates across the tree of life. Furthermore, our results suggest that selection for certain LHTs may be associated with deterministic molecular evolutionary outcomes.

Protein evolution depends on multiple distinct population size parameters

That population size affects the fate of new mutations arising in genomes, modulating both how frequently they arise and how efficiently natural selection is able to filter them, is well established. It is therefore clear that these distinct roles for population size that characterize different processes should affect the evolution of proteins and need to be carefully defined. Empirical evidence is consistent with a role for demography in influencing protein evolution, supporting the idea that functional constraints alone do not determine the composition of coding sequences.

Given that the relationship between population size, mutant fitness and fixation probability has been well characterized, estimating fitness from observed substitutions is well within reach with well-formulated models. Molecular evolution research has, therefore, increasingly begun to leverage concepts from population genetics to quantify the selective effects associated with different classes of mutation. However, in order for this type of analysis to provide meaningful information about the intra- and inter-specific evolution of coding sequences, a clear definition of concepts of population size, what they influence, and how they are best parameterized is essential.

Here, we present an overview of the many distinct concepts that “population size” and “effective population size” may refer to, what they represent for studying proteins, and how this knowledge can be harnessed to produce better specified models of protein evolution.

A mixed relaxed clock model

Over recent years, several alternative relaxed clock models have been proposed in the context of Bayesian dating. These models fall in two distinct categories: uncorrelated and autocorrelated across branches. The choice between these two classes of relaxed clocks is still an open question. More fundamentally, the true process of rate variation may have both long-term trends and short-term fluctuations, suggesting that more sophisticated clock models unfolding over multiple time scales should ultimately be developed. Here, a mixed relaxed clock model is introduced, which can be mechanistically interpreted as a rate variation process undergoing short-term fluctuations on the top of Brownian long-term trends. Statistically, this mixed clock represents an alternative solution to the problem of choosing between autocorrelated and uncorrelated relaxed clocks, by proposing instead to combine their respective merits. Fitting this model on a dataset of 105 placental mammals, using both node-dating and tip-dating approaches, suggests that the two pure clocks, Brownian and white noise, are rejected in favour of a mixed model with approximately equal contributions for its uncorrelated and autocorrelated components. The tip-dating analysis is particularly sensitive to the choice of the relaxed clock model. In this context, the classical pure Brownian relaxed clock appears to be overly rigid, leading to biases in divergence time estimation. By contrast, the use of a mixed clock leads to more recent and more reasonable estimates for the crown ages of placental orders and superorders. Altogether, the mixed clock introduced here represents a first step towards empirically more adequate models of the patterns of rate variation across phylogenetic trees.This article is part of the themed issue 'Dating species divergences using rocks and clocks'.

The influence of molecular markers and methods on inferring the phylogenetic relationships between the representatives of the Arini (parrots, Psittaciformes), determined on the basis of their complete mitochondrial genomes

BACKGROUND

Conures are a morphologically diverse group of Neotropical parrots classified as members of the tribe Arini, which has recently been subjected to a taxonomic revision. The previously broadly defined Aratinga genus of this tribe has been split into the 'true' Aratinga and three additional genera, Eupsittula, Psittacara and Thectocercus. Popular markers used in the reconstruction of the parrots' phylogenies derive from mitochondrial DNA. However, current phylogenetic analyses seem to indicate conflicting relationships between Aratinga and other conures, and also among other Arini members. Therefore, it is not clear if the mtDNA phylogenies can reliably define the species tree. The inconsistencies may result from the variable evolution rate of the markers used or their weak phylogenetic signal. To resolve these controversies and to assess to what extent the phylogenetic relationships in the tribe Arini can be inferred from mitochondrial genomes, we compared representative Arini mitogenomes as well as examined the usefulness of the individual mitochondrial markers and the efficiency of various phylogenetic methods.

RESULTS

Single molecular markers produced inconsistent tree topologies, while different methods offered various topologies even for the same marker. A significant disagreement in these tree topologies occurred for cytb, nd2 and nd6 genes, which are commonly used in parrot phylogenies. The strongest phylogenetic signal was found in the control region and RNA genes. However, these markers cannot be used alone in inferring Arini phylogenies because they do not provide fully resolved trees. The most reliable phylogeny of the parrots under study is obtained only on the concatenated set of all mitochondrial markers. The analyses established significantly resolved relationships within the former Aratinga representatives and the main genera of the tribe Arini. Such mtDNA phylogeny can be in agreement with the species tree, owing to its match with synapomorphic features in plumage colouration.

CONCLUSIONS

Phylogenetic relationships inferred from single mitochondrial markers can be incorrect and contradictory. Therefore, such phylogenies should be considered with caution. Reliable results can be produced by concatenated sets of all or at least the majority of mitochondrial genes and the control region. The results advance a new view on the relationships among the main genera of Arini and resolve the inconsistencies between the taxa that were previously classified as the broadly defined genus Aratinga. Although gene and species trees do not always have to be consistent, the mtDNA phylogenies for Arini can reflect the species tree.

Hierarchical tissue organization as a general mechanism to limit the accumulation of somatic mutations

How can tissues generate large numbers of cells, yet keep the divisional load (the number of divisions along cell lineages) low in order to curtail the accumulation of somatic mutations and reduce the risk of cancer? To answer the question we consider a general model of hierarchically organized self-renewing tissues and show that the lifetime divisional load of such a tissue is independent of the details of the cell differentiation processes, and depends only on two structural and two dynamical parameters. Our results demonstrate that a strict analytical relationship exists between two seemingly disparate characteristics of self-renewing tissues: divisional load and tissue organization. Most remarkably, we find that a sufficient number of progressively slower dividing cell types can be almost as efficient in minimizing the divisional load, as non-renewing tissues. We argue that one of the main functions of tissue-specific stem cells and differentiation hierarchies is the prevention of cancer.

To limit the accumulation of somatic mutations, renewing tissues must minimize the number of times each cell divides during differentiation. Here, the authors analytically derive the lower limit of lifetime divisional load of a tissue, show that hierarchically differentiating tissues can approach this limit, and that this depends on uneven divisional rates across the hierarchy.

Darwinism for the Genomic Age: Connecting Mutation to Diversification

A growing body of evidence suggests that rates of diversification of biological lineages are correlated with differences in genome-wide mutation rate. Given that most research into differential patterns of diversification rate have focused on species traits or ecological parameters, a connection to the biochemical processes of genome change is an unexpected observation. While the empirical evidence for a significant association between mutation rate and diversification rate is mounting, there has been less effort in explaining the factors that mediate this connection between genetic change and species richness. Here we draw together empirical studies and theoretical concepts that may help to build links in the explanatory chain that connects mutation to diversification. First we consider the way that mutation rates vary between species. We then explore how differences in mutation rates have flow-through effects to the rate at which populations acquire substitutions, which in turn influences the speed at which populations become reproductively isolated from each other due to the acquisition of genomic incompatibilities. Since diversification rate is commonly measured from phylogenetic analyses, we propose a conceptual approach for relating events of reproductive isolation to bifurcations on molecular phylogenies. As we examine each of these relationships, we consider theoretical models that might shine a light on the observed association between rate of molecular evolution and diversification rate, and critically evaluate the empirical evidence for these links, focusing on phylogenetic comparative studies. Finally, we ask whether we are getting closer to a real understanding of the way that the processes of molecular evolution connect to the observable patterns of diversification.

Complete mitochondrial genome of the Italian slow-worm Anguis veronensis Pollini, 1818, and its comparison with mitogenomes of other Anguis species

In this paper, we present complete mitochondrial genome of the Italian legless lizard species Anguis veronensis Pollini, 1818. The complete mtDNA consisted of 13 protein-coding genes, 22 tRNAs, and two rRNA genes which in total formed a DNA strand of 17,322 bp. Anguis veronensis mitogenome had the same gene order as two other compared Anguis spp., i.e. A. cephallonica and A. fragilis. The base composition of A. veronensis mitochondrial genome was A - 30.8%, T - 24.9%, C - 29.9%, G - 14.4%, with an A + T bias (55.7%). The newly described genome provides valuable data for future comparative mitogenomic analysis within Anguis genus.

Complete mitochondrial genome of golden conure (Guaruba guarouba)

Arini tribe with 19 genera is the most diversified tribe of neotropical parrots. Six of them are classified as macaws and nine as conures. The presence of bare facial area distinguishes macaws from conures and other members of this tribe. However, such morphological division seems to be disputable as the smallest macaw (monotypic Diopsittaca genus) turned out to be more closely related to three monotypic conures genera (Guaruba, Leptosittaca, Thectocercus) than to other macaws. We sequenced the complete mitochondrial genome of Guaruba guarouba to enrich the resource of molecular markers for examination of phylogenetic relationships between macaws and conures.

Complete mitochondrial genome of white-eyed parakeet (Psittacara leucophthalmus): the basal species to other Psittacara

Recently resurrected Psittacara genus is one of the 19 recognized in parrot tribe Arini. The status of taxa within Psittacara remains controversial because some forms are treated as species or subspecies depending on authorities. Evolutionary history of Psittacara is also unclear because related phylogenetic clades contain taxa from distant and non-overlapping regions. However, the basal placement of Psittacara leucophthalmus with wide South American distribution suggests that other taxa with restricted range could emerge by a local split of larger population. We sequenced P. leucophthalmus mitogenome to increase the set of sequences required to determine taxonomic level and phylogeny of Psittacara taxa.

The first complete mitochondrial genome sequence from the blue-headed parrot (Pionus menstruus menstruus): a representative for the genus

Androglossini is one of four tribes recognized within a neotropical parrot subfamily Arinae. The tribe includes 10 genera of which Pionus is represented by eight species. However, its evolutionary diversification and relationship with other Androglossini members are still unclear. Depending on studied molecular markers, Pionus is closely related with Amazona genus or two monotypic genera Alipiopsitta and Graydidascalus or the clade in which Amazona genus is sister to Alipiopsitta and Graydidascalus. Therefore, we sequenced Pionus menstruus menstruus mitogenome to gain molecular data appropriate for future studies to resolve these discrepancies obtained in various phylogenetic analyses published so far.

Bone-associated gene evolution and the origin of flight in birds

BACKGROUND

Bones have been subjected to considerable selective pressure throughout vertebrate evolution, such as occurred during the adaptations associated with the development of powered flight. Powered flight evolved independently in two extant clades of vertebrates, birds and bats. While this trait provided advantages such as in aerial foraging habits, escape from predators or long-distance travels, it also imposed great challenges, namely in the bone structure.

RESULTS

We performed comparative genomic analyses of 89 bone-associated genes from 47 avian genomes (including 45 new), 39 mammalian, and 20 reptilian genomes, and demonstrate that birds, after correcting for multiple testing, have an almost two-fold increase in the number of bone-associated genes with evidence of positive selection (~52.8 %) compared with mammals (~30.3 %). Most of the positive-selected genes in birds are linked with bone regulation and remodeling and thirteen have been linked with functional pathways relevant to powered flight, including bone metabolism, bone fusion, muscle development and hyperglycemia levels. Genes encoding proteins involved in bone resorption, such as TPP1, had a high number of sites under Darwinian selection in birds.

CONCLUSIONS

Patterns of positive selection observed in bird ossification genes suggest that there was a period of intense selective pressure to improve flight efficiency that was closely linked with constraints on body size.

Complete mitochondrial genome of Blue-headed Macaw (Primolius couloni): its comparison with mitogenome of Blue-throated Macaw (Ara glaucogularis)

Primolius is a genus of midsized Macaws comprising three species. Blue-headed Macaws (Primolius couloni) are native to eastern Peru, extreme western Brazil and north-western Bolivia. In this study, full mitochondrial genome of considered species was sequenced. It is 16,995 bp long and contains 13 protein-coding genes, 2 rRNAs, 22 tRNAs and a control region. Its comparison with published Blue-throated Macaw (Ara glaucogularis) mitogenome revealed their high degree of identity. Primolius couloni mitogenome is the first complete genomic sequence of this genus. It will be indispensable to refine the phylogenetic relationships within the tribe Arini and will enrich the resource of markers for systematic, phylogenetic and population genetic studies.

Kr/Kc but not dN/dS correlates positively with body mass in birds, raising implications for inferring lineage-specific selection

Background

The ratio of the rates of non-synonymous and synonymous substitution (d_N/d_S) is commonly used to estimate selection in coding sequences. It is often suggested that, all else being equal, d_N/d_S should be lower in populations with large effective size (N_e) due to increased efficacy of purifying selection. As N_e is difficult to measure directly, life history traits such as body mass, which is typically negatively associated with population size, have commonly been used as proxies in empirical tests of this hypothesis. However, evidence of whether the expected positive correlation between body mass and d_N/d_S is consistently observed is conflicting.

Results

Employing whole genome sequence data from 48 avian species, we assess the relationship between rates of molecular evolution and life history in birds. We find a negative correlation between d_N/d_S and body mass, contrary to nearly neutral expectation. This raises the question whether the correlation might be a method artefact. We therefore in turn consider non-stationary base composition, divergence time and saturation as possible explanations, but find no clear patterns. However, in striking contrast to d_N/d_S, the ratio of radical to conservative amino acid substitutions (K_r/K_c) correlates positively with body mass.

Conclusions

Our results in principle accord with the notion that non-synonymous substitutions causing radical amino acid changes are more efficiently removed by selection in large populations, consistent with nearly neutral theory. These findings have implications for the use of d_N/d_S and suggest that caution is warranted when drawing conclusions about lineage-specific modes of protein evolution using this metric.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0542-8) contains supplementary material, which is available to authorized users.

Complete mitochondrial genome of Mitred Conure (Psittacara mitratus): its comparison with mitogenome of Socorro Conure (Psittacara brevipes)

The Mitred Conure (Psittacara mitratus) is a species native to the forests and woodlands in the Andes in central Peru, south through west-central Bolivia, to northwestern Argentina. The genus Psittacara have recently been distinguished on the basis of molecular revision of the Aratinga species. In consequence, it became one of the nine genera, which form morphologically diverse group termed as Conures. Although, the number of species/subspecies belonging to Psittacara remains controversial, it is now believed that the genus is represented by 11 species. Taxonomic position of Mitred Conure was determined by molecular research and makes the species crucial for examination of evolutionary diversification of the genus. Therefore, we sequenced complete mitochondrial genome of P. mitratus mitogenome to gain a source of comprehensive molecular data. It will be indispensable to refine the phylogenetic relationships within "Psittacara group" as well as within the tribe Arini.

Complete mitochondrial genome of Red-throated Conure (Psittacara rubritorquis): its comparison with mitogenome of Socorro Conure (Psittacara brevipes)

According to some taxonomists the Red-throated Conure (Psittacara rubritorquis) is considered a subspecies of Green Conure (Psittacara holochlora). Some other classifications treat rubritorquis as a separate species based on relatively minor morphological differences between both species/subspecies. So far, taxonomic position of P. rubritorquis was determined by molecular researches using only ND2 gene sequence or incomplete combined mitochondrial ND2, COI and CYTB gene sequences. Obtained outcomes found that P. rubritorquis should be treated as a subspecies of P. holochlora. However, the lack of P. h. brewsterii and P. h. strenua samples as well as incompleteness of combined mitochondrial sequence do not exclude opposite scenario. Therefore, we sequenced P. rubritorquis mitogenome to gain a source of molecular data appropriate for future examination of evolutionary diversification of the P. holochlora group.

Complete mitochondrial genome of Red-bellied Macaw (Orthopsittaca manilata): its comparison with mitogenome of Blue-throated Macaw (Ara glaucogularis)

The Red-bellied Macaw (Orthopsittaca manilata) is a species of the monotypic genus Orthopsittaca. The genus is one of the six genera, which form morphologically diverse group termed as Macaws. Individuals of Orthopsittaca manilata species are found in extremely large Amazonian area of South America. In this study, full mitochondrial genome of considered species was sequenced. It is 16,985 bp long and contains 13 protein-coding genes, 2 rRNAs, 22 tRNAs and a control region. Its comparison with published Blue-throated Macaw (Ara glaucogularis) mitogenome revealed their high degree of identity. Presented Orthopsittaca manilata mitogenome is the first complete genomic sequence of this genus. It will enrich the resource of molecular markers for future examination of evolutionary diversification of Macaws. It will be also indispensable to refine the phylogenetic relationships within the tribe Arini.

Monte Carlo algorithms for Brownian phylogenetic models

MOTIVATION

Brownian models have been introduced in phylogenetics for describing variation in substitution rates through time, with applications to molecular dating or to the comparative analysis of variation in substitution patterns among lineages. Thus far, however, the Monte Carlo implementations of these models have relied on crude approximations, in which the Brownian process is sampled only at the internal nodes of the phylogeny or at the midpoints along each branch, and the unknown trajectory between these sampled points is summarized by simple branchwise average substitution rates.

RESULTS

A more accurate Monte Carlo approach is introduced, explicitly sampling a fine-grained discretization of the trajectory of the (potentially multivariate) Brownian process along the phylogeny. Generic Monte Carlo resampling algorithms are proposed for updating the Brownian paths along and across branches. Specific computational strategies are developed for efficient integration of the finite-time substitution probabilities across branches induced by the Brownian trajectory. The mixing properties and the computational complexity of the resulting Markov chain Monte Carlo sampler scale reasonably with the discretization level, allowing practical applications with up to a few hundred discretization points along the entire depth of the tree. The method can be generalized to other Markovian stochastic processes, making it possible to implement a wide range of time-dependent substitution models with well-controlled computational precision.

AVAILABILITY

The program is freely available at www.phylobayes.org.

Mitochondrial DNA as a tool for reconstructing past life-history traits in mammals

Reconstructing the ancestral characteristics of species is a major goal in evolutionary and comparative biology. Unfortunately, fossils are not always available and sufficiently informative, and phylogenetic methods based on models of character evolution can be unsatisfactory. Genomic data offer a new opportunity to estimate ancestral character states, through (i) the correlation between DNA evolutionary processes and species life-history traits and (ii) available reliable methods for ancestral sequence inference. Here, we assess the relevance of mitochondrial DNA--the most popular molecular marker in animals--as a predictor of ancestral life-history traits in mammals, using the order of Cetartiodactyla as a benchmark. Using the complete set of 13 mitochondrial protein-coding genes, we show that the lineage-specific nonsynonymous over synonymous substitution rate ratio (dN/dS) is closely correlated with the species body mass, longevity and age of sexual maturity in Cetartiodactyla and can be used as a marker of ancestral traits provided that the noise introduced by short branches is appropriately dealt with. Based on ancestral dN/dS estimates, we predict that the first cetartiodactyls were relatively small animals (around 20 kg). This finding is in accordance with Cope's rule and the fossil record but could not be recovered via continuous character evolution methods.

Complete mitochondrial genome of Critically Endangered Blue-throated Macaw (Ara glaucogularis): its comparison with partial mitogenome of Scarlet Macaw (Ara macao)

Blue-throated Macaw (Ara glaucogularis) is Critically Endangered species of parrot endemic to small Bolivian area. In this study, full mitochondrial genome of considered species was sequenced. It is 16,983 bp long and contains 13 protein-coding genes, 2 rRNAs, 22 tRNAs and a control region. It`s comparison with published Ara macao mitogenome revealed their high degree of identity. On the other hand, analysis of both the genome compositions showed incompleteness of Ara macao CYTB gene. Hence, mitogenome of Ara macao species occurred only partial sequence. In consequence, Ara glaucogularis mitogenome is the first complete Macaw sequence, which will be indispensable to refine the phylogenetic relationships within the tribe Arini and will enrich the resource of markers for systematic, phylogenetic and population genetic studies.

Evidence for GC-biased gene conversion as a driver of between-lineage differences in avian base composition

BACKGROUND

While effective population size (Ne) and life history traits such as generation time are known to impact substitution rates, their potential effects on base composition evolution are less well understood. GC content increases with decreasing body mass in mammals, consistent with recombination-associated GC biased gene conversion (gBGC) more strongly impacting these lineages. However, shifts in chromosomal architecture and recombination landscapes between species may complicate the interpretation of these results. In birds, interchromosomal rearrangements are rare and the recombination landscape is conserved, suggesting that this group is well suited to assess the impact of life history on base composition.

RESULTS

Employing data from 45 newly and 3 previously sequenced avian genomes covering a broad range of taxa, we found that lineages with large populations and short generations exhibit higher GC content. The effect extends to both coding and non-coding sites, indicating that it is not due to selection on codon usage. Consistent with recombination driving base composition, GC content and heterogeneity were positively correlated with the rate of recombination. Moreover, we observed ongoing increases in GC in the majority of lineages.

CONCLUSIONS

Our results provide evidence that gBGC may drive patterns of nucleotide composition in avian genomes and are consistent with more effective gBGC in large populations and a greater number of meioses per unit time; that is, a shorter generation time. Thus, in accord with theoretical predictions, base composition evolution is substantially modulated by species life history.

Domestication and the mitochondrial genome: comparing patterns and rates of molecular evolution in domesticated mammals and birds and their wild relatives

Studies of domesticated animals have led to the suggestion that domestication could have significant effects on patterns of molecular evolution. In particular, analyses of mitochondrial genome sequences from domestic dogs and yaks have yielded higher ratios of non-synonymous to synonymous substitutions in the domesticated lineages than in their wild relatives. These results are important because they imply that changes to selection or population size operating over a short timescale can cause significant changes to the patterns of mitochondrial molecular evolution. In this study, our aim is to test whether the impact on mitochondrial genome evolution is a general feature of domestication or whether it is specific to particular examples. We test whether domesticated mammals and birds have consistently different patterns of molecular evolution than their wild relatives for 16 phylogenetically independent comparisons of mitochondrial genome sequences. We find no consistent difference in branch lengths or dN/dS between domesticated and wild lineages. We also find no evidence that our failure to detect a consistent pattern is due to the short timescales involved or low genetic distance between domesticated lineages and their wild relatives. However, removing comparisons where the wild relative may also have undergone a bottleneck does reveal a pattern consistent with reduced effective population size in domesticated lineages. Our results suggest that, although some domesticated lineages may have undergone changes to selective regime or effective population size that could have affected mitochondrial evolution, it is not possible to generalize these patterns over all domesticated mammals and birds.

Eight new mtDNA sequences of glass sponges reveal an extensive usage of +1 frameshifting in mitochondrial translation

Three previously studied mitochondrial genomes of glass sponges (phylum Porifera, class Hexactinellida) contained single nucleotide insertions in protein coding genes inferred as sites of +1 translational frameshifting. To investigate the distribution and evolution of these sites and to help elucidate the mechanism of frameshifting, we determined eight new complete or nearly complete mtDNA sequences from glass sponges and examined individual mitochondrial genes from three others. We found nine new instances of single nucleotide insertions in these sequences and analyzed them both comparatively and phylogenetically. The base insertions appear to have been gained and lost repeatedly in hexactinellid mt protein genes, suggesting no functional significance for the frameshifting sites. A high degree of sequence conservation, the presence of unusual tRNAs, and a distinct pattern of codon usage suggest the "out-of-frame pairing" model of translational frameshifting. Additionally, we provide evidence that relaxed selection pressure on glass sponge mtDNA - possibly a result of their low growth rates and deep-water lifestyle - has allowed frameshift insertions to be tolerated for hundreds of millions of years. Our study provides the first example of a phylogenetically diverse and extensive usage of translational frameshifting in animal mitochondrial coding sequences.