Body mass is less important than bird order in determining the molecular rate for bird mitochondrial DNA

A negative relationship between body mass and molecular evolution rates has been suggested, and recently a correlation equation has been published based on mitochondrial genomic data of 475 bird species and their body masses. Here, we re‐analysed these data and show that the bird order as a proxy of monophyletic groups was a stronger predictor of the molecular rate than the body mass. We provide order‐specific molecular substitution rates. Only three orders (Galliformes, Gruiformes, Pelecaniformes) showed a very clear negative correlation, and specific correlation equations are given for these. The molecular rates of bird orders differed strongly at similar mean body masses, and we suggest that the previously described trend across all birds may arise as smaller species also tend to have characteristic life histories, namely faster turnover of generations, higher fecundity and shorter lifespans.     

Evolutionary rates of mitochondrial genomes correspond to diversification rates and to contemporary species richness in birds and reptiles

Rates of biological diversification should ultimately correspond to rates of genome evolution. Recent studies have compared diversification rates with phylogenetic branch lengths, but incomplete phylogenies hamper such analyses for many taxa. Herein, we use pairwise comparisons of confamilial sauropsid (bird and reptile) mitochondrial DNA (mtDNA) genome sequences to estimate substitution rates. These molecular evolutionary rates are considered in light of the age and species richness of each taxonomic family, using a random-walk speciation–extinction process to estimate rates of diversification. We find the molecular clock ticks at disparate rates in different families and at different genes. For example, evolutionary rates are relatively fast in snakes and lizards, intermediate in crocodilians and slow in turtles and birds. There was also rate variation across genes, where non-synonymous substitution rates were fastest at ATP8 and slowest at CO3. Family-by-gene interactions were significant, indicating that local clocks vary substantially among sauropsids. Most importantly, we find evidence that mitochondrial genome evolutionary rates are positively correlated with speciation rates and with contemporary species richness. Nuclear sequences are poorly represented among reptiles, but the correlation between rates of molecular evolution and species diversification also extends to 18 avian nuclear genes we tested. Thus, the nuclear data buttress our mtDNA findings.     

Mitochondrial rate variation among lineages of passerine birds

The order Passeriformes comprises the majority of extant avian species. Analyses of molecular data have provided important insights into the evolution of this diverse order. However, molecular estimates of the evolutionary and demographic timescales of passerine species have been hindered by a lack of reliable calibrations. This has led to a reliance on the application of standard substitution rates to mitochondrial DNA data, particularly rates estimated from analyses of the gene encoding cytochrome b (CYTB). To investigate patterns of rate variation across passerine lineages, we used a Bayesian phylogenetic approach to analyse the protein‐coding genes of 183 mitochondrial genomes. We found that the most commonly used mitochondrial marker, CYTB, has low variation in rates across passerine lineages. This lends support to its widespread use as a molecular clock in birds. However, we also found that the patterns of among‐lineage rate variation in CYTB are only weakly related to the evolutionary rate of the mitochondrial genome as a whole. Our analyses confirmed the presence of mutational saturation at third codon positions across the protein‐coding genes of the mitochondrial genome, reinforcing the view that these sites should be excluded in studies of deep passerine relationships. The results of our analyses have provided information that will be useful for molecular‐clock studies of passerine evolution.     

Evolution of modern birds revealed by mitogenomics: timing the radiation and origin of major orders

Mitochondrial (mt) genes and genomes are among the major sources of data for evolutionary studies in birds. This places mitogenomic studies in birds at the core of intense debates in avian evolutionary biology. Indeed, complete mt genomes are actively been used to unveil the phylogenetic relationships among major orders, whereas single genes (e.g., cytochrome c oxidase I [COX1]) are considered standard for species identification and defining species boundaries (DNA barcoding). In this investigation, we study the time of origin and evolutionary relationships among Neoaves orders using complete mt genomes. First, we were able to solve polytomies previously observed at the deep nodes of the Neoaves phylogeny by analyzing 80 mt genomes, including 17 new sequences reported in this investigation. As an example, we found evidence indicating that columbiforms and charadriforms are sister groups. Overall, our analyses indicate that by improving the taxonomic sampling, complete mt genomes can solve the evolutionary relationships among major bird groups. Second, we used our phylogenetic hypotheses to estimate the time of origin of major avian orders as a way to test if their diversification took place prior to the Cretaceous/Tertiary (K/T) boundary. Such timetrees were estimated using several molecular dating approaches and conservative calibration points. Whereas we found time estimates slightly younger than those reported by others, most of the major orders originated prior to the K/T boundary. Finally, we used our timetrees to estimate the rate of evolution of each mt gene. We found great variation on the mutation rates among mt genes and within different bird groups. COX1 was the gene with less variation among Neoaves orders and the one with the least amount of rate heterogeneity across lineages. Such findings support the choice of COX 1 among mt genes as target for developing DNA barcoding approaches in birds.     

