Age at first reproduction explains rate variation in the strepsirrhine molecular clock

Genomic basis of deep-water adaptation in Arctic Charr (Salvelinus alpinus) morphs

The post-glacial colonization of Gander Lake in Newfoundland, Canada, by Arctic Charr (Salvelinus alpinus) provides the opportunity to study the genomic basis of adaptation to extreme deep-water environments. Colonization of deep-water (>50 m) habitats often requires extensive adaptation to cope with novel environmental challenges from high hydrostatic pressure, low temperature, and low light, but the genomic mechanisms underlying evolution in these environments are rarely known. Here, we compare genomic divergence between a deep-water morph adapted to depths of up to 288 m and a larger, piscivorous pelagic morph occupying shallower depths. Using both a SNP array and resequencing of whole nuclear and mitochondrial genomes, we find clear genetic divergence (F_ST  = 0.11-0.15) between deep and shallow water morphs, despite an absence of morph divergence across the mitochondrial genome. Outlier analyses identified many diverged genomic regions containing genes enriched for processes such as gene expression and DNA repair, cardiac function, and membrane transport. Detection of putative copy number variants (CNVs) uncovered 385 genes with CNVs distinct to piscivorous morphs, and 275 genes with CNVs distinct to deep-water morphs, enriched for processes associated with synapse assembly. Demographic analyses identified evidence for recent and local morph divergence, and ongoing reductions in diversity consistent with postglacial colonization. Together, these results show that Arctic Charr morph divergence has occurred through genome-wide differentiation and elevated divergence of genes underlying multiple cellular and physiological processes, providing insight into the genomic basis of adaptation in a deep-water habitat following postglacial recolonization.

An empirical estimate of the generation time of mouse lemurs

The generation time of organisms drives the rate of change in populations and across evolutionary times. In long-lived species, generation time should also account for overlapping generations, and the average age of parents has been proposed as a best approximation under these conditions. This study uses this definition to estimate the generation time of a widely studied small primate, Microcebus murinus, based on parentage data generated for a free-living population over a 6-year period in northwestern Madagascar. The average age of parents was calculated separately for mothers and fathers of three different offspring cohorts that differed in the degree of demographic uncertainty. In addition, adult survival rates were calculated for males and females based on long-term capture data from the same population to estimate the possible upper limits of generation time. Adult survival was low with only 44% of adult females and 38% of adult males being recaptured at the beginning of their second breeding season. The average age of mothers was 1.56-1.91 years, pointing toward a 2-year female generation time due to the high proportion of 1-year old mothers in all three cohorts. Female generation time estimates were fairly stable across the three offspring cohorts. In contrast, the average age of fathers differed by more than 1 year from the first to the third offspring cohort (1.71-2.83 years) pointing toward a 3-year generation time, but also suggesting a higher degree of demographic uncertainty in the early years of the study. For future modeling purposes, we, therefore, propose to use the average, 2.5 years, of male and female values as new estimate for the generation time of mouse lemurs.

Support for the evolutionary speed hypothesis from intraspecific population genetic data in the non-biting midge Chironomus riparius

The evolutionary speed hypothesis (ESH) proposes a causal mechanism for the latitudinal diversity gradient. The central idea of the ESH is that warmer temperatures lead to shorter generation times and increased mutation rates. On an absolute time scale, both should lead to an acceleration of selection and drift. Based on the ESH, we developed predictions regarding the distribution of intraspecific genetic diversity: populations of ectothermic species with more generations per year owing to warmer ambient temperatures should be more differentiated from each other, accumulate more mutations and show evidence for increased mutation rates compared with populations in colder regions. We used the multivoltine insect species Chironomus riparius to test these predictions with cytochrome oxidase I (COI) sequence data and found that populations from warmer regions are indeed significantly more differentiated and have significantly more derived haplotypes than populations from colder regions. We also found a significant correlation of the annual mean temperature with the population mutation parameter θ that serves as a proxy for the per generation mutation rate under certain assumptions. This pattern could be corroborated with two nuclear loci. Overall, our results support the ESH and indicate that the thermal regime experienced may be crucially driving the evolution of ectotherms and may thus ultimately govern their speciation rate.

Patterns of Substitution Rate Variation at Many Nuclear Loci in Two Species Trios in the Brassicaceae Partitioned with ANOVA

