The determinants of the molecular substitution process in turtles

Taking advantage of reference-guided assembly in a slowly-evolving lineage: Application to Testudo graeca

BACKGROUND

Obtaining de novo chromosome-level genome assemblies greatly enhances conservation and evolutionary biology studies. For many research teams, long-read sequencing technologies (that produce highly contiguous assemblies) remain unaffordable or unpractical. For the groups that display high synteny conservation, these limitations can be overcome by a reference-guided assembly using a close relative genome. Among chelonians, tortoises (Testudinidae) are considered one of the most endangered taxa, which calls for more genomic resources. Here we make the most of high synteny conservation in chelonians to produce the first chromosome-level genome assembly of the genus Testudo with one of the most iconic tortoise species in the Mediterranean basin: Testudo graeca.

RESULTS

We used high-quality, paired-end Illumina sequences to build a reference-guided assembly with the chromosome-level reference of Gopherus evgoodei. We reconstructed a 2.29 Gb haploid genome with a scaffold N50 of 107.598 Mb and 5.37% gaps. We sequenced 25,998 protein-coding genes, and identified 41.2% of the assembly as repeats. Demographic history reconstruction based on the genome revealed two events (population decline and recovery) that were consistent with previously suggested phylogeographic patterns for the species. This outlines the value of such reference-guided assemblies for phylogeographic studies.

CONCLUSIONS

Our results highlight the value of using close relatives to produce de novo draft assemblies in species where such resources are unavailable. Our annotated genome of T. graeca paves the way to delve deeper into the species' evolutionary history and provides a valuable resource to enhance direct conservation efforts on their threatened populations.

Genomic estimates of mutation and substitution rates contradict the evolutionary speed hypothesis of the latitudinal diversity gradient

The latitudinal diversity gradient (LDG) refers to a decrease in biodiversity from the equator to the poles. The evolutionary speed hypothesis, backed by the metabolic theory of ecology, asserts that nucleotide mutation and substitution rates per site per year are higher and thereby speciation rates are higher at higher temperatures, generating the LDG. However, prior empirical investigations of the relationship between the temperature and mutation or substitution rate were based on a few genes and the results were mixed. We here revisit this relationship using genomic data. No significant correlation between the temperature and mutation rate is found in 13 prokaryotes or in 107 eukaryotes. An analysis of 234 diverse trios of bacterial taxa indicates that the synonymous substitution rate is not significantly associated with the growth temperature. The same data, however, reveal a significant negative association between the nonsynonymous substitution rate and temperature, which is explainable by a larger fraction of detrimental nonsynonymous mutations at higher temperatures due to a stronger demand for protein stability. We conclude that the evolutionary speed hypothesis of the LDG is unsupported by genomic data and advise that future mechanistic studies of the LDG should focus on other hypotheses.

Genomic analyses of the scorpion mud turtle (Kinosternon scorpioides) (Linnaeus, 1766) in insular and continental Colombia: Evidence for multiple conservation and taxonomic units

The turtle genus Kinosternon is widespread with at least 25 species distributed from Mexico to northern Argentina. The taxonomy of this genus is controversial and requires a full revision using both morphological and molecular approaches. In this study, we did a genomic analysis on the species Kinosternon scorpioides distributed in insular and continental Colombia in order to define conservation units. Total DNA was extracted from 24 tissue samples and RADseq genotyping analysis was done. In addition, the intron R35 was amplified and sequenced for a subset of samples. A total of 35,507 SNPs combined with 1,047 bp of the intron were used for spatiotemporal colonization pattern reconstruction and phylogenetic analyses. In addition, SNPs were used for population structure inferences and allele frequency-based analyses. Reciprocal monophyly, significant differences in allele frequencies (Fst = 0.32 - 0.78), and evidence of reproductive isolation (no admixture/geneflow), indicate long-term divergence between groups (2-8 MYA), possibly due to geographical barriers. Four Evolutionarily Significant Units (ESUs) were defined within our samples. One ESU was represented by the insular subspecies K. scorpioides albogulare, found on San Andrés island, and three ESUs were defined for the subspecies K. s. scorpioides in continental Colombia: one trans-Andean, found in northwestern Colombia (Caribbean region) and two cis-Andean, found in eastern and southeastern Colombia in the Orinoco and Amazon regions, respectively. Colonization of this species occurred from an ancestral area on South of Central America region (~ 8.43 MYA), followed by an establishing of current populations on San Andrés Island and then, in the continent. First, in the Colombian Caribbean, next, in the Orinoco, and more recently, in the Amazon. We hypothesize that the emergence of the Panamá Isthmus, as well as the final uplift of the North Eastern Andes and Vaupes Arch, were key event leading to the differentiation of these ESUs. For management and conservation purposes, each of these ESUs should be considered as a separate management unit. A full revision of the taxonomy of the genus Kinosternon is warranted.

