Modelling heterotachy in phylogenetic inference by reversible-jump Markov chain Monte Carlo

The rate at which a given site in a gene sequence alignment evolves over time may vary. This phenomenon--known as heterotachy--can bias or distort phylogenetic trees inferred from models of sequence evolution that assume rates of evolution are constant. Here, we describe a phylogenetic mixture model designed to accommodate heterotachy. The method sums the likelihood of the data at each site over more than one set of branch lengths on the same tree topology. A branch-length set that is best for one site may differ from the branch-length set that is best for some other site, thereby allowing different sites to have different rates of change throughout the tree. Because rate variation may not be present in all branches, we use a reversible-jump Markov chain Monte Carlo algorithm to identify those branches in which reliable amounts of heterotachy occur. We implement the method in combination with our 'pattern-heterogeneity' mixture model, applying it to simulated data and five published datasets. We find that complex evolutionary signals of heterotachy are routinely present over and above variation in the rate or pattern of evolution across sites, that the reversible-jump method requires far fewer parameters than conventional mixture models to describe it, and serves to identify the regions of the tree in which heterotachy is most pronounced. The reversible-jump procedure also removes the need for a posteriori tests of 'significance' such as the Akaike or Bayesian information criterion tests, or Bayes factors. Heterotachy has important consequences for the correct reconstruction of phylogenies as well as for tests of hypotheses that rely on accurate branch-length information. These include molecular clocks, analyses of tempo and mode of evolution, comparative studies and ancestral state reconstruction. The model is available from the authors' website, and can be used for the analysis of both nucleotide and morphological data.     

Artifactual phylogenies caused by correlated distribution of substitution rates among sites and lineages: the good, the bad, and the ugly

Despite the advances in understanding molecular evolution, current phylogenetic methods barely take account of a fraction of the complexity of evolution. We are chiefly constrained by our incomplete knowledge of molecular evolutionary processes and the limits of computational power. These limitations lead to the establishment of either biologically simplistic models that rarely account for a fraction of the complexity involved or overfitting models that add little resolution to the problem. Such oversimplified models may lead us to assign high confidence to an incorrect tree (inconsistency). Rate-across-site (RAS) models are commonly used evolutionary models in phylogenetic studies. These account for heterogeneity in the evolutionary rates among sites but do not account for changing within-site rates across lineages (heterotachy). If heterotachy is common, using RAS models may lead to systematic errors in tree inference. In this work we show possible misleading effects in tree inference when the assumption of constant within-site rates across lineages is violated using maximum likelihood. Using a simulation study, we explore the ways in which gamma stationary models can lead to wrong topology or to deceptive bootstrap support values when the within-site rates change across lineages. More precisely, we show that different degrees of heterotachy mislead phylogenetic inference when the model assumed is stationary. Finally, we propose a geometry-based approach to visualize and to test for the possible existence of bias due to heterotachy.     

Phylogeny and biogeography of the American live oaks (Quercus subsection Virentes): a genomic and population genetics approach

The nature and timing of evolution of niche differentiation among closely related species remains an important question in ecology and evolution. The American live oak clade, Virentes, which spans the unglaciated temperate and tropical regions of North America and Mesoamerica, provides an instructive system in which to examine speciation and niche evolution. We generated a fossil‐calibrated phylogeny of Virentes using RADseq data to estimate divergence times and used nuclear microsatellites, chloroplast sequences and an intron region of nitrate reductase (NIA‐i3) to examine genetic diversity within species, rates of gene flow among species and ancestral population size of disjunct sister species. Transitions in functional and morphological traits associated with ecological and climatic niche axes were examined across the phylogeny. We found the Virentes to be monophyletic with three subclades, including a southwest clade, a southeastern US clade and a Central American/Cuban clade. Despite high leaf morphological variation within species and transpecific chloroplast haplotypes, RADseq and nuclear SSR data showed genetic coherence of species. We estimated a crown date for Virentes of 11 Ma and implicated the formation of the Sea of Cortés in a speciation event ~5 Ma. Tree height at maturity, associated with fire tolerance, differs among the sympatric species, while freezing tolerance appears to have diverged repeatedly across the tropical–temperate divide. Sympatric species thus show evidence of ecological niche differentiation but share climatic niches, while allopatric and parapatric species conserve ecological niches, but diverge in climatic niches. The mode of speciation and/or degree of co‐occurrence may thus influence which niche axis plants diverge along.     

Molecular evolution and phylogeny of the atpB-rbcL spacer of chloroplast DNA in the true mosses.

