Speciation Times Research Articles

From gene trees to a dated allopolyploid network: insights from the angiosperm genus Viola (Violaceae).

Allopolyploidization accounts for a significant fraction of speciation events in many eukaryotic lineages. However, existing phylogenetic and dating methods require tree-like topologies and are unable to handle the network-like phylogenetic relationships of lineages containing allopolyploids. No explicit framework has so far been established for evaluating competing network topologies, and few attempts have been made to date phylogenetic networks. We used a four-step approach to generate a dated polyploid species network for the cosmopolitan angiosperm genus Viola L. (Violaceae Batch.). The genus contains ca 600 species and both recent (neo-) and more ancient (meso-) polyploid lineages distributed over 16 sections. First, we obtained DNA sequences of three low-copy nuclear genes and one chloroplast region, from 42 species representing all 16 sections. Second, we obtained fossil-calibrated chronograms for each nuclear gene marker. Third, we determined the most parsimonious multilabeled genome tree and its corresponding network, resolved at the section (not the species) level. Reconstructing the “correct” network for a set of polyploids depends on recovering all homoeologs, i.e., all subgenomes, in these polyploids. Assuming the presence of Viola subgenome lineages that were not detected by the nuclear gene phylogenies (“ghost subgenome lineages”) significantly reduced the number of inferred polyploidization events. We identified the most parsimonious network topology from a set of five competing scenarios differing in the interpretation of homoeolog extinctions and lineage sorting, based on (i) fewest possible ghost subgenome lineages, (ii) fewest possible polyploidization events, and (iii) least possible deviation from expected ploidy as inferred from available chromosome counts of the involved polyploid taxa. Finally, we estimated the homoploid and polyploid speciation times of the most parsimonious network. Homoploid speciation times were estimated by coalescent analysis of gene tree node ages. Polyploid speciation times were estimated by comparing branch lengths and speciation rates of lineages with and without ploidy shifts. Our analyses recognize Viola as an old genus (crown age 31 Ma) whose evolutionary history has been profoundly affected by allopolyploidy. Between 16 and 21 allopolyploidizations are necessary to explain the diversification of the 16 major lineages (sections) of Viola, suggesting that allopolyploidy has accounted for a high percentage—between 67% and 88%—of the speciation events at this level. The theoretical and methodological approaches presented here for (i) constructing networks and (ii) dating speciation events within a network, have general applicability for phylogenetic studies of groups where allopolyploidization has occurred. They make explicit use of a hitherto underexplored source of ploidy information from chromosome counts to help resolve phylogenetic cases where incomplete sequence data hampers network inference. Importantly, the coalescent-based method used herein circumvents the assumption of tree-like evolution required by most techniques for dating speciation events.

PyRate: a new program to estimate speciation and extinction rates from incomplete fossil data

Summary Despite the advancement of phylogenetic methods to estimate speciation and extinction rates, their power can be limited under variable rates, in particular for clades with high extinction rates and small number of extant species. Fossil data can provide a powerful alternative source of information to investigate diversification processes. Here, we present PyRate, a computer program to estimate speciation and extinction rates and their temporal dynamics from fossil occurrence data. The rates are inferred in a Bayesian framework and are comparable to those estimated from phylogenetic trees. We describe how PyRate can be used to explore different models of diversification. In addition to the diversification rates, it provides estimates of the parameters of the preservation process (fossilization and sampling) and the times of speciation and extinction of each species in the data set. Moreover, we develop a new birth–death model to correlate the variation of speciation/extinction rates with changes of a continuous trait. Finally, we demonstrate the use of Bayes factors for model selection and show how the posterior estimates of a PyRate analysis can be used to generate calibration densities for Bayesian molecular clock analysis. PyRate is an open‐source command‐line Python program available at http://sourceforge.net/projects/pyrate/.

Integrative taxonomy and evolutionary history of a newly revealed spiderDysdera ninniicomplex (Araneae: Dysderidae)