Temperature-dependent evolutionary speed shapes the evolution of biodiversity patterns across tetrapod radiations

Biodiversity varies predictably with environmental energy around the globe, but the underlaying mechanisms remain incompletely understood. The evolutionary speed hypothesis predicts that environmental kinetic energy shapes variation in speciation rates through temperature- or life history-dependent rates of evolution. To test whether variation in evolutionary speed can explain the relationship between energy and biodiversity in birds, mammals, amphibians, and reptiles, we simulated diversification over 65 million years of geological and climatic change with a spatially explicit eco-evolutionary simulation model. We modelled four distinct evolutionary scenarios in which speciation-completion rates were dependent on temperature (M1), life history (M2), temperature and life history (M3), or were independent of temperature and life-history (M0). To assess the agreement between simulated and empirical data, we performed model selection by fitting supervised machine learning models to multidimensional biodiversity patterns. We show that a model with temperature-dependent rates of speciation (M1) consistently had the strongest support. In contrast to statistical inferences, which showed no general relationships between temperature and speciation rates in tetrapods, we demonstrate how process-based modelling can disentangle the causes behind empirical biodiversity patterns. Our study highlights how environmental energy has played a fundamental role in the evolution of biodiversity over deep time.

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.

Temperature predicts the rate of molecular evolution in Australian Eugongylinae skinks

Temperature differences over time and space have been hypothesized to cause variation in the rate of molecular evolution of species, but empirical evidence is mixed. To further test this hypothesis, we utilized a large exon-capture sequence data of Australian Eugongylinae skinks, exemplifying a radiation of temperature-sensitive ectotherms spanning a large latitudinal gradient. The association between temperature (and other species traits) and long-term substitution rate was assessed based on 1268 sequenced exons of 44 species pairs from the Eugongylinae subfamily using regression analyses. Temperature is the strongest, positively correlated predictor of variation in substitution rate across the Australian Eugongylinae. It explains 45% of variation in synonymous substitution rate, and 11% after controlling for all the other factors. Synonymous substitution rate is also negatively associated with body size, with a 6% variation explained by body size after controlling for the effects of temperature. Other factors are not associated with synonymous substitution rate after controlling for temperature. Overall, this study points to temperature as a strong predictor of the molecular evolution rate in the Eugongylinae subfamily, and demonstrates the power of large-scale exonic data to identify correlates of the rate of molecular evolution.

Why Do Some Vertebrates Have Microchromosomes?

With more than 70,000 living species, vertebrates have a huge impact on the field of biology and research, including karyotype evolution. One prominent aspect of many vertebrate karyotypes is the enigmatic occurrence of tiny and often cytogenetically indistinguishable microchromosomes, which possess distinctive features compared to macrochromosomes. Why certain vertebrate species carry these microchromosomes in some lineages while others do not, and how they evolve remain open questions. New studies have shown that microchromosomes exhibit certain unique characteristics of genome structure and organization, such as high gene densities, low heterochromatin levels, and high rates of recombination. Our review focuses on recent concepts to expand current knowledge on the dynamic nature of karyotype evolution in vertebrates, raising important questions regarding the evolutionary origins and ramifications of microchromosomes. We introduce the basic karyotypic features to clarify the size, shape, and morphology of macro- and microchromosomes and report their distribution across different lineages. Finally, we characterize the mechanisms of different evolutionary forces underlying the origin and evolution of microchromosomes.

