A time-calibrated salamander phylogeny including 765 species and 503 genes

Recent time-calibrated amphibian phylogenies agree on the family-level relationships among extant salamanders but had disparate sampling regimes and inferred very different divergence times. For example, a recent phylogenomic study based on 220 nuclear loci had limited taxon sampling (41 species) and estimated relatively young divergence dates, whereas a more extensive supermatrix study based on 15 genes and 481 species estimated dates that were 22-45 million years older for major clades. Here, we combined phylogenomic and supermatrix approaches to estimate the largest salamander phylogeny to date based on molecular markers. Our matrix contained 765 salamander species and 503 genes (with 92.3% missing data overall). We included 284 more species than the previous largest salamander phylogeny (59% increase) and sampled approximately 93% of all currently described salamander species. Our dating analyses incorporated more than twice as many fossil calibration points within salamanders as previous studies. Maximum-likelihood estimates of tree topology yielded family-level relationships that were consistent with earlier studies. Nearly all species were placed in the expected genera, despite extensive missing data in many species. Bootstrap support was generally high across the tree but was poor in some clades where sampling of genes was limited (e.g., among some bolitoglossine salamanders). The dating analyses yielded age estimates for major clades that were generally intermediate between those from the previous phylogenomic and supermatrix analyses. We also provide a set of 200 time-calibrated trees for use in comparative analyses.

Revisiting phylogeny, systematics, and biogeography of a Pleistocene radiation

PREMISE

Pilosocereus (Cactaceae) is an important dry forest element in all subregions and transitional zones of the neotropics, with the highest diversity in eastern Brazil. The genus is subdivided into informal taxonomic groups; however, most of these are not supported by recent molecular phylogenetic inferences. This lack of confidence is probably due to the use of an insufficient number of loci and the complexity of cactus diversification. Here, we explored the species relationships in Pilosocereus in more detail, integrating multilocus phylogenetic approaches with the assessment of the ancestral range and the effect of geography on diversification shifts.

METHODS

We used 28 nuclear, plastid, and mitochondrial loci from 54 plant samples of 31 Pilosocereus species for phylogenetic analyses. We used concatenated and coalescent phylogenetic trees and Bayesian models to estimate the most likely ancestral range and diversification shifts.

RESULTS

All Pilosocereus species were clustered in the same branch, except P. bohlei. The phylogenetic relationships were more associated with the geographic distribution than taxonomic affinities among taxa. The genus began diversifying during the Plio-Pleistocene transition in the Caatinga domain and experienced an increased diversification rate during the Calabrian age.

CONCLUSIONS

We recovered a well-supported multispecies coalescent phylogeny. Our results refine the pattern of rapid diversification of Pilosocereus species across neotropical drylands during the Pleistocene and highlight the need for taxonomic rearrangements in the genus. We recovered a pulse of diversification during the Pleistocene that was likely driven by multiple dispersal and vicariance events within and among the Caatinga, Cerrado, and Atlantic Forest domains.

Target enrichment and extensive population sampling help untangle the recent, rapid radiation of Oenothera sect. Calylophus

Oenothera sect. Calylophus is a North American group of 13 recognized taxa in the evening primrose family (Onagraceae) with an evolutionary history that may include independent origins of bee pollination, edaphic endemism, and permanent translocation heterozygosity. Like other groups that radiated relatively recently and rapidly, taxon boundaries within Oenothera sect. Calylophus have remained challenging to circumscribe. In this study, we used target enrichment, flanking non-coding regions, gene tree/species tree methods, tests for gene flow modified for target-enrichment data, and morphometric analysis to reconstruct phylogenetic hypotheses, evaluate current taxon circumscriptions, and examine character evolution in Oenothera sect. Calylophus. Because sect. Calylophus comprises a clade with a relatively restricted geographic range, we were able to extensively sample across the range of geographic, edaphic and morphological diversity in the group. We found that the combination of exons and flanking non-coding regions led to improved support for species relationships. We reconstructed potential hybrid origins of some accessions and note that if processes such as hybridization are not taken into account, the number of inferred evolutionary transitions may be artificially inflated. We recovered strong evidence for multiple evolutionary origins of bee pollination from ancestral hawkmoth pollination, edaphic specialization on gypsum, and permanent translocation heterozygosity. This study applies newly emerging techniques alongside dense infraspecific sampling and morphological analyses to effectively reconstruct the recalcitrant history of a rapid radiation.

Slender salamanders (genus Batrachoseps) reveal Southern California to be a center for the diversification, persistence, and introduction of salamander lineages

Background

