Phylogeographic analysis and environmental niche modeling of the plain-bellied watersnake (Nerodia erythrogaster) reveals low levels of genetic and ecological differentiation

The Phylogenomic and Biogeographic History of the Gartersnakes, Watersnakes, and Allies (Natricidae: Thamnophiini)

North American Thamnophiini (gartersnakes, watersnakes, brownsnakes, and swampsnakes) are an ecologically and phenotypically diverse temperate clade of snakes representing 61 species across 10 genera. In this study, we estimate phylogenetic trees using ∼3,700 ultraconserved elements (UCEs) for 76 specimens representing 75% of all Thamnophiini species. We infer phylogenies using multispecies coalescent methods and time calibrate them using the fossil record. We also conducted ancestral area estimation to identify how major biogeographic boundaries in North America affect broadscale diversification in the group. While most nodes exhibited strong statistical support, analysis of concordant data across gene trees reveals substantial heterogeneity. Ancestral area estimation demonstrated that the genus Thamnophis was the only taxon in this subfamily to cross the Western Continental Divide, even as other taxa dispersed southward toward the tropics. Additionally, levels of gene tree discordance are overall higher in transition zones between bioregions, including the Rocky Mountains. Therefore, the Western Continental Divide may be a significant transition zone structuring the diversification of Thamnophiini during the Neogene and Pleistocene. Here we show that despite high levels of discordance across gene trees, we were able to infer a highly resolved and well-supported phylogeny for Thamnophiini, which allows us to understand broadscale patterns of diversity and biogeography.

Existence of Bov-B LINE Retrotransposons in Snake Lineages Reveals Recent Multiple Horizontal Gene Transfers with Copy Number Variation

Transposable elements (TEs) are dynamic elements present in all eukaryotic genomes. They can “jump” and amplify within the genome and promote segmental genome rearrangements on both autosomes and sex chromosomes by disruption of gene structures. The Bovine-B long interspersed nuclear element (Bov-B LINE) is among the most abundant TE-retrotransposon families in vertebrates due to horizontal transfer (HT) among vertebrate lineages. Recent studies have shown multiple HTs or the presence of diverse Bov-B LINE groups in the snake lineage. It is hypothesized that Bov-B LINEs are highly dynamic and that the diversity reflects multiple HTs in snake lineages. Partial sequences of Bov-B LINE from 23 snake species were characterized. Phylogenetic analysis resolved at least two Bov-B LINE groups that might correspond to henophidian and caenophidian snakes; however, the tree topology differed from that based on functional nuclear and mitochondrial gene sequences. Several Bov-B LINEs of snakes showed greater than 80% similarity to sequences obtained from insects, whereas the two Bov-B LINE groups as well as sequences from the same snake species classified in different Bov-B LINE groups showed sequence similarities of less than 80%. Calculation of estimated divergence time and pairwise divergence between all individual Bov-B LINE copies suggest invasion times ranging from 79.19 to 98.8 million years ago in snakes. Accumulation of elements in a lineage-specific fashion ranged from 9 × 10⁻⁶% to 5.63 × 10⁻²% per genome. The genomic proportion of Bov-B LINEs varied among snake species but was not directly associated with genome size or invasion time. No differentiation in Bov-B LINE copy number between males and females was observed in any of the snake species examined. Incongruence in tree topology between Bov-B LINEs and other snake phylogenies may reflect past HT events. Sequence divergence of Bov-B LINEs between copies suggests that recent multiple HTs occurred within the same evolutionary timeframe in the snake lineage. The proportion of Bov-B LINEs varies among species, reflecting species specificity in TE invasion. The rapid speciation of snakes, coinciding with Bov-B LINE invasion in snake genomes, leads us to better understand the effect of Bov-B LINEs on snake genome evolution.

Regional patterns of genetic structure and environmental differentiation in willow populations (Salix humboldtiana Willd.) from Central Mexico

AIM

To infer the geological and climatic factors that have shaped the genetic diversity and structure of a willow species (Salix humboldtiana) in three basins of Central Mexico.

LOCATION

Central Mexico.

METHODS

We collected samples from 11 populations across two hydrological basins (Balsas and Lerma) and one population from another basin (Ameca) within the Mexican Central Plateau (MCP). Individuals were analyzed using sequences of two chloroplast DNA (cpDNA) regions and eight nuclear simple sequence repeats (nSSR). Population genetic diversity and structure were determined from these data. To evaluate whether genetic structure was associated with ecological niche differentiation, we determined whether there is niche equivalence, overlap, or divergence between the Balsas and Lerma basins. Also, we evaluated the relative contributions of geographic distribution and climatic variation on population genetic structuring through redundancy analysis (RDA) and partial RDA.

RESULTS

