CALIBRATING DIVERGENCE TIMES ON SPECIES TREES VERSUS GENE TREES: IMPLICATIONS FOR SPECIATION HISTORY OF APHELOCOMA JAYS

Convergence in reduced body size, head size, and blood glucose in three island reptiles

Many oceanic islands harbor diverse species that differ markedly from their mainland relatives with respect to morphology, behavior, and physiology. A particularly common morphological change exhibited by a wide range of species on islands worldwide involves either a reduction in body size, termed island dwarfism, or an increase in body size, termed island gigantism. While numerous instances of dwarfism and gigantism have been well documented, documentation of other morphological changes on islands remains limited. Furthermore, we lack a basic understanding of the physiological mechanisms that underlie these changes, and whether they are convergent. A major hypothesis for the repeated evolution of dwarfism posits selection for smaller, more efficient body sizes in the context of low resource availability. Under this hypothesis, we would expect the physiological mechanisms known to be downregulated in model organisms exhibiting small body sizes due to dietary restriction or artificial selection would also be downregulated in wild species exhibiting dwarfism on islands. We measured body size, relative head size, and circulating blood glucose in three species of reptiles-two snakes and one lizard-in the California Channel Islands relative to mainland populations. Collating data from 6 years of study, we found that relative to mainland population the island populations had smaller body size (i.e., island dwarfism), smaller head sizes relative to body size, and lower levels of blood glucose, although with some variation by sex and year. These findings suggest that the island populations of these three species have independently evolved convergent physiological changes (lower glucose set point) corresponding to convergent changes in morphology that are consistent with a scenario of reduced resource availability and/or changes in prey size on the islands. This provides a powerful system to further investigate ecological, physiological, and genetic variables to elucidate the mechanisms underlying convergent changes in life history on islands.

Phylogenetic surveys on the newt genus Tylototriton sensu lato (Salamandridae, Caudata) reveal cryptic diversity and novel diversification promoted by historical climatic shifts

Global climatic transitions and Tibetan Plateau uplifts are hypothesized to have profoundly impacted biodiversity in southeastern Asia. To further test the hypotheses related to the impacts of these incidents, we investigated the diversification patterns of the newt genus Tylototriton sensu lato, distributed across the mountain ranges of southeastern Asia. Gene-tree and species-tree analyses of two mitochondrial genes and two nuclear genes revealed five major clades in the genus, and suggested several cryptic species. Dating estimates suggested that the genus originated in the early-to-middle Miocene. Under different species delimitating scenarios, diversification analyses with birth-death likelihood tests indicated that the genus held a higher diversification rate in the late Miocene-to-Pliocene era than that in the Pleistocene. Ancestral area reconstructions indicated that the genus originated from the northern Indochina Peninsula. Accordingly, we hypothesized that the Miocene Climatic Transition triggered the diversification of the genus, and the reinforcement of East Asian monsoons associated with the stepwise uplifts of the Tibetan Plateau promoted the radiation of the genus in southeastern Asia during the Miocene-to-Pliocene period. Quaternary glacial cycles likely had limited effects on speciation events in the genus, but mainly had contributions on their intraspecific differentiations.

Bayesian Phylogenetic Estimation of Clade Ages Supports Trans-Atlantic Dispersal of Cichlid Fishes

Divergence-time estimation based on molecular phylogenies and the fossil record has provided insights into fundamental questions of evolutionary biology. In Bayesian node dating, phylogenies are commonly time calibrated through the specification of calibration densities on nodes representing clades with known fossil occurrences. Unfortunately, the optimal shape of these calibration densities is usually unknown and they are therefore often chosen arbitrarily, which directly impacts the reliability of the resulting age estimates. As possible solutions to this problem, two nonexclusive alternative approaches have recently been developed, the “fossilized birth–death” (FBD) model and “total-evidence dating.” While these approaches have been shown to perform well under certain conditions, they require including all (or a random subset) of the fossils of each clade in the analysis, rather than just relying on the oldest fossils of clades. In addition, both approaches assume that fossil records of different clades in the phylogeny are all the product of the same underlying fossil sampling rate, even though this rate has been shown to differ strongly between higher level taxa. We here develop a flexible new approach to Bayesian age estimation that combines advantages of node dating and the FBD model. In our new approach, calibration densities are defined on the basis of first fossil occurrences and sampling rate estimates that can be specified separately for all clades. We verify our approach with a large number of simulated data sets, and compare its performance to that of the FBD model. We find that our approach produces reliable age estimates that are robust to model violation, on par with the FBD model. By applying our approach to a large data set including sequence data from over 1000 species of teleost fishes as well as 147 carefully selected fossil constraints, we recover a timeline of teleost diversification that is incompatible with previously assumed vicariant divergences of freshwater fishes. Our results instead provide strong evidence for transoceanic dispersal of cichlids and other groups of teleost fishes.

Disentangling Incomplete Lineage Sorting and Introgression to Refine Species-Tree Estimates for Lake Tanganyika Cichlid Fishes

