Molecular phylogenies link rates of evolution and speciation

Combining Molecular, Macroevolutionary, and Macroecological Perspectives on the Generation of Diversity

Charles Darwin presented a unified process of diversification driven by the gradual accumulation of heritable variation. The growth in DNA databases and the increase in genomic sequencing, combined with advances in molecular phylogenetic analyses, gives us an opportunity to realize Darwin's vision, connecting the generation of variation to the diversification of lineages. The rate of molecular evolution is correlated with the rate of diversification across animals and plants, but the relationship between genome change and speciation is complex: Mutation rates evolve in response to life history and niche; substitution rates are influenced by mutation, selection, and population size; rates of acquisition of reproductive isolation vary between populations; and traits, niches, and distribution can influence diversification rates. The connection between mutation rate and diversification rate is one part of the complex and varied story of speciation, which has theoretical importance for understanding the generation of biodiversity and also practical impacts on the use of DNA to understand the dynamics of speciation over macroevolutionary timescales.

Punctuational evolution is pervasive in distal site metastatic colonization

The evolution of metastasis represents a lethal stage of cancer progression. Yet, the evolutionary kinetics of metastatic disease remain unresolved. Here, using single cell CRISPR-Cas9 lineage tracing data, we show that in metastatic disease, gradual molecular evolution is punctuated by episodes of rapid evolutionary change associated with lineage divergence. By measuring punctuational effects across the metastatic cascade, we show that punctuational effects contribute more to the molecular diversity at distal site metastases compared to the paired primary tumor, suggesting qualitatively different modes of evolution may drive primary and metastatic tumor progression. This is the first empirical evidence for distinct patterns of molecular evolution at early and late stages of metastasis and demonstrates the complex interplay of cell intrinsic and extrinsic factors that shape lethal cancer.

HOST SWITCHING IN DICYEMIDS (PHYLUM DICYEMIDA)

Dicyemids (phylum Dicyemida) are the most common and most characteristic endosymbionts in the renal sacs of benthic cephalopod molluscs: octopuses and cuttlefishes. Typically, 2 or 3 dicyemid species are found in a single specimen of the host, and most dicyemids have high host specificity. Host-specific parasites are restricted to a limited range of host species by ecological barriers that impede dispersal and successful establishment; therefore, phylogenies of interacting groups are often congruent due to repeated co-speciation. Most frequently, however, host and parasite phylogenies are not congruent, which can be explained by processes such as host switching and other macro-evolutionary events. Here, the history of dicyemids and their host cephalopod associations were studied by comparing their phylogenies. Dicyemid species were collected from 8 decapodiform species and 12 octopodiform species in Japanese waters. Using whole mitochondrial cytochrome c oxidase subunit 1 (COI) sequences, a phylogeny of 37 dicyemid species, including 4 genera representing the family Dicyemidae, was reconstructed. Phylogenetic trees derived from analyses of COI genes consistently suggested that dicyemid species should be separated into 3 major clades and that the most common genera, Dicyema and Dicyemennea, are not monophyletic. Thus, morphological classification does not reflect the phylogenetic relationships of these 2 genera. Divergence (speciation) of dicyemid species seems to have occurred within a single host species. Possible host-switching events may have occurred between the Octopodiformes and Decapodiformes or within the Octopodiformes or the Decapodiformes. Therefore, the mechanism of dicyemid speciation may be a mixture of host switching and intra-host speciation. This is the first study in which the process of dicyemid diversification involving cephalopod hosts has been evaluated with a large number of dicyemid species and genera.

Stability and Assembly Mechanisms of Butterfly Communities across Environmental Gradients of a Subtropical Mountain

Mountain ecosystems harbor evolutionarily unique and exceptionally rich biodiversity, particularly in insects. In this study, we characterized the diversity, community stability, and assembly mechanisms of butterflies on a subtropical mountain in the Chebaling National Nature Reserve, Guangdong Province, China, using grid-based monitoring across the entire region for two years. The results showed that species richness, abundance, and Faith's phylogenetic diversity decreased with increasing elevation; taxonomic diversity played a considerable role in mediating the effects of environmental changes on stability. Moreover, our results showed that stochastic processes are dominant in governing the assembly of butterfly communities across all elevational gradients, with habitats at an elevation of 416-580 m subjected to the strongest stochastic processes, whereas heterogeneous selection processes displayed stronger effects on the assembly of butterfly communities at 744-908 m, 580-744 m, and 908-1072 m, with abiotic factors inferred as the main driving forces. In addition, significant differences were detected between the barcode tree and the placement tree for the calculated β-NTI values at 416-580 m. Overall, this study provides new insights into the effects of environmental change on the stability and assembly of butterflies in Chebaling, which will be beneficial for biodiversity conservation and policy development.

Detecting punctuated evolution in SARS-CoV-2 over the first year of the pandemic

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) evolved slowly over the first year of the Coronavirus Disease 19 (COVID-19) pandemic with differential mutation rates across lineages. Here, we explore how this variation arose. Whether evolutionary change accumulated gradually within lineages or during viral lineage branching is unclear. Using phylogenetic regression models, we show that ~13% of SARS-CoV-2 genomic divergence up to May 2020 is attributable to lineage branching events (punctuated evolution). The net number of branching events along lineages predicts ~5% of the deviation from the strict molecular clock. We did not detect punctuated evolution in SARS-CoV-1, possibly due to the small sample size, and in sarbecovirus broadly, likely due to a different evolutionary process altogether. Punctuation in SARS-CoV-2 is probably neutral because most mutations were not positively selected and because the strength of the punctuational effect remained constant over time, at least until May 2020, and across continents. However, the small punctuational contribution to SARS-CoV-2 diversity is consistent with the founder effect arising from narrow transmission bottlenecks. Therefore, punctuation in SARS-CoV-2 may represent the macroevolutionary consequence (rate variation) of a microevolutionary process (transmission bottleneck).