The spider genus Dysdera is a species-rich clade of specialized woodlice predators, composed typically of complexes of sibling species. Here, we analyse the Dysdera ninnii complex, distinguishing three species that exhibit slight but constant differences in the morphology of their copulatory organs, and in their genetic background. We designate a neotype for D. ninnii and redescribe it. We consider Dysdera pavesii Thorell, 1873 to be a junior synonym of Dysdera ninnii Canestrini, 1868. In addition, we describe two new species (Dysdera moravica sp. nov. and Dysdera microdonta sp. nov.). All three species occur in the region of north-eastern Italy, Slovenia, and north-western Croatia. Dysdera moravica sp. nov. expanded to central Europe. The species occur allopatrically or parapatrically. All three species possess the same diploid number and X0 sex chromosome determination. In some individuals we found chromosome fusions, and such polymorphism is common in spiders with holokinetic chromosomes. The analysis of mitochondrial (cytochrome c oxidase subunit I, COI) and nuclear ribosomal (internal transcribed spacer 2, ITS2) DNA markers revealed two clades, one formed by D. ninnii and D. microdonta sp. nov., and a second by D. moravica sp. nov. Species of the first clade are not well defined by DNA markers. We noticed only weak separation of maternally inherited COI, and even overlap of autosomally inherited ITS2 sequences. We suggest that either short speciation time, unfinished lineage sorting, or rare hybridization events caused this pattern. In one sample of D. microdonta sp. nov. we detected the coxA gene of a Rickettsia species, which is the first record of this parasitic bacteria from the spider family Dysderidae. Dysdera microdonta sp. nov. occurs at higher altitudes than D. ninnii, and their distribution ranges form a long contact zone. Remarkably, we did not record any overlap of the two distribution ranges, suggesting that the lack of a precopulatory interspecific barrier causes a loss of reproduction potential. We hypothesize that because of the unsolved interspecific barrier together with only tiny differences in morphology and COI sequences, and no differences in karyotypes and ITS2 sequences, the D. ninnii species complex is evolutionarily young. © 2014 The Linnean Society of London

Unsorted homology within locus and species trees.