The nucleotide variation of a noncoding region between the atpB and rbcL genes of the chloroplast genome was used to estimate the phylogeny of 11 species of true mosses (subclass Bryidae). The A+T rich (82.6%) spacer sequence is conserved with 48% of bases showing no variation between the ingroup and outgroup. Rooted at liverworts, Marchantia and Bazzania, the monophyly of true mosses was supported cladistically and statistically. A nonparametric Wilcoxon Signed-Ranks test Ts statistic for testing the taxonomic congruence showed no significant differences between gene trees and organism trees as well as between parsimony trees and neighbor-joining trees. The reconstructed phylogeny based on the atpB-rbcL spacer sequences indicated the validity of the division of acrocarpous and pleurocarpous mosses. The size of the chloroplast spacer in mosses fits into an evolutionary trend of increasing spacer length from liverworts through ferns to seed plants. According to the relative rate tests, the hypothesis of a molecular clock was supported in all species except for Thuidium, which evolved relatively fast. The evolutionary rate of the chloroplast DNA spacer in mosses was estimated to be (1.12 +/- 0.019) x 10(-10) nucleotides per site per year, which is close to the nonsynonymous substitution rates of the rbcL gene in the vascular plants. The constrained molecular evolution (total nucleotide substitutions, K approximately 0.0248) of the chloroplast DNA spacer is consistent with the slow evolution in morphological traits of mosses. Based on the calibrated evolutionary rate, the time of the divergence of true mosses was estimated to have been as early as 220 million years ago.     

Molecular phylogeny ofIva (Asteraceae,Heliantheae) based on chloroplast DNA restriction site variation

Iva s.str. (comprising ten species) was examined by cpDNA restriction site variation to determine phyletic relationships within the group. The results were compared with relationships proposed from other data. A total of 86 restriction site mutations was detected, 47 of which proved phylogenetically informative. A single most parsimonious tree was obtained using both Wagner and Dollo parsimony. The tree revealed three main lineages that are congruent with the three chromosome lineages (base numbers of x = 16, 17, 18). The monophyly of the x = 16 and 18 groups was supported strongly by molecular data, while the monophyly of x = 17 lineage was only supported moderately. Relationships among the three lineages indicate that the sect.Iva is paraphyletic because sect.Linearbractea is nested within it. Both morphological data and the secondary chemical data are in agreement with the proposed cpDNA phylogeny. Because of this agreement, sect.Iva is revised such that,I. axillaris was excluded and positioned within the newly proposed sect.Rhizoma. Patterns and rates of cpDNA evolution were also examined. The results indicated an uneven evolution in the chloroplast genome with different rates of cpDNA evolution in at least a few species ofIva. However, the evolutionary clock hypothesis can not be rejected within most of the lineages inIva.     

Late Pleistocene speciation of three closely related tree peonies endemic to the Qinling–Daba Mountains,a major glacial refugium in Central China

Determining the factors promoting speciation is a major task in ecological and evolutionary research and can be aided by phylogeographic analysis. The Qinling–Daba Mountains (QDM) located in central China form an important geographic barrier between southern subtropical and northern temperate regions, and exhibit complex topography, climatic, and ecological diversity. Surprisingly, few phylogeographic analyses and studies of plant speciation in this region have been conducted. To address this issue, we investigated the genetic divergence and evolutionary histories of three closely related tree peony species (Paeonia qiui, P. jishanensis, and P. rockii) endemic to the QDM. Forty populations of the three tree peony species were genotyped using 22 nuclear simple sequence repeat markers (nSSRs) and three chloroplast DNA sequences to assess genetic structure and phylogenetic relationships, supplemented by morphological characterization and ecological niche modeling (ENM). Morphological and molecular genetic analyses showed the three species to be clearly differentiated from each other. In addition, coalescent analyses using DIYABC conducted on nSSR variation indicated that the species diverged from each other in the late Pleistocene, while ecological niche modeling (ENM) suggested they occupied a larger area during the Last Glacial Maximum (LGM) than at present. The combined genetic evidence from nuclear and chloroplast DNA and the results of ENM indicate that each species persisted through the late Pleistocene in multiple refugia in the Qinling, Daba, and Taihang Mountains with divergence favored by restricted gene flow caused by geographic isolation, ecological divergence, and limited pollen and seed dispersal. Our study contributes to a growing understanding of the origin and population structure of tree peonies and provides insights into the high level of plant endemism present in the Qinling–Daba Mountains of Central China.     

Effects of genotype on rates of substitution during experimental evolution