Nucleotide Substitutions during Speciation may Explain Substitution Rate Variation

Although molecular mechanisms associated with the generation of mutations are highly conserved across taxa, there is widespread variation in mutation rates between evolutionary lineages. When phylogenies are reconstructed based on nucleotide sequences, such variation is typically accounted for by the assumption of a relaxed molecular clock, which is a statistical distribution of mutation rates without much underlying biological mechanism. Here, we propose that variation in accumulated mutations may be partly explained by an elevated mutation rate during speciation. Using simulations, we show how shifting mutations from branches to speciation events impacts inference of branching times in phylogenetic reconstruction. Furthermore, the resulting nucleotide alignments are better described by a relaxed than by a strict molecular clock. Thus, elevated mutation rates during speciation potentially explain part of the variation in substitution rates that is observed across the tree of life. [Molecular clock; phylogenetic reconstruction; speciation; substitution rate variation.].

Investigating the reliability of molecular estimates of evolutionary time when substitution rates and speciation rates vary

Background

An accurate timescale of evolutionary history is essential to testing hypotheses about the influence of historical events and processes, and the timescale for evolution is increasingly derived from analysis of DNA sequences. But variation in the rate of molecular evolution complicates the inference of time from DNA. Evidence is growing for numerous factors, such as life history and habitat, that are linked both to the molecular processes of mutation and fixation and to rates of macroevolutionary diversification. However, the most widely used methods rely on idealised models of rate variation, such as the uncorrelated and autocorrelated clocks, and molecular dating methods are rarely tested against complex models of rate change. One relationship that is not accounted for in molecular dating is the potential for interaction between molecular substitution rates and speciation, a relationship that has been supported by empirical studies in a growing number of taxa. If these relationships are as widespread as current evidence suggests, they may have a significant influence on molecular dates.

Results

We simulate phylogenies and molecular sequences under three different realistic rate variation models—one in which speciation rates and substitution rates both vary but are unlinked, one in which they covary continuously and one punctuated model in which molecular change is concentrated in speciation events, using empirical case studies to parameterise realistic simulations. We test three commonly used “relaxed clock” molecular dating methods against these realistic simulations to explore the degree of error in molecular dates under each model. We find average divergence time inference errors ranging from 12% of node age for the unlinked model when reconstructed under an uncorrelated rate prior using BEAST 2, to up to 91% when sequences evolved under the punctuated model are reconstructed under an autocorrelated prior using PAML.

Conclusions

We demonstrate the potential for substantial errors in molecular dates when both speciation rates and substitution rates vary between lineages. This study highlights the need for tests of molecular dating methods against realistic models of rate variation generated from empirical parameters and known relationships.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02015-8.

Diversification Rate is Associated with Rate of Molecular Evolution in Ray-Finned Fish (Actinopterygii)

Understanding the factors that drive diversification of taxa across the tree of life is a key focus of macroevolutionary research. While the effects of life history, ecology, climate and geography on diversity have been studied for many taxa, the relationship between molecular evolution and diversification has received less attention. However, correlations between rates of molecular evolution and diversification rate have been detected in a range of taxa, including reptiles, plants and birds. A correlation between rates of molecular evolution and diversification rate is a prediction of several evolutionary theories, including the evolutionary speed hypothesis which links variation in mutation rates to differences in speciation rates. If it is widespread, such correlations could also have significant practical impacts, if they are not adequately accounted for in phylogenetic inference of evolutionary rates and timescales. Ray-finned fish (Actinopterygii) offer a prime target to test for this relationship due to their extreme variation in clade size suggesting a wide range of diversification rates. We employ both a sister-pairs approach and a whole-tree approach to test for correlations between substitution rate and net diversification. We also collect life history and ecological trait data and account for potential confounding factors including body size, latitude, max depth and reef association. We find evidence to support a relationship between diversification and synonymous rates of nuclear evolution across two published backbone phylogenies, as well as weak evidence for a relationship between mitochondrial nonsynonymous rates and diversification at the genus level.

Comparative genomics of tadpole shrimps (Crustacea, Branchiopoda, Notostraca): Dynamic genome evolution against the backdrop of morphological stasis

This analysis presents five genome assemblies of four Notostraca taxa. Notostraca origin dates to the Permian/Upper Devonian and the extant forms show a striking morphological similarity to fossil taxa. The comparison of sequenced genomes with other Branchiopoda genomes shows that, despite the morphological stasis, Notostraca share a dynamic genome evolution with high turnover for gene families' expansion/contraction and a transposable elements content comparable to other branchiopods. While Notostraca substitutions rate appears similar or lower in comparison to other branchiopods, a subset of genes shows a faster evolutionary pace, highlighting the difficulty of generalizing about genomic stasis versus dynamism. Moreover, we found that the variation of Triops cancriformis transposable elements content appeared linked to reproductive strategies, in line with theoretical expectations. Overall, besides providing new genomic resources for the study of these organisms, which appear relevant for their ecology and evolution, we also confirmed the decoupling of morphological and molecular evolution.

Hb adaptation to hypoxia in high-altitude fishes: Fresh evidence from schizothoracinae fishes in the Qinghai-Tibetan Plateau

