Nuclear–mitochondrial discordance and gene flow in a recent radiation of toads

Morphology and vocalization comparison of the Houston Toad and the Dwarf American Toad: implications for their historic range

Documenting changes in the distribution and abundance of a given taxon requires historical data. In the absence of long-term monitoring data collected throughout the range of a taxon, conservation biologists often rely on preserved museum specimens to determine the past or present, putative geographic distribution. Distributional data for the Houston Toad (Anaxyrus houstonensis) has consistently been confounded by similarities with a sympatric congener, the Dwarf American Toad (A. americanus charlesmithi), both in monitoring data derived from chorusing surveys, and in historical data via museum specimens. In this case, misidentification can have unintended impacts on conservation efforts, where the Houston Toad is federally endangered, and the Dwarf American Toad is of least concern. Previously published reports have compared these two taxon on the basis of their male advertisement call and morphological appearance, often with the goal of using these characters to substantiate their taxonomic status prior to the advent of DNA sequencing technology. However, numerous studies report findings that contradict one another, and no consensus on the true differences or similarities can be drawn. Here, we use contemporary recordings of wild populations of each taxon to test for quantifiable differences in male advertisement call. Additionally, we quantitatively examine a subset of vouchered museum specimens representing each taxon to test previously reported differentiating morphometric characters used to distinguish among other Bufonids of East-Central Texas, USA. Finally, we assemble and qualitatively evaluate a database of photographs representing catalogued museum vouchers for each taxon to determine if their previously documented historic ranges may be larger than are currently accepted. Our findings reveal quantifiable differences between two allopatric congeners with respect to their male advertisement call, whereas we found similarities among their detailed morphology. Additionally, we report on the existence of additional, historically overlooked, museum records for the Houston Toad in the context of its putative historic range, and discuss errors associated with the curation of these specimens whose identity and nomenclature have not been consistent through time. These results bookend decades of disagreement regarding the morphology, voice, and historic distribution of these taxa, and alert practitioners of conservation efforts for the Houston Toad to previously unreported locations of occurrence.

Using genome-wide data to ascertain taxonomic status and assess population genetic structure for Houston toads (Bufo [= Anaxyrus] houstonensis)

The Houston toad (Bufo [= Anaxyrus] houstonensis) is an endangered amphibian with a small geographic range. Land-use changes have primarily driven decline in B. houstonensis with population supplementation predominant among efforts to reduce its current extinction risk. However, there has been historic uncertainty regarding the evolutionary and conservation significance of B. houstonensis. To this end, we used 1170 genome-wide nuclear DNA markers to examine phylogenetic relationships between our focal taxon, representatives of the Nearctic B. americanus group, and B. nebulifer, a sympatric Middle American species. Phylogenetic analyses indicate B. houstonensis is a taxon that is distinct from B. americanus. We corroborated such genetic distinctiveness with an admixture analysis that provided support for recent reproductive isolation between B. americanus and B. houstonensis. However, ABBA-BABA tests for ancient admixture indicated historic gene flow between Nearctic species while no signal of historic gene flow was detected between Nearctic and Middle-American species. We used an admixture analysis to recognize four Management Units (MU) based on observed genetic differentiation within B. houstonensis and recommend captive propagation, population supplementation, and habitat restoration efforts specific to each MU. Our results re-affirm the evolutionary novelty of an endangered relict.

A New Assessment of Robust Capuchin Monkey (Sapajus) Evolutionary History Using Genome-Wide SNP Marker Data and a Bayesian Approach to Species Delimitation

Robust capuchin monkeys, Sapajus genus, are among the most phenotypically diverse and widespread groups of primates in South America, with one of the most confusing and often shifting taxonomies. We used a ddRADseq approach to generate genome-wide SNP markers for 171 individuals from all putative extant species of Sapajus to access their evolutionary history. Using maximum likelihood, multispecies coalescent phylogenetic inference, and a Bayes Factor method to test for alternative hypotheses of species delimitation, we inferred the phylogenetic history of the Sapajus radiation, evaluating the number of discrete species supported. Our results support the recognition of three species from the Atlantic Forest south of the São Francisco River, with these species being the first splits in the robust capuchin radiation. Our results were congruent in recovering the Pantanal and Amazonian Sapajus as structured into three monophyletic clades, though new morphological assessments are necessary, as the Amazonian clades do not agree with previous morphology-based taxonomic distributions. Phylogenetic reconstructions for Sapajus occurring in the Cerrado, Caatinga, and northeastern Atlantic Forest were less congruent with morphology-based phylogenetic reconstructions, as the bearded capuchin was recovered as a paraphyletic clade, with samples from the Caatinga biome being either a monophyletic clade or nested with the blond capuchin monkey.