Evolutionary Rates are Correlated Between Buchnera Endosymbionts and the Mitochondrial Genomes of Their Aphid Hosts

The evolution of bacterial endosymbiont genomes is strongly influenced by host-driven selection. Factors affecting host genome evolution will potentially affect endosymbiont genomes in similar ways. One potential outcome is correlations in molecular rates between the genomes of the symbiotic partners. Recently, we presented the first evidence of such correlations between the mitochondrial genomes of cockroaches and the genomes of their endosymbiont (Blattabacterium cuenoti). Here we investigate whether similar patterns are found in additional host-symbiont partners. We use partial genome data from multiple strains of the bacterial endosymbionts Buchnera aphidicola and Sulcia muelleri, and the mitochondrial genomes of their sap-feeding insect hosts. Both endosymbionts show phylogenetic congruence with the mitochondria of their hosts, a result that is expected due to their identical mode of inheritance. We compared root-to-tip distances and branch lengths of phylogenetically independent species pairs. Both analyses showed a highly significant correlation of molecular rates between the genomes of Buchnera and the mitochondrial genomes of their hosts. A similar correlation was detected between Sulcia and their hosts, but was not statistically significant. Our results indicate that evolutionary rate correlations between hosts and long-term symbionts may be a widespread phenomenon.

The Effects of Ecological Traits on the Rate of Molecular Evolution in Ray-Finned Fishes: A Multivariable Approach

Myriad environmental and biological traits have been investigated for their roles in influencing the rate of molecular evolution across various taxonomic groups. However, most studies have focused on a single trait, while controlling for additional factors in an informal way, generally by excluding taxa. This study utilized a dataset of cytochrome c oxidase subunit I (COI) barcode sequences from over 7000 ray-finned fish species to test the effects of 27 traits on molecular evolutionary rates. Environmental traits such as temperature were considered, as were traits associated with effective population size including body size and age at maturity. It was hypothesized that these traits would demonstrate significant correlations with substitution rate in a multivariable analysis due to their associations with mutation and fixation rates, respectively. A bioinformatics pipeline was developed to assemble and analyze sequence data retrieved from the Barcode of Life Data System (BOLD) and trait data obtained from FishBase. For use in phylogenetic regression analyses, a maximum likelihood tree was constructed from the COI sequence data using a multi-gene backbone constraint tree covering 71% of the species. A variable selection method that included both single- and multivariable analyses was used to identify traits that contribute to rate heterogeneity estimated from different codon positions. Our analyses revealed that molecular rates associated most significantly with latitude, body size, and habitat type. Overall, this study presents a novel and systematic approach for integrative data assembly and variable selection methodology in a phylogenetic framework.

Interstitial Telomeric Repeats Are Rare in Turtles

Telomeres are nucleoprotein complexes protecting chromosome ends in most eukaryotic organisms. In addition to chromosome ends, telomeric-like motifs can be accumulated in centromeric, pericentromeric and intermediate (i.e., between centromeres and telomeres) positions as so-called interstitial telomeric repeats (ITRs). We mapped the distribution of (TTAGGG)_n repeats in the karyotypes of 30 species from nine families of turtles using fluorescence in situ hybridization. All examined species showed the expected terminal topology of telomeric motifs at the edges of chromosomes. We detected ITRs in only five species from three families. Combining our and literature data, we inferred seven independent origins of ITRs among turtles. ITRs occurred in turtles in centromeric positions, often in several chromosomal pairs, in a given species. Their distribution does not correspond directly to interchromosomal rearrangements. Our findings support that centromeres and non-recombining parts of sex chromosomes are very dynamic genomic regions, even in turtles, a group generally thought to be slowly evolving. However, in contrast to squamate reptiles (lizards and snakes), where ITRs were found in more than half of the examined species, and birds, the presence of ITRs is generally rare in turtles, which agrees with the expected low rates of chromosomal rearrangements and rather slow karyotype evolution in this group.