The concept of homology lies at the root of evolutionary biology. Since the seminal work of Fitch (1970), three main categories of homology relationships have been defined at the molecular level: orthology, paralogy, and xenology. In brief, if two gene copies arose by duplication they are paralogs, whereas if they arose through speciation they are orthologs. If one of them was transferred from a contemporaneous species, we call them xenologs (Supplementary Fig. S1 in Supplementary Material online, available at http://dx.doi.org/10.5061/dryad.87k57; see Gray and Fitch (1983); Fitch (2000)). Indeed, these terms were coined under a phylogenetic framework in which species were represented by single individuals, and as such they have remained very much intact during the last four decades—although particular cases within these categories have received specific names (Mindell and Meyer 2001). However, advances in sequencing technology have changed the field, and it is now very common to collect data sets containing multiple gene loci and/or multiple individuals per species. In general, such genome-wide data sets not only have unveiled extensive phylogenomic incongruence (Jeffroy et al. 2006; Salichos and Rokas 2013) but have brought back to the spotlight the consideration of how ancestral polymorphisms sort within populations (Edwards 2009). Altogether, phylogenomic data make imperative the explicit distinction between organismal and gene histories. Let us consider phylogenetic relationships at three different levels: species, loci, and gene copies (Fig. 1). The distinction between species/population trees and gene trees has been known for decades (Goodman et al. 1979; Pamilo and Nei 1988; Takahata 1989), whereas the introduction of locus trees into these models is very recent (Rasmussen and Kellis 2012). In brief, a species tree depicts the evolutionary history of the sampled organisms. In this case, the nodes represent speciation events, connected by branches that reflect the population history along these periods, and where their widths represent effective population size (Ne) and their lengths represent time (usually in years or number of generations). Apart from speciations, only evolutionary processes that affect species as a whole are represented at this level, like hybridization. Note that species trees are equivalent to population trees when the organismal units of interest are conspecific populations. In this case, the nodes of the population trees represent isolation events. In general, we will refer to “species” as any diverging, interbreeding group of individuals regardless of its taxonomic rank. On the other hand, a locus tree represents the evolutionary history of the sampled loci for a given gene family (see Rasmussen and Kellis 2012). Since the loci exist inside individuals evolving as part of a population, the locus tree is embedded within the species tree. In a locus tree, the nodes depict either genetic divergence due to speciation in the embedding species tree or locus-level events such as duplication, losses, or horizontal gene transfers, whereas the branch lengths and widths represent time and Ne, respectively. Here, we assume that the locuslevel events get immediately fixed in the population, so these Ne are equivalent to those in the species tree and are the same for every locus. Finally, a gene tree represents the evolutionary history of the sampled gene copies that evolve inside the locus tree. Gene tree nodes indicate coalescent events, which looking forward in time correspond to the process of DNA replication and divergence, and that can occur around the speciation time, well before (deep coalescence) or afterwards (migration in population trees). The branches of the gene tree usually represent amount of substitutions per site, and can also represent number of generations or other measures of time. Importantly, these three historical layers do not necessarily coincide. True species/population trees can differ from true locus trees due to gene duplications, losses, and/or horizontal gene transfers, whereas true gene trees can differ from their embedding locus and species trees if there is incomplete lineage sorting (ILS) (Maddison 1997; Page and Charleston 1997) (and migration in the case of population trees). In this regard,

The fossilized birth-death process for coherent calibration of divergence-time estimates.

Time-calibrated species phylogenies are critical for addressing a wide range of questions in evolutionary biology, such as those that elucidate historical biogeography or uncover patterns of coevolution and diversification. Because molecular sequence data are not informative on absolute time, external data--most commonly, fossil age estimates--are required to calibrate estimates of species divergence dates. For Bayesian divergence time methods, the common practice for calibration using fossil information involves placing arbitrarily chosen parametric distributions on internal nodes, often disregarding most of the information in the fossil record. We introduce the "fossilized birth-death" (FBD) process--a model for calibrating divergence time estimates in a Bayesian framework, explicitly acknowledging that extant species and fossils are part of the same macroevolutionary process. Under this model, absolute node age estimates are calibrated by a single diversification model and arbitrary calibration densities are not necessary. Moreover, the FBD model allows for inclusion of all available fossils. We performed analyses of simulated data and show that node age estimation under the FBD model results in robust and accurate estimates of species divergence times with realistic measures of statistical uncertainty, overcoming major limitations of standard divergence time estimation methods. We used this model to estimate the speciation times for a dataset composed of all living bears, indicating that the genus Ursus diversified in the Late Miocene to Middle Pliocene.

Laminar flame speeds, counterflow ignition, and kinetic modeling of the butene isomers

Laminar flame speeds and ignition temperatures of nonpremixed counterflow were measured experimentally for the butene isomers at normal and elevated pressures. Results show that the flame speed increases in the order of isobutene, trans-2-butene, and 1-butene. Furthermore, isobutene has the highest ignition temperature, while those of trans-2-butene and cis-2-butene are quite similar to each other and are slightly higher than that of 1-butene. These results indicate that the reactivities of the butene isomers increase in the order of isobutene, 2-butene, and 1-butene. Furthermore, the critical reaction paths for butene isomer oxidation during ignition were analyzed by sensitivity analysis and a new kinetic model was developed with updated rate coefficients from ab initio calculation and kinetic theories. Simulations with the present kinetic model agree well with the measured laminar flame speeds, counterflow ignition temperatures as well as recent laminar flame speciation and ignition delay times for all butene isomers.

Quantifying the effects of anagenetic and cladogenetic evolution

An ongoing debate in evolutionary biology is whether phenotypic change occurs predominantly around the time of speciation or whether it instead accumulates gradually over time. In this work I propose a general framework incorporating both types of change, quantify the effects of speciational change via the correlation between species and attribute the proportion of change to each type. I discuss results of parameter estimation of Hominoid body size in this light. I derive mathematical formulae related to this problem, the probability generating functions of the number of speciation events along a randomly drawn lineage and from the most recent common ancestor of two randomly chosen tip species for a conditioned Yule tree. Additionally I obtain in closed form the variance of the distance from the root to the most recent common ancestor of two randomly chosen tip species.

Climatic niche breadth and species richness in temperate treefrogs

AbstractAimPatterns of species richness are often closely linked with climate, but the specific mechanisms by which species' climatic niches underlie large‐scale richness patterns remain poorly understood. It has been hypothesized that reduced temperature seasonality in the tropics promotes the evolution of species with narrow temperature niche breadths, and that this hypothesis helps explain high tropical richness. However, the relationship between species' climatic niche breadths and species richness has yet to be tested. We have addressed this issue using treefrogs (Hylidae) in eastern North America.LocationEastern North America.MethodsWe characterized climatic niches and niche breadths for all 24 hylid species in eastern North America using temperature and precipitation variables. We then examined the relationships between species richness, climatic niche positions and climatic niche breadths using phylogenetic comparative methods.ResultsSpecies richness was negatively associated with mean climatic niche breadth, such that high‐richness climates had species with narrower climatic niches. Our results also supported the roles of niche conservatism and the time‐for‐speciation effect in generating the relationship between climate and species richness in the region (more species in warm, wet regions that have been inhabited longer). Importantly, we show that the invasion of low‐richness climates has occurred primarily through recent intraspecific niche expansion into these climates rather than evolution of species that are narrowly specialized for these conditions (although the two hylid clades studied showed somewhat different patterns).Main conclusionsWe found that climatic zones with high species richness contain more species with narrower climatic niche breadths. Our results suggest that this pattern arose because narrow climatic niche breadths restricted the dispersal of most hylid species out of the ancestral, warm, moist climatic zones, allowing more time for speciation to build up higher species richness in these zones.

Speciation, population structure, and demographic history of the Mojave Fringe-toed Lizard (Uma scoparia), a species of conservation concern.

The North American deserts were impacted by both Neogene plate tectonics and Quaternary climatic fluctuations, yet it remains unclear how these events influenced speciation in this region. We tested published hypotheses regarding the timing and mode of speciation, population structure, and demographic history of the Mojave Fringe-toed Lizard (Uma scoparia), a sand dune specialist endemic to the Mojave Desert of California and Arizona. We sampled 109 individual lizards representing 22 insular dune localities, obtained DNA sequences for 14 nuclear loci, and found that U. scoparia has low genetic diversity relative to the U. notata species complex, comparable to that of chimpanzees and southern elephant seals. Analyses of genotypes using Bayesian clustering algorithms did not identify discrete populations within U. scoparia. Using isolation-with-migration (IM) models and a novel coalescent-based hypothesis testing approach, we estimated that U. scoparia diverged from U. notata in the Pleistocene epoch. The likelihood ratio test and the Akaike Information Criterion consistently rejected nested speciation models that included parameters for migration and population growth of U. scoparia. We reject the Neogene vicariance hypothesis for the speciation of U. scoparia and define this species as a single evolutionarily significant unit for conservation purposes.

Vicariant versus dispersal processes in the settlement of <italic>Pseudamnicola</italic> (Caenogastropoda, Hydrobiidae) in the Mediterranean Balearic Islands

The genus Pseudamnicola comprises a group of tiny springsnails inhabiting several continental and insular regions of the Mediterranean basin. Given the limited dispersal capabilities of these animals, it is difficult to explain the wide distribution range of the genus and, more specifically, its presence in isolated habitats, such as on islands. Thus, to investigate the process(es) that may explain these distribution patterns, we morphologically re-described and genetically analysed the six Pseudamnicola (Pseudamnicola) species occurring in the Iberian Peninsula and the nearby Balearic Islands. Genetic relationships were explored by sequencing two mitochondrial (cytochrome c oxidase subunit I and 16S rRNA) and one nuclear (28S rRNA) gene in 19 populations. Our morphological study confirmed the presence of previously described species, whereas our phylogenetic results revealed three lineages within the subgenus: one clade grouping the species from Minorca Island with an Iberian Peninsula species, a second clade grouping the three species from Majorca Island, and a third clade that consists of a single species, which occurs in both the Iberian Peninsula and Ibiza Island. Calculated speciation times show that the cladogenetic events involving the insular species seem to have occurred after the current conformation of the Balearic Islands (c. 20 Mya). Therefore, the speciation process may have been related to subsequent transmarine colonizations, probably during the Messinian Salinity Crisis, and the Pleistocene glaciations when landmass corridors connected the islands with the continent. © 2014 The Linnean Society of London

Vicariant versus dispersal processes in the settlement ofPseudamnicola(Caenogastropoda, Hydrobiidae) in the Mediterranean Balearic Islands

The genus Pseudamnicola comprises a group of tiny springsnails inhabiting several continental and insular regions of the Mediterranean basin. Given the limited dispersal capabilities of these animals, it is difficult to explain the wide distribution range of the genus and, more specifically, its presence in isolated habitats, such as on islands. Thus, to investigate the process(es) that may explain these distribution patterns, we morphologically re-described and genetically analysed the six Pseudamnicola (Pseudamnicola) species occurring in the Iberian Peninsula and the nearby Balearic Islands. Genetic relationships were explored by sequencing two mitochondrial (cytochrome c oxidase subunit I and 16S rRNA) and one nuclear (28S rRNA) gene in 19 populations. Our morphological study confirmed the presence of previously described species, whereas our phylogenetic results revealed three lineages within the subgenus: one clade grouping the species from Minorca Island with an Iberian Peninsula species, a second clade grouping the three species from Majorca Island, and a third clade that consists of a single species, which occurs in both the Iberian Peninsula and Ibiza Island. Calculated speciation times show that the cladogenetic events involving the insular species seem to have occurred after the current conformation of the Balearic Islands (c. 20 Mya). Therefore, the speciation process may have been related to subsequent transmarine colonizations, probably during the Messinian Salinity Crisis, and the Pleistocene glaciations when landmass corridors connected the islands with the continent. © 2014 The Linnean Society of London

NICHE AND RANGE SIZE PATTERNS SUGGEST THAT SPECIATION BEGINS IN SMALL, ECOLOGICALLY DIVERGED POPULATIONS IN NORTH AMERICAN MONKEYFLOWERS (MIMULUS SPP.)

Closely related species (e.g., sister taxa) often occupy very different ecological niches and can exhibit large differences in geographic distributions despite their shared evolutionary history. Budding speciation is one process that may partially explain how differences in niche and distribution characteristics may rapidly evolve. Budding speciation is the process through which new species form as initially small colonizing populations that acquire reproductive isolation. This mode of species formation predicts that, at the time of speciation, sister species should have highly asymmetrical distributions. We tested this hypothesis in North American monkeyflowers, a diverse clade with a robust phylogeny, using data on geographical ranges, climate, and plant community attributes. We found that recently diverged sister pairs have highly asymmetrical ranges and niche breadths, relative to older sister pairs. Additionally, we found that sister species occupy distinct environmental niche positions, and that 80% of sister species have completely or partially overlapping distributions (i.e., are broadly sympatric). Together, these results suggest that budding speciation has occurred frequently in Mimulus, that it has likely taken place both inside the range and on the range periphery, and that observed divergences in habitat and resource use could be associated with speciation in small populations.

Bayesian Estimation of Speciation and Extinction from Incomplete Fossil Occurrence Data

The temporal dynamics of species diversity are shaped by variations in the rates of speciation and extinction, and there is a long history of inferring these rates using first and last appearances of taxa in the fossil record. Understanding diversity dynamics critically depends on unbiased estimates of the unobserved times of speciation and extinction for all lineages, but the inference of these parameters is challenging due to the complex nature of the available data. Here, we present a new probabilistic framework to jointly estimate species-specific times of speciation and extinction and the rates of the underlying birth-death process based on the fossil record. The rates are allowed to vary through time independently of each other, and the probability of preservation and sampling is explicitly incorporated in the model to estimate the true lifespan of each lineage. We implement a Bayesian algorithm to assess the presence of rate shifts by exploring alternative diversification models. Tests on a range of simulated data sets reveal the accuracy and robustness of our approach against violations of the underlying assumptions and various degrees of data incompleteness. Finally, we demonstrate the application of our method with the diversification of the mammal family Rhinocerotidae and reveal a complex history of repeated and independent temporal shifts of both speciation and extinction rates, leading to the expansion and subsequent decline of the group. The estimated parameters of the birth-death process implemented here are directly comparable with those obtained from dated molecular phylogenies. Thus, our model represents a step towards integrating phylogenetic and fossil information to infer macroevolutionary processes.

Erosive processes after tectonic uplift stimulate vicariant and adaptive speciation: evolution in an Afrotemperate-endemic paper daisy genus

BackgroundThe role of tectonic uplift in stimulating speciation in South Africa’s only alpine zone, the Drakensberg, has not been explicitly examined. Tectonic processes may influence speciation both through the creation of novel habitats and by physically isolating plant populations. We use the Afrotemperate endemic daisy genus Macowania to explore the timing and mode (geographic versus adaptive) of speciation in this region. Between sister species pairs we expect high morphological divergence where speciation has happened in sympatry (adaptive) while with geographic (vicariant) speciation we may expect to find less morphological divergence and a greater degree of allopatry. A dated molecular phylogenetic hypothesis for Macowania elucidates species’ relationships and is used to address the potential impact of uplift on diversification. Morphological divergence of a small sample of reproductive and vegetative characters, used as a proxy for adaptive divergence, is measured against species’ range distributions to estimate mode of speciation across two subclades in the genus.ResultsThe Macowania crown age is consistent with the hypothesis of post-uplift diversification, and we find evidence for both vicariant and adaptive speciation between the two subclades within Macowania. Both subclades exhibit strong signals of range allopatry, suggesting that geographic isolation was important in speciation. One subclade, associated with dry, rocky environments at high altitudes, shows very little morphological and ecological differentiation but high range allopatry. The other subclade occupies a greater variety of habitats and exhibits far greater morphological differentiation, but contains species with overlapping distribution ranges.ConclusionsSpecies in Macowania are likely to have diversified in response to tectonic uplift, and we invoke uplift and uplift-mediated erosion as the main drivers of speciation. The greater relative morphological divergence in sympatric species of Macowania indicates that speciation in the non-sympatric taxa may not have required obvious adaptive differences, implying that simple geographic isolation was the driving force for speciation (‘neutral speciation’).

Time evolution of Pb(II) speciation in Pampa soil fractions

Soil pollution by heavy metals, particularly lead, is an important environmental concern; the bioavailability of such pollutants is strongly dependent on their chemical form. Here, the speciation of Pb(II) in soil fractions as a function of time shortly after its incorporation is studied, using a selective sequential extraction method. The sample came from an Argentinean Pampas region and was extensively characterised, including Rietveld analysis of the silt+clay fraction XRD pattern to find the major mineral components. Experiments were run twice, once in the winter and once in the summer. The results show different speciation time profiles in both cases, showing faster changes in winter due to the higher water content. The summer experiment corresponds to an earlier stage in the speciation profile evolution compared with winter. The soluble/exchangeable fraction decreases with time in summer but shows a lower and constant value in the winter. A high proportion is found to be adsorbed onto the stable (aluminosilicates+quartz) mineral fraction. The results strongly suggest that, even at a short time following soil pollution with Pb, a high proportion is adsorbed onto the mineral fraction, with only a low fraction being bioavailable. The most stable (mineral incorporated) form is observed to increase with time. Soil water content appears to be more important than temperature in determining the differences between the two.

Robust estimates of divergence times and selection with a poisson random field model: a case study of comparative phylogeographic data.

Mutation frequencies can be modeled as a Poisson random field (PRF) to estimate speciation times and the degree of selection on newly arisen mutations. This approach provides a quantitative theory for comparing intraspecific polymorphism with interspecific divergence in the presence of selection and can be used to estimate population genetic parameters. Although the original PRF model has been extended to more general biological settings to make statistical inference about selection and divergence among model organisms, it has not been incorporated into phylogeographic studies that focus on estimating population genetic parameters for nonmodel organisms. Here, we modified a recently developed time-dependent PRF model to independently estimate genetic parameters from a nuclear and mitochondrial DNA data set of 22 sister pairs of birds that have diverged across a biogeographic barrier. We found that species that inhabit humid habitats had more recent divergence times and larger effective population sizes than those that inhabit drier habitats, and divergence time estimated from the PRF model were similar to estimates from a coalescent species-tree approach. Selection coefficients were higher in sister pairs that inhabited drier habitats than in those in humid habitats, but overall the mitochondrial DNA was under weak selection. Our study indicates that PRF models are useful for estimating various population genetic parameters and serve as a framework for incorporating estimates of selection into comparative phylogeographic studies.

Birth–death models and coalescent point processes: The shape and probability of reconstructed phylogenies

Forward-in-time models of diversification (i.e., speciation and extinction) produce phylogenetic trees that grow “vertically” as time goes by. Pruning the extinct lineages out of such trees leads to natural models for reconstructed trees (i.e., phylogenies of extant species). Alternatively, reconstructed trees can be modelled by coalescent point processes (CPPs), where trees grow “horizontally” by the sequential addition of vertical edges. Each new edge starts at some random speciation time and ends at the present time; speciation times are drawn from the same distribution independently. CPPs lead to extremely fast computation of tree likelihoods and simulation of reconstructed trees. Their topology always follows the uniform distribution on ranked tree shapes (URT).We characterize which forward-in-time models lead to URT reconstructed trees and among these, which lead to CPP reconstructed trees. We show that for any “asymmetric” diversification model in which speciation rates only depend on time and extinction rates only depend on time and on a non-heritable trait (e.g., age), the reconstructed tree is CPP, even if extant species are incompletely sampled. If rates additionally depend on the number of species, the reconstructed tree is (only) URT (but not CPP). We characterize the common distribution of speciation times in the CPP description, and discuss incomplete species sampling as well as three special model cases in detail: (1) the extinction rate does not depend on a trait; (2) rates do not depend on time; (3) mass extinctions may happen additionally at certain points in the past.

Robustness to Divergence Time Underestimation When Inferring Species Trees from Estimated Gene Trees

To infer species trees from gene trees estimated from phylogenomic data sets, tractable methods are needed that can handle dozens to hundreds of loci. We examine several computationally efficient approaches-MP-EST, STAR, STEAC, STELLS, and STEM-for inferring species trees from gene trees estimated using maximum likelihood (ML) and Bayesian approaches. Among the methods examined, we found that topology-based methods often performed better using ML gene trees and methods employing coalescent times typically performed better using Bayesian gene trees, with MP-EST, STAR, STEAC, and STELLS outperforming STEM under most conditions. We examine why the STEM tree (also called GLASS or Maximum Tree) is less accurate on estimated gene trees by comparing estimated and true coalescence times, performing species tree inference using simulations, and analyzing a great ape data set keeping track of false positive and false negative rates for inferred clades. We find that although true coalescence times are more ancient than speciation times under the multispecies coalescent model, estimated coalescence times are often more recent than speciation times. This underestimation can lead to increased bias and lack of resolution with increased sampling (either alleles or loci) when gene trees are estimated with ML. The problem appears to be less severe using Bayesian gene-tree estimates.

Evolution of the Sex Chromosomes in Salmonid Fishes

Most of the information on sex chromosomes in salmonid fishes is for species in the 3 genera of the subfamily Salmoninae found in North America: Salvelinus, Salmo and Oncorhynchus. All of the species are male heterogametic with XY sex determination. Morphologically distinguishable sex chromosomes are found only in Salvelinus namaycush,S. fontinalis and Oncorhynchus mykiss. Sex chromosomes have been identified in the other species using a combination of chromosome mapping and fluorescence in situ hybridization with probes containing sex-linked markers. Although all species share conserved linkage groups, the major sex-determining locus (SD) is found at the telomere of a different linkage group in almost every species, suggesting that the SD often transposes to a new location at the time of speciation. In a couple of species, intraspecific variation has been found in the chromosomal location of the SD. Recently, sdY has been identified as the major sex-determining gene in rainbow trout, and it maps to the sex linkage group in all of these species. BACs containing sdY have been isolated and sequenced in O.mykiss, and the genetic markers adjacent to sdY are not sex-linked in the other Oncorhynchus species, suggesting that the transposed region is very small. Possible explanations for the frequent occurrence of transposition of the SD are discussed.

DrML: Probabilistic Modeling of Gene Duplications

DrML is a software program for inferring evolutionary scenarios from a gene tree and a species tree with speciation time estimates that is based on a general maximum likelihood model. The program implements novel algorithms that efficiently infer most likely scenarios of gene duplication and loss events. Our comparative studies suggest that the general maximum likelihood model provides more credible estimates than standard parsimony reconciliation, especially when speciation times differ significantly. DrML is an open source project written in Python, and along with an on-line manual and sample data sets publicly available.