Adaptive radiation is thought to be responsible for the evolution of a great portion of the past and present diversity of life. Instances of adaptive radiation, characterized by the rapid emergence of an array of species as a consequence to their adaptation to distinct ecological niches, are important study systems in evolutionary biology. However, because of the rapid lineage formation in these groups, and occasional gene flow between the participating species, it is often difficult to reconstruct the phylogenetic history of species that underwent an adaptive radiation. In this study, we present a novel approach for species-tree estimation in rapidly diversifying lineages, where introgression is known to occur, and apply it to a multimarker data set containing up to 16 specimens per species for a set of 45 species of East African cichlid fishes (522 individuals in total), with a main focus on the cichlid species flock of Lake Tanganyika. We first identified, using age distributions of most recent common ancestors in individual gene trees, those lineages in our data set that show strong signatures of past introgression. This led us to formulate three hypotheses of introgression between different lineages of Tanganyika cichlids: the ancestor of Boulengerochromini (or of Boulengerochromini and Bathybatini) received genomic material from the derived H-lineage; the common ancestor of Cyprichromini and Perissodini experienced, in turn, introgression from Boulengerochromini and/or Bathybatini; and the Lake Tanganyika Haplochromini and closely related riverine lineages received genetic material from Cyphotilapiini. We then applied the multispecies coalescent model to estimate the species tree of Lake Tanganyika cichlids, but excluded the lineages involved in these introgression events, as the multispecies coalescent model does not incorporate introgression. This resulted in a robust species tree, in which the Lamprologini were placed as sister lineage to the H-lineage (including the Eretmodini), and we identify a series of rapid splitting events at the base of the H-lineage. Divergence ages estimated with the multispecies coalescent model were substantially younger than age estimates based on concatenation, and agree with the geological history of the Great Lakes of East Africa. Finally, we formally tested the three hypotheses of introgression using a likelihood framework, and find strong support for introgression between some of the cichlid tribes of Lake Tanganyika.

Different diversity-dependent declines in speciation rate unbalances species richness in terrestrial slugs

Two genera of terrestrial slugs (Arion and Geomalacus) display a striking disproportion in species richness in the Iberian Peninsula. While there are 17 Iberian endemic species in Arion, morphological criteria only recognize four species within Geomalacus. Sequence data were used to test whether these differences could result from: (1) cryptic diversity within Geomalacus; (2) an earlier origin for Arion (older clades are expected to accumulate more species); (3) distinct patterns of diversification rates (higher initial speciation rates in Arion), and (4) some combination of the above factors (e.g., an older clade with higher speciation rates). Species delimitation tests based on mitochondrial and nuclear data revealed eight cryptic lineages within Geomalacus that lessened the asymmetry; nevertheless, the disparity required further investigation. No meaningful differences in crown group ages of each recovered clade were found. Regardless the different premises of the two equally plausible diversification models (similar initial speciation rates vs. higher initial speciation rates in Geomalacus), both coincide on diversity-dependent diversification for the two groups but weaker rate declines in Arion best explains the observed asymmetry in species richness. Also, the broader environmental tolerance combined with a faster dispersal and wider distribution may have represented an evolutionary advantage for Arion.

How and Why to Build a Unified Tree of Life

Phylogenetic trees are a crucial backbone for a wide breadth of biological research spanning systematics, organismal biology, ecology, and medicine. In 2015, the Open Tree of Life project published a first draft of a comprehensive tree of life, summarizing digitally available taxonomic and phylogenetic knowledge. This paper reviews, investigates, and addresses the following questions as a follow-up to that paper, from the perspective of researchers involved in building this summary of the tree of life: Is there a tree of life and should we reconstruct it? Is available data sufficient to reconstruct the tree of life? Do we have access to phylogenetic inferences in usable form? Can we combine different phylogenetic estimates across the tree of life? And finally, what is the future of understanding the tree of life?

Geographical features are the predominant driver of molecular diversification in widely distributed North American whipsnakes

Allopatric divergence following the formation of geographical features has been implicated as a major driver of evolutionary diversification. Widespread species complexes provide opportunities to examine allopatric divergence across varying degrees of isolation in both time and space. In North America, several geographical features may play such a role in diversification, including the Mississippi River, Pecos River, Rocky Mountains, Cochise Filter Barrier, Gulf of California and Isthmus of Tehuantepec. We used thousands of nuclear single nucleotide polymorphisms (SNPs) and mitochondrial DNA from several species of whipsnakes (genera Masticophis and Coluber) distributed across North and Central America to investigate the role that these geographical features have played on lineage divergence. We hypothesize that these features restrict gene flow and separate whipsnakes into diagnosable genomic clusters. We performed genomic clustering and phylogenetic reconstructions at the species and population levels using Bayesian and likelihood analyses and quantified migration levels across geographical features to assess the degree of genetic isolation due to allopatry. Our analyses suggest that (i) major genetic divisions are often consistent with isolation by geographical features, (ii) migration rates between clusters are asymmetrical across major geographical features, and (iii) areas that receive proportionally more migrants possess higher levels of genetic diversity. Collectively, our findings suggest that multiple features of the North American landscape contributed to allopatric divergence in this widely distributed snake group.

Phylogeographic investigation and ecological niche modelling of the endemic frog species Nanorana pleskei revealed multiple refugia in the eastern Tibetan Plateau

