Coalescent-based Approaches Research Articles

Comparison of colonial volvocine algae based on phylotranscriptomic analysis of gene family evolution and natural selection

ABSTRACT This study is the first to determine the transcriptomes of eight colonial volvocine algae: Pandorina morum, P. colemaniae, Volvulina compacta, Eudorina elegans, E. cylindrica, Yamagishiella unicocca, Colemanosphaera charkowiensis and C. angeleri. By analysing these data together with the whole-genome data of Chlamydomonas reinhardtii and other colonial volvocine algae, we determined 26 979 gene families and 105 single-copy orthogroups. Using the single-copy orthogroups, we adopted a supermatrix and coalescent-based approach to perform phylotranscriptomic analysis of colonial volvocine algae. The phylogenetic relationships of most species were consistent with those of a previous study, and the results showed paraphyly of the genus Pandorina. Then, we analysed gene family expansion and contraction in colonial volvocine algae and found that they exhibited a large number of expanded and contracted gene families. In particular, the common ancestor of the Volvocaceae had the largest number of expanded gene families, which may be associated with the increased volume of extracellular matrix. Finally, using the ML method to analyse the substitution rates among colonial volvocine algae based on single-copy orthogroups, we showed a significantly higher non-synonymous substitution (dN) rate in Tetrabaena socialis. By analysing the differences in the dN/dS ratio, we found that more orthogroups exhibited a significantly higher dN/dS in T. socialis; we speculate that the higher dN/dS could indicate relaxed selection pressure and that this lower constraint allows the accumulation of more non-synonymous mutations in T. socialis.

Influence of directional positive Darwinian selection-driven evolution on arboviruses Dengue and Zika virulence and pathogenesis

Dengue (DENV) and Zika (ZIKV) viruses are antigenically and evolutionarily related; immunological cross-reactions between them have been associated to both cross-protection and infection-enhanced mechanisms. Here, DENV-1-4 and ZIKV were investigated through Bayesian coalescent-based approaches and selection-driven Darwinian evolution methods using robust datasets. Our findings show that both DENV and ZIKV, driven essentially by directional positive selection, have undergone evolution and diversification and that their entire polyproteins are subject to an intense directional evolution. Interestingly, positively selected codons mapped here are directly associated to DENV-1-2 virulence as well as the ZIKV burgeoning 2015–16 outbreak in the Americas, therefore, having impact on the pathogenesis of these viruses. Biochemical prediction analysis focusing on markers involved in virulence and viral transmission dynamics identified alterations in N-Glycosylation-, Phosphorylation- and Palmitoylation-sites in ZIKV sampled from different countries, hosts and isolation sources. Taking into account both DENV-ZIKV co-circulation either into and/or out of flavivirus-endemic regions, as well as recombination and quasispecies scenarios, these results indicate the action of a selection-driven evolution affecting the biology, virulence and pathogenesis of these pathogens in a non-randomized environment.

Coalescent-based species delimitation in the sand lizards of the Liolaemus wiegmannii complex (Squamata: Liolaemidae)

Coalescent-based algorithms coupled with the access to genome-wide data have become powerful tools for assessing questions on recent or rapid diversification, as well as delineating species boundaries in the absence of reciprocal monophyly. In southern South America, the diversification of Liolaemus lizards during the Pleistocene is well documented and has been attributed to the climatic changes that characterized this recent period of time. Past climatic changes had harsh effects at extreme latitudes, including Patagonia, but habitat changes at intermediate latitudes of South America have also been recorded, including expansion of sand fields over northern Patagonia and Pampas). In this work, we apply a coalescent-based approach to study the diversification of the Liolaemus wiegmannii species complex, a morphologically conservative clade that inhabits sandy soils across northwest and south-central Argentina, and the south shores of Uruguay. Using four standard sequence markers (mitochondrial DNA and three nuclear loci) along with ddRADseq data we inferred species limits and a time-calibrated species tree for the L. wiegmannii complex in order to evaluate the influence of Quaternary sand expansion/retraction cycles on diversification. We also evaluated the evolutionary independence of the recently described L. gardeli and inferred its phylogenetic position relative to L. wiegmannii. We find strong evidence for six allopatric candidate species within L. wiegmannii, which diversified during the Pleistocene. The Great Patagonian Glaciation (∼1 million years before present) likely split the species complex into two main groups: one composed of lineages associated with sub-Andean sedimentary formations, and the other mostly related to sand fields in the Pampas and northern Patagonia. We hypothesize that early speciation within L. wiegmannii was influenced by the expansion of sand dunes throughout central Argentina and Pampas. Finally, L. gardeli is supported as a distinct lineage nested within the L. wiegmannii complex.