The southern California biodiversity hotspot has had a complex geological history, with both plate tectonic forces and sea level changes repeatedly reconfiguring the region, and likely driving both lineage splittings and extinctions. Here we investigate patterns of genetic divergence in two species of slender salamanders (Plethodontidae: Batrachoseps) in this region. The complex geological history in combination with several organismal traits led us to predict that these species harbor multiple ancient mitochondrial lineages endemic to southern California. These species belong to a clade characterized by fine-scale mitochondrial structure, which has been shown to track ancient splits. Both focal species, Batrachoseps major and B. nigriventris, are relatively widely distributed in southern California, and estimated to have persisted there across millions of years. Recently several extralimital populations of Batrachoseps were found in the San Joaquin Valley of California, a former desert area that has been extensively modified for agriculture. The origins of these populations are unknown, but based on morphology, they are hypothesized to result from human-mediated introductions of B. major.

Methods

We sequenced the mitochondrial gene cytochrome b from a geographically comprehensive sampling of the mitochondrial lineages of B. major and B. nigriventris that are endemic to southern California. We used phylogenetic analyses to characterize phylogeographic structure and identify mitochondrial contact zones. We also included the San Joaquin Valley samples to test whether they resulted from introductions. We used a bootstrap resampling approach to compare the strength of isolation-by-distance in both Batrachoseps species and four other salamander species with which they co-occur in southern California.

Results

The northern lineage of B. major harbors at least eight deeply differentiated, geographically cohesive mitochondrial subclades. We identify geographic contact between many of these mtDNA lineages and some biogeographic features that are concordant with lineage boundaries. Batrachoseps nigriventris also has multiple deeply differentiated clades within the region. Comparative analyses highlight the smaller spatial scales over which mitochondrial divergence accumulates in Batrachoseps relative to most other salamander species in southern California. The extralimital populations of Batrachoseps from the San Joaquin Valley are assigned to B. major and are shown to result from at least two independent introductions from different source populations. We also suggest that B. major on Catalina Island, where it is considered native, may be the result of an introduction. Some of the same traits that facilitate the build-up of deep phylogeographic structure in Batrachoseps likely also contribute to its propensity for introductions, and we anticipate that additional introduced populations will be discovered.

Integrative taxonomy by molecular species delimitation: multi-locus data corroborate a new species of Balkan Drusinae micro-endemics

BACKGROUND

Taxonomy offers precise species identification and delimitation and thus provides basic information for biological research, e.g. through assessment of species richness. The importance of molecular taxonomy, i.e., the identification and delimitation of taxa based on molecular markers, has increased in the past decade. Recently developed exploratory tools now allow estimating species-level diversity in multi-locus molecular datasets.

RESULTS

Here we use molecular species delimitation tools that either quantify differences in intra- and interspecific variability of loci, or divergence times within and between species, or perform coalescent species tree inference to estimate species-level entities in molecular genetic datasets. We benchmark results from these methods against 14 morphologically readily differentiable species of a well-defined subgroup of the diverse Drusinae subfamily (Trichoptera, Limnephilidae). Using a 3798 bp (6 loci) molecular data set we aim to corroborate a geographically isolated new species by integrating comparative morphological studies and molecular taxonomy.

CONCLUSIONS

Our results indicate that only multi-locus species delimitation provides taxonomically relevant information. The data further corroborate the new species Drusus zivici sp. nov. We provide differential diagnostic characters and describe the male, female and larva of this new species and discuss diversity patterns of Drusinae in the Balkans. We further discuss potential and significance of molecular species delimitation. Finally we argue that enhancing collaborative integrative taxonomy will accelerate assessment of global diversity and completion of reference libraries for applied fields, e.g., conservation and biomonitoring.

New species in the Sitalcina sura species group (Opiliones, Laniatores, Phalangodidae), with evidence for a biogeographic link between California desert canyons and Arizona sky islands