The largest plateau Tibetan Plateau supplied an excellent opportunity to investigate the influence of the Pleistocene events on the high-elevation species. To test for the alternative hypotheses of Pleistocene glacial refugia, we used partial sequences of two mitochondrial genes and one nuclear gene to examine the phylogeographic patterns of the endemic frog species Nanorana pleskei across its known range in the eastern Tibetan Plateau, and conducted species distribution modelling (SDM) to explore changes of its distribution range through current and paleo periods. In all data sets, the species was divided into lineage north occupying open plateau platform and lineage south colonizing the mountainous plateau. The divergence of two major clades was estimated at the early Pleistocene. In mtDNA, lineage north contained northeastern and northwestern sublineages, and lineage south had two overlapping-distributed sublineages. Different lineages possessed distinct demographic characteristics, i.e., subdivision in the northeastern sublineage, historical bottleneck effects and recent expansions in the northwestern sublineage and the southeastern sublineage. SDMs depicted that stable suitable habitats had existed in the upper-middle streams of the Yellow River, Dadu River, Jinsha River and Yalong River. These regions were also recognized as the ancestral areas of different lineages. In conclusion, Nanorana pleskei lineages have probably experienced long-term separations. Stable suitable habitats existing in upper-middle streams of major rivers on the eastern Tibetan Plateau and distinct demographic dynamics of different lineages indicated that the lineages possessed independent evolutionary processes in multiple glacial refugia. The findings verified the profound effects of Pleistocene climatic fluctuations on the plateau endemic species.

Epigean gammarids survived millions of years of severe climatic fluctuations in high latitude refugia throughout the Western Carpathians

Isolated glacial refugia have been documented in Central Europe for a number of taxa, but conclusive evidence for epigean aquatic species has remained elusive. Using molecular data (mitochondrial and nuclear markers), we compared the spatial patterns of lineage diversity of the widely distributed Gammarus fossarum species complex between two adjacent biogeographically and geomorphologically distinct Central European regions: the Bohemian Massif and the Western Carpathians. We investigated if the observed patterns of spatial diversity are more likely to stem from historical or present-day factors. Phylogenetic and phylogeographic analyses revealed eight phylogenetically diverse lineages: two exhibiting local signatures of recent demographic expansion inhabit both regions, while the other six display a relict distributional pattern and are found only in the Western Carpathians. Molecular dating indicates that these lineages are old and probably diverged throughout the Miocene (7-18Ma). Furthermore, their distribution does not seem to be constrained by the present boundaries of river catchments or topography. The contrasting spatial patterns of diversity observed between the two regions thus more likely result from historical rather than contemporaneous or recent factors. Our results indicate that despite the high latitude and proximity to the Pleistocene ice sheets, the Western Carpathians functioned as long-term glacial refugia for permanent freshwater fauna, allowing the uninterrupted survival of ancient lineages through millions of years of drastic climatic fluctuations.

Partial support for the central-marginal hypothesis within a population: reduced genetic diversity but not increased differentiation at the range edge of an island endemic bird

Large-scale population comparisons have contributed to our understanding of the evolution of geographic range limits and species boundaries, as well as the conservation value of populations at range margins. The central-marginal hypothesis (CMH) predicts a decline in genetic diversity and an increase in genetic differentiation toward the periphery of species' ranges due to spatial variation in genetic drift and gene flow. Empirical studies on a diverse array of taxa have demonstrated support for the CMH. However, nearly all such studies come from widely distributed species, and have not considered if the same processes can be scaled down to single populations. Here, we test the CMH on a species composed of a single population: the Island Scrub-Jay (Aphelocoma insularis), endemic to a 250 km² island. We examined microsatellite data from a quarter of the total population and found that homozygosity increased toward the island's periphery. However, peripheral portions of the island did not exhibit higher genetic differentiation. Simulations revealed that highly localized dispersal and small total population size, but not spatial variation in population density, were critical for generating fine-scale variation in homozygosity. Collectively, these results demonstrate that microevolutionary processes driving spatial variation in genetic diversity among populations can also be important for generating spatial variation in genetic diversity within populations.

Around the World in Eight Million Years: Historical Biogeography and Evolution of the Spray Zone Spider Amaurobioides (Araneae: Anyphaenidae)

Closely related organisms with transoceanic distributions have long been the focus of historical biogeography, prompting the question of whether long-distance dispersal, or tectonic-driven vicariance shaped their current distribution. Regarding the Southern Hemisphere continents, this question deals with the break-up of the Gondwanan landmass, which has also affected global wind and oceanic current patterns since the Miocene. With the advent of phylogenetic node age estimation and parametric bioinformatic advances, researchers have been able to disentangle historical evolutionary processes of taxa with greater accuracy. In this study, we used the coastal spider genus Amaurobioides to investigate the historical biogeographical and evolutionary processes that shaped the modern-day distribution of species of this exceptional genus of spiders. As the only genus of the subfamily Amaurobioidinae found on three Southern Hemisphere continents, its distribution is well-suited to study in the context of Gondwanic vicariance versus long-distance, transoceanic dispersal. Ancestral species of the genus Amaurobioides appear to have undergone several long-distance dispersal events followed by successful establishments and speciation, starting from the mid-Miocene through to the Pleistocene. The most recent common ancestor of all present-day Amaurobioides species is estimated to have originated in Africa after arriving from South America during the Miocene. From Africa the subsequent dispersals are likely to have taken place predominantly in an eastward direction. The long-distance dispersal events by Amaurobioides mostly involved transoceanic crossings, which we propose occurred by rafting, aided by the Antarctic Circumpolar Current and the West Wind Drift.