A congruent topology for deep gastropod relationships.

Gastropod molluscs are among the most diverse and abundant animals in the oceans, and are successful colonizers of terrestrial and freshwater environments. Past phylogenetic efforts to resolve gastropod relationships resulted in a range of conflicting hypotheses. Here, we use phylogenomics to address deep relationships between the five major gastropod lineages—Caenogastropoda, Heterobranchia, Neritimorpha, Patellogastropoda and Vetigastropoda—and provide one congruent and well-supported topology. We substantially expand taxon sampling for outgroups and for previously underrepresented gastropod lineages, presenting new transcriptomes for neritimorphs and patellogastropods. We conduct analyses under maximum-likelihood, Bayesian inference and a coalescent-based approach, accounting for the most pervasive sources of systematic errors in large datasets: compositional heterogeneity, site heterogeneity, heterotachy, variation in evolutionary rates among genes, matrix completeness, outgroup choice and gene tree conflict. We find that vetigastropods and patellogastropods are sister taxa, and that neritimorphs are the sister group to caenogastropods and heterobranchs. We name these two major unranked clades Psilogastropoda and Angiogastropoda, respectively. We additionally provide the first genomic-scale data for internal relationships of neritimorphs and patellogastropods. Our results highlight the need for reinterpreting the evolution of morphological and developmental characters in gastropods, especially for inferring their ancestral states.

Gene Selection and Evolutionary Modeling Affect Phylogenomic Inference of Neuropterida Based on Transcriptome Data.

Neuropterida is a super order of Holometabola that consists of the orders Megaloptera (dobsonflies, fishflies, and alderflies), Neuroptera (lacewings) and Raphidioptera (snakeflies). Several proposed higher-level relationships within Neuropterida, such as the relationships between the orders or between the families, have been extensively debated. To further understand the evolutionary history of Neuropterida, we conducted phylogenomic analyses of all 13 published transcriptomes of the neuropterid species, as well as of a new transcriptome of the fishfly species Ctenochauliodes similis of Liu and Yang, 2006 (Megaloptera: Corydalidae: Chauliodinae) that we sequenced. Our phylogenomic data matrix contained 1392 ortholog genes from 22 holometabolan species representing six families from Neuroptera, two families from Raphidioptera, and two families from Megaloptera as the ingroup taxa, and nine orders of Holometabola as outgroups. Phylogenetic reconstruction was performed using both concatenation and coalescent-based approaches under a site-homogeneous model as well as under a site-heterogeneous model. Surprisingly, analyses using the site-homogeneous model strongly supported a paraphyletic Neuroptera, with Coniopterygidae assigned as the sister group of all other Neuropterida. In contrast, analyses using the site-heterogeneous model recovered Neuroptera as monophyletic. The monophyly of Neuroptera was also recovered in concatenation and coalescent-based analyses using genes with stronger phylogenetic signals [i.e., higher average bootstrap support (ABS) values and higher relative tree certainty including all conflicting bipartitions (RTCA) values] under the site-homogeneous model. The present study illustrated how both data selection and model selection influence phylogenomic analyses of large-scale data matrices comprehensively.

Large ancestral effective population size explains the difficult phylogenetic placement of owl monkeys.