Uncovering the genetic basis of hypoxic adaptation is one of the most active research areas in evolutionary biology. Among air-breathing vertebrates, modifications of hemoglobin (Hb) play a pivotal role in mediating an adaptive response to high-altitude hypoxia. However, the relative contributions in water-breathing organisms are still unclear. Here, we tested the Hb concentration of fish at different altitudes. All species showed species-specific Hb concentration, which has a non-positive correlation with altitude. Moreover, we investigated the expression of Hb genes by the RNA-seq and quantitative real-time PCR (qRT-PCR), and Hb composition by two-dimensional electrophoresis (2-DE). The results showed that the multiple Hb genes and isoforms are co-expressed in schizothoracinae fishes endemic to the Qinghai-Tibetan Plateau (QTP). Phylogenetic analyses of Hb genes indicated that the evolutionary relationships are not easily reconciled with the organismal phylogeny. Furthermore, evidence of positive selection was found in the Hb genes of schizothoracinae fishes through the selection pressure analysis. We demonstrated that positively selected sites likely facilitated the functional divergence of Hb isoforms. Taken together, this study indicated that the long-term maintenance of high Hb concentration may be a disadvantage for physiologically acclimating to high altitude hypoxia. Meanwhile, the genetically based modification of Hb-O₂ affinity in schizothoracinae fishes might facilitate the evolutionary adaptation to Tibetan aqueous environments.

Why are there not more herbivorous insect species?

Insect species richness is estimated to exceed three million species, of which roughly half is herbivorous. Despite the vast number of species and varied life histories, the proportion of herbivorous species among plant-consuming organisms is lower than it could be due to constraints that impose limits to their diversification. These include ecological factors, such as vague interspecific competition; anatomical and physiological limits, such as neural limits and inability of handling a wide range of plant allelochemicals; phylogenetic constraints, like niche conservatism; and most importantly, a low level of concerted genetic variation necessary to a phyletic conversion. It is suggested that diversification ultimately depends on what we call the intrinsic trend of diversification of the insect genome. In support of the above, we survey the major types of host-specificity, the mechanisms and constraints of host specialization, possible pathways of speciation, and hypotheses concerning insect diversification.

Genetic insights into the evolution of genera with the eastern Asia-eastern North America floristic disjunction: a transcriptomics analysis

PREMISE

Large disjunctions in species distributions provide excellent opportunities to study processes that shape biogeographic patterns. One such disjunction is the eastern Asia-eastern North America (EA-ENA) floristic disjunction. For many genera with this disjunction, species richness is greater in EA than in ENA; this pattern has been attributed, in part, to higher rates of molecular evolution and speciation in EA. Longer branch lengths have been found in some EA clades, relative to their ENA sister clades, suggesting that the EA lineages have evolved at a higher rate, possibly due to environmental heterogeneity, potentially contributing to the species richness anomaly.

METHODS

To evaluate whether rates of molecular evolution are elevated in EA relative to ENA, we used transcriptomes from species in 11 genera displaying this disjunction. Rates of molecular evolution were estimated for up to 385 orthologous nuclear loci per genus.

RESULTS

No statistically significant differences were identified in pairwise comparisons between EA and ENA sister species, suggesting equal rates of molecular evolution for both species; the data also suggest similar selection pressures in both regions. For larger genera, evidence likewise argues against more species-rich clades having higher molecular evolutionary rates, regardless of region. Our results suggest that genes across multiple gene ontology categories are evolving at similar rates under purifying selection in species in both regions.

CONCLUSIONS

Our data support the hypothesis that greater species richness in EA than ENA is due to factors other than an overall increase in rates of molecular evolution in EA.

Model-Based Inference of Punctuated Molecular Evolution

In standard models of molecular evolution, DNA sequences evolve through asynchronous substitutions according to Poisson processes with a constant rate (called the molecular clock) or a rate that can vary (relaxed clock). However, DNA sequences can also undergo episodes of fast divergence that will appear as synchronous substitutions affecting several sites simultaneously at the macroevolutionary timescale. Here, we develop a model, which we call the Relaxed Clock with Spikes model, combining basal, clock-like molecular substitutions with episodes of fast divergence called spikes arising at speciation events. Given a multiple sequence alignment and its time-calibrated species phylogeny, our model is able to detect speciation events (including hidden ones) cooccurring with spike events and to estimate the probability and amplitude of these spikes on the phylogeny. We identify the conditions under which spikes can be distinguished from the natural variance of the clock-like component of molecular substitutions and from variations of the clock. We apply the method to genes underlying snake venom proteins and identify several spikes at gene-specific locations in the phylogeny. This work should pave the way for analyses relying on whole genomes to inform on modes of species diversification.

Genomic insights into adaptation to heterogeneous environments for the ancient relictual Circaeaster agrestis (Circaeasteraceae, Ranunculales)

Investigating the interaction between environmental heterogeneity and local adaptation is critical for understanding the evolutionary history of a species, providing the premise for studying the response of organisms to rapid climate change. However, for most species how exactly the spatial heterogeneity promotes population divergence and how genomic variations contribute to adaptive evolution remain poorly understood. We examine the contributions of geographical and environmental variables to population divergence of the relictual, alpine herb Circaeaster agrestis, as well as the genetic basis of local adaptation using RAD-seq and plastome data. We detected significant genetic structure with an extraordinary disequilibrium of genetic diversity among regions, and signals of isolation-by-distance along with isolation-by-resistance. The populations were estimated to begin diverging in the late Miocene, along with a possible ancestral distribution of the Hengduan Mountains and adjacent regions. Both environmental gradient and redundancy analyses revealed significant association between genetic variation and temperature variables. Genome-environment association analyses identified 16 putatively adaptive loci related mainly to biotic and abiotic stress resistance. Our genome-wide data provide new insights into the important role of environmental heterogeneity in shaping genetic structure, and access the footprints of local adaptation in an ancient relictual species, informing future conservation efforts.

The tuatara genome reveals ancient features of amniote evolution

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

Patterns, Mechanisms and Genetics of Speciation in Reptiles and Amphibians