Hidden Genetic Diversity in an Asexually Reproducing Lichen Forming Fungal Group

Asexual species with vegetative propagation of both symbiont partners (soredia) in lichens may harbor lower species diversity because they may indeed represent evolutionary dead ends or clones. In this study we aim to critically examine species boundaries in the sorediate lichen forming fungi Parmotrema reticulatum–Parmotrema pseudoreticulatum complex applying coalescent-based approaches and other recently developed DNA-based methods. To this end, we gathered 180 samples from Africa, Asia, Australasia, Europe, North and South America and generated sequences of internal transcribed spacer of nuclear ribosomal DNA (ITS) and DNA replication licensing factor MCM7 (MCM7). The dataset was analysed using different approaches such as traditional phylogeny–maximum likelihood and Bayesian–genetic distances, automatic barcode gap discovery and coalescent-based methods–PTP, GMYC, spedeSTEM and *Beast–in order to test congruence among results. Additionally, the divergence times were also estimated to elucidate diversification events. Delimitations inferred from the different analyses are comparable with only minor differences, and following a conservative approach we propose that the sampled specimens of the P. reticulatum–P. pseudoreticulatum complex belong to at least eight distinct species-level lineages. Seven are currently classified under P. reticulatum and one as P. pseudoreticulatum. In this work we discuss one of only few examples of cryptic species that have so far been found in sorediate reproducing lichen forming fungi. Additionally our estimates suggest a recent origin of the species complex–during the Miocene. Consequently, the wide distribution of several of the cryptic species has to be explained by intercontinental long-distance dispersal events.

Revision of widespread red squirrels (genus: Tamiasciurus) highlights the complexity of speciation within North American forests

Integration of molecular methods, ecological modeling, and statistical hypothesis testing are increasing our understanding of differentiation within species and phylogenetic relationships among species by revealing environmental connections to evolutionary processes. Within mammals, novel diversity is being discovered and characterized as more complete geographic sampling is coupled with newer multi-disciplinary approaches. North American red squirrels exemplify a forest obligate genus whose species are monitored as indicators of forest ecosystem condition, yet phylogenetic relationships reflecting evolutionary history within this genus remain tentative. Through testing of competing systematic and niche-based divergence hypotheses, we recognize three species, Tamiasciurus douglasii, T. hudsonicus, and T. fremonti. Our data provide evidence of regional differences in evolutionary dynamics and continental gradients of complexity that are important both for future management and for investigating multiple pathways that can lead to the formation of new species.

Biogeographic barriers drive co-diversification within associated eukaryotes of the Sarracenia alata pitcher plant system

Understanding if the members of an ecological community have co-diversified is a central concern of evolutionary biology, as co-diversification suggests prolonged association and possible coevolution. By sampling associated species from an ecosystem, researchers can better understand how abiotic and biotic factors influence diversification in a region. In particular, studies of co-distributed species that interact ecologically can allow us to disentangle the effect of how historical processes have helped shape community level structure and interactions. Here we investigate the Sarracenia alata pitcher plant system, an ecological community where many species from disparate taxonomic groups live inside the fluid-filled pitcher leaves. Direct sequencing of the eukaryotes present in the pitcher plant fluid enables us to better understand how a host plant can shape and contribute to the genetic structure of its associated inquilines, and to ask whether genetic variation in the taxa are structured in a similar manner to the host plant. We used 454 amplicon-based metagenomics to demonstrate that the pattern of genetic diversity in many, but not all, of the eukaryotic community is similar to that of S. alata, providing evidence that associated eukaryotes share an evolutionary history with the host pitcher plant. Our work provides further evidence that a host plant can influence the evolution of its associated commensals.

Multilocus coalescent species delimitation reveals widespread cryptic differentiation among Drakensberg mountain-living freshwater crabs (Decapoda : Potamonautes)

Cryptic lineages present major challenges for evolutionary and conservation studies, particularly where these lineages remain undiscovered. Freshwater crabs are known to harbour cryptic diversity, in most cases with limited morphological differences. During the present study, we used a multilocus (12S rRNA, 16S rRNA, COI, 28S rRNA, DecapANT and PEPCK) Bayesian species delimitation to examine cryptic diversity within a freshwater crab species complex (Potamonautes clarus/P. depressus). We sampled 25 highland rivers in the Tugela and uMkomazi River drainage systems of the Drakensberg Mountain range, in the KwaZulu–Natal province of South Africa. Our results showed there to be at least eight lineages: six novel potamonautid freshwater crabs, and two described taxa P. clarus and P. depressus. Divergence from the most recent common ancestor occurred between the mid- and late Miocene (12.1 Mya), while divergence within the species complex occurred ~10.3 Mya up until the Holocene (0.11 Mya). The discovery of six novel lineages of freshwater crabs from a seemingly restricted distribution range has conservation implications, but to date most conservation planning strategies have focussed on freshwater vertebrates. By conducting a fine-scale phylogenetic survey using invertebrates, this study provides a platform for the inclusion of freshwater invertebrates in future conservation assessments.