Evolutionary rates are correlated between cockroach symbionts and mitochondrial genomes

Bacterial endosymbionts evolve under strong host-driven selection. Factors influencing host evolution might affect symbionts in similar ways, potentially leading to correlations between the molecular evolutionary rates of hosts and symbionts. Although there is evidence of rate correlations between mitochondrial and nuclear genes, similar investigations of hosts and symbionts are lacking. Here, we demonstrate a correlation in molecular rates between the genomes of an endosymbiont (Blattabacterium cuenoti) and the mitochondrial genomes of their hosts (cockroaches). We used partial genome data for multiple strains of B. cuenoti to compare phylogenetic relationships and evolutionary rates for 55 cockroach/symbiont pairs. The phylogenies inferred for B. cuenoti and the mitochondrial genomes of their hosts were largely congruent, as expected from their identical maternal and cytoplasmic mode of inheritance. We found a correlation between evolutionary rates of the two genomes, based on comparisons of root-to-tip distances and on comparisons of the branch lengths of phylogenetically independent species pairs. Our results underscore the profound effects that long-term symbiosis can have on the biology of each symbiotic partner.

Evolutionary Rates of Bumblebee Genomes Are Faster at Lower Elevations

The importance of climate in determining biodiversity patterns has been well documented. However, the relationship between climate and rates of genetic evolution remains controversial. Latitude and elevation have been associated with rates of change in genetic markers such as cytochrome b. What is not known, however, is the strength of such associations and whether patterns found among these genes apply across entire genomes. Here, using bumblebee genetic data from seven subgenera of Bombus, we demonstrate that all species occupying warmer elevations have undergone faster genome-wide evolution than those in the same subgenera occupying cooler elevations. Our findings point to a critical biogeographic role in the relative rates of whole species evolution, potentially influencing global biodiversity patterns.

Visualization of Genetic Drift Processes Using the Conserved Collagen 1α1 GXY Domain

Speciation proceeds by the accumulation of DNA differences in time. The genetic code changes as a result of genetic drift and by selective pressure. In variable domains, exposure to high selective pressure obscures the view on background mutations. Therefore, we characterized and visualized background mutations using the highly conserved collagen 1α1 GXY domain. Typical change routes were identified and the data set showed several indications that changes in the collagen 1α1 GXY domain have taken place randomly within a functionally restricted space. The types of nucleotide and codon group differences are similar across the vertebrate subphylum and gradually become less functionally neutral with increasing distance between species, which offers the opportunity for rapid visualization of evolutionary relations using a single domain. It was concluded that the findings and approach of the study could be important for analytical method development in authenticity research, especially when conserved domains are targeted.

Turtles of the genera Geoemyda and Pangshura (Testudines: Geoemydidae) lack differentiated sex chromosomes: the end of a 40-year error cascade for Pangshura

For a long time, turtles of the family Geoemydidae have been considered exceptional because representatives of this family were thought to possess a wide variety of sex determination systems. In the present study, we cytogenetically studied Geoemyda spengleri and G. japonica and re-examined the putative presence of sex chromosomes in Pangshura smithii. Karyotypes were examined by assessing the occurrence of constitutive heterochromatin, by comparative genome hybridization and in situ hybridization with repetitive motifs, which are often accumulated on differentiated sex chromosomes in reptiles. We found similar karyotypes, similar distributions of constitutive heterochromatin and a similar topology of tested repetitive motifs for all three species. We did not detect differentiated sex chromosomes in any of the species. For P. smithii, a ZZ/ZW sex determination system, with differentiated sex chromosomes, was described more than 40 years ago, but this finding has never been re-examined and was cited in all reviews of sex determination in reptiles. Here, we show that the identification of sex chromosomes in the original report was based on the erroneous pairing of chromosomes in the karyogram, causing over decades an error cascade regarding the inferences derived from the putative existence of female heterogamety in geoemydid turtles.

Global Rate Variation in Bony Vertebrates