In this contribution, the aspects of reptile and amphibian speciation that emerged from research performed over the past decade are reviewed. First, this study assesses how patterns and processes of speciation depend on knowing the taxonomy of the group in question, and discuss how integrative taxonomy has contributed to speciation research in these groups. This study then reviews the research on different aspects of speciation in reptiles and amphibians, including biogeography and climatic niches, ecological speciation, the relationship between speciation rates and phenotypic traits, and genetics and genomics. Further, several case studies of speciation in reptiles and amphibians that exemplify many of these themes are discussed. These include studies of integrative taxonomy and biogeography in South American lizards, ecological speciation in European salamanders, speciation and phenotypic evolution in frogs and lizards. The final case study combines genomics and biogeography in tortoises. The field of amphibian and reptile speciation research has steadily moved forward from the assessment of geographic and ecological aspects, to incorporating other dimensions of speciation, such as genetic mechanisms and evolutionary forces. A higher degree of integration among all these dimensions emerges as a goal for future research.

Mitochondrial genome diversity and evolution in Branchiopoda (Crustacea)

BACKGROUND

The crustacean class Branchiopoda includes fairy shrimps, clam shrimps, tadpole shrimps, and water fleas. Branchiopods, which are well known for their great variety of reproductive strategies, date back to the Cambrian and extant taxa can be mainly found in freshwater habitats, also including ephemeral ponds. Mitochondrial genomes of the notostracan taxa Lepidurus apus lubbocki (Italy), L. arcticus (Iceland) and Triops cancriformis (an Italian and a Spanish population) are here characterized for the first time and analyzed together with available branchiopod mitogenomes.

RESULTS

Overall, branchiopod mitogenomes share the basic structure congruent with the ancestral Pancrustacea model. On the other hand, rearrangements involving tRNAs and the control region are observed among analyzed taxa. Remarkably, an unassigned region in the L. apus lubbocki mitogenome showed a chimeric structure, likely resulting from a non-homologous recombination event between the two flanking trnC and trnY genes. Notably, Anostraca and Onychocaudata mitogenomes showed increased GC content compared to both Notostraca and the common ancestor, and a significantly higher substitution rate, which does not correlate with selective pressures, as suggested by dN/dS values.

CONCLUSIONS

Branchiopod mitogenomes appear rather well-conserved, although gene rearrangements have occurred. For the first time, it is reported a putative non-homologous recombination event involving a mitogenome, which produced a pseudogenic tRNA sequence. In addition, in line with data in the literature, we explain the higher substitution rate of Anostraca and Onychocaudata with the inferred GC substitution bias that occurred during their evolution.

The use of DNA barcodes to estimate phylogenetic diversity in forest communities of southern China

To elucidate potential ecological and evolutionary processes associated with the assembly of plant communities, there is now widespread use of estimates of phylogenetic diversity that are based on a variety of DNA barcode regions and phylogenetic construction methods. However, relatively few studies consider how estimates of phylogenetic diversity may be influenced by single DNA barcodes incorporated into a sequence matrix (conservative regions vs. hypervariable regions) and the use of a backbone family-level phylogeny. Here, we use general linear mixed-effects models to examine the influence of different combinations of core DNA barcodes (rbcL, matK, ITS, and ITS2) and phylogeny construction methods on a series of estimates of community phylogenetic diversity for two subtropical forest plots in Guangdong, southern China. We ask: (a) What are the relative influences of single DNA barcodes on estimates phylogenetic diversity metrics? and (b) What is the effect of using a backbone family-level phylogeny to estimate topology-based phylogenetic diversity metrics? The combination of more than one barcode (i.e., rbcL + matK + ITS) and the use of a backbone family-level phylogeny provided the most parsimonious explanation of variation in estimates of phylogenetic diversity. The use of a backbone family-level phylogeny showed a stronger effect on phylogenetic diversity metrics that are based on tree topology compared to those that are based on branch lengths. In addition, the variation in the estimates of phylogenetic diversity that was explained by the top-rank models ranged from 0.1% to 31% and was dependent on the type of phylogenetic community structure metric. Our study underscores the importance of incorporating a multilocus DNA barcode and the use of a backbone family-level phylogeny to infer phylogenetic diversity, where the type of DNA barcode employed and the phylogenetic construction method used can serve as a significant source of variation in estimates of phylogenetic community structure.

When Darwin's Special Difficulty Promotes Diversification in Insects

Eusociality, Darwin's special difficulty, has been widely investigated but remains a topic of great debate in organismal biology. Eusocial species challenge existing theories, and the impact of highly integrated societies on diversification dynamics is controversial with opposing assertions and hypotheses in the literature. Here, using phylogenetic approaches in termites-the first group that has evolved eusociality-we assessed the fundamental prediction that eusocial lineages have higher diversification rates than non-eusocial clades. We found multiple lines of evidence that eusociality provided higher diversification as compared to non-eusociality. This is particularly exacerbated for eusocial species with "true" workers as compared to species with "false" workers. Because most species with "true" workers have an entirely prokaryotic microbiota, the latter feature is also related to higher diversification rates, but it should be investigated further, notably in relation to angiosperm diversification. Overall, this study suggests that societies with "true" workers are not only more successful at ecological timescales but also over millions of years, which further implies that both organism- and species-level traits act on species selection.

Comparative Genomics Reveals Accelerated Evolution in Conserved Pathways during the Diversification of Anole Lizards