A Bayesian evaluation of human mitochondrial substitution rates

Accurate estimates of mitochondrial substitution rates are central to molecular studies of human evolution, but meaningful comparisons of published studies are problematic because of the wide range of methodologies and data sets employed. These differences are nowhere more pronounced than among rates estimated from phylogenies, genealogies, and pedigrees. By using a data set comprising mitochondrial genomes from 177 humans, we estimate substitution rates for various data partitions by using Bayesian phylogenetic analysis with a relaxed molecular clock. We compare the effect of multiple internal calibrations with the customary human-chimpanzee split. The analyses reveal wide variation among estimated substitution rates and divergence times made with different partitions and calibrations, with evidence of substitutional saturation, natural selection, and significant rate heterogeneity among lineages and among sites. Collectively, the results support dates for migration out of Africa and the common mitochondrial ancestor of humans that are considerably more recent than most previous estimates. Our results also demonstrate that human mitochondrial genomes exhibit a number of molecular evolutionary complexities that necessitate the use of sophisticated analytical models for genetic analyses.     

Additive Uncorrelated Relaxed Clock Models for the Dating of Genomic Epidemiology Phylogenies

Phylogenetic dating is one of the most powerful and commonly used methods of drawing epidemiological interpretations from pathogen genomic data. Building such trees requires considering a molecular clock model which represents the rate at which substitutions accumulate on genomes. When the molecular clock rate is constant throughout the tree then the clock is said to be strict, but this is often not an acceptable assumption. Alternatively, relaxed clock models consider variations in the clock rate, often based on a distribution of rates for each branch. However, we show here that the distributions of rates across branches in commonly used relaxed clock models are incompatible with the biological expectation that the sum of the numbers of substitutions on two neighboring branches should be distributed as the substitution number on a single branch of equivalent length. We call this expectation the additivity property. We further show how assumptions of commonly used relaxed clock models can lead to estimates of evolutionary rates and dates with low precision and biased confidence intervals. We therefore propose a new additive relaxed clock model where the additivity property is satisfied. We illustrate the use of our new additive relaxed clock model on a range of simulated and real data sets, and we show that using this new model leads to more accurate estimates of mean evolutionary rates and ancestral dates.     

Relaxed clocks and inferences of heterogeneous patterns of nucleotide substitution and divergence time estimates across whales and dolphins (Mammalia: Cetacea)

Various nucleotide substitution models have been developed to accommodate among lineage rate heterogeneity, thereby relaxing the assumptions of the strict molecular clock. Recently developed "uncorrelated relaxed clock" and "random local clock" (RLC) models allow decoupling of nucleotide substitution rates between descendant lineages and are thus predicted to perform better in the presence of lineage-specific rate heterogeneity. However, it is uncertain how these models perform in the presence of punctuated shifts in substitution rate, especially between closely related clades. Using cetaceans (whales and dolphins) as a case study, we test the performance of these two substitution models in estimating both molecular rates and divergence times in the presence of substantial lineage-specific rate heterogeneity. Our RLC analyses of whole mitochondrial genome alignments find evidence for up to ten clade-specific nucleotide substitution rate shifts in cetaceans. We provide evidence that in the uncorrelated relaxed clock framework, a punctuated shift in the rate of molecular evolution within a subclade results in posterior rate estimates that are either misled or intermediate between the disparate rate classes present in baleen and toothed whales. Using simulations, we demonstrate abrupt changes in rate isolated to one or a few lineages in the phylogeny can mislead rate and age estimation, even when the node of interest is calibrated. We further demonstrate how increasing prior age uncertainty can bias rate and age estimates, even while the 95% highest posterior density around age estimates decreases; in other words, increased precision for an inaccurate estimate. We interpret the use of external calibrations in divergence time studies in light of these results, suggesting that rate shifts at deep time scales may mislead inferences of absolute molecular rates and ages.     