This study investigated long-term substitution rate differences using three calibration points, divergences between lobe-finned vertebrates and ray-finned fish, between mammals and sauropsids, and between holosteans (gar and bowfin) and teleost fish with amino acid sequence data of 625 genes for 25 bony vertebrates. The result showed that the substitution rate was two to three times higher in the stem branches of lobe-finned vertebrates before the mammal-sauropsid divergence than in amniotes. The rate in the stem branch of ray-finned fish before the holostean-teleost fish divergence was also a few times higher than the holostean rate, whereas it was similar to or somewhat slower than the teleost fish rate. The phylogenetic relationship of coelacanth and lungfish with tetrapod was difficult to determine because of the short interval of the divergences. Considering the high rate in the stem branches, the divergences of coelacanth and lungfish from the stem branch were estimated as 408–427 Ma and 399–414 Ma, respectively, with the interval of 9–13 Myr. With the external calibration of the mammal-sauropsid split, the estimated times for ordinal divergences within eutherian mammals tend to be smaller than those in previous studies that used the calibration points within the lineage, with deeper divergences before the Cretaceous–Paleogene boundary and shallower ones after the boundary. In contrast the estimated times within birds were larger than those of previous studies, with the divergence between Galliformes and Anseriformes ∼80 Ma and that between Galloanserae and Neoaves 110 Ma.

Is molecular evolution faster in the tropics?

The evolutionary speed hypothesis (ESH) suggests that molecular evolutionary rates are higher among species inhabiting warmer environments. Previously, the ESH has been investigated using small numbers of latitudinally-separated sister lineages; in animals, these studies typically focused on subsets of Chordata and yielded mixed support for the ESH. This study analyzed public DNA barcode sequences from the cytochrome c oxidase subunit I (COI) gene for six of the largest animal phyla (Arthropoda, Chordata, Mollusca, Annelida, Echinodermata, and Cnidaria) and paired latitudinally-separated taxa together informatically. Of 8037 lineage pairs, just over half (51.6%) displayed a higher molecular rate in the lineage inhabiting latitudes closer to the equator, while the remainder (48.4%) displayed a higher rate in the higher-latitude lineage. To date, this study represents the most comprehensive analysis of latitude-related molecular rate differences across animals. While a statistically-significant pattern was detected from our large sample size, our findings suggest that the EHS may not serve as a strong universal mechanism underlying the latitudinal diversity gradient and that COI molecular clocks may generally be applied across latitudes. This study also highlights the merits of using automation to analyze large DNA barcode datasets.

Molecular phylogeny and divergence of the map turtles (Emydidae: Graptemys)

The map turtles (genus Graptemys) comprise a morphologically diverse clade that forms a major component of the southeastern US hotspot of chelonian diversity. Map turtles have experienced both recent and rapid diversification resulting in long-standing uncertainty regarding species boundaries and phylogenetic relationships within the genus as well as timing of their divergence. We present a phylogeny for the group that includes geographically representative sampling for all described species and subspecies. We make use of an empirical prior on rates of molecular evolution to estimate divergence times with a molecular clock under a coalescent framework. Together, the phylogeny and divergence time estimates suggest that diversification has been both more recent and more rapid than has so far been suspected. We provide a well-supported evolutionary framework for Graptemys that is necessary for understanding map turtle diversity, biogeography, and for conservation of this threatened clade of turtles.

Adaptive Patterns of Mitogenome Evolution Are Associated with the Loss of Shell Scutes in Turtles