Anonymous nuclear markers reveal taxonomic incongruence and long-term disjunction in a cactus species complex with continental-island distribution in South America

The Pilosocereus aurisetus complex consists of eight cactus species with a fragmented distribution associated to xeric enclaves within the Cerrado biome in eastern South America. The phylogeny of these species is incompletely resolved, and this instability complicates evolutionary analyses. Previous analyses based on both plastid and microsatellite markers suggested that this complex contained species with inherent phylogeographic structure, which was attributed to recent diversification and recurring range shifts. However, limitations of the molecular markers used in these analyses prevented some questions from being properly addressed. In order to better understand the relationship among these species and make a preliminary assessment of the genetic structure within them, we developed anonymous nuclear loci from pyrosequencing data of 40 individuals from four species in the P. aurisetus complex. The data obtained from these loci were used to identify genetic clusters within species, and to investigate the phylogenetic relationship among these inferred clusters using a species tree methodology. Coupled with a palaeodistributional modelling, our results reveal a deep phylogenetic and climatic disjunction between two geographic lineages. Our results highlight the importance of sampling more regions from the genome to gain better insights on the evolution of species with an intricate evolutionary history. The methodology used here provides a feasible approach to develop numerous genealogical molecular markers throughout the genome for non-model species. These data provide a more robust hypothesis for the relationship among the lineages of the P. aurisetus complex.

Sequence capture of ultraconserved elements from bird museum specimens

New DNA sequencing technologies are allowing researchers to explore the genomes of the millions of natural history specimens collected prior to the molecular era. Yet, we know little about how well specific next-generation sequencing (NGS) techniques work with the degraded DNA typically extracted from museum specimens. Here, we use one type of NGS approach, sequence capture of ultraconserved elements (UCEs), to collect data from bird museum specimens as old as 120 years. We targeted 5060 UCE loci in 27 western scrub-jays (Aphelocoma californica) representing three evolutionary lineages that could be species, and we collected an average of 3749 UCE loci containing 4460 single nucleotide polymorphisms (SNPs). Despite older specimens producing fewer and shorter loci in general, we collected thousands of markers from even the oldest specimens. More sequencing reads per individual helped to boost the number of UCE loci we recovered from older specimens, but more sequencing was not as successful at increasing the length of loci. We detected contamination in some samples and determined that contamination was more prevalent in older samples that were subject to less sequencing. For the phylogeny generated from concatenated UCE loci, contamination led to incorrect placement of some individuals. In contrast, a species tree constructed from SNPs called within UCE loci correctly placed individuals into three monophyletic groups, perhaps because of the stricter analytical procedures used for SNP calling. This study and other recent studies on the genomics of museum specimens have profound implications for natural history collections, where millions of older specimens should now be considered genomic resources.

Multilocus species trees and species delimitation in a temporal context: application to the water shrews of the genus Neomys

Background

Multilocus data are becoming increasingly important in determining the phylogeny of closely related species and delimiting species. In species complexes where unequivocal fossil calibrations are not available, rigorous dating of the coalescence-based species trees requires accurate mutation rates of the loci under study but, generally, these rates are unknown. Here, we obtained lineage-specific mutation rates of these loci from a higher-level phylogeny with a reliable fossil record and investigated how different choices of mutation rates and species tree models affected the split time estimates. We implemented this strategy with a genus of water shrews, Neomys, whose taxonomy has been contentious over the last century.

Results

We sequenced 13 introns and cytochrome b from specimens of the three species currently recognized in this genus including two subspecies of N. anomalus that were originally described as species. A Bayesian multilocus species delimitation method and estimation of gene flow supported that these subspecies are distinct evolutionary lineages that should be treated as distinct species: N. anomalus (sensu stricto), limited to part of the Iberian Peninsula, and N. milleri, with a larger Eurasian range. We then estimated mutation rates from a Bayesian relaxed clock analysis of the mammalian orthologues with several fossil calibrations. Next, using the estimated Neomys-specific rates for each locus in an isolation-with-migration model, the split time for these sister taxa was dated at 0.40 Myr ago (with a 95 % confidence interval of 0.26 – 0.86 Myr), likely coinciding with one of the major glaciations of the Middle Pleistocene. We also showed that the extrapolation of non-specific rates or the use of simpler models would lead to very different split time estimates.

Conclusions

We showed that the estimation of rigorous lineage-specific mutation rates for each locus allows the inference of robust split times in a species tree framework. These times, in turn, afford a better understanding of the timeframe required to achieve isolation and, eventually, speciation in sister lineages. The application of species delimitation methods and an accurate dating strategy to the genus Neomys helped to clarify its controversial taxonomy.

Electronic supplementary material

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

Formal Comment to Pettengill: The Time to Most Recent Common Ancestor Does Not (Usually) Approximate the Date of Divergence