Complex cycles of divergence and migration shape lineage structure in the common kingsnake species complex

Aim

The Nearctic is a complex patchwork of habitats and geologic features that form barriers to gene flow resulting in phylogeographic structure and speciation in many lineages. Habitats are rarely stable over geologic time, and the Nearctic has undergone major climatic changes in the past few million years. We use the common kingsnake species complex to study how climate, geography, and history influence lineage formation over a large, complex landscape.

Location

Nearctic/North America.

Taxon

Common kingsnake, Lampropeltis getula, species complex.

Methods

We analyzed genome-wide sequence data from 51 snakes spanning the majority of the species complex's range. We used population clustering, generalized dissimilarity modeling, and coalescent methods to identify the number of genetic clusters within the L. getula complex, infer the environmental correlates of genetic differentiation, and estimate models of divergence and gene flow among lineages.

Results

We identified three major lineages within the L. getula complex and further continuous spatial structure within lineages. The most important ecological correlates of genetic distance in the complex are related to aridity and precipitation, consistent with lineage breaks at the Great Plains/Desert ecotone and the Cochise Filter Barrier. Lineages are estimated to have undergone multiple rounds of isolation and secondary contact, with highly asymmetric migration occurring at present.

Main conclusions

Changing climates combined with a large and geologically complex landscape have resulted in a mosaic of discrete and spatially continuous genetic structure. Multiple rounds of isolation and secondary contact as climate fluctuated over the past ~4.4 My have likely driven the evolution of discrete lineages that maintain high levels of gene flow. Continuous structure is strongly shaped by aridity and precipitation, suggesting roles for major precipitation gradients in helping to maintain lineage identity in the face of gene flow when lineages are in geographic contact.

Phylogenomics, introgression, and demographic history of South American true toads (Rhinella)

The effects of genetic introgression on species boundaries and how they affect species' integrity and persistence over evolutionary time have received increased attention. The increasing availability of genomic data has revealed contrasting patterns of gene flow across genomic regions, which impose challenges to inferences of evolutionary relationships and of patterns of genetic admixture across lineages. By characterizing patterns of variation across thousands of genomic loci in a widespread complex of true toads (Rhinella), we assess the true extent of genetic introgression across species thought to hybridize to extreme degrees based on natural history observations and multi-locus analyses. Comprehensive geographic sampling of five large-ranged Neotropical taxa revealed multiple distinct evolutionary lineages that span large geographic areas and, at times, distinct biomes. The inferred major clades and genetic clusters largely correspond to currently recognized taxa; however, we also found evidence of cryptic diversity within taxa. While previous phylogenetic studies revealed extensive mito-nuclear discordance, our genetic clustering analyses uncovered several admixed individuals within major genetic groups. Accordingly, historical demographic analyses supported that the evolutionary history of these toads involved cross-taxon gene flow both at ancient and recent times. Lastly, ABBA-BABA tests revealed widespread allele sharing across species boundaries, a pattern that can be confidently attributed to genetic introgression as opposed to incomplete lineage sorting. These results confirm previous assertions that the evolutionary history of Rhinella was characterized by various levels of hybridization even across environmentally heterogeneous regions, posing exciting questions about what factors prevent complete fusion of diverging yet highly interdependent evolutionary lineages.

Cytonuclear discordance, reticulation and cryptic diversity in one of North America's most common frogs