Complete mitochondrial DNA genome sequences show that modern birds are not descended from transitional shorebirds

To test the hypothesis put forward by Feduccia of the origin of modern birds from transitional birds, we sequenced the first two complete mitochondrial genomes of shorebirds (ruddy turnstone and blackish oystercatcher) and compared their sequences with those of already published avian genomes. When corrected for rate heterogeneity across sites and non-homogeneous nucleotide compositions among lineages in maximum likelihood (ML), the optimal tree places palaeognath birds as sister to the neognaths including shorebirds. This optimal topology is a re-rooting of recently published ordinal-level avian trees derived from mitochondrial sequences. Using a penalized likelihood (PL) rate-smoothing process in conjunction with dates estimated from fossils, we show that the basal splits in the bird tree are much older than the Cretaceous-Tertiary (K-T) boundary, reinforcing previous molecular studies that rejected the derivation of modern birds from transitional shorebirds. Our mean estimate for the origin of modern birds at about 123 million years ago (Myr ago) is quite close to recent estimates using both nuclear and mitochondrial genes, and supports theories of continental break-up as a driving force in avian diversification. Not only did many modern orders of birds originate well before the K-T boundary, but the radiation of major clades occurred over an extended period of at least 40 Myr ago, thus also falsifying Feduccia's rapid radiation scenario following a K-T bottleneck.     

A bi-organellar phylogenomic study of Pandanales: inference of higher-order relationships and unusual rate-variation patterns

We used a bi-organellar phylogenomic approach to address higher-order relationships in Pandanales, including the first molecular phylogenetic study of the panama-hat family, Cyclanthaceae. Our genus-level study of plastid and mitochondrial gene sets includes a comprehensive sampling of photosynthetic lineages across the order, and provides a framework for investigating clade ages, biogeographic hypotheses and organellar molecular evolution. Using multiple inference methods and both organellar genomes, we recovered mostly congruent and strongly supported relationships within and between families, including the placement of fully mycoheterotrophic Triuridaceae. Cyclanthaceae and Pandanaceae plastomes have slow substitution rates, contributing to weakly supported plastid-based relationships in Cyclanthaceae. While generally slowly evolving, mitochondrial genomes exhibit sporadic rate elevation across the order. However, we infer well-supported relationships even for slower evolving mitochondrial lineages in Cyclanthaceae. Clade age estimates across photosynthetic lineages are largely consistent with previous studies, are well correlated between the two organellar genomes (with slightly younger inferences from mitochondrial data), and support several biogeographic hypotheses. We show that rapidly evolving non-photosynthetic lineages may bias age estimates upwards at neighbouring photosynthetic nodes, even using a relaxed clock model. Finally, we uncovered new genome structural variants in photosynthetic taxa at plastid inverted repeat boundaries that show promise as interfamilial phylogenetic markers.     

Inconsistencies in estimating the age of HIV-1 subtypes due to heterotachy

Rate heterogeneity among lineages is a common feature of molecular evolution, and it has long impeded our ability to accurately estimate the age of evolutionary divergence events. The development of relaxed molecular clocks, which model variable substitution rates among lineages, was intended to rectify this problem. Major subtypes of pandemic HIV-1 group M are thought to exemplify closely related lineages with different substitution rates. Here, we report that inferring the time of most recent common ancestor of all these subtypes in a single phylogeny under a single (relaxed) molecular clock produces significantly different dates for many of the subtypes than does analysis of each subtype on its own. We explore various methods to ameliorate this problem. We conclude that current molecular dating methods are inadequate for dealing with this type of substitution rate variation in HIV-1. Through simulation, we show that heterotachy causes root ages to be overestimated.     

Overdispersion of the molecular clock: temporal variation of gene-specific substitution rates in Drosophila