The phylogenetic position of owl monkeys, grouped in the genus Aotus, has been a controversial issue for understanding Neotropical primate evolution. Explanations of the difficult phylogenetic assignment of owl monkeys have been elusive, frequently relying on insufficient data (stochastic error) or scenarios of rapid speciation (adaptive radiation) events. Using a coalescent-based approach, we explored the population-level mechanisms likely explaining these topological discrepancies. We examined the topological variance of 2,192 orthologous genes shared between representatives of the three major Cebidae lineages and the outgroup. By employing a methodological framework that allows for reticulated tree topologies, our analysis explicitly tested for non-dichotomous evolutionary processes impacting the finding of the position of owl monkeys in the cebid phylogeny. Our findings indicated that Aotus is a sister lineage of the callitrichines. Most gene trees (>50%) failed to recover the species tree topology, although the distribution of gene trees mismatching the true species topology followed the standard expectation of the multispecies coalescent without reticulation. We showed that the large effective population size of the common ancestor of Aotus and callitrichines was the most likely factor responsible for generating phylogenetic uncertainty. On the other hand, fast speciation scenarios or introgression played minor roles. We propose that the difficult phylogenetic placement of Aotus is explained by population-level processes associated with the large ancestral effective size. These results shed light on the biogeography of the early cebid diversification in the Miocene, highlighting the relevance of evaluating phylogenetic relationships employing population-aware approaches.

Phylogeny and species delimitation of near Eastern Neurergus newts (Salamandridae) based on genome-wide RADseq data analysis

We reconstruct the molecular phylogeny of Near Eastern mountain brook newts of the genus Neurergus (family Salamandridae) based on newly determined RADseq data, and compare the outcomes of concatenation-based phylogenetic reconstruction with species-tree inference. Furthermore, we test the current taxonomy of Neurergus (with four species: Neurergus strauchii, N. crocatus, N. kaiseri, and N. derjugini) against coalescent-based species-delimitation approaches of our genome-wide genetic data set. While the position of N. strauchii as sister species to all other Neurergus species was consistent in all of our analyses, the phylogenetic relationships between the three remaining species changed depending on the applied method. The concatenation approach, as well as quartet-based species-tree inference, supported a topology with N. kaiseri as the closest relative to N. derjugini, while full-coalescent species-tree inference approaches supported N. crocatus as sister species of N. derjugini. Investigating the individual signal of gene trees highlighted an extensive variation among gene histories, most likely resulting from incomplete lineage sorting. Coalescent-based species-delimitation models suggest that the current taxonomy might underestimate the species richness within Neurergus and supports seven species. Based on the current sampling, our analysis suggests that N. strauchii, N. derjugini and N. kaiseri might each be subdivided into further species. However, as amphibian species are known to be composed of deep conspecific lineages that do not always warrant species status, these results need to be cautiously interpreted in an integrative taxonomic framework. We hypothesize that the rather shallow divergences detected within N. kaiseri and N. derjugini likely reflect an ongoing speciation process and thus require further investigation. On the contrary, the much deeper genetic divergence found between the two morphologically and geographically differentiated subspecies of N. strauchii leads us to propose that N. s. barani should be considered a distinct species, Neurergus barani Öz, 1994.

Anchored hybrid enrichment phylogenomics resolves the backbone of erebine moths