Complicated phylogenetic histories benefit from diverse sources of inference. Pseudacris crucifer (spring peeper) spans most of eastern North America and comprises six mtDNA lineages that form multiple contact zones. The putative Miocene or early Pliocene origins of the oldest lineages within Pseudacris crucifer imply sufficient time for species-level divergence. To understand why this species appears unified while congeners have radiated, we analyze and compare male advertisement calls, mitochondrial, and nuclear markers and speak to the complex processes that have potentially influenced its contemporary patterns. We find extensive geographic and topological mitonuclear discordance, with three nuclear lineages containing 6 more-structured mtDNA lineages, and nuclear introgression at some contact zones. Male advertisement call differentiation is incongruent with the genetic structure as only one lineage appears differentiated. Occupying the Interior Highlands of the central United States, this Western lineage also has the most concordant mitochondrial and nuclear geographic patterns. Based on our findings we suggest that the antiquity of common ancestors was not as important as the maintenance of allopatry in the divergence in P. crucifer genetic lineages. We use multiple lines of evidence to generate hypotheses of isolation, reticulation, and discordance within this species and to expand our understanding of the early stages of speciation.

DNA barcoding reveals cryptic diversity in the underestimated genus Triplophysa (Cypriniformes: Cobitidae, Nemacheilinae) from the northeastern Qinghai-Tibet Plateau

Background

The northeastern part of the Qinghai-Tibet Plateau (QTP) presents a high number of plateau loach species. As one of the three major groups of fishes distributed on the QTP, plateau loach has high ecological value. However, the taxonomy and systematics of these fish are still controversial, and a large number of new species have been reported. The reason for this phenomenon is that the degree of morphological variation is low, the phylogenetic information provided by morphological and anatomical features used for species identification is relatively poor, and many cryptic species are observed. Based on the high-density sampling points from the biodiversity hotspots surveyed, this study aims to evaluate the biodiversity of plateau loach in the northeastern part of the QTP and reveal the hidden diversity by comparing morphological species with molecular operational taxonomic units (MOTUs).

Results

After careful identification and comparison of the morphology and DNA barcoding of 1630 specimens, 22 species were identified, with 20 considered valid local species and two identified as new species that had not been previously described. Based on the combination of morphological and molecular methods, a total of 24 native species were found, two of which were cryptic species: Triplophysa robusta sp1 and Triplophysa minxianensis sp1. Fourteen of the 24 species form clusters of barcodes that allow them to be reliably identified. The remaining cases involved 10 closely related species, including rapidly differentiated species and species that seemed to have experienced incomplete lineage sorting or showed introgressions.

Conclusions

The results highlight the need to combine traditional taxonomies with molecular methods to correctly identify species, especially closely related species, such as the plateau loach. This study provides a basis for protecting the biodiversity of plateau loach.

Population Analysis and Evolution of Saccharomyces cerevisiae Mitogenomes

The study of mitogenomes allows the unraveling of some paths of yeast evolution that are often not exposed when analyzing the nuclear genome. Although both nuclear and mitochondrial genomes are known to determine phenotypic diversity and fitness, no concordance has yet established between the two, mainly regarding strains' technological uses and/or geographical distribution. In the current work, we proposed a new method to align and analyze yeast mitogenomes, overcoming current difficulties that make it impossible to obtain comparable mitogenomes for a large number of isolates. To this end, 12,016 mitogenomes were considered, and we developed a novel approach consisting of the design of a reference sequence intended to be comparable between all mitogenomes. Subsequently, the population structure of 6646 Saccharomyces cerevisiae mitogenomes was assessed. Results revealed the existence of particular clusters associated with the technological use of the strains, in particular regarding clinical isolates, laboratory strains, and yeasts used for wine-associated activities. As far as we know, this is the first time that a positive concordance between nuclear and mitogenomes has been reported for S. cerevisiae, in terms of strains' technological applications. The results obtained highlighted the importance of including the mtDNA genome in evolutionary analysis, in order to clarify the origin and history of yeast species.

Quaternary land bridges have not been universal conduits of gene flow