There are marked variations among loci and among lineages in rates of nucleotide substitution. The generation time hypothesis (GTH) is a neutral explanation for substitution rate heterogeneity that has genomewide application, predicting that species with shorter generation times accumulate DNA sequence substitutions faster than species with longer generation times do since faster genome replication provides more opportunities for mutations to occur and reach fixation by genetic drift. Relatively few studies have rigorously evaluated the GTH in plants, and there are numerous alternative hypotheses for plant substitution rate variation. One major challenge has been finding pairs of closely related plant species with contrasting generation times and appropriate outgroup taxa that all also have DNA sequence data for numerous loci. To test for causes of rate variation, we obtained sequence data for 256 genes for Arabidopsis thaliana, normally reproducing every year, and the biennial Arabidopsis lyrata with three closely related outgroup taxa (Brassica rapa, Capsella grandiflora, and Neslia paniculata) as well as the biennial Brassica oleracea and the annual B. rapa lineage with the outgroup N. paniculata. A sign test indicated that more loci than expected by chance have faster rates of substitution on the branch leading to the annual than to the perennial for one three-species trio but not another. Tajima's 1D and 2D tests, and a likelihood ratio test that incorporated saturation correction, rejected rate homogeneity for up to 26 genes (up to 14 genes when correcting for multiple tests), consistently showing faster rates for the annual lineage in the Arabidopsis species trio. ANOVA showed significant rate heterogeneity between the Arabidopsis and Brassica species trios (about 6 % of rate variation) and among loci (about 26-32 % of rate variation). The lineage-by-locus interaction which would be caused by locus- and lineage-specific natural selection explained about 13 % of substitution rate variation in one ANOVA model using substitution rates from genes partitioned into odd and even codons but was not a significant effect without partitioned genes. Annual/perennial lineage and species trio by annual/perennial lineage each explained about 1 % of substitution rate variation.

Inflation of Molecular Clock Rates and Dates: Molecular Phylogenetics, Biogeography, and Diversification of a Global Cicada Radiation from Australasia (Hemiptera: Cicadidae: Cicadettini)

Dated phylogenetic trees are important for studying mechanisms of diversification, and molecular clocks are important tools for studies of organisms lacking good fossil records. However, studies have begun to identify problems in molecular clock dates caused by uncertainty of the modeled molecular substitution process. Here we explore Bayesian relaxed-clock molecular dating while studying the biogeography of ca. 200 species from the global cicada tribe Cicadettini. Because the available fossils are few and uninformative, we calibrate our trees in part with a cytochrome oxidase I (COI) clock prior encompassing a range of literature estimates for arthropods. We show that tribe-level analyses calibrated solely with the COI clock recover extremely old dates that conflict with published estimates for two well-studied New Zealand subclades within Cicadettini. Additional subclade analyses suggest that COI relaxed-clock rates and maximum-likelihood branch lengths become inflated relative to EF-1[Formula: see text] intron and exon rates and branch lengths as clade age increases. We present corrected estimates derived from: (i) an extrapolated EF-1[Formula: see text] exon clock derived from COI-calibrated analysis within the largest New Zealand subclade; (ii) post hoc scaling of the tribe-level chronogram using results from subclade analyses; and (iii) exploitation of a geological calibration point associated with New Caledonia. We caution that considerable uncertainty is generated due to dependence of substitution estimates on both the taxon sample and the choice of model, including gamma category number and the choice of empirical versus estimated base frequencies. Our results suggest that diversification of the tribe Cicadettini commenced in the early- to mid-Cenozoic and continued with the development of open, arid habitats in Australia and worldwide. We find that Cicadettini is a rare example of a global terrestrial animal group with an Australasian origin, with all non-Australasian genera belonging to two distal clades. Within Australia, we show that Cicadettini is more widely distributed than any other cicada tribe, diverse in temperate, arid and monsoonal habitats, and nearly absent from rainforests. We comment on the taxonomic implications of our findings for thirteen cicada genera.

Generation time, life history and the substitution rate of neutral mutations

Our understanding of molecular evolution is hampered by a lack of quantitative predictions about how life-history (LH) traits should correlate with substitution rates. Comparative studies have shown that neutral substitution rates vary substantially between species, and evidence shows that much of this diversity is associated with variation in LH traits. However, while these studies often agree, some unexplained and contradictory results have emerged. Explaining these results is difficult without a clear theoretical understanding of the problem. In this study, we derive predictions for the relationships between LH traits and substitution rates in iteroparous species by using demographic theory to relate commonly measured life-history traits to genetic generation time, and by implication to neutral substitution rates. This provides some surprisingly simple explanations for otherwise confusing patterns, such as the association between fecundity and substitution rates. The same framework can be applied to more complex life histories if full life-tables are available.