Squamates include all lizards and snakes, and display some of the most diverse and extreme morphological adaptations among vertebrates. However, compared with birds and mammals, relatively few resources exist for comparative genomic analyses of squamates, hampering efforts to understand the molecular bases of phenotypic diversification in such a speciose clade. In particular, the ∼400 species of anole lizard represent an extensive squamate radiation. Here, we sequence and assemble the draft genomes of three anole species-Anolis frenatus, Anolis auratus, and Anolis apletophallus-for comparison with the available reference genome of Anolis carolinensis. Comparative analyses reveal a rapid background rate of molecular evolution consistent with a model of punctuated equilibrium, and strong purifying selection on functional genomic elements in anoles. We find evidence for accelerated evolution in genes involved in behavior, sensory perception, and reproduction, as well as in genes regulating limb bud development and hindlimb specification. Morphometric analyses of anole fore and hindlimbs corroborated these findings. We detect signatures of positive selection across several genes related to the development and regulation of the forebrain, hormones, and the iguanian lizard dewlap, suggesting molecular changes underlying behavioral adaptations known to reinforce species boundaries were a key component in the diversification of anole lizards.

Fast Genes and Slow Clades: Comparative Rates of Molecular Evolution in Mammals

Although interest in the rate of molecular evolution and the molecular clock remains high, our knowledge for most groups in these areas is derived largely from a patchwork of studies limited in both their taxon coverage and the number of genes examined. Using a comprehensive molecular data set of 44 genes (18 nDNA, 11 tRNA and 15 additional mtDNA genes) together with a virtually complete and dated phylogeny of extant mammals, I 1) describe differences in the rate of molecular evolution (i.e. substitution rate) within this group in an explicit phylogenetic and quantitative framework and 2) present the first attempt to localize the phylogenetic positions of any rate shifts. Significant rate differences were few and confirmed several long-held trends, including a progressive rate slowdown within hominids and a reduced substitution rate within Cetacea. However, many new patterns were also uncovered, including the mammalian orders being characterized generally by basal rate slowdowns. A link between substitution rate and the size of a clade (which derives from its net speciation rate) is also suggested, with the species-poor major clades (“orders”) showing more decreased rates that often extend throughout the entire clade. Significant rate increases were rare, with the rates within (murid) rodents being fast, but not significantly so with respect to other mammals as a whole. Despite clear lineage-specific differences, rates generally change gradually along these lineages, supporting the potential existence of a local molecular clock in mammals. Together, these results will lay the foundation for a broad-scale analysis to establish the correlates and causes of the rate of molecular evolution in mammals.

Evolutionary dynamics of language systems

Understanding how and why language subsystems differ in their evolutionary dynamics is a fundamental question for historical and comparative linguistics. One key dynamic is the rate of language change. While it is commonly thought that the rapid rate of change hampers the reconstruction of deep language relationships beyond 6,000-10,000 y, there are suggestions that grammatical structures might retain more signal over time than other subsystems, such as basic vocabulary. In this study, we use a Dirichlet process mixture model to infer the rates of change in lexical and grammatical data from 81 Austronesian languages. We show that, on average, most grammatical features actually change faster than items of basic vocabulary. The grammatical data show less schismogenesis, higher rates of homoplasy, and more bursts of contact-induced change than the basic vocabulary data. However, there is a core of grammatical and lexical features that are highly stable. These findings suggest that different subsystems of language have differing dynamics and that careful, nuanced models of language change will be needed to extract deeper signal from the noise of parallel evolution, areal readaptation, and contact.

Phylogenetic estimates of diversification rate are affected by molecular rate variation

Molecular phylogenies are increasingly being used to investigate the patterns and mechanisms of macroevolution. In particular, node heights in a phylogeny can be used to detect changes in rates of diversification over time. Such analyses rest on the assumption that node heights in a phylogeny represent the timing of diversification events, which in turn rests on the assumption that evolutionary time can be accurately predicted from DNA sequence divergence. But there are many influences on the rate of molecular evolution, which might also influence node heights in molecular phylogenies, and thus affect estimates of diversification rate. In particular, a growing number of studies have revealed an association between the net diversification rate estimated from phylogenies and the rate of molecular evolution. Such an association might, by influencing the relative position of node heights, systematically bias estimates of diversification time. We simulated the evolution of DNA sequences under several scenarios where rates of diversification and molecular evolution vary through time, including models where diversification and molecular evolutionary rates are linked. We show that commonly used methods, including metric-based, likelihood and Bayesian approaches, can have a low power to identify changes in diversification rate when molecular substitution rates vary. Furthermore, the association between the rates of speciation and molecular evolution rate can cause the signature of a slowdown or speedup in speciation rates to be lost or misidentified. These results suggest that the multiple sources of variation in molecular evolutionary rates need to be considered when inferring macroevolutionary processes from phylogenies.

Correlated evolutionary rates across genomic compartments in Annonaceae

The molecular clock hypothesis is an important concept in biology. Deviations from a constant rate of nucleotide substitution have been found widely among lineages, genomes, genes and individual sites. Phylogenetic research can accommodate for these differences in applying specific models of evolution. Lineage-specific rate heterogeneity however can generate bi- or multimodal distributions of substitution rates across the branches of a tree and this may mislead phylogenetic inferences with currently available models. The plant family Annonaceae is an excellent case to study lineage-specific rate heterogeneity. The two major sister subfamilies, Annonoideae and Malmeoideae, have shown great discrepancies in branch lengths. We used high-throughput sequencing data of 72 genes, 99 spacers and 16 introns from 24 chloroplast genomes and nuclear ribosomal DNA of 23 species to study the molecular rate of evolution in Annonaceae. In all analyses, longer branch lengths and/or higher substitution rates were found for the Annonoideae compared to the Malmeoideae. The Annonaceae had wide variability in chloroplast length, ranging from minimal 175,684bp to 201,723 for Annonoideae and minimal 152,357 to 170,985bp in Malmeoideae, mostly reflecting variation in inverted-repeat length. The Annonoideae showed a higher GC-content in the conserved parts of the chloroplast genome and higher omega (d_N/d_S)-ratios than the Malmeoideae, which could indicate less stringent purifying selection, a pattern that has been found in groups with small population sizes. This study generates new insights into the processes causing lineage-specific rate heterogeneity, which could lead to improved phylogenetic methods.