In 2013 Zhou et al. concluded that Salmonella enterica serovar Agona represents a genetically monomorphic lineage of recent ancestry, whose most recent common ancestor existed in 1932, or earlier. The Abstract stated 'Agona consists of three lineages with minimal mutational diversity: only 846 single nucleotide polymorphisms (SNPs) have accumulated in the non-repetitive, core genome since Agona evolved in 1932 and subsequently underwent a major population expansion in the 1960s.' These conclusions have now been criticized by Pettengill, who claims that the evolutionary models used to date Agona may not have been appropriate, the dating estimates were inaccurate, and the age of emergence of Agona should have been qualified by an upper limit reflecting the date of its divergence from an outgroup, serovar Soerenga. We dispute these claims. Firstly, Pettengill's analysis of Agona is not justifiable on technical grounds. Secondly, an upper limit for divergence from an outgroup would only be meaningful if the outgroup were closely related to Agona, but close relatives of Agona are yet to be identified. Thirdly, it is not possible to reliably date the time of divergence between Agona and Soerenga. We conclude that Pettengill's criticism is comparable to a tempest in a teapot.

The effects of recombination, mutation and selection on the evolution of the Rp1 resistance genes in grasses

Plant immune genes, or resistance genes, are involved in a co-evolutionary arms race with a diverse range of pathogens. In agronomically important grasses, such R genes have been extensively studied because of their role in pathogen resistance and in the breeding of resistant cultivars. In this study, we evaluate the importance of recombination, mutation and selection on the evolution of the R gene complex Rp1 of Sorghum, Triticum, Brachypodium, Oryza and Zea. Analyses show that recombination is widespread, and we detected 73 independent instances of sequence exchange, involving on average 1567 of 4692 nucleotides analysed (33.4%). We were able to date 24 interspecific recombination events and found that four occurred postspeciation, which suggests that genetic introgression took place between different grass species. Other interspecific events seemed to have been maintained over long evolutionary time, suggesting the presence of balancing selection. Significant positive selection (i.e. a relative excess of nonsynonymous substitutions (dN /dS >1)) was detected in 17-95 codons (0.42-2.02%). Recombination was significantly associated with areas with high levels of polymorphism but not with an elevated dN /dS ratio. Finally, phylogenetic analyses show that recombination results in a general overestimation of the divergence time (mean = 14.3%) and an alteration of the gene tree topology if the tree is not calibrated. Given that the statistical power to detect recombination is determined by the level of polymorphism of the amplicon as well as the number of sequences analysed, it is likely that many studies have underestimated the importance of recombination relative to the mutation rate.

Multilocus phylogenetic analyses of Hispaniolan and Bahamian trunk anoles (distichus species group)

The distichus species group includes six species and 21 subspecies of trunk ecomorph anoles distributed across Hispaniola and its satellite islands as well as the northern Bahamas. Although this group has long served as a model system for studies of reproductive character displacement, adaptation, behavior and speciation, it has never been the subject of a comprehensive phylogenetic analysis. Our goal here is to generate a multilocus phylogenetic dataset (one mitochondrial and seven nuclear loci) and to use this dataset to infer phylogenetic relationships among the majority of the taxa assigned to the distichus species group. We use these phylogenetic trees to address three topics about the group's evolution. First, we consider longstanding taxonomic controversies about the status of several species and subspecies assigned to the distichus species group. Second, we investigate the biogeographic history of the group and specifically test the hypotheses that historical division of Hispaniola into two paleo-islands contributed to the group's diversification and that Bahamian and Hispaniolan satellite island populations are derived from colonists from the main Hispaniolan landmass. Finally, third, we use comparative phylogenetic analyses to test the hypothesis that divergence between pale yellow and darkly pigmented orange or red dewlap coloration has occurred repeatedly across the distichus species group.

Islands within an island: repeated adaptive divergence in a single population

Physical barriers to gene flow were once viewed as prerequisites for adaptive evolutionary divergence. However, a growing body of theoretical and empirical work suggests that divergence can proceed within a single population. Here we document genetic structure and spatially replicated patterns of phenotypic divergence within a bird species endemic to 250 km(2) Santa Cruz Island, California, USA. Island scrub-jays (Aphelocoma insularis) in three separate stands of pine habitat had longer, shallower bills than jays in oak habitat, a pattern that mirrors adaptive differences between allopatric populations of the species' mainland congener. Variation in both bill measurements was heritable, and island scrub-jays mated nonrandomly with respect to bill morphology. The population was not panmictic; instead, we found a continuous pattern of isolation by distance across the east-west axis of the island, as well as a subtle genetic discontinuity across the boundary between the largest pine stand and adjacent oak habitat. The ecological factors that appear to have facilitated adaptive differentiation at such a fine scale--environmental heterogeneity and localized dispersal--are ubiquitous in nature. These findings support recent arguments that microgeographic patterns of adaptive divergence may be more common than currently appreciated, even in mobile taxonomic groups like birds.

The phylogeny of pikas (Ochotona) inferred from a multilocus coalescent approach