Quaternary climate oscillations are a well-known driver of animal diversification, but their effects are most well studied in areas where glaciations lead to habitat fragmentation. In large areas of the planet, however, glaciations have had the opposite effect, but here their impacts are much less well understood. This is especially true in Southeast Asia, where cyclical changes in land distribution have generated enormous land expansions during glacial periods. In this study, we selected a panel of five songbird species complexes covering a range of ecological specificities to investigate the effects Quaternary land bridges have had on the connectivity of Southeast Asian forest biota. Specifically, we combined morphological and bioacoustic analysis with an arsenal of population genomic and modelling approaches applied to thousands of genome-wide DNA markers across a total of more than 100 individuals. Our analyses show that species dependent on forest understorey exhibit deep differentiation between Borneo and western Sundaland, with no evidence of gene flow during the land bridges accompanying the last 1-2 ice ages. In contrast, dispersive canopy species and habitat generalists have experienced more recent gene flow. Our results argue that there remains much cryptic species-level diversity to be discovered in Southeast Asia even in well-known animal groups such as birds, especially in nondispersive forest understorey inhabitants. We also demonstrate that Quaternary land bridges have not been equally suitable conduits of gene flow for all species complexes and that life history is a major factor in predicting relative population divergence time across Quaternary climate fluctuations.

Discordant evolution of mitochondrial and nuclear yeast genomes at population level

Background

Mitochondria are essential organelles partially regulated by their own genomes. The mitochondrial genome maintenance and inheritance differ from the nuclear genome, potentially uncoupling their evolutionary trajectories. Here, we analysed mitochondrial sequences obtained from the 1011 Saccharomyces cerevisiae strain collection and identified pronounced differences with their nuclear genome counterparts.

Results

In contrast with pre-whole genome duplication fungal species, S. cerevisiae mitochondrial genomes show higher genetic diversity compared to the nuclear genomes. Strikingly, mitochondrial genomes appear to be highly admixed, resulting in a complex interconnected phylogeny with a weak grouping of isolates, whereas interspecies introgressions are very rare. Complete genome assemblies revealed that structural rearrangements are nearly absent with rare inversions detected. We tracked intron variation in COX1 and COB to infer gain and loss events throughout the species evolutionary history. Mitochondrial genome copy number is connected with the nuclear genome and linearly scale up with ploidy. We observed rare cases of naturally occurring mitochondrial DNA loss, petite, with a subset of them that do not suffer the expected growth defect in fermentable rich media.

Conclusions

Overall, our results illustrate how differences in the biology of two genomes coexisting in the same cells can lead to discordant evolutionary histories.

Phylogeographic diversification and postglacial range dynamics shed light on the conservation of the kelp Saccharina japonica

Studies of postglacial range shifts could enhance our understanding of seaweed species' responses to climate change and hence facilitate the conservation of natural resources. However, the distribution dynamics and phylogeographic diversification of the commercially and ecologically important kelp Saccharina japonica in the Northwest Pacific (NWP) are still poorly surveyed. In this study, we analyzed the evolutionary history of S. japonica using two mitochondrial markers and 24 nuclear microsatellites. A STRUCTURE analysis revealed two partially isolated lineages: lineage H, which is scattered along the coast of Japan; and lineage P, which occurs along the west coast of the Japan Sea. Ecological niche modeling projections to the Last Glacial Maximum (LGM) revealed that the southern coasts of the Japan Sea and the Pacific side of the Oshima and Honshu Peninsulas provided the most suitable habitats for S. japonica, implying that these regions served as ancient refugia during the LGM. Ancient isolation in different refugia may explain the observed divergence between lineages P and H. An approximate Bayesian computation analysis indicated that the two lineages experienced post-LGM range expansion and that postglacial secondary contact occurred in Sakhalin. Model projections into the year 2,100 predicted that S. japonica will shift northwards and lose its genetic diversity center on the Oshima Peninsula in Hokkaido and Shimokita Peninsula in Honshu. The range shifts and evolutionary history of S. japonica improve our understanding of how climate change impacted the distribution range and diversity of this species and provide useful information for the conservation of natural resources under ongoing environmental change in the NWP.

Molecular systematics of the Triplophysa robusta (Cobitoidea) complex: Extensive gene flow in a depauperate lineage