The subfamily Erebinae (Lepidoptera, Erebidae) includes approximately 10,000 species with many still undescribed. It is one of the most diverse clades within the moth superfamily Noctuoidea and encompasses a diversity of ecological habits. Erebine caterpillars feed on a broad range of host plants including several economically important crops. Adults possess a unique array of adaptations for predator defense, including some of the most sensitive hearing organs (tympana) across the Lepidoptera and striking wing coloration to startle visual predators. Despite the relevance of these moths to agriculture and ecological research, a robust phylogenetic framework is lacking. Here we used anchored hybrid enrichment, a relatively new approach in phylogenomics, to resolve relationships among the subfamily. Using the recently developed Lep1 anchored hybrid enrichment probe set, 658 gene fragments with an average length of 320 bp were captured from an exemplar set of 75 erebine species, representing 73 genera and 23 tribes. While the total number of erebine tribes is not firmly established, this represents at least 75% of known tribal level diversity. Anchored hybrid enrichment data were partitioned by locus and by codon position for maximum likelihood phylogenetic analysis and coalescent-based species-tree approaches. Results from our study provided strong nodal support (BP ≥ 95) for nearly all nodes in the partitioned ML tree, solidifying many relationships that were previously uncertain or moderately supported based on morphology or a smaller number of gene fragments. Likelihood analyses confidently resolved the placement of Acantholipini as a sister tribe to Sypnini and all other Erebinae. The remaining tribes were placed in a single, strongly supported clade split into two major subclades. Additionally, 25 tropical species that did not have previous tribal assignments are confidently placed on the phylogeny. Statistical comparisons with Shimodaira-Hasegawa (SH) tests found that our maximum likelihood trees were significantly more likely than alternative hypotheses. This study demonstrates the utility of anchored phylogenomics for resolving relationships within subfamilies of Lepidoptera.

Anchored phylogenomics unravels the evolution of spider flies (Diptera, Acroceridae) and reveals discordance between nucleotides and amino acids

The onset of phylogenomics has contributed to the resolution of numerous challenging evolutionary questions while offering new perspectives regarding biodiversity. However, in some instances, analyses of large genomic datasets can also result in conflicting estimates of phylogeny. Here, we present the first phylogenomic scale study of a dipteran parasitoid family, built upon anchored hybrid enrichment and transcriptomic data of 240 loci of 43 ingroup acrocerid taxa. A new hypothesis for the timing of spider fly evolution is proposed, wielding recent advances in divergence time dating, including the fossilized birth-death process to show that the origin of Acroceridae is younger than previously proposed. To test the robustness of our phylogenetic inferences, we analyzed our datasets using different phylogenetic estimation criteria, including supermatrix and coalescent-based approaches, maximum-likelihood and Bayesian methods, combined with other approaches such as permutations of the data, homogeneous versus heterogeneous models, and alternative data and taxon sets. Resulting topologies based on amino acids and nucleotides are both strongly supported but critically discordant, primarily in terms of the monophyly of Panopinae. Conflict was not resolved by controlling for compositional heterogeneity and saturation in third codon positions, which highlights the need for a better understanding of how different biases affect different data sources. In our study, results based on nucleotides were both more robust to alterations of the data and different analytical methods and more compatible with our current understanding of acrocerid morphology and patterns of host usage.

Phylogenomic Analysis of the Explosive Adaptive Radiation of the Espeletia Complex (Asteraceae) in the Tropical Andes.

The subtribe Espeletiinae (Asteraceae), endemic to the high-elevations in the Northern Andes, exhibits an exceptional diversity of species, growth-forms, and reproductive strategies. This complex of 140 species includes large trees, dichotomous trees, shrubs and the extraordinary giant caulescent rosettes, considered as a classic example of adaptation in tropical high-elevation ecosystems. The subtribe has also long been recognized as a prominent case of adaptive radiation, but the understanding of its evolution has been hampered by a lack of phylogenetic resolution. Herein, we produce the first fully resolved phylogeny of all morphological groups of Espeletiinae, using whole plastomes and about a million nuclear nucleotides obtained with an original de novo assembly procedure without reference genome, and analyzed with traditional and coalescent-based approaches that consider the possible impact of incomplete lineage sorting and hybridization on phylogenetic inference. We show that the diversification of Espeletiinae started from a rosette ancestor about 2.3 Ma, after the final uplift of the Northern Andes. This was followed by two independent radiations in the Colombian and Venezuelan Andes, with a few trans-cordilleran dispersal events among low-elevation tree lineages but none among high-elevation rosettes. We demonstrate complex scenarios of morphological change in Espeletiinae, usually implying the convergent evolution of growth-forms with frequent loss/gains of various traits. For instance, caulescent rosettes evolved independently in both countries, likely as convergent adaptations to life in tropical high-elevation habitats. Tree growth-forms evolved independently three times from the repeated colonization of lower elevations by high-elevation rosette ancestors. The rate of morphological diversification increased during the early phase of the radiation, after which it decreased steadily towards the present. On the other hand, the rate of species diversification in the best-sampled Venezuelan radiation was on average very high (3.1 spp/My), with significant rate variation among growth-forms (much higher in polycarpic caulescent rosettes). Our results point out a scenario where both adaptive morphological evolution and geographical isolation due to Pleistocene climatic oscillations triggered an exceptionally rapid radiation for a continental plant group.