Evolutionary Modeling of Genotype-Phenotype Associations, and Application to Primate Coding and Non-coding mtDNA Rate Variation

Variation in substitution rates across a phylogeny can be indicative of shifts in the evolutionary dynamics of a protein or non-protein coding regions. One way to understand these signals is to seek the phenotypic correlates of rate variation. Here, we extended a previously published likelihood method designed to detect evolutionary associations between genotypic evolutionary rate and phenotype over a phylogeny. In simulation with two discrete categories of phenotype, the method has a low false-positive rate and detects greater than 80% of true-positives with a tree length of three or greater and a three-fold or greater change in substitution rate given the phenotype. In addition, we successfully extend the test from two to four phenotype categories and evaluated its performance. We then applied the method to two major hypotheses for rate variation in the mitochondrial genome of primates-longevity and generation time as well as body mass which is correlated with many aspects of life history-using three categories of phenotype through discretization of continuous values. Similar to previous results for mammals, we find that the majority of mitochondrial protein-coding genes show associations consistent with the longevity and body mass predictions and that the predominant signal of association comes from the third codon position. We also found a significant association between maximum lifespan and the evolutionary rate of the control region of the mtDNA. In contrast, 24 protein-coding genes from the nuclear genome do not show a consistent pattern of association, which is inconsistent with the generation time hypothesis. These results show the extended method can robustly identify genotype-phenotype associations up to at least four phenotypic categories, and demonstrate the successful application of the method to study factors affecting neutral evolutionary rate in protein-coding and non-coding loci.

Exploring the correlations between sequence evolution rate and phenotypic divergence across the Mammalian tree provides insights into adaptive evolution

Sequence evolution behaves in a relatively consistent manner, leading to one of the fundamental paradigms in biology, the existence of a 'molecular clock'. The molecular clock can be distilled to the concept of accumulation of substitutions, through time yielding a stable rate from which we can estimate lineage divergence. Over the last 50 years, evolutionary biologists have obtained an in-depth understanding of this clock's nuances. It has been fine-tuned by taking into account the vast heterogeneity in rates across lineages and genes, leading to 'relaxed' molecular clock methods for timetree reconstruction. Sequence rate varies with life history traits including body size, generation time and metabolic rate, and we review recent studies on this topic. However, few studies have explicitly examined correlates between molecular evolution and morphological evolution. The patterns observed across diverse lineages suggest that rates of molecular and morphological evolution are largely decoupled. We discuss how identifying the molecular mechanisms behind rapid functional radiations are central to understanding evolution. The vast functional divergence within mammalian lineages that have relatively 'slow' sequence evolution refutes the hypotheses that pulses in diversification yielding major phenotypic change are the result of steady accumulation of substitutions. Patterns rather suggest phenotypic divergence is likely caused by regulatory alterations mediated through mechanisms such as insertions/deletions in functional regions. These can rapidly arise and sweep to fixation faster than predicted from a lineage's sequence neutral substitution rate, enabling species to leapfrog between phenotypic 'islands'. We suggest research directions that could illuminate mechanisms behind the functional diversity we see today.

Phylogenetic analyses uncover a novel clade of transferrin in nonmammalian vertebrates

Transferrin is a protein super-family involved in iron transport, a central process in cellular homeostasis. Throughout the evolution of vertebrates, transferrin members have diversified into distinct subfamilies including serotransferrin, ovotransferrin, lactoferrin, melanotransferrin, the inhibitor of carbonic anhydrase, pacifastin, and the major yolk protein in sea urchin. Previous phylogenetic analyses have established the branching order of the diverse transferrin subfamilies but were mostly focused on the transferrin repertoire present in mammals. Here, we conduct a comprehensive phylogenetic analysis of transferrin protein sequences in sequenced vertebrates, placing a special focus on the less-studied nonmammalian vertebrates. Our analyses uncover a novel transferrin clade present across fish, sauropsid, and amphibian genomes but strikingly absent from mammals. Our reconstructed scenario implies that this novel class emerged through a duplication event at the vertebrate ancestor, and that it was subsequently lost in the lineage leading to mammals. We detect footprints of accelerated evolution following the duplication event, which suggest positive selection and early functional divergence of this novel clade. Interestingly, the loss of this novel class of transferrin in mammals coincided with the divergence by duplication of lactoferrin and serotransferrin in this lineage. Altogether, our results provide novel insights on the evolution of iron-binding proteins in the various vertebrate groups.