Gene flow between populations assumed to be isolated frequently leads to incorrect inferences of evolutionary history. Understanding gene flow and its causes has long been a key topic in evolutionary biology. In this study, we explored the evolutionary history of the Triplophysa robusta complex, using a combination of multilocus analyses and coalescent simulation. Our multilocus approach detected conspicuous mitonuclear discordances in the T. robusta complex. Mitochondrial results showed reticular clades, whereas the nuclear results corresponded with the morphological data. Coalescent simulation indicated that gene flow was the source of these discordances. Molecular clock analysis combined with geological processes suggest that intense geological upheavals have shaped a complicated evolutionary history for the T. robusta complex since the late Miocene, causing extensive gene flow which has distorted the molecular systematics of the T. robusta complex. We suggest that frequent gene flow may restrict speciation in the T. robusta complex, leading to such a depauperate lineage. Based on this comprehensive understanding, we provide our proposals for taxonomic revision of the T. robusta complex.

Hidden diversity within the depauperate genera of the snake tribe Lampropeltini (Serpentes, Colubridae)

Accurate representation of lineage diversity through complete taxon sampling is crucial to understanding the evolution of biodiversity, particularly when using molecular phylogenetics to estimate evolutionary relationships. In this interest, taxonomic diversity is often used as a proxy for lineage diversity even though the two concepts are not synonymous. We explore this within the snake tribe Lampropeltini which includes some of the most conspicuous and heavily studied snakes in North America. Both the taxonomy and hypothesized relationships within this tribe have been in flux. The number of species has increased from 23 to 51 over the last thirty years, predominately within three of the nine genera (Lampropeltis, Pantherophis, Pituophis). The remaining six depauperate genera (Arizona, Bogertophis, Cemophora, Pseudelaphe, Rhinocheilus, and Senticolis) have been poorly represented in phylogenetic studies. To estimate evolutionary relationships and determine if the dichotomy in depauperate and speciose genera within Lampropeltini is a function of taxon sampling or truly represents the lineage diversity, we estimated the phylogeny of this group using nuclear and mitochondrial loci in a concatenated and coalescent framework with the largest sampling of the six depauperate genera to date. In addition, we estimated the divergence dates among the genera to assess whether the instability of Lampropeltini phylogenetic relationships is due to an adaptive radiation. While some nodes still remain unresolved, the generic-level relationships we recovered agree with those of a recent next-generation study that used a much larger set of loci for fewer individuals. We also tested two putative species, Arizona pacata and Pseudelaphe phaescens, for the first time phylogenetically and find evidence that they are distinct lineages. Overall, we find that the taxonomic and genetic diversity are not correlated in Lampropeltini and that representing putative diversity in phylogenies will lead to a better estimate of evolutionary histories, especially in groups with complex radiations.

A novel mitochondrial genome of Arborophila and new insight into Arborophila evolutionary history

The lineage of the Bar-backed Partridge (Arborophila brunneopectus) was investigated to determine the phylogenetic relationships within Arborophila as the species is centrally distributed within an area covered by the distributions of 22 South-east Asian hill partridge species. The complete mitochondrial genome (mitogenome) of A. brunneopectus was determined and compared with four other hill partridge species mitogenomes. NADH subunit genes are radical in hill partridge mitogenomes and contain the most potential positive selective sites around where variable sites are abundant. Together with 44 other mitogenomes of closely related species, we reconstructed highly resolved phylogenetic trees using maximum likelihood (ML) and Bayesian inference (BI) analyses and calculated the divergence and dispersal history of Arborophila using combined datasets composed of their 13-protein coding sequences. Arborophila is reportedly be the oldest group in Phasianidae whose ancestors probably originated in Asia. A. rufipectus shares a closer relationship with A. ardens and A. brunneopectus compared to A. gingica and A. rufogularis, and such relationships were supported and profiled by NADH dehydrogenase subunit 5 (ND5). The intragenus divergence of all five Arborophila species occurred in the Miocene (16.84~5.69 Mya) when there were periods of climate cooling. We propose that these cooling events in the Miocene forced hill partridges from higher to lower altitudes, which led to geographic isolation and speciation. We demonstrated that the apparently deleterious +1 frameshift mutation in NADH dehydrogenase subunit 3 (ND3) found in all Arborophila is an ancient trait that has been eliminated in some younger lineages, such as Passeriformes. It is unclear of the biological advantages of this elimination for the relevant taxa and this requires further investigation.

Phylogenetic relationships of Rhizoctonia fungi within the Cantharellales