Rates of molecular evolution may vary widely between populations, yet the causes of this variation are still incompletely understood. Genetic differences between populations may make an important contribution to variation in rates of evolution, owing to differences in fitness, population size, mutation rates, or in the distribution of fitness effects (DFEs) of available beneficial mutations. By whole genome sequencing of Escherichia coli populations experimentally evolved in the presence of a quinolone antibiotic, we found that rates of substitution varied by genotype, with evidence for a contribution from a genotype's starting fitness. Subsequent targeted sequencing showed that genotypes with high average substitution rates were more likely to undergo the simultaneous fixation of several mutations, consistent with theoretical models of multiple mutation dynamics. Moreover, patterns of substitution were indicative of epistatic relationships between known resistance mutations.     

Slow molecular evolution in 18S rDNA, rbcL and nad5 genes of mosses compared with higher plants

The evolutionary potential of bryophytes (mosses, liverworts and hornworts) has been debated for decades. Fossil record and biogeographical distribution patterns suggest very slow morphological evolution and the retainment of several ancient traits since the split with vascular plants some 450 million years ago. Many have argued that bryophytes may evolve as rapidly as higher plants on the molecular level, but this hypothesis has not been tested so far. Here, it is shown that mosses have experienced significantly lower rates of molecular evolution than higher plants within 18S rDNA (nuclear), rbcL (chloroplast) and nad5 (mitochondrial) genes. Mosses are on an average evolving 2-3 times slower than ferns, gymnosperms and angiosperms; and also green algae seem to be evolving faster than nonvascular plants. These results support the observation of a general correlation between morphological and molecular evolutionary rates in plants and also show that mosses are 'evolutionary sphinxes' regarding both morphological and molecular evolutionary potential.     

Macroecology and macroevolution of body size in Anolis lizards

Body size is one of the most influential traits affecting many ecological and physiological processes across animal and plant taxa. Studies of the environmental factors shaping body size patterns may evaluate either temporal or spatial dimensions. Here, we analyzed body size evolution in the radiation of Anolis lizards across both geographical and temporal dimensions. We used a set of macroecological and macroevolutionary methods to test current and past environmental effects on geographical gradients of body size and its evolutionary rates. First, we test whether a set of current ecological/physiological hypotheses (heat balance, productivity and seasonality) explains spatial body size gradients. Second, we evaluate how tempo (i.e. evolutionary rates) and mode (i.e. evolutionary process) of body size evolution changed through time and the role of paleo-temperatures on rates of body size evolution during the Cenozoic. We did not find a signature of current environmental variables driving spatial body size gradients. By contrast, we found strong support for a correlation between temperature changes (i.e. climate cooling) during the Cenozoic and rates of body size evolution (i.e. body size diversification). We suggest that patterns of body size evolution in Anolis lizards might be influenced by thermoregulatory behavior across clades and regions.     

Evolutionary and biogeographical implications of variation in skull morphology of raccoon dogs (Nyctereutes procyonoides,Mammalia: Carnivora)

Morphological characteristics reflect geographical variation resulting from adaptation to varying environmental conditions. Carnivore species distributed over a wide geographical range generally have highly polymorphic morphological variation. The raccoon dog (Nyctereutes procyonoides) has a longitudinal distribution restricted to East Asia and the northern Indochina Peninsula. Its unique geographical range makes it an appropriate model to examine how morphological differences are influenced by geography. To demonstrate morphological evolution of Russian, Chinese, Korean and Japanese raccoon dogs predicted by geographical differences, we tested the island rule and Bergmann's rule. We compared craniodental variation among populations and examined morphological implications for intraspecific taxonomic status. Insular raccoon dogs possessed substantially smaller body size than those from the mainland. Moreover, different island effects among Japanese islands were demonstrated by markedly larger occipital condyle breath in the Hokkaido population. Larger skull size in Russian and Hokkaido raccoon dogs could be explained by Bergmann's rule. Based on previous chromosomal and molecular studies and results of our morphological analyses, we suggest Japanese raccoon dogs are a distinct species from the mainland N. procyonoides.     

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

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

Complete chloroplast genome sequence of a tree fern Alsophila spinulosa : insights into evolutionary changes in fern chloroplast genomes

Background  

Ferns have generally been neglected in studies of chloroplast genomics. Before this study, only one polypod and two basal ferns had their complete chloroplast (cp) genome reported. Tree ferns represent an ancient fern lineage that first occurred in the Late Triassic. In recent phylogenetic analyses, tree ferns were shown to be the sister group of polypods, the most diverse group of living ferns. Availability of cp genome sequence from a tree fern will facilitate interpretation of the evolutionary changes of fern cp genomes. Here we have sequenced the complete cp genome of a scaly tree fern Alsophila spinulosa (Cyatheaceae).     