The determinants of the molecular substitution process in turtles

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

The lemur revolution starts now: the genomic coming of age for a non-model organism

Morris Goodman was a revolutionary. Together with a mere handful of like-minded scientists, Morris established himself as a leader in the molecular phylogenetic revolution of the 1960s. The effects of this revolution are most evident in this journal, which he founded in 1992. Happily for lemur biologists, one of Morris Goodman's primary interests was in reconstructing the phylogeny of the primates, including the tooth-combed Lorisifomes of Africa and Asia, and the Lemuriformes of Madagascar (collectively referred to as the suborder Strepsirrhini). This paper traces the development of molecular phylogenetic and evolutionary genetic trends and methods over the 50-year expanse of Morris Goodman's career, particularly as they apply to our understanding of lemuriform phylogeny, biogeography, and biology. Notably, this perspective reveals that the lemuriform genome is sufficiently rich in phylogenetic signal such that the very earliest molecular phylogenetic studies - many of which were conducted by Goodman himself - have been validated by contemporary studies that have exploited advanced computational methods applied to phylogenomic scale data; studies that were beyond imagining in the earliest days of phylogeny reconstruction. Nonetheless, the frontier still beckons. New technologies for gathering and analyzing genomic data will allow investigators to build upon what can now be considered a nearly-known phylogeny of the Lemuriformes in order to ask innovative questions about the evolutionary mechanisms that generate and maintain the extraordinary breadth and depth of biological diversity within this remarkable clade of primates.

Morris Goodman's hominoid rate slowdown: the importance of being neutral

Half a century ago, when the field of molecular evolution did not even exist, Morris Goodman analyzed profiles of immunological interactions between species and reached the following two remarkable conclusions: first, protein evolution slowed down in the human lineage compared to other primate lineages; second, this slowdown was more pronounced for proteins whose functions were likely to be neutral. It took several decades of research to fully grasp these ideas and document the pattern of hominoid rate slowdown. Along the way, studies of hominoid rate slowdown led to major progresses in understanding determinants of neutral molecular evolution, which in turn is used to calibrate rates of adaptive evolution. Furthermore, the growing knowledge on the origin of mutations provides a basis for understanding differential evolutionary rates between sex chromosomes and autosomes, which has deep implications for inferring human evolutionary histories, and other aspects of molecular evolution. Primate genomics in particular stand to provide critical information in these pursuits, due to the abundance of genomic data, relatively rich documentation of life history traits, and several model systems, including our own species.

Evidence for a convergent slowdown in primate molecular rates and its implications for the timing of early primate evolution

A long-standing problem in primate evolution is the discord between paleontological and molecular clock estimates for the time of crown primate origins: the earliest crown primate fossils are ~56 million y (Ma) old, whereas molecular estimates for the haplorhine-strepsirrhine split are often deep in the Late Cretaceous. One explanation for this phenomenon is that crown primates existed in the Cretaceous but that their fossil remains have not yet been found. Here we provide strong evidence that this discordance is better-explained by a convergent molecular rate slowdown in early primate evolution. We show that molecular rates in primates are strongly and inversely related to three life-history correlates: body size (BS), absolute endocranial volume (EV), and relative endocranial volume (REV). Critically, these traits can be reconstructed from fossils, allowing molecular rates to be predicted for extinct primates. To this end, we modeled the evolutionary history of BS, EV, and REV using data from both extinct and extant primates. We show that the primate last common ancestor had a very small BS, EV, and REV. There has been a subsequent convergent increase in BS, EV, and REV, indicating that there has also been a convergent molecular rate slowdown over primate evolution. We generated a unique timescale for primates by predicting molecular rates from the reconstructed phenotypic values for a large phylogeny of living and extinct primates. This analysis suggests that crown primates originated close to the K-Pg boundary and possibly in the Paleocene, largely reconciling the molecular and fossil timescales of primate evolution.

Evolution of the couple cytochrome c and cytochrome c oxidase in primates