SFREEMAP - A simulation-free tool for stochastic mapping

BACKGROUND

Stochastic mapping is frequently used in comparative biology to simulate character evolution, enabling the probabilistic computation of statistics such as number of state transitions along a tree and distribution of states in its internal nodes. Common implementations rely on Continuous-time Markov Chain simulations whose parameters are difficult to adjust and subjected to inherent inaccuracy. Thus, researchers must run a large number of simulations in order to obtain adequate estimates. Although execution time tends to be relatively small when simulations are performed on a single tree assumed to be the "true" topology, it may become an issue if analyses are conducted on several trees, such as the ones that make up posterior distributions obtained via Bayesian phylogenetic inference. Working with such distributions is preferable to working with a single tree, for they allow the integration of phylogenetic uncertainty into parameter estimation. In such cases, detailed character mapping becomes less important than parameter integration across topologies. Here, we present an R-based implementation (SFREEMAP) of an analytical approach to obtain accurate, per-branch expectations of numbers of state transitions and dwelling times. We also introduce an intuitive way of visualizing the results by integrating over the posterior distribution and summarizing the parameters onto a target reference topology (such as a consensus or MAP tree) provided by the user.

RESULTS

We benchmarked SFREEMAP's performance against make.simmap, a popular R-based implementation of stochastic mapping. SFREEMAP confirmed theoretical expectations outperforming make.simmap in every experiment and reducing computation time of relatively modest datasets from hours to minutes. We have also demonstrated that SFREEMAP returns estimates which were not only similar to the ones obtained by averaging across make.simmap mappings, but also more accurate, according to simulated data. We illustrate our visualization strategy using previously published data on the evolution of coloniality in scleractinian corals.

CONCLUSION

SFREEMAP is an accurate and fast alternative to ancestral state reconstruction via simulation-based stochastic mapping.

Ecology has contrasting effects on genetic variation within species versus rates of molecular evolution across species in water beetles

Comparative analysis is a potentially powerful approach to study the effects of ecological traits on genetic variation and rate of evolution across species. However, the lack of suitable datasets means that comparative studies of correlates of genetic traits across an entire clade have been rare. Here, we use a large DNA-barcode dataset (5062 sequences) of water beetles to test the effects of species ecology and geographical distribution on genetic variation within species and rates of molecular evolution across species. We investigated species traits predicted to influence their genetic characteristics, such as surrogate measures of species population size, latitudinal distribution and habitat types, taking phylogeny into account. Genetic variation of cytochrome oxidase I in water beetles was positively correlated with occupancy (numbers of sites of species presence) and negatively with latitude, whereas substitution rates across species depended mainly on habitat types, and running water specialists had the highest rate. These results are consistent with theoretical predictions from nearly-neutral theories of evolution, and suggest that the comparative analysis using large databases can give insights into correlates of genetic variation and molecular evolution.

Exploring the Relationships between Mutation Rates, Life History, Genome Size, Environment, and Species Richness in Flowering Plants

A new view is emerging of the interplay between mutation at the genomic level, substitution at the population level, and diversification at the lineage level. Many studies have suggested that rate of molecular evolution is linked to rate of diversification, but few have evaluated competing hypotheses. By analyzing sequences from 130 families of angiosperms, we show that variation in the synonymous substitution rate is correlated among genes from the mitochondrial, chloroplast, and nuclear genomes and linked to differences in traits among families (average height and genome size). Within each genome, synonymous rates are correlated to nonsynonymous substitution rates, suggesting that increasing the mutation rate results in a faster rate of genome evolution. Substitution rates are correlated with species richness in protein-coding sequences from the chloroplast and nuclear genomes. These data suggest that species traits contribute to lineage-specific differences in the mutation rate that drive both synonymous and nonsynonymous rates of change across all three genomes, which in turn contribute to greater rates of divergence between populations, generating higher rates of diversification. These observations link mutation in individuals to population-level processes and to patterns of lineage divergence.

Clock model makes a large difference to age estimates of long-stemmed clades with no internal calibration: a test using Australian grasstrees

BACKGROUND

Estimating divergence times in phylogenies using a molecular clock depends on accurate modeling of nucleotide substitution rates in DNA sequences. Rate heterogeneity among lineages is likely to affect estimates, especially in lineages with long stems and short crowns ("broom" clades) and no internal calibration. We evaluate the performance of the random local clocks model (RLC) and the more routinely employed uncorrelated lognormal relaxed clock model (UCLN) in situations in which a significant rate shift occurs on the stem branch of a broom clade. We compare the results of simulations to empirical results from analyses of a real rate-heterogeneous taxon - Australian grass trees (Xanthorrhoea) - whose substitution rate is slower than in its sister groups, as determined by relative rate tests.

RESULTS

In the simulated datasets, the RLC model performed much better than UCLN: RLC correctly estimated the age of the crown node of slow-rate broom clades, whereas UCLN estimates were consistently too young. Similarly, in the Xanthorrhoea dataset, UCLN returned significantly younger crown ages than RLC (mean estimates respectively 3-6 Ma versus 25-35 Ma). In both real and simulated datasets, Bayes Factor tests strongly favored the RLC model over the UCLN model.

CONCLUSIONS

The choice of an unsuitable molecular clock model can strongly bias divergence time estimates. In particular, for data predicted to have more rate variation among than within clades, dating with RLC is much more likely to be accurate than with UCLN. The choice of clocks should be informed by the biology of the study group (e.g., life-form) or assessed with relative rate tests and post-hoc model comparisons.

Dispersal capacity predicts both population genetic structure and species richness in reef fishes