The western United States is home to numerous narrowly endemic harvestman taxa (Arachnida, Opiliones), including members of the genus Sitalcina Banks, 1911. Sitalcina is comprised of three species groups, including the monospecific Sitalcina californica and Sitalcina lobata groups, and the Sitalcina sura group with eight described species. All species in the Sitalcina sura group have very small geographic distributions, with group members distributed like disjunct "beads on a string" from Monterey south to southern California and southeast to the sky-island mountain ranges of southern Arizona. Here, molecular phylogenetic and species delimitation analyses were conducted for all described species in the Sitalcina sura group, plus several newly discovered populations. Species trees were reconstructed using multispecies coalescent methods implemented in *BEAST, and species delimitation was accomplished using Bayes Factor Delimitation (BFD). Based on quantitative species delimitation results supported by consideration of morphological characters, two new species (Sitalcina oasiensis sp. n., Sitalcina ubicki sp. n.) are described. We also provide a description of the previously unknown male of Sitalcina borregoensis Briggs, 1968. Molecular phylogenetic evidence strongly supports distinctive desert versus coastal clades, with desert canyon taxa from southern California more closely related to Arizona taxa than to geographically proximate California coastal taxa. We hypothesize that southern ancestry and plate tectonics have played a role in the diversification history of this animal lineage, similar to sclerophyllous plant taxa of the Madro-Tertiary Geoflora. Molecular clock analyses for the Sitalcina sura group are generally consistent with these hypotheses. We also propose that additional Sitalcina species await discovery in the desert canyons of southern California and northern Baja, and the mountains of northwestern mainland Mexico.

An integrative approach to phylogeography: investigating the effects of ancient seaways, climate, and historical geology on multi-locus phylogeographic boundaries of the Arboreal Salamander (Aneides lugubris)

BACKGROUND

Phylogeography is an important tool that can be used to reveal cryptic biodiversity and to better understand the processes that promote lineage diversification. We studied the phylogeographic history of the Arboreal Salamander (Aneides lugubris), a wide-ranging species endemic to the California floristic province. We used multi-locus data to reconstruct the evolutionary history of A. lugubris and to discover the geographic location of major genetic breaks within the species. We also used species distribution modeling and comparative phylogeography to better understand the environmental factors that have shaped the genetic history of A. lugubris.

RESULTS

We found six major mitochondrial clades in A. lugubris. Nuclear loci supported the existence of at least three genetically distinct groups, corresponding to populations north of the San Francisco Bay and in the Sierra Nevada, in the Santa Cruz Mountains, and in the central coast and southern California. All of the genetic breaks in mitochondrial and nuclear loci corresponded to regions where historical barriers to dispersal have been observed in other species. Geologic or water barriers likely were the most important factors restricting gene flow among clades. Climatic unsuitability during glacial maximum may have contributed to the isolation of the mitochondrial clades in the central coast and southern California. A projection of our species distribution model to a future scenario with a moderate amount of climate change suggests that most of the range of A. lugubris will remain climatically suitable, but climatic conditions in the Sierra Nevada and low elevation areas in Southern California are likely to deteriorate.

CONCLUSIONS

Aneides lugubris contains substantial cryptic genetic diversity as a result of historical isolation of populations. At least two (and perhaps three) evolutionarily significant units in A. lugubris merit protection; all six mitochondrial clades should be considered as management units within the species.

Multilocus phylogeography of the common midwife toad, Alytes obstetricans (Anura, Alytidae): Contrasting patterns of lineage diversification and genetic structure in the Iberian refugium

Recent investigations on the evolutionary history of the common midwife toad (Alytes obstetricans) revealed high levels of geographically structured genetic diversity but also a situation where delineation of major historical lineages and resolution of their relationships are much more complex than previously thought. We studied sequence variation in one mitochondrial and four nuclear genes throughout the entire distribution range of all recognized A. obstetricans subspecies to infer the evolutionary processes that shaped current patterns of genetic diversity and population subdivision. We found six divergent, geographically structured mtDNA haplogroups diagnosing population lineages, and varying levels of admixture in nuclear markers. Given the timeframe inferred for the splits between major lineages, the climatic and environmental changes that occurred during the Pleistocene seem to have shaped the diversification history of A. obstetricans. Survival of populations in allopatric refugia through the Ice Ages supports the generality of the "refugia-within-refugia" scenario for the Iberian Peninsula. However, lineages corresponding to subspecies A. o. almogavarii, A. o. pertinax, A. o. obstetricans, and A. o. boscai responded differently to Pleistocene climatic oscillations after diverging from a common ancestor. Alytes o. obstetricans expanded northward from a northern Iberian refugium through the western Pyrenees, leaving a signal of contrasting patterns of genetic diversity, with a single mtDNA haplotype north of the Pyrenees from SW France to Germany. Both A. o. pertinax and A. o. boscai are widespread and genetically diverse in Iberia, the latter comprising two divergent lineages with a long independent history. Finally, A. o. almogavarii is mostly restricted to the north-eastern corner of Iberia north of the Ebro river, with additional populations in a small region in south-eastern France. This taxon exhibits unparalleled levels of genetic diversity and little haplotype sharing with other lineages, suggesting a process of incipient speciation.