Both cpDNA and nSSRs data indicated the presence of three highly differentiated genetic groups, mostly geographically congruent with the three main hydrological basins. According to nSSRs, the three genetic groups can be further subdivided into eight subgroups corresponding to different rivers within the main basins. The niche equivalency test showed that the niches of the species in the Balsas and Lerma basins are significantly nonequivalent. The RDA indicated a significant association of genetic variation among populations with climate variables (particularly those related to the precipitation regime), while controlling for geographic distribution.

MAIN CONCLUSIONS

The genetic structure of S. humboldtiana is strongly associated with the historical and current geological configuration of the basins and the rivers within basins. The observed hierarchical genetic differentiation can be due to gene flow limitation resulting from physical barriers to the dispersal of S. humboldtiana, but also to some degree of isolation by environment, as suggested by the significant association between genetic variation among populations and precipitation regime.

Geographic patterns of phenotypic diversity in incipient species of North American blister beetles (Coleoptera: Meloidae) are not determined by species niches, but driven by demography along the speciation process

The Epicauta stigmata complex is a group of blister beetles composed of three parapatric or sympatric species that occur in central Mexico to southern USA: E. stigmata, E. uniforma and E. melanochroa. These species are morphologically very similar, and are mainly distinguished by body colour differences. Here we assessed whether phenotypic divergence in coloration patterns define evolutionary units within the complex. We studied the phylogenetic relationships, demographic history and concordances between morphological and ecological traits in the group. The complex apparently had a demographic history of recent population expansion during the last glaciation period 75000 to 9500 years ago. The three species show no reciprocal monophyly, and thus their allospecificity was not confirmed. The current distribution of haplotypes and the genetic divergences in these taxa can be explained by either recent mitochondrial introgression events caused by hybridisation or by incomplete lineage sorting. Colour pattern differences in the complex are not likely a product of local selection acting over a common genetic background. We suggest that phenotypic divergence in colour patterns during an incipient speciation process might be seen as an enhancing factor of cohesion within each of the three evolutionary units.

Assessment of snake DNA barcodes based on mitochondrial COI and Cytb genes revealed multiple putative cryptic species in Thailand

DNA barcodes of mitochondrial cytochrome c oxidase I (COI), cytochrome b (Cytb) genes, and their combined data sets were constructed from 35 snake species in Thailand. No barcoding gap was detected in either of the two genes from the observed intra- and interspecific sequence divergences. Intra- and interspecific sequence divergences of the COI gene differed 14 times, with barcode cut-off scores ranging over 2%-4% for threshold values differentiated among most of the different species; the Cytb gene differed 6 times with cut-off scores ranging over 2%-6%. Thirty-five specific nucleotide mutations were also found at interspecific level in the COI gene, identifying 18 snake species, but no specific nucleotide mutation was observed for Cytb in any single species. This suggests that COI barcoding was a better marker than Cytb. Phylogenetic clustering analysis indicated that most species were represented by monophyletic clusters, suggesting that these snake species could be clearly differentiated using COI barcodes. However, the two-marker combination of both COI and Cytb was more effective, differentiating snake species by over 2%-4%, and reducing species numbers in the overlap value between intra- and interspecific divergences. Three species delimitation algorithms (general mixed Yule-coalescent, automatic barcoding gap detection, and statistical parsimony network analysis) were extensively applied to a wide range of snakes based on both barcodes. This revealed cryptic diversity for eleven snake species in Thailand. In addition, eleven accessions from the database previously grouped under the same species were represented at different species level, suggesting either high genetic diversity, or the misidentification of these sequences in the database as a consequence of cryptic species.

Phylogeography and evolution of the Red Salamander (Pseudotriton ruber)

Phylogeographic studies frequently result in the elevation of subspecific taxa to species given monophyly, or the synonymy of subspecies that are not monophyletic. However, given limited or incongruent datasets, retention of subspecies can be useful to describe hypothesized incipient species or to illustrate interesting biological phenomena driving morphological diversity. Four subspecific taxa have been used to describe largely allopatric geographic variation within the species Pseudotriton ruber, a plethodontid salamander occupying stream and spring habitats across eastern North America: P. r. vioscai occurs in lowland Coastal Plain habitats, while P. r. ruber, P. r. nitidus, and P. r. schencki occupy upland regions in and around the Appalachian Mountains. Pseudotriton ruber co-occurs through its distribution with the aposematic newt Notophthalmus viridescens, and both species are hypothesized to be part of a Müllerian mimicry complex. In this study, we sequenced regions of two mitochondrial (cytochrome b, NADH dehydrogenase subunit 2) and one single copy nuclear protein-coding gene (pro-opiomelanocortin) from individuals sampled across much of the distribution of P. ruber and then used maximum-likelihood and Bayesian phylogenetic inference to test the monophyly of subspecies, reconstruct biogeographic history, and make inferences about morphological evolution. Phylogeographic hypotheses from mitochondrial and nuclear datasets described structure among populations of P. ruber which separated Coastal Plain and upland Appalachian populations, but subspecies were not monophyletic. Biogeographic reconstruction estimated the ancestor of all populations to have occupied and initially diverged in the Coastal Plain during the Pliocene (∼3.6mya), before one lineage subsequently invaded upland areas of Appalachia. Bold bright coloration of high elevation subspecies P. r. nitidus and P. r. schencki appears to have evolved twice. We hypothesize that the Müllerian mimicry complex with N. viridescens and P. ruber may provide a selective mechanism driving the co-evolution of striking bright and dull morphological variation among populations of both species. While P. ruber subspecies were not consistent with our criteria for diagnosing species (monophyly) and therefore could not be elevated to species, we advocate for the retention of subspecies because they describe hypotheses about an incipient species (P. r. vioscai) and how Müllerian mimicry may shape morphological diversity of species.