The clarification of the systematics of pikas (genus Ochotona) has been hindered by largely overlapping morphological characters among species and the lack of a comprehensive molecular phylogeny. Here we estimate the first multilocus phylogeny of the genus to date, by analysing 12 nuclear DNA markers (total of 7.5Kb) in 11 species of pikas from the four classified subgenera (Pika, Ochotona, Lagotona and Conothoa) using a multispecies coalescent-based framework. The species-tree confirmed the subgeneric classification by retrieving as monophyletic the subgenera represented here by more than one species. Contrary to previous phylogenies based on mtDNA alone, Lagotona was found to be sister to Pika. Also, support for the monophyly of the alpina group was not strong, thus caution should be used in future analyses of this group. A relaxed molecular clock calibrated using the Ochotonidae-Leporidae divergence resulted in more recent estimates of divergence times relative to previous studies. Strong concordance with inferences based on fossil records was found, suggesting that the initial diversification of the genus took place by the end of late Miocene. Finally, this work sets up methodologies and gathers molecular markers that can be used to extend the understanding of the evolutionary history of the genus.

The roles of history and ecology in chloroplast phylogeographic patterns of the bird-dispersed plant parasite Phoradendron californicum (Viscaceae) in the Sonoran Desert

•

PREMISE OF THE STUDY

A recurrent explanation for phylogeographic discontinuities in the Baja California Peninsula and the Sonoran Desert Region has been the association of vicariant events with Pliocene and Pleistocene seaway breaks. Nevertheless, despite its relevance for plant dispersal, other explanations such as ecological and paleoclimatic factors have received little attention. Here, we analyzed the role of several of these factors to describe the phylogeographic patterns of the desert mistletoe, Phoradendron californicum.•

METHODS

Using noncoding chloroplast regions, we assess the marginal probability of 19 a priori hypotheses related to geological and ecological factors to predict the cpDNA variation in P. californicum using a Bayesian coalescent framework. Complementarily, we used the macrofossil record and niche model projections on Last Glacial Maximum climatic conditions for hosts, mistletoe, and a bird specialist to interpret phylogeographic patterns.•

KEY RESULTS

Genealogical reconstructions revealed five clades, which suggest a combination of cryptic divergence, long-distance seed dispersal, and isolating postdivergence events. Bayesian hypothesis test favored a series of Pliocene and Pleistocene geological events related to the formation of the Baja California Peninsula and seaways across the peninsula as the most supported explanation for this genealogical pattern. However, age estimates, niche projections, and fossil records show dynamic host-mistletoe interactions and evidence of host races, indicating that ecological and geological factors have been interacting during the formation and structuring of phylogeographic divergence.•

CONCLUSIONS

Variation in cpDNA across the species range results from the interplay of vicariant events, past climatic oscillations, and more dynamic factors related to ecological processes at finer temporal and spatial scales.

Spiders on a Hot Volcanic Roof: Colonisation Pathways and Phylogeography of the Canary Islands Endemic Trap-Door Spider Titanidiops canariensis (Araneae, Idiopidae)

Studies conducted on volcanic islands have greatly contributed to our current understanding of how organisms diversify. The Canary Islands archipelago, located northwest of the coast of northern Africa, harbours a large number of endemic taxa. Because of their low vagility, mygalomorph spiders are usually absent from oceanic islands. The spider Titanidiops canariensis, which inhabits the easternmost islands of the archipelago, constitutes an exception to this rule. Here, we use a multi-locus approach that combines three mitochondrial and four nuclear genes to investigate the origins and phylogeography of this remarkable trap-door spider. We provide a timeframe for the colonisation of the Canary Islands using two alternative approaches: concatenation and species tree inference in a Bayesian relaxed clock framework. Additionally, we investigate the existence of cryptic species on the islands by means of a Bayesian multi-locus species delimitation method. Our results indicate that T. canariensis colonised the Canary Islands once, most likely during the Miocene, although discrepancies between the timeframes from different approaches make the exact timing uncertain. A complex evolutionary history for the species in the archipelago is revealed, which involves two independent colonisations of Fuerteventura from the ancestral range of T. canariensis in northern Lanzarote and a possible back colonisation of southern Lanzarote. The data further corroborate a previously proposed volcanic refugium, highlighting the impact of the dynamic volcanic history of the island on the phylogeographic patterns of the endemic taxa. T. canariensis includes at least two different species, one inhabiting the Jandia peninsula and central Fuerteventura and one spanning from central Fuerteventura to Lanzarote. Our data suggest that the extant northern African Titanidiops lineages may have expanded to the region after the islands were colonised and, hence, are not the source of colonisation. In addition, T. maroccanus may harbour several cryptic species.

Influence of gene flow on divergence dating - implications for the speciation history of Takydromus grass lizards