Phylogenetic relationships of Rhizoctonia fungi within the order Cantharellales were studied using sequence data from portions of the ribosomal DNA cluster regions ITS-LSU, rpb2, tef1, and atp6 for 50 taxa, and public sequence data from the rpb2 locus for 165 taxa. Data sets were analysed individually and combined using Maximum Parsimony, Maximum Likelihood, and Bayesian Phylogenetic Inference methods. All analyses supported the monophyly of the family Ceratobasidiaceae, which comprises the genera Ceratobasidium and Thanatephorus. Multi-locus analysis revealed 10 well-supported monophyletic groups that were consistent with previous separation into anastomosis groups based on hyphal fusion criteria. This analysis coupled with analyses of a larger sample of 165 rpb2 sequences of fungi in the Cantharellales supported a sister relationship between the Botryobasidiaceae and Ceratobasidiaceae and a sister relationship of the Tulasnellaceae with the rest of the Cantharellales. The inclusion of additional sequence data did not clarify incongruences observed in previous studies of Rhizoctonia fungi in the Cantharellales based on analyses of a single or multiple genes. The diversity of ecological and morphological characters associated with these fungi requires further investigation on character evolution for re-evaluating homologous and homoplasious characters.

Phylogenetic relationships of toads of the Rhinella granulosa group (Anura: Bufonidae): a molecular perspective with comments on hybridization and introgression

The Rhinella granulosa group consists of 13 species of toads distributed throughout open areas of South America and Panama. In this paper we perform a phylogenetic analysis considering all but one species of the group, employing five nuclear and four mitochondrial genes, for up to 7910 bp per specimen. Separate phylogenetic analyses under direct optimization (DO) of nuclear and mitochondrial sequences recovered the R. granulosa group as monophyletic and revealed topological incongruence that can be explained mainly by multiple events of hybridization and introgression, both mitochondrial and nuclear. The DO combined analysis, after the exclusion of putatively introgressed or heterozygous genomes, resulted in a phylogenetic hypothesis for the R. granulosa group in which most of the species are recovered as monophyletic, but with interspecific relationships poorly supported. The optimization of morphological (adult and larval), chromosomal, and behavioural characters resulted in 12 putative phenotypic synapomorphies for this species group and some other synapomorphies for internal clades. Our results indicate the need for additional population genetic studies on R. dorbignyi and R. fernandezae to corroborate the taxonomic status of both taxa. Finally, we discuss biological and genetic characteristics of Bufonidae, as possible explanations for the common occurrence of hybridization and introgression observed in some lineages of this family.

Semi-permeable species boundaries in Iberian barbels (Barbus and Luciobarbus, Cyprinidae)

BACKGROUND

The evolution of species boundaries and the relative impact of selection and gene flow on genomic divergence are best studied in populations and species pairs exhibiting various levels of divergence along the speciation continuum. We studied species boundaries in Iberian barbels, Barbus and Luciobarbus, a system of populations and species spanning a wide degree of genetic relatedness, as well as geographic distribution and range overlap. We jointly analyze multiple types of molecular markers and morphological traits to gain a comprehensive perspective on the nature of species boundaries in these cyprinid fishes.

RESULTS

Intraspecific molecular and morphological differentiation is visible among many populations. Genomes of all sympatric species studied are porous to gene flow, even if they are not sister species. Compared to their allopatric counterparts, sympatric representatives of different species share alleles and show an increase in all measures of nucleotide polymorphism (S, Hd, K, π and θ). High molecular diversity is particularly striking in L. steindachneri from the Tejo and Guadiana rivers, which co-varies with other sympatric species. Interestingly, different nuclear markers introgress across species boundaries at various levels, with distinct impacts on population trees. As such, some loci exhibit limited introgression and population trees resemble the presumed species tree, while alleles at other loci introgress more freely and population trees reflect geographic affinities and interspecific gene flow. Additionally, extent of introgression decreases with increasing genetic divergence in hybridizing species pairs.

CONCLUSIONS