Bayesian tests of topology hypotheses with an example from diving beetles

We review Bayesian approaches to model testing in general and to the assessment of topological hypotheses in particular. We show that the standard way of setting up Bayes factor tests of the monophyly of a group, or the placement of a sample sequence in a known reference tree, can be misleading. The reason for this is related to the well-known dependency of Bayes factors on model-specific priors. Specifically, when testing tree hypotheses it is important that each hypothesis is associated with an appropriate tree space in the prior. This can be achieved by using appropriately constrained searches or by filtering trees in the posterior sample, but in a more elaborate way than typically implemented. If it is difficult to find the appropriate tree sets to be contrasted, then the posterior model odds may be more informative than the Bayes factor. We illustrate the recommended techniques using an empirical test case addressing the issue of whether two genera of diving beetles (Coleoptera: Dytiscidae), Suphrodytes and Hydroporus, should be synonymized. Our refined Bayes factor tests, in contrast to standard analyses, show that there is strong support for Suphrodytes nesting inside Hydroporus, and the genera are therefore synonymized.

Reevaluation of a classic phylogeographic barrier: new techniques reveal the influence of microgeographic climate variation on population divergence

We evaluated the mtDNA divergence and relationships within Geomys pinetis to assess the status of formerly recognized Geomys taxa. Additionally, we integrated new hypothesis-based tests in ecological niche models (ENM) to provide greater insight into causes for divergence and potential barriers to gene flow in Southeastern United States (Alabama, Florida, and Georgia). Our DNA sequence dataset confirmed and strongly supported two distinct lineages within G. pinetis occurring east and west of the ARD. Divergence date estimates showed that eastern and western lineages diverged about 1.37 Ma (1.9 Ma-830 ka). Predicted distributions from ENMs were consistent with molecular data and defined each population east and west of the ARD with little overlap. Niche identity and background similarity tests were statistically significant suggesting that ENMs from eastern and western lineages are not identical or more similar than expected based on random localities drawn from the environmental background. ENMs also support the hypothesis that the ARD represents a ribbon of unsuitable climate between more suitable areas where these populations are distributed. The estimated age of divergence between eastern and western lineages of G. pinetis suggests that the divergence was driven by climatic conditions during Pleistocene glacial-interglacial cycles. The ARD at the contact zone of eastern and western lineages of G. pinetis forms a significant barrier promoting microgeographic isolation that helps maintain ecological and genetic divergence.

DNA barcoding amphibians and reptiles

Only a few major research programs are currently targeting COI barcoding of amphibians and reptiles (including chelonians and crocodiles), two major groups of tetrapods. Amphibian and reptile species are typically old, strongly divergent, and contain deep conspecific lineages which might lead to problems in species assignment with incomplete reference databases. As far as known, there is no single pair of COI primers that will guarantee a sufficient rate of success across all amphibian and reptile taxa, or within major subclades of amphibians and reptiles, which means that the PCR amplification strategy needs to be adjusted depending on the specific research question. In general, many more amphibian and reptile taxa have been sequenced for 16S rDNA, which for some purposes may be a suitable complementary marker, at least until a more comprehensive COI reference database becomes available. DNA barcoding has successfully been used to identify amphibian larval stages (tadpoles) in species-rich tropical assemblages. Tissue sampling, DNA extraction, and amplification of COI is straightforward in amphibians and reptiles. Single primer pairs are likely to have a failure rate between 5 and 50% if taxa of a wide taxonomic range are targeted; in such cases the use of primer cocktails or subsequent hierarchical usage of different primer pairs is necessary. If the target group is taxonomically limited, many studies have followed a strategy of designing specific primers which then allow an easy and reliable amplification of all samples.