Mitochondrial energy metabolism has been affected by a broad set of ancient and recent evolutionary events. The oldest example is the endosymbiosis theory that led to mitochondria and a recently proposed example is adaptation to cold climate by anatomically modern human lineages. Mitochondrial energy metabolism has also been associated with an important area in anthropology and evolutionary biology, brain enlargement in human evolution. Indeed, several studies have pointed to the need for a major metabolic rearrangement to supply a sufficient amount of energy for brain development in primates.The genes encoding for the coupled cytochrome c (Cyt c) and cytochrome c oxidase (COX, complex IV, EC 1.9.3.1) seem to have an exceptional pattern of evolution in the anthropoid lineage. It has been proposed that this evolution was linked to the rearrangement of energy metabolism needed for brain enlargement. This hypothesis is reinforced by the fact that the COX enzyme was proposed to have a large role in control of the respiratory chain and thereby global energy production.After summarizing major events that occurred during the evolution of COX and cytochrome c on the primate lineage, we review the different evolutionary forces that could have influenced primate COX evolution and discuss the probable causes and consequences of this evolution. Finally, we discuss and review the co-occurring primate phenotypic evolution.

A genome sequence resource for the aye-aye (Daubentonia madagascariensis), a nocturnal lemur from Madagascar

We present a high-coverage draft genome assembly of the aye-aye (Daubentonia madagascariensis), a highly unusual nocturnal primate from Madagascar. Our assembly totals ~3.0 billion bp (3.0 Gb), roughly the size of the human genome, comprised of ~2.6 million scaffolds (N50 scaffold size = 13,597 bp) based on short paired-end sequencing reads. We compared the aye-aye genome sequence data with four other published primate genomes (human, chimpanzee, orangutan, and rhesus macaque) as well as with the mouse and dog genomes as nonprimate outgroups. Unexpectedly, we observed strong evidence for a relatively slow substitution rate in the aye-aye lineage compared with these and other primates. In fact, the aye-aye branch length is estimated to be ~10% shorter than that of the human lineage, which is known for its low substitution rate. This finding may be explained, in part, by the protracted aye-aye life-history pattern, including late weaning and age of first reproduction relative to other lemurs. Additionally, the availability of this draft lemur genome sequence allowed us to polarize nucleotide and protein sequence changes to the ancestral primate lineage-a critical period in primate evolution, for which the relevant fossil record is sparse. Finally, we identified 293,800 high-confidence single nucleotide polymorphisms in the donor individual for our aye-aye genome sequence, a captive-born individual from two wild-born parents. The resulting heterozygosity estimate of 0.051% is the lowest of any primate studied to date, which is understandable considering the aye-aye's extensive home-range size and relatively low population densities. Yet this level of genetic diversity also suggests that conservation efforts benefiting this unusual species should be prioritized, especially in the face of the accelerating degradation and fragmentation of Madagascar's forests.

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

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

The neutral theory of molecular evolution in the genomic era

The neutral theory of molecular evolution has been widely accepted and is the guiding principle for studying evolutionary genomics and the molecular basis of phenotypic evolution. Recent data on genomic evolution are generally consistent with the neutral theory. However, many recently published papers claim the detection of positive Darwinian selection via the use of new statistical methods. Examination of these methods has shown that their theoretical bases are not well established and often result in high rates of false-positive and false-negative results. When the deficiencies of these statistical methods are rectified, the results become largely consistent with the neutral theory. At present, genome-wide analyses of natural selection consist of collections of single-locus analyses. However, because phenotypic evolution is controlled by the interaction of many genes, the study of natural selection ought to take such interactions into account. Experimental studies of evolution will also be crucial.

Evolutionary rates in Veronica L. (Plantaginaceae): disentangling the influence of life history and breeding system

The evolutionary rate at which DNA sequences evolve is known to differ between different groups of organisms. However, the reasons for these different rates are seldom known. Among plants, the generation-time hypothesis, which states that organisms that reproduce faster also have more DNA substitutions per time, has gained most popularity. We evaluate the generation-time hypothesis using 131 DNA sequences from the plastid trnLF region and the nuclear ribosomal ITS region of the genus Veronica (Plantaginaceae). We also examine the alternative hypothesis that a higher substitution rate is correlated with selfing breeding system. Selfing is associated with annual life history in many organisms and may thus often be the underlying reason for observed correlations of annual life history with other characters. We provide evidence that annual life history is more likely to be the responsible factor for higher substitution rates in Veronica than a selfing breeding system. Nevertheless, the way in which annual life history may influence substitution rate in detail remains unknown, and some possibilities are discussed.