Molecular evolution of key metabolic genes during transitions to C4 and CAM photosynthesis

Next-generation sequencing facilitates rapid production of well-sampled phylogenies built from very large genetic data sets, which can then be subsequently exploited to examine the molecular evolution of the genes themselves. We present an evolutionary analysis of 83 gene families (19 containing carbon-concentrating mechanism (CCM) genes, 64 containing non-CCM genes) in the portullugo clade (Caryophyllales), a diverse lineage of mostly arid-adapted plants that contains multiple evolutionary origins of all known photosynthesis types in land plants (C3 , C4 , CAM, C4 -CAM, and various intermediates). We inferred a phylogeny of 197 individuals from 167 taxa using coalescent-based approaches and individual gene family trees using maximum likelihood. Positive selection analyses were conducted on individual gene family trees with a mixed effects model of evolution (MEME). We devised new indices to compare levels of convergence and prevalence of particular residues between CCM and non-CCM genes and between species with different photosynthetic pathways. Contrary to expectations, there were no significant differences in the levels of positive selection detected in CCM versus non-CCM genes. However, we documented a significantly higher level of convergent amino acid substitutions in CCM genes, especially inC4 taxa. Our analyses reveal a new suite of amino acid residues putatively important for C4 and CAM function. We discuss both the advantages and challenges of using targeted enrichment sequence data for exploratory studies of molecular evolution.

Inverse European Latitudinal Cline at the timeless Locus of Drosophila melanogaster Reveals Selection on a Clock Gene: Population Genetics of ls-tim.

The spread of adaptive genetic variants in populations is a cornerstone of evolutionary theory but with relatively few biologically well-understood examples. Previous work on the ls-tim variant of timeless, which encodes the light-sensitive circadian regulator in Drosophila melanogaster, suggests that it may have originated in southeastern Italy. Flies characterized by the new allele show photoperiod-related phenotypes likely to be adaptive in seasonal environments. ls-tim may be spreading from its point of origin in Italy by directional selection, but there are alternative explanations for its observed clinal geographical distribution, including balancing selection and demography. From population analyses of ls-tim frequencies collected on the eastern side of the Iberian Peninsula, we show that ls-tim frequencies are inverted compared with those in Italy. This pattern is consistent with a scenario of directional selection rather than latitude-associated balancing selection. Neutrality tests further reveal the signature of directional selection at the ls-tim site, which is reduced a few kb pairs either side of ls-tim. A reanalysis of allele frequencies from a large number of microsatellite loci do not demonstrate any frequent ls-tim-like spatial patterns, so a general demographic effect or population expansion from southeastern Italy cannot readily explain current ls-tim frequencies. Finally, a revised estimate of the age of ls-tim allele using linkage disequilibrium and coalescent-based approaches reveals that it may be only 300 to 3000 years old, perhaps explaining why it has not yet gone to fixation. ls-tim thus provides a rare temporal snapshot of a new allele that has come under selection before it reaches equilibrium.

Coalescent-based delimitation outperforms distance-based methods for delineating less divergent species: the case of Kurixalus odontotarsus species group