Dispersal is a fundamental species characteristic that should directly affect both rates of gene flow among spatially distributed populations and opportunities for speciation. Yet no single trait associated with dispersal has been demonstrated to affect both micro- and macroevolutionary patterns of diversity across a diverse biological assemblage. Here, we examine patterns of genetic differentiation and species richness in reef fishes, an assemblage of over 7,000 species comprising approximately one-third of the extant bony fishes and over one-tenth of living vertebrates. In reef fishes, dispersal occurs primarily during a planktonic larval stage. There are two major reproductive and parental investment syndromes among reef fishes, and the differences between them have implications for dispersal: (1) benthic guarding fishes lay negatively buoyant eggs, typically guarded by the male parent, and from these eggs hatch large, strongly swimming larvae; in contrast, (2) pelagic spawning fishes release small floating eggs directly into the water column, which drift unprotected before small weakly swimming larvae hatch. Using phylogenetic comparative methods, we show that benthic guarders have significantly greater population structure than pelagic spawners and additionally that taxonomic families of benthic guarders are more species rich than families of pelagic spawners. Our findings provide a compelling case for the continuity between micro- and macroevolutionary processes of biological diversification and underscore the importance of dispersal-related traits in influencing the mode and tempo of evolution.

Structural phylogenomics retrodicts the origin of the genetic code and uncovers the evolutionary impact of protein flexibility

The genetic code shapes the genetic repository. Its origin has puzzled molecular scientists for over half a century and remains a long-standing mystery. Here we show that the origin of the genetic code is tightly coupled to the history of aminoacyl-tRNA synthetase enzymes and their interactions with tRNA. A timeline of evolutionary appearance of protein domain families derived from a structural census in hundreds of genomes reveals the early emergence of the 'operational' RNA code and the late implementation of the standard genetic code. The emergence of codon specificities and amino acid charging involved tight coevolution of aminoacyl-tRNA synthetases and tRNA structures as well as episodes of structural recruitment. Remarkably, amino acid and dipeptide compositions of single-domain proteins appearing before the standard code suggest archaic synthetases with structures homologous to catalytic domains of tyrosyl-tRNA and seryl-tRNA synthetases were capable of peptide bond formation and aminoacylation. Results reveal that genetics arose through coevolutionary interactions between polypeptides and nucleic acid cofactors as an exacting mechanism that favored flexibility and folding of the emergent proteins. These enhancements of phenotypic robustness were likely internalized into the emerging genetic system with the early rise of modern protein structure.

Molecular evolution and the latitudinal biodiversity gradient

Species density is higher in the tropics (low latitude) than in temperate regions (high latitude) resulting in a latitudinal biodiversity gradient (LBG). The LBG must be generated by differential rates of speciation and/or extinction and/or immigration among regions, but the role of each of these processes is still unclear. Recent studies examining differences in rates of molecular evolution have inferred a direct link between rate of molecular evolution and rate of speciation, and postulated these as important drivers of the LBG. Here we review the molecular genetic evidence and examine the factors that might be responsible for differences in rates of molecular evolution. Critical to this is the directionality of the relationship between speciation rates and rates of molecular evolution.

Molecular hyperdiversity and evolution in very large populations

The genomic density of sequence polymorphisms critically affects the sensitivity of inferences about ongoing sequence evolution, function and demographic history. Most animal and plant genomes have relatively low densities of polymorphisms, but some species are hyperdiverse with neutral nucleotide heterozygosity exceeding 5%. Eukaryotes with extremely large populations, mimicking bacterial and viral populations, present novel opportunities for studying molecular evolution in sexually reproducing taxa with complex development. In particular, hyperdiverse species can help answer controversial questions about the evolution of genome complexity, the limits of natural selection, modes of adaptation and subtleties of the mutation process. However, such systems have some inherent complications and here we identify topics in need of theoretical developments. Close relatives of the model organisms Caenorhabditis elegans and Drosophila melanogaster provide known examples of hyperdiverse eukaryotes, encouraging functional dissection of resulting molecular evolutionary patterns. We recommend how best to exploit hyperdiverse populations for analysis, for example, in quantifying the impact of noncrossover recombination in genomes and for determining the identity and micro-evolutionary selective pressures on noncoding regulatory elements.

Rates of molecular evolution and diversification in plants: chloroplast substitution rates correlate with species-richness in the Proteaceae

BACKGROUND

Many factors have been identified as correlates of the rate of molecular evolution, such as body size and generation length. Analysis of many molecular phylogenies has also revealed correlations between substitution rates and clade size, suggesting a link between rates of molecular evolution and the process of diversification. However, it is not known whether this relationship applies to all lineages and all sequences. Here, in order to investigate how widespread this phenomenon is, we investigate patterns of substitution in chloroplast genomes of the diverse angiosperm family Proteaceae. We used DNA sequences from six chloroplast genes (6278bp alignment with 62 taxa) to test for a correlation between diversification and the rate of substitutions.

RESULTS

Using phylogenetically-independent sister pairs, we show that species-rich lineages of Proteaceae tend to have significantly higher chloroplast substitution rates, for both synonymous and non-synonymous substitutions.

CONCLUSIONS

We show that the rate of molecular evolution in chloroplast genomes is correlated with net diversification rates in this large plant family. We discuss the possible causes of this relationship, including molecular evolution driving diversification, speciation increasing the rate of substitutions, or a third factor causing an indirect link between molecular and diversification rates. The link between the synonymous substitution rate and clade size is consistent with a role for the mutation rate of chloroplasts driving the speed of reproductive isolation. We find no significant differences in the ratio of non-synonymous to synonymous substitutions between lineages differing in net diversification rate, therefore we detect no signal of population size changes or alteration in selection pressures that might be causing this relationship.

A universal trend among proteomes indicates an oily last common ancestor