Simple models of molecular evolution assume that sequences evolve by a Poisson process in which nucleotide or amino acid substitutions occur as rare independent events. In these models, the expected ratio of the variance to the mean of substitution counts equals 1, and substitution processes with a ratio greater than 1 are called overdispersed. Comparing the genomes of 10 closely related species of Drosophila, we extend earlier evidence for overdispersion in amino acid replacements as well as in four-fold synonymous substitutions. The observed deviation from the Poisson expectation can be described as a linear function of the rate at which substitutions occur on a phylogeny, which implies that deviations from the Poisson expectation arise from gene-specific temporal variation in substitution rates. Amino acid sequences show greater temporal variation in substitution rates than do four-fold synonymous sequences. Our findings provide a general phenomenological framework for understanding overdispersion in the molecular clock. Also, the presence of substantial variation in gene-specific substitution rates has broad implications for work in phylogeny reconstruction and evolutionary rate estimation.     

The determinants of the molecular substitution process in turtles

Neutral rates of molecular evolution vary across species, and this variation has been shown to be related to biological traits. One of the first patterns to be observed in vertebrates has been an inverse relationship between body mass (BM) and substitution rates. The effects of three major life‐history traits (LHT) that covary with BM – metabolic rate, generation time and longevity (LON) – have been invoked to explain this relationship. However, most of the theoretical and empirical evidence supporting this relationship comes from endothermic vertebrates, that is, mammals and birds, in which the environmental conditions, especially temperature, do not have a direct impact on cellular and molecular biology. We analysed the variations in mitochondrial and nuclear rates of synonymous substitution across 224 turtle species and examined their correlation with two LHT (LON and BM) and two environmental variables [latitude (LAT) and habitat]. Our analyses indicate that in turtles, neutral rates of molecular evolution are hardly correlated with LON or BM. Rather, both the mitochondrial and nuclear substitution rates are significantly correlated with LAT – faster evolution in the tropics – and especially so for aquatic species. These results question the generality of the relationships reported in mammals and birds and suggest that environmental factors might be the strongest determinants of the mutation rate in ectotherms.     

Exploring patterns and extent of bias in estimating divergence time from mitochondrial DNA sequence data in a particular lineage: a case study of salamanders (order Caudata)

In the practice of molecular dating, substitution saturation will bias the results if not properly modeled. Date estimates based on commonly used mitochondrial DNA sequences likely suffer from this problem because of their high substitution rate. Nevertheless, the patterns and extent of such expected bias remain unknown for many major evolutionary lineages, which often differ in ages, available calibrations, and substitution rates of their mitochondrial genome. In this case study of salamanders, we used estimates based on multiple nuclear exons to assess the effects of saturation on dating divergences using mitochondrial genome sequences on a timescale of ~200-300 My. The results indicated that, due to saturation for older divergences and in the absence of younger effective calibration points, dates derived from the mitochondrial data were considerably overestimated and systematically biased toward the calibration point for the ingroup root. The overestimate might be as great as 3-10 times (about 20 My) older than actual divergence dates for recent splitting events and 40 My older for events that are more ancient. For deep divergences, dates estimated were strongly compressed together. Furthermore, excluding the third codon positions of protein-coding genes or only using the RNA genes or second codon positions did not considerably improve the performance. In the order Caudata, slowly evolving markers such as nuclear exons are preferred for dating a phylogeny covering a relatively wide time span. Dates estimated from these markers can be used as secondary calibrations for dating recent events based on rapidly evolving markers for which mitochondrial DNA sequences are attractive candidates due to their short coalescent time. In other groups, similar evaluation should be performed to facilitate the choice of markers for molecular dating and making inferences from the results.     

A mitogenomic timescale for birds detects variable phylogenetic rates of molecular evolution and refutes the standard molecular clock