Few empirical studies have compared coalescent-based methods to distance-based methods for delimitation of less divergent species. In this study, we used two coalescent-based (BFD and BPP) and two distance-based barcoding (ABGD and jMOTU) methods to delimit closely related species in the Kurixalus odontotarsus species group. Phylogenetic analyses revealed that the K. odontotarsus species group comprises 11 distinct maternal clades with strong support values. Based on the genetic and morphological evidences, we consider that species diversity in the K. odontotarsus species group was underestimated and the 11 clades represent 11 species, of which six are unnamed. The coalescent-based delimitations decisively supported the scenario of 11-species corresponding to the 11 clades. However, the distance-based ABGD only obtained 3–6 candidate species, which is not consistent with morphological evidence. These results indicate that BFD and BPP are more conservative than ABGD to false negatives (lumping). Method of fixed threshold (jMOTU) may obtain a resolution similar to that inferred by BFD and BPP, but it severely relies on subjective choice of the threshold and lacks statistical support. We consider that coalescent-based BFD and BPP approaches outperform distance-based methods for delineation of less divergent species.

Targeted Enrichment of Large Gene Families for Phylogenetic Inference: Phylogeny and Molecular Evolution of Photosynthesis Genes in the Portullugo Clade (Caryophyllales).

Hybrid enrichment is an increasingly popular approach for obtaining hundreds of loci for phylogenetic analysis across many taxa quickly and cheaply. The genes targeted for sequencing are typically single-copy loci, which facilitate a more straightforward sequence assembly and homology assignment process. However, this approach limits the inclusion of most genes of functional interest, which often belong to multi-gene families. Here, we demonstrate the feasibility of including large gene families in hybrid enrichment protocols for phylogeny reconstruction and subsequent analyses of molecular evolution, using a new set of bait sequences designed for the "portullugo" (Caryophyllales), a moderately sized lineage of flowering plants (~ 2200 species) that includes the cacti and harbors many evolutionary transitions to C$_{\mathrm{4}}$ and CAM photosynthesis. Including multi-gene families allowed us to simultaneously infer a robust phylogeny and construct a dense sampling of sequences for a major enzyme of C$_{\mathrm{4}}$ and CAM photosynthesis, which revealed the accumulation of adaptive amino acid substitutions associated with C$_{\mathrm{4}}$ and CAM origins in particular paralogs. Our final set of matrices for phylogenetic analyses included 75-218 loci across 74 taxa, with ~ 50% matrix completeness across data sets. Phylogenetic resolution was greatly improved across the tree, at both shallow and deep levels. Concatenation and coalescent-based approaches both resolve the sister lineage of the cacti with strong support: Anacampserotaceae $+$ Portulacaceae, two lineages of mostly diminutive succulent herbs of warm, arid regions. In spite of this congruence, BUCKy concordance analyses demonstrated strong and conflicting signals across gene trees. Our results add to the growing number of examples illustrating the complexity of phylogenetic signals in genomic-scale data.

Genetic diversity in the Andes: variation within and between the South American species of Oreobolus R. Br. (Cyperaceae)

This study examines genetic relationships among and within the South American species of Oreobolus that span the temperate and tropical Andes hotspots and represent a good case study to investigate diversification in the Paramo. A total of 197 individuals covering the distributional range of most of these species were sequenced for the nuclear ribosomal internal transcribed spacer (ITS) and 118 individuals for three chloroplast DNA regions (trnL-F, trnH-psbA and rpl32-trnL). Haplotype networks and measures of genetic diversity were calculated at different taxonomic and geographic levels. To test for possible geographic structure, a spatial analysis of molecular variance (SAMOVA) was undertaken and species relationships were recovered using a coalescent-based approach. Results indicate complex relationships among the five South American species of Oreobolus, which are likely to have been confounded by incomplete lineage sorting, though hybridization cannot be completely discarded as an influence on genetic patterns, particularly among the northern populations of O. obtusangulus and O. cleefii. We report a case of cryptic speciation in O. obtusangulus where northern and southern populations of morphologically similar individuals are genetically distinct in all analyses. At the population level, the genetic evidence is consistent with contraction and expansion of islands of Paramo vegetation during the climatic fluctuations of the Quaternary, highlighting the role of these processes in shaping modern diversity in that ecosystem.