Despite progresses in ancestral protein sequence reconstruction, much needs to be unraveled about the nature of the putative last common ancestral proteome that served as the prototype of all extant lifeforms. Here, we present data that indicate a steady decline (oil escape) in proteome hydrophobicity over species evolvedness (node number) evident in 272 diverse proteomes, which indicates a highly hydrophobic (oily) last common ancestor (LCA). This trend, obtained from simple considerations (free from sequence reconstruction methods), was corroborated by regression studies within homologous and orthologous protein clusters as well as phylogenetic estimates of the ancestral oil content. While indicating an inherent irreversibility in molecular evolution, oil escape also serves as a rare and universal reaction-coordinate for evolution (reinforcing Darwin's principle of Common Descent), and may prove important in matters such as (i) explaining the emergence of intrinsically disordered proteins, (ii) developing composition- and speciation-based "global" molecular clocks, and (iii) improving the statistical methods for ancestral sequence reconstruction.

A total-evidence approach to dating with fossils, applied to the early radiation of the hymenoptera

Phylogenies are usually dated by calibrating interior nodes against the fossil record. This relies on indirect methods that, in the worst case, misrepresent the fossil information. Here, we contrast such node dating with an approach that includes fossils along with the extant taxa in a Bayesian total-evidence analysis. As a test case, we focus on the early radiation of the Hymenoptera, mostly documented by poorly preserved impression fossils that are difficult to place phylogenetically. Specifically, we compare node dating using nine calibration points derived from the fossil record with total-evidence dating based on 343 morphological characters scored for 45 fossil (4--20 complete) and 68 extant taxa. In both cases we use molecular data from seven markers (∼5 kb) for the extant taxa. Because it is difficult to model speciation, extinction, sampling, and fossil preservation realistically, we develop a simple uniform prior for clock trees with fossils, and we use relaxed clock models to accommodate rate variation across the tree. Despite considerable uncertainty in the placement of most fossils, we find that they contribute significantly to the estimation of divergence times in the total-evidence analysis. In particular, the posterior distributions on divergence times are less sensitive to prior assumptions and tend to be more precise than in node dating. The total-evidence analysis also shows that four of the seven Hymenoptera calibration points used in node dating are likely to be based on erroneous or doubtful assumptions about the fossil placement. With respect to the early radiation of Hymenoptera, our results suggest that the crown group dates back to the Carboniferous, ∼309 Ma (95% interval: 291--347 Ma), and diversified into major extant lineages much earlier than previously thought, well before the Triassic. [Bayesian inference; fossil dating; morphological evolution; relaxed clock; statistical phylogenetics.].

The origin of modern frogs (Neobatrachia) was accompanied by acceleration in mitochondrial and nuclear substitution rates

Background

Understanding the causes underlying heterogeneity of molecular evolutionary rates among lineages is a long-standing and central question in evolutionary biology. Although several earlier studies showed that modern frogs (Neobatrachia) experienced an acceleration of mitochondrial gene substitution rates compared to non-neobatrachian relatives, no further characterization of this phenomenon was attempted. To gain new insights on this topic, we sequenced the complete mitochondrial genomes and nine nuclear loci of one pelobatoid (Pelodytes punctatus) and five neobatrachians, Heleophryne regis (Heleophrynidae), Lechriodus melanopyga (Limnodynastidae), Calyptocephalella gayi (Calyptocephalellidae), Telmatobius bolivianus (Ceratophryidae), and Sooglossus thomasseti (Sooglossidae). These represent major clades not included in previous mitogenomic analyses, and most of them are remarkably species-poor compared to other neobatrachians.

Results

We reconstructed a fully resolved and robust phylogeny of extant frogs based on the new mitochondrial and nuclear sequence data, and dated major cladogenetic events. The reconstructed tree recovered Heleophryne as sister group to all other neobatrachians, the Australasian Lechriodus and the South American Calyptocephalella formed a clade that was the sister group to Nobleobatrachia, and the Seychellois Sooglossus was recovered as the sister group of Ranoides. We used relative-rate tests and direct comparison of branch lengths from mitochondrial and nuclear-based trees to demonstrate that both mitochondrial and nuclear evolutionary rates are significantly higher in all neobatrachians compared to their non-neobatrachian relatives, and that such rate acceleration started at the origin of Neobatrachia.

Conclusions

Through the analysis of the selection coefficient (ω) in different branches of the tree, we found compelling evidence of relaxation of purifying selection in neobatrachians, which could (at least in part) explain the observed higher mitochondrial and nuclear substitution rates in this clade. Our analyses allowed us to discard that changes in substitution rates could be correlated with increased mitochondrial genome rearrangement or diversification rates observed in different lineages of neobatrachians.

Diversification and the rate of molecular evolution: no evidence of a link in mammals

BACKGROUND

Recent research has indicated a positive association between rates of molecular evolution and diversification in a number of taxa. However debate continues concerning the universality and cause of this relationship. Here, we present the first systematic investigation of this relationship within the mammals. We use phylogenetically independent sister-pair comparisons to test for a relationship between substitution rates and clade size at a number of taxonomic levels. Total, non-synonymous and synonymous substitution rates were estimated from mitochondrial and nuclear DNA sequences.

RESULTS

We found no evidence for an association between clade size and substitution rates in mammals, for either the nuclear or the mitochondrial sequences. We found significant associations between body size and substitution rates, as previously reported.

CONCLUSIONS

Our results present a contrast to previous research, which has reported significant positive associations between substitution rates and diversification for birds, angiosperms and reptiles. There are three possible reasons for the differences between the observed results in mammals versus other clades. First, there may be no link between substitution rates and diversification in mammals. Second, this link may exist, but may be much weaker in mammals than in other clades. Third, the link between substitution rates and diversification may exist in mammals, but may be confounded by other variables.