Current understanding of the diversification of birds is hindered by their incomplete fossil record and uncertainty in phylogenetic relationships and phylogenetic rates of molecular evolution. Here we performed the first comprehensive analysis of mitogenomic data of 48 vertebrates, including 35 birds, to derive a Bayesian timescale for avian evolution and to estimate rates of DNA evolution. Our approach used multiple fossil time constraints scattered throughout the phylogenetic tree and accounts for uncertainties in time constraints, branch lengths, and heterogeneity of rates of DNA evolution. We estimated that the major vertebrate lineages originated in the Permian; the 95% credible intervals of our estimated ages of the origin of archosaurs (258 MYA), the amniote-amphibian split (356 MYA), and the archosaur-lizard divergence (278 MYA) bracket estimates from the fossil record. The origin of modern orders of birds was estimated to have occurred throughout the Cretaceous beginning about 139 MYA, arguing against a cataclysmic extinction of lineages at the Cretaceous/Tertiary boundary. We identified fossils that are useful as time constraints within vertebrates. Our timescale reveals that rates of molecular evolution vary across genes and among taxa through time, thereby refuting the widely used mitogenomic or cytochrome b molecular clock in birds. Moreover, the 5-Myr divergence time assumed between 2 genera of geese (Branta and Anser) to originally calibrate the standard mitochondrial clock rate of 0.01 substitutions per site per lineage per Myr (s/s/l/Myr) in birds was shown to be underestimated by about 9.5 Myr. Phylogenetic rates in birds vary between 0.0009 and 0.012 s/s/l/Myr, indicating that many phylogenetic splits among avian taxa also have been underestimated and need to be revised. We found no support for the hypothesis that the molecular clock in birds "ticks" according to a constant rate of substitution per unit of mass-specific metabolic energy rather than per unit of time, as recently suggested. Our analysis advances knowledge of rates of DNA evolution across birds and other vertebrates and will, therefore, aid comparative biology studies that seek to infer the origin and timing of major adaptive shifts in vertebrates.     

鞘翅目昆虫线粒体基因组研究进展

鞘翅目（Coleoptera）是世界上最具多样性的类群,具有很高的生态和形态多样性,这些多样性吸引了很多进化生物学家和分类学家的关注。随着分子生物学的发展,分子生物学技术广泛应用于鞘翅目系统学的研究,但随着研究的深入,简单的分子片段已经不能满足研究的需求,需要发掘更新的分子标记。近年来,线粒体全基因组已经成为鞘翅目分子系统学研究中很重要的分子标记之一,并广泛地应用于鞘翅目昆虫各个阶元的研究中。本文就鞘翅目线粒体全基因组的概况、研究进展及存在问题进行了总结和讨论。目前,鞘翅目线粒体基因组的研究主要包括物种线粒体基因组组成与结构、分子系统学和分子进化等方面。线粒体基因组在解决系统发育和进化方面表现出了很多的优越性,然而也存在着一些缺点,如序列难获得、基因类型单一、各基因进化速率不同、应用较局限等。     

Testicular melanization has evolved in birds with high mtDNA mutation rates

Melanin is mainly found in the integument of animals, but it also appears in several extracutaneous tissues. The presence of melanin in testes has been anecdotally reported in all vertebrate groups, but the causes and functions of this melanin remain unknown. Similar to other extracutaneous melanins, testicular melanin may protect male germ cells from oxidative stress. Given the high respiratory activity of spermatozoa, oxidative stress generated by mitochondrial dysfunction as a consequence of mtDNA mutations directly affects sperm viability. Thus, natural selection may favour testicular melanization in males of species with high historical mutation rates in the mitochondrial genome. Here, we tested this hypothesis using information on occurrence of testicular melanization and mutation accumulation as reflected by cytochrome b mtDNA base pair substitution rates in a large set of 134 species of birds, controlling for the confounding effects of body mass, reproductive activity and phylogeny. We found that testicular melanization has evolved in species with high rates of accumulated mitochondrial mutations and propose that this is an adaptive response related to the protective capacity of melanin against oxidative stress. In support of this hypothesis, testicular melanization was more frequently observed during the breeding season of birds (i.e. when spermatogenesis is likely to occur) than during reproductive inactivity. In contrast to other extracutaneous melanins whose abundance seems to reflect skin and coat colour, we did not find a correlation between the proportion of plumage coloured by melanins and occurrence of testicular melanization. Whereas future experimental studies should test these hypotheses, our study highlights for the first time that melanization patterns in animals may evolve as a response to historical mutation rates.     

Bayesian random local clocks,or one rate to rule them all

Background  

Relaxed molecular clock models allow divergence time dating and "relaxed phylogenetic" inference, in which a time tree is estimated in the face of unequal rates across lineages. We present a new method for relaxing the assumption of a strict molecular clock using Markov chain Monte Carlo to implement Bayesian modeling averaging over random local molecular clocks. The new method approaches the problem of rate variation among lineages by proposing a series of local molecular clocks, each extending over a subregion of the full phylogeny. Each branch in a phylogeny (subtending a clade) is a possible location for a change of rate from one local clock to a new one. Thus, including both the global molecular clock and the unconstrained model results, there are a total of 2^2n-2 possible rate models available for averaging with 1, 2, ..., 2n - 2 different rate categories.     