Allopatric divergence, demographic history, and conservation implications of an endangered conifer Cupressus chengiana in the eastern Qinghai-Tibet Plateau

Isolation and demographic history are key factors that affect lineage divergence of tree species in topographic complex areas, such as the Qinghai-Tibet Plateau (QTP), yet few studies have evaluated these factors in a coalescent-based modeling framework. In the present study, we surveyed ten nuclear DNA sequence loci (nDNA) and six nuclear microsatellite loci (nSSR) for an endangered conifer, Cupressus chengiana, throughout its natural range in the eastern QTP. BARRIER analyses revealed a strong genetic barrier between Gansu and Sichuan populations of C. chengiana, and isolation with migration models detected limited gene flow between them, supporting the division of this species into two evolutionary significant units (ESUs). Two independent coalescent-based approaches suggest a Quaternary divergence between ESUs, their consensus age range ((0.09–) 0.59–1.53 (–2.71) Mya) largely overlaps with the time period when the largest glaciation occurred in the QTP. Both demographic inferences, IMa2 and DIYABC, suggest that both ESUs may have experienced a bottleneck or population contraction event during the late Quaternary. A documented massive recent anthropogenic habitat loss and fragmentation may have led to further decrease of the natural distribution of this conifer. We propose that the conservation and management of both natural stands and plantations of C. chengiana should be reconsidered in the light of our findings.

Evolutionary analysis of Chironius snakes unveils cryptic diversity and provides clues to diversification in the Neotropics

Recent hypotheses to explain tropical diversity involves the Neogene and Quaternary geoclimatic dynamics, but the absence of unambiguous data permitting the choice between alternative hypotheses makes a general theory for the origin of tropical biodiversity far to be achieved. The occurrence of Chironius snakes in well-defined biogeographical regions led us to adopt Chironius as a model to unveil patterns of vertebrate diversification in the Neotropics. Here, we used molecular markers and records on geographic distribution to investigate Chironius evolution and, subsequently, providing hints on diversification in the Neotropics. To avoid analyzing nominal species that do not constitute exclusive evolutionary lineages, we firstly conducted a species delimitation study prior to carrying out the species distribution modeling analysis. We generated 161 sequences of 12S, 16S, c-mos and rag2 for 15 species and 50 specimens, and included additional data from GenBank yielding a matrix of 137 terminals, and performed the following evolutionary analyses: inference of a concatenated gene tree, estimation of gene divergence times, inference of the coalescent-based phylogeny of Chironius, estimation of effective population sizes and modeling potential distribution of species across the last millennia. We tested for species boundaries within Chironius by implementing a coalescent-based Bayesian species delimitation approach. Our analyses supported the monophyly of Chironius, although our findings underscored cryptic candidate species in C. flavolineatus and C. exoletus. The inferred timetree suggested that Chironius snakes have evolved in the early Miocene (ca. 20.2Mya) and began to diversify from the late Miocene to the early Pliocene, values that are much older than previously reported. Following genetic divergence of virtually all extant Chironius species investigated, the effective sizes of the populations have expanded when compared to their MRCAs. The evolutionary and SDM data from C. brazili and C. diamantina provided additional evidence of the origin of species in the Neotropics. We argue that temperature, instead of altitude, has been the major driving factor in the evolution of both species, and thus we present a case for the consequences of global warming.

Multilocus coalescent species delimitation to evaluate traditionally defined morphotypes in Hydrangea sect. Asperae (Hydrangeaceae)