The mitochondrial genome encodes several protein components of the oxidative phosphorylation (OXPHOS) pathway and is critical for aerobic respiration. These proteins have evolved adaptively in many taxa, but linking molecular-level patterns with higher-level attributes (e.g., morphology, physiology) remains a challenge. Turtles are a promising system for exploring mitochondrial genome evolution as different species face distinct respiratory challenges and employ multiple strategies for ensuring efficient respiration. One prominent adaptation to a highly aquatic lifestyle in turtles is the secondary loss of keratenized shell scutes (i.e., soft-shells), which is associated with enhanced swimming ability and, in some species, cutaneous respiration. We used codon models to examine patterns of selection on mitochondrial protein-coding genes along the three turtle lineages that independently evolved soft-shells. We found strong evidence for positive selection along the branches leading to the pig-nosed turtle (Carettochelys insculpta) and the softshells clade (Trionychidae), but only weak evidence for the leatherback (Dermochelys coriacea) branch. Positively selected sites were found to be particularly prevalent in OXPHOS Complex I proteins, especially subunit ND2, along both positively selected lineages, consistent with convergent adaptive evolution. Structural analysis showed that many of the identified sites are within key regions or near residues involved in proton transport, indicating that positive selection may have precipitated substantial changes in mitochondrial function. Overall, our study provides evidence that physiological challenges associated with adaptation to a highly aquatic lifestyle have shaped the evolution of the turtle mitochondrial genome in a lineage-specific manner.

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'.

Body mass-corrected molecular rate for bird mitochondrial DNA

Mitochondrial DNA remains one of the most widely used molecular markers to reconstruct the phylogeny and phylogeography of closely related birds. It has been proposed that bird mitochondrial genomes evolve at a constant rate of ~0.01 substitution per site per million years, that is that they evolve according to a strict molecular clock. This molecular clock is often used in studies of bird mitochondrial phylogeny and molecular dating. However, rates of mitochondrial genome evolution vary among bird species and correlate with life history traits such as body mass and generation time. These correlations could cause systematic biases in molecular dating studies that assume a strict molecular clock. In this study, we overcome this issue by estimating corrected molecular rates for birds. Using complete or nearly complete mitochondrial genomes of 475 species, we show that there are strong relationships between body mass and substitution rates across birds. We use this information to build models that use bird species' body mass to estimate their substitution rates across a wide range of common mitochondrial markers. We demonstrate the use of these corrected molecular rates on two recently published data sets. In one case, we obtained molecular dates that are twice as old as the estimates obtained using the strict molecular clock. We hope that this method to estimate molecular rates will increase the accuracy of future molecular dating studies in birds.

Transcriptome sequencing reveals genome-wide variation in molecular evolutionary rate among ferns

BACKGROUND

Transcriptomics in non-model plant systems has recently reached a point where the examination of nuclear genome-wide patterns in understudied groups is an achievable reality. This progress is especially notable in evolutionary studies of ferns, for which molecular resources to date have been derived primarily from the plastid genome. Here, we utilize transcriptome data in the first genome-wide comparative study of molecular evolutionary rate in ferns. We focus on the ecologically diverse family Pteridaceae, which comprises about 10 % of fern diversity and includes the enigmatic vittarioid ferns-an epiphytic, tropical lineage known for dramatically reduced morphologies and radically elongated phylogenetic branch lengths. Using expressed sequence data for 2091 loci, we perform pairwise comparisons of molecular evolutionary rate among 12 species spanning the three largest clades in the family and ask whether previously documented heterogeneity in plastid substitution rates is reflected in their nuclear genomes. We then inquire whether variation in evolutionary rate is being shaped by genes belonging to specific functional categories and test for differential patterns of selection.

RESULTS

We find significant, genome-wide differences in evolutionary rate for vittarioid ferns relative to all other lineages within the Pteridaceae, but we recover few significant correlations between faster/slower vittarioid loci and known functional gene categories. We demonstrate that the faster rates characteristic of the vittarioid ferns are likely not driven by positive selection, nor are they unique to any particular type of nucleotide substitution.

CONCLUSIONS

Our results reinforce recently reviewed mechanisms hypothesized to shape molecular evolutionary rates in vittarioid ferns and provide novel insight into substitution rate variation both within and among fern nuclear genomes.

Molecular evolutionary rates are not correlated with temperature and latitude in Squamata: an exception to the metabolic theory of ecology?

Background

The metabolic theory of ecology stipulates that molecular evolutionary rates should correlate with temperature and latitude in ectothermic organisms. Previous studies have shown that most groups of vertebrates, such as amphibians, turtles and even endothermic mammals, have higher molecular evolutionary rates in regions where temperature is high. However, the association between molecular evolutionary rates and temperature or latitude has never been tested in Squamata.