Dating the time of divergence and understanding speciation processes are central to the study of the evolutionary history of organisms but are notoriously difficult. The difficulty is largely rooted in variations in the ancestral population size or in the genealogy variation across loci. To depict the speciation processes and divergence histories of three monophyletic Takydromus species endemic to Taiwan, we sequenced 20 nuclear loci and combined with one mitochondrial locus published in GenBank. They were analysed by a multispecies coalescent approach within a Bayesian framework. Divergence dating based on the gene tree approach showed high variation among loci, and the divergence was estimated at an earlier date than when derived by the species-tree approach. To test whether variations in the ancestral population size accounted for the majority of this variation, we conducted computer inferences using isolation-with-migration (IM) and approximate Bayesian computation (ABC) frameworks. The results revealed that gene flow during the early stage of speciation was strongly favoured over the isolation model, and the initiation of the speciation process was far earlier than the dates estimated by gene- and species-based divergence dating. Due to their limited dispersal ability, it is suggested that geographical isolation may have played a major role in the divergence of these Takydromus species. Nevertheless, this study reveals a more complex situation and demonstrates that gene flow during the speciation process cannot be overlooked and may have a great impact on divergence dating. By using multilocus data and incorporating Bayesian coalescence approaches, we provide a more biologically realistic framework for delineating the divergence history of Takydromus.

Evolution at a different pace: distinctive phylogenetic patterns of cone snails from two ancient oceanic archipelagos

Ancient oceanic archipelagos of similar geological age are expected to accrue comparable numbers of endemic lineages with identical life history strategies, especially if the islands exhibit analogous habitats. We tested this hypothesis using marine snails of the genus Conus from the Atlantic archipelagos of Cape Verde and Canary Islands. Together with Azores and Madeira, these archipelagos comprise the Macaronesia biogeographic region and differ remarkably in the diversity of this group. More than 50 endemic Conus species have been described from Cape Verde, whereas prior to this study, only two nonendemic species, including a putative species complex, were thought to occur in the Canary Islands. We combined molecular phylogenetic data and geometric morphometrics with bathymetric and paleoclimatic reconstructions to understand the contrasting diversification patterns found in these regions. Our results suggest that species diversity is even lower than previously thought in the Canary Islands, with the putative species complex corresponding to a single species, Conus guanche. One explanation for the enormous disparity in Conus diversity is that the amount of available habitat may differ, or may have differed in the past due to eustatic (global) sea level changes. Historical bathymetric data, however, indicated that sea level fluctuations since the Miocene have had a similar impact on the available habitat area in both Cape Verde and Canary archipelagos and therefore do not explain this disparity. We suggest that recurrent gene flow between the Canary Islands and West Africa, habitat losses due to intense volcanic activity in combination with unsuccessful colonization of new Conus species from more diverse regions, were all determinant in shaping diversity patterns within the Canarian archipelago. Worldwide Conus species diversity follows the well-established pattern of latitudinal increase of species richness from the poles towards the tropics. However, the eastern Atlantic revealed a striking pattern with two main peaks of Conus species richness in the subtropical area and decreasing diversities toward the tropical western African coast. A Random Forests model using 12 oceanographic variables suggested that sea surface temperature is the main determinant of Conus diversity either at continental scales (eastern Atlantic coast) or in a broader context (worldwide). Other factors such as availability of suitable habitat and reduced salinity due to the influx of large rivers in the tropical area also play an important role in shaping Conus diversity patterns in the western coast of Africa.

Evaluating evolutionary history in the face of high gene tree discordance in Australian Gehyra (Reptilia: Gekkonidae)

Species tree methods have provided improvements for estimating species relationships and the timing of diversification in recent radiations by allowing for gene tree discordance. Although gene tree discordance is often observed, most discordance is attributed to incomplete lineage sorting rather than other biological phenomena, and the causes of discordance are rarely investigated. We use species trees from multi-locus data to estimate the species relationships, evolutionary history and timing of diversification among Australian Gehyra-a group renowned for taxonomic uncertainty and showing a large degree of gene tree discordance. We find support for a recent Asian origin and two major clades: a tropically adapted clade and an arid adapted clade, with some exceptions, but no support for allopatric speciation driven by chromosomal rearrangement in the group. Bayesian concordance analysis revealed high gene tree discordance and comparisons of Robinson-Foulds distances showed that discordance between gene trees was significantly higher than that generated by topological uncertainty within each gene. Analysis of gene tree discordance and incomplete taxon sampling revealed that gene tree discordance was high whether terminal taxon or gene sampling was maximized, indicating discordance is due to biological processes, which may be important in contributing to gene tree discordance in many recently diversified organisms.

Late Miocene lineage divergence and ecological differentiation of rare endemic Juniperus blancoi: clues for the diversification of North American conifers

Western North America and Mexico contain a large number of conifer species. This diversity could be the product of orographic and climate changes of the late Tertiary and Quaternary. In this study, we focus on the evolutionary history of Juniperus blancoi, in order to determine the impact of climate change and environmental heterogeneity on population differentiation. We estimated the population structure, phylogenetic relationships and historical demography of J. blancoi populations using nuclear genes. We correlated genetic structure with ecological differentiation, divergence times and changes in population size. Populations of J. blancoi are differentiated into three lineages that correspond to low-, mid- and high-altitude populations. The three groups diversified in the late Miocene, early Pliocene, with only a few events of gene flow since then. Two lineages in the north exhibited a pattern of population growth during the Pleistocene that could be linked to climate changes. Populations of J. blancoi experienced significant ecological differentiation and early divergence events, which correspond to periods of global cooling and mountain uplift during the Miocene. This suggests that mountain ranges in tropical and subtropical latitudes play an important role in the speciation and persistence of conifer taxa in diversity hotspots, by providing diverse environmental conditions.