We show that reproductive isolation in Iberian Barbus and Luciobarbus is not complete and species boundaries are semi-permeable to (some) gene flow, as different species (including non-sister) are exchanging genes in areas of sympatry. Our results support a speciation-with-gene-flow scenario with heterogeneous barriers to gene flow across the genome, strengthening with genetic divergence. This is consistent with observations coming from other systems and supports the notion that speciation is not instantaneous but a gradual process, during which different species are still able to exchange some genes, while selection prevents gene flow at other loci. We also provide evidence for a hybrid origin of a barbel ecotype, L. steindachneri, suggesting that ecology plays a key role in species coexistence and hybridization in Iberian barbels. This ecotype with intermediate, yet variable, molecular, morphological, trophic and ecological characteristics is the local product of introgressive hybridization of L. comizo with up to three different species (with L. bocagei in the Tejo, with L. microcephalus and L. sclateri in the Guadiana). In spite of the homogenizing effects of ongoing gene flow, species can still be discriminated using a combination of morphological and molecular markers. Iberian barbels are thus an ideal system for the study of species boundaries, since they span a wide range of genetic divergences, with diverse ecologies and degrees of sympatry.

Cryptic genetic diversity and complex phylogeography of the boreal North American scorpion, Paruroctonus boreus (Vaejovidae)

Diverse studies in western North America have revealed the role of topography for dynamically shaping genetic diversity within species though vicariance, dispersal and range expansion. We examined patterns of phylogeographical diversity in the widespread but poorly studied North American vaejovid scorpion, Paruroctonus boreus Girard 1854. We used mitochondrial sequence data and parsimony, likelihood, and Bayesian inference to reconstruct phylogenetic relationships across the distributional range of P. boreus, focusing on intermontane western North America. Additionally, we developed a species distribution model to predict its present and historical distributions during the Last Glacial Maximum and the Last Interglacial Maximum. Our results documented complex phylogeographic relationships within P. boreus, with multiple, well-supported crown clades that are either geographically-circumscribed or widespread and separated by short, poorly supported internodes. We also observed subtle variation in predicted habitat suitability, especially at the northern, eastern and southern edges of the predicted distributional range under past climatic conditions. The complex phylogenetic relationships of P. boreus suggests that historical isolation and expansion of populations may have occurred. Variation in the predicted distributional range over time may implicate past climatic fluctuations in generating the patterns of genetic diversity observed in P. boreus. These findings highlight both the potential for cryptic biodiversity in widespread North American scorpion species and the importance of phylogeographical studies for understanding the factors responsible for generating the biodiversity of western North America.

Delimiting genetic units in Neotropical toads under incomplete lineage sorting and hybridization

Background

Delimiting genetic units is useful to enhance taxonomic discovery and is often the first step toward understanding evolutionary mechanisms generating diversification. The six species within the Rhinella crucifer group of toads were defined under morphological criteria alone. Previous data suggest limited correspondence of these species to mitochondrial lineages, and morphological intergradation at transitions between forms suggests hybridization. Here we extensively sampled populations throughout the geographic distribution of the group and analyzed mitochondrial and nuclear sequence data to delimit genetic units using tree–based and allele frequency–based approaches.

Results

These approaches yielded complementary results, with allele frequency-based methods performing unexpectedly well given the limited number of loci examined. Both mitochondrial and nuclear markers supported a genetic structure of five units within the group, with three of the inferred units distributed within its main range, while two other units occur in separate isolates. The inferred units are mostly discordant with currently described forms: unequivocal association exists for only two of the six species in the group. Genetic evidence for hybridization exists for two pairs of units, with clear cyto–nuclear allele mixing observed in one case.

Conclusions

Our results confirmed that current taxonomy does not represent evolutionary units in the Rhinella crucifer group. Correspondence between genetically distinguishable units and the currently recognized species is only possible for Rhinella henseli and R. inopina. The recognition of other species relies on the reassessment of the geographic range of R. crucifer, the examination of the type series of R. ornata for hybrids, and on the use of additional markers to verify the genetic distinctiveness of R. abei. We state that R. pombali should not remain a valid species since its description appears to be based on hybrids, and that the name R. pombali should be considered a synonym of both R. crucifer and R. ornata. The fifth inferred but undescribed genetic unit may represent a new species. Our results underscore the potential of the R. crucifer species group to contribute to a better understanding of diversification processes and hybridization patterns in the Neotropics, and provide the basis for future evolutionary and taxonomic studies.