Results

We used a large dataset including the spatial distributions and environmental variables for 1,651 species of Squamata and compared the contrast of the rates of molecular evolution with the contrast of temperature and latitude between sister species. Using major axis regressions and a new algorithm to choose independent sister species pairs, we found that temperature and absolute latitude were not associated with molecular evolutionary rates.

Conclusions

This absence of association in such a diverse ectothermic group questions the mechanisms explaining current pattern of species diversity in Squamata and challenges the presupposed universality of the metabolic theory of ecology.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0666-4) contains supplementary material, which is available to authorized users.

Experimental Crossing of Two Distinct Species of Leopard Geckos, Eublepharis angramainyu and E. macularius: Viability, Fertility and Phenotypic Variation of the Hybrids

Hybridization between distinct species of animals and subsequent genetic introgression plays a considerable role in the speciation process and the emergence of adaptive characters. Fitness of between-species hybrids usually sharply decreases with the divergence time of the concerned species and the divergence depth, which still allows for a successful crossing differs among principal clades of vertebrates. Recently, a review of hybridization events among distinct lizard species revealed that lizards belong to vertebrates with a highly developed ability to hybridize. In spite of this, reliable reports of experimental hybridizations between genetically fairly divergent species are only exceptional. Here, we show the results of the crossing of two distinct allopatric species of eyelid geckos possessing temperature sex determination and lacking sex chromosomes: Eublepharis macularius distributed in Pakistan/Afghanistan area and E. angramainyu, which inhabits Mesopotamia and adjacent areas. We demonstrated that F₁ hybrids were viable and fertile, and the introgression of E. angramainyu genes into the E. macularius genome can be enabled via a backcrossing. The examined hybrids (except those of the F₂ generation) displayed neither malformations nor a reduced survival. Analyses of morphometric and coloration traits confirmed phenotypic distinctness of both parental species and their F₁ hybrids. These findings contrast with long-term geographic and an evolutionary separation of the studied species. Thus, the occurrence of fertile hybrids of comparably divergent species, such as E. angramainyu and E. macularius, may also be expected in other taxa of squamates. This would violate the current estimates of species diversity in lizards.

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.

Population genomics of the endangered giant Galápagos tortoise

BACKGROUND

The giant Galápagos tortoise, Chelonoidis nigra, is a large-sized terrestrial chelonian of high patrimonial interest. The species recently colonized a small continental archipelago, the Galápagos Islands, where it has been facing novel environmental conditions and limited resource availability. To explore the genomic consequences of this ecological shift, we analyze the transcriptomic variability of five individuals of C. nigra, and compare it to similar data obtained from several continental species of turtles.

RESULTS

Having clarified the timing of divergence in the Chelonoidis genus, we report in C. nigra a very low level of genetic polymorphism, signatures of a weakened efficacy of purifying selection, and an elevated mutation load in coding and regulatory sequences. These results are consistent with the hypothesis of an extremely low long-term effective population size in this insular species. Functional evolutionary analyses reveal a reduced diversity of immunity genes in C. nigra, in line with the hypothesis of attenuated pathogen diversity in islands, and an increased selective pressure on genes involved in response to stress, potentially related to the climatic instability of its environment and its elongated lifespan. Finally, we detect no population structure or homozygosity excess in our five-individual sample.

CONCLUSIONS

These results enlighten the molecular evolution of an endangered taxon in a stressful environment and point to island endemic species as a promising model for the study of the deleterious effects on genome evolution of a reduced long-term population size.

The human mutation rate is increasing, even as it slows

Substitution rates vary between species, and many explanations regarding the causes of this variation have been proposed. Here we consider how new genomic data on the per-generation mutation rate impinge on proposed hypotheses for substitution rate variation in primates. We propose that the generation-time effect as it is usually understood cannot explain the observed rate variation, but instead that selection for decreased somatic mutation rates can. By considering the disparate causes underlying mutation rate changes in recent human history, we also show that the per-generation mutation rate is increasing even as the per-cell-division rate is decreasing.