The tempo of genetic evolution in birds: body mass and climate effects

Aim Negative relationships between body mass and substitution rates have previously been reported. However, most of these studies have involved contrasted taxa that, due to their highly divergent phylogenetic histories, also differ in many additional characteristics other than mass. In particular, there has been little examination of the potentially confounding effects of climate or population size. Here we test for differences in rates of microevolution among bird species that, although differing in mass, are nonetheless very closely related phylogenetic pairs. We additionally tested for latitudinal/elevational and population size effects across these contrasts. Location Global. Methods The tempo of microevolution within the cytochrome b gene of mitochondrial DNA was compared between closely related bird species that differed in body mass, using 130 phylogenetically independent species pairs. In order to minimize climate effects, pairs not having overlapping latitudinal ranges were discarded. In addition, a subset of pairs was identified and analysed that involved comparisons between species that have different latitudinal or elevational midpoints. Results Species with smaller mass had substitution rates marginally faster than those with larger mass (small : large median ratio = 1.05). However, this result was only statistically significant when data were pruned to eliminate comparisons in which population or range size also varied substantially between contrasted species. Latitude and elevation had a much stronger association with substitution rates than body mass within the subset of pairs (n = 30) that also differed in their spatial distributions: lower elevation or latitude species had substantially more substitutions than those at higher latitudes or elevations (low : high ratio = 1.35). Furthermore, when the dataset was pruned of pairs in which body mass was confounded by latitude or elevation, the body mass effect was eliminated. Main conclusions Body mass is known to correlate with latitude, so that the latitudinal/elevational association with microevolution we found might either be additive to, or causal of, the body mass effect. These results are consistent with the evolutionary speed hypothesis, which suggests that latitudinal diversity gradients derive from variation in the rate of microevolution. Our findings also serve to raise concerns about biogeographical studies that use genetic distances between taxa to estimate time since divergence.     

The complete mitochondrial genome of the bird schistosome Trichobilharzia regenti (Platyhelminthes: Digenea), causative agent of cercarial dermatitis

The complete mitochondrial (mt) genome of the neuropathogenic bird schistosome Trichobilharzia regenti was fully sequenced in order to develop molecular markers for future diagnostic, molecular ecological, population, and phylogenetic studies. The genome was 14,838 bp in length, with a 68.4% AT bias in protein coding regions. A repeat element (3 x 184 bp) between trnV and trnW distinguished a single short noncoding region. As 9 of 14 genera of schistosomes parasitize birds, future characterization of their mt genomes is desirable for species-specific and strain- or population-specific diagnostic markers; this concerns not only the nasal representatives, e.g., T. regenti characterized in this study, but also numerous species within the predominant group of visceral (blood dwelling) bird schistosomes.     

Phylogenomic datasets provide both precision and accuracy in estimating the timescale of placental mammal phylogeny

The fossil record suggests a rapid radiation of placental mammals following the Cretaceous-Paleogene (K-Pg) mass extinction 65 million years ago (Ma); nevertheless, molecular time estimates, while highly variable, are generally much older. Early molecular studies suffer from inadequate dating methods, reliance on the molecular clock, and simplistic and over-confident interpretations of the fossil record. More recent studies have used Bayesian dating methods that circumvent those issues, but the use of limited data has led to large estimation uncertainties, precluding a decisive conclusion on the timing of mammalian diversifications. Here we use a powerful Bayesian method to analyse 36 nuclear genomes and 274 mitochondrial genomes (20.6 million base pairs), combined with robust but flexible fossil calibrations. Our posterior time estimates suggest that marsupials diverged from eutherians 168-178 Ma, and crown Marsupialia diverged 64-84 Ma. Placentalia diverged 88-90 Ma, and present-day placental orders (except Primates and Xenarthra) originated in a ～20 Myr window (45-65 Ma) after the K-Pg extinction. Therefore we reject a pre K-Pg model of placental ordinal diversification. We suggest other infamous instances of mismatch between molecular and palaeontological divergence time estimates will be resolved with this same approach.