Phylogenetic Relationships and Morphological Character Evolution of Photosynthetic Euglenids (Excavata) Inferred from Taxon‐rich Analyses of Five Genes

Photosynthetic euglenids acquired chloroplasts by secondary endosymbiosis, which resulted in changes to their mode of nutrition and affected the evolution of their morphological characters. Mapping morphological characters onto a reliable molecular tree could elucidate major trends of those changes. We analyzed nucleotide sequence data from regions of three nuclear‐encoded genes (nSSU, nLSU, hsp90), one chloroplast‐encoded gene (cpSSU) and one nuclear‐encoded chloroplast gene (psbO) to estimate phylogenetic relationships among 59 photosynthetic euglenid species. Our results were consistent with previous works; most genera were monophyletic, except for the polyphyletic genus Euglena, and the paraphyletic genus Phacus. We also analyzed character evolution in photosynthetic euglenids using our phylogenetic tree and eight morphological traits commonly used for generic and species diagnoses, including: characters corresponding to well‐defined clades, apomorphies like presence of lorica and mucilaginous stalks, and homoplastic characters like rigid cells and presence of large paramylon grains. This research indicated that pyrenoids were lost twice during the evolution of phototrophic euglenids, and that mucocysts, which only occur in the genus Euglena, evolved independently at least twice. In contrast, the evolution of cell shape and chloroplast morphology was difficult to elucidate, and could not be unambiguously reconstructed in our analyses.     

Evidence for Maternal Transmission of Chloroplast DNA in the Genus Asplenium (Aspleniaceae,Pteridophyta)

Little is known about the mode of transmission of chloroplasts in ferns, despite the importance of such knowledge for molecular phylogenetic and biosystematic studies. Andersson-Kötto (1930, 1931) inferred from crossing experiments that chloroplasts are inherited biparentally in Asplenium scolopendrium L. Here we present evidence from artificial hybridisation that demonstrates maternal inheritance of chloroplast DNA in the genus Asplenium (Aspleniaceae, Pteridophyta) using length variation in a non-coding spacer (trnL_UAA-trnF_GAA) in the chloroplast DNA.     

Chloroplast DNA variation and population structure in the widespread forest tree, Eucalyptus grandis

Recognition of genetic structure of populations and the ability to identify vulnerable populations is useful for the formation of conservation management strategies for plants. Eucalyptus grandis is a tall forest tree that has a major area of occurrence in subtropical eastern Australia, with smaller populations located in the east coast tropics. Many widespread forest species exhibit population differentiation that corresponds to geographic regions. However, Eucalyptus grandis appears to be an exception based on isozyme and morphological data. This is intriguing given a large discontinuity between northern populations and those in the southern part of the species range. In this study, the distribution of a maternally inherited chloroplast locus was examined because it was more likely to reveal genetic structure due to the slower evolution of the chloroplast genome and limited dispersal of seed in eucalypts. As expected, the G _ST for chloroplast DNA was higher than that for nuclear DNA but indicated low population differentiation for a forest tree species. Phylogeographic analysis indicated that the 15 populations grouped into three broad geographical regions; however, overall population structure was weak suggesting that the large geographical disjunction in the distribution of E. grandis may be relatively recent. A paradigm for conservation management of E. grandis based on chloroplast DNA haplotype distribution would take into account the low differentiation among populations.     

Effective population size and the rate and pattern of nucleotide substitutions

Both the overall rate of nucleotide substitution and the relative proportions of synonymous and non-synonymous substitutions are predicted to vary between species that differ in effective population size (N_e). Our understanding of the genetic processes underlying these lineage-specific differences in molecular evolution is still developing. Empirical analyses indicate that variation in substitution rates and patterns caused by differences in N_e is often substantial, however, and must be accounted for in analyses of molecular evolution.     

Scombroid Fishes Provide Novel Insights into the Trait/Rate Associations of Molecular Evolution