The number of species recognized in section Asperae of the flowering plant genus Hydrangea differs widely between subsequent revisions. This variation is largely centered around the H. aspera species complex, with numbers of recognized species varying from one to nearly a dozen. Despite indications of molecular variation in this complex, no sequence-based species delimitation methods have been employed to evaluate the primarily morphology-based species boundaries. In the present study, a multi-locus coalescent-based approach to species delimitation is employed in order to identify separate evolutionary lines within H. sect. Asperae, using four chloroplast and four nuclear molecular markers. Eight lineages were recovered within the focal group, of which five correspond with named morphotypes. The other three lineages illustrate types of conflict between molecular species delimitation and traditional morphology-based taxonomy. One molecular lineage comprises two named morphotypes, which possibly diverged recently enough to not have developed sufficient molecular divergence. A second conflict is found in H. strigosa. This morphotype is recovered as a separate lineage when occurring in geographic isolation, but when occurring in sympatry with two other morphotypes (H. aspera and H. robusta), the coalescent species delimitation lumps these taxa into a single putative species.

Carex sect. Rhynchocystis (Cyperaceae): a Miocene subtropical relict in the Western Palaearctic showing a dispersal‐derived Rand Flora pattern

Aim To evaluate how Cenozoic climate changes shaped the evolution and distribution of Carex section Rhynchocystis. Location Western Palaearctic and Afrotropical regions (Rand Flora pattern). Methods DNA regions ITS, ETS (nuclear), matK and rpl32-trnLUAG (plastid) were amplified for 86 samples of species from section Rhynchocystis. Phylogenetic and phylogeographical relationships were inferred using maximum parsimony, Bayesian inference and coalescent-based species tree approaches. Divergence times and ancestral areas were also inferred. Results Carex section Rhynchocystis is a clade that diversified during the middle Miocene in Europe. Most cladogenesis events date to the middle and late Miocene. The Afrotropical group seems to have originated from a colonization event from Europe that occurred in the late Miocene. Main conclusions Species of the section Rhynchocystis in the Western Palaearctic are Miocene relicts. Late Miocene-Pliocene aridification of the Mediterranean rather than the more commonly reported Pleistocene glaciations seems to have shaped the phylogeography of the group. Putative Miocene-Pliocene refugia were probably located in the Mediterranean peninsulas and islands, as well as in the eastern shores of the ancient Paratethys Sea. The colonization of Africa could have been facilitated by Miocene-Pliocene global cooling.

Species trees, temporal divergence and historical biogeography of coastal rove beetles (Coleoptera: Staphylinidae) reveal their early Miocene origin and show that most divergence events occurred in the early Pliocene along the Pacific coasts.

We investigated the species trees, temporal divergence and historical biogeography of Emplenota Casey and Triochara Bernhauer (ET clade hereafter), subgenera of the Aleochara coastal staphylinid beetles, to determine their place of origin. We explored a multilocus dataset to infer gene trees and species trees based on traditional concatenated and multispecies coalescent-based approaches using both model- (maximum likelihood and Bayesian) and parsimony-based methods. The multilocus dataset comprised DNA sequences from five nuclear genes (ArgK, CAD, EF1-α, wg and 28S) and three mitochondrial genes (COI, COII and 16S). *BEAST analysis according to a coalescent-based approach resolved the following phylogenetic relationships: ((A.trisulcata (A.nubis+A.zerchei)) (A.obscurella ((A.pacifica (A.curtidens + A.litoralis)) (A.hayamai (A.yamato (A.fucicola (A.segregata+A.puetzi))))))). Using a relaxed molecular clock, we reconstructed a time-calibrated phylogeny for this group. Furthermore, to account for the historical biogeography of the ET clade at the ages of major divergences, we investigated the ancestral area based on a time-calibrated phylogeny. Biogeographical analyses suggested that the ancestor of the ET clade was widely distributed along the eastern and western Palearctic and the western Nearctic coasts in the early Miocene. According to reconstruction of the ancestral area, one dispersal and three vicariance events were required, and the analyses indicated that vicariance events were important in shaping its current distribution patterns. Most of the divergence events occurred in the late Miocene and early Pliocene along the Pacific coasts, and the East Asian seacoasts harboured the most species-rich fauna, including eight of the 13 species in the ET clade.