The study of which life history traits primarily affect molecular evolutionary rates is often confounded by the covariance of these traits. Scombroid fishes (billfishes, tunas, barracudas, and their relatives) are unusual in that their mass-specific metabolic rate is positively associated with body size. This study exploits this atypical pattern of trait variation, which allows for direct tests of whether mass-specific metabolic rate or body size is the more important factor of molecular evolutionary rates. We inferred a phylogeny for scombroids from a supermatrix of molecular and morphological characters and used new phylogenetic comparative approaches to assess the associations of body size and mass-specific metabolic rate with substitution rate. As predicted by the body size hypothesis, there is a negative correlation between body size and substitution rate. However, unexpectedly, we also find a negative association between mass-specific metabolic and substitution rates. These relationships are supported by analyses of the total molecular data, separate mitochondrial and nuclear genes, and individual loci, and they are robust to phylogenetic uncertainty. The molecular evolutionary rates of scombroids are primarily tied to body size. This study demonstrates that groups with novel patterns of trait variation can be particularly informative for identifying which life history traits are the primary factors of molecular evolutionary rates.     

Characterizing morphological (co)variation using structural equation models: Body size,allometric relationships and evolvability in a house sparrow metapopulation

Body size plays a key role in the ecology and evolution of all organisms. Therefore, quantifying the sources of morphological (co)variation, dependent and independent of body size, is of key importance when trying to understand and predict responses to selection. We combine structural equation modeling with quantitative genetics analyses to study morphological (co)variation in a meta‐population of house sparrows (Passer domesticus). As expected, we found evidence of a latent variable “body size,” causing genetic and environmental covariation between morphological traits. Estimates of conditional evolvability show that allometric relationships constrain the independent evolution of house sparrow morphology. We also found spatial differences in general body size and its allometric relationships. On islands where birds are more dispersive and mobile, individuals were smaller and had proportionally longer wings for their body size. Although on islands where sparrows are more sedentary and nest in dense colonies, individuals were larger and had proportionally longer tarsi for their body size. We corroborated these results using simulations and show that our analyses produce unbiased allometric slope estimates. This study highlights that in the short term allometric relationships may constrain phenotypic evolution, but that in the long term selection pressures can also shape allometric relationships.     

Flight loss linked to faster molecular evolution in insects

The loss of flight ability has occurred thousands of times independently during insect evolution. Flight loss may be linked to higher molecular evolutionary rates because of reductions in effective population sizes (N_e) and relaxed selective constraints. Reduced dispersal ability increases population subdivision, may decrease geographical range size and increases (sub)population extinction risk, thus leading to an expected reduction in N_e. Additionally, flight loss in birds has been linked to higher molecular rates of energy-related genes, probably owing to relaxed selective constraints on energy metabolism. We tested for an association between insect flight loss and molecular rates through comparative analysis in 49 phylogenetically independent transitions spanning multiple taxa, including moths, flies, beetles, mayflies, stick insects, stoneflies, scorpionflies and caddisflies, using available nuclear and mitochondrial protein-coding DNA sequences. We estimated the rate of molecular evolution of flightless (FL) and related flight-capable lineages by ratios of non-synonymous-to-synonymous substitutions (dN/dS) and overall substitution rates (OSRs). Across multiple instances of flight loss, we show a significant pattern of higher dN/dS ratios and OSRs in FL lineages in mitochondrial but not nuclear genes. These patterns may be explained by relaxed selective constraints in FL ectotherms relating to energy metabolism, possibly in combination with reduced N_e.     

Variable extent of parallelism in respiratory,circulatory, and neurological traits across lake whitefish species pairs

Parallel adaptive radiation events provide a powerful framework for investigations of ecology's contribution to phenotypic diversification. Ecologically driven divergence has been invoked to explain the repeated evolution of sympatric dwarf and normal lake whitefish (Coregonus clupeaformis) species in multiple lakes in eastern North America. Nevertheless, links between most putatively adaptive traits and ecological variation remain poorly defined within and among whitefish species pairs. Here, we examine four species pairs for variation in gill, heart, and brain size; three traits predicted to show strong phenotypic responses to ecological divergence. In each of the species pairs, normals exhibited larger body size standardized gills compared to dwarfs – a pattern that is suggestive of a common ecological driver of gill size divergence. Within lakes, the seasonal hypoxia experienced in the benthic environment is a likely factor leading to the requirement for larger gills in normals. Interestingly, the morphological pathways used to achieve larger gills varied between species pairs from Québec and Maine, which may imply subtle non‐parallelism in gill size divergence related to differences in genetic background. There was also a non‐significant trend toward larger hearts in dwarfs, the more active species of the two, whereas brain size varied exclusively among the lake populations. Taken together, our results suggest that the diversification of whitefish has been driven by parallel and non‐parallel ecological conditions across lakes. Furthermore, the phenotypic response to ecological variation may depend on genetic background of each population.