Mitochondrial DNA-Based phylogeography of North American rubber boas, Charina bottae (Serpentes: Boidae)

Reference genome of the rubber boa, Charina bottae (Serpentes: Boidae)

The rubber boa, Charina bottae is a semi-fossorial, cold-temperature adapted snake that ranges across the wetter and cooler ecoregions of the California Floristic Province. The rubber boa is 1 of 2 species in the family Boidae native to California and currently has 2 recognized subspecies, the Northern rubber boa C. bottae bottae and the Southern rubber boa C. bottae umbratica. Recent genomic work on C. bottae indicates that these 2 subspecies are collectively composed of 4 divergent lineages that separated during the late Miocene. Analysis of habitat suitability indicates that C. bottae umbratica montane sky-island populations from southern California will lose the majority of their habit over the next 70 yr, and is listed as Threatened under the California Endangered Species Act. Here, we report a new, chromosome-level assembly of C. bottae bottae as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genome strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises 289 scaffolds covering 1,804,944,895 bp, has a contig N50 of 37.3 Mb, a scaffold N50 of 97 Mb, and BUSCO completeness score of 96.3%, and represents the first reference genome for the Boidae snake family. This genome will enable studies of genetic differentiation and connectivity among C. bottae bottae and C. bottae umbratica populations across California and help manage locally endemic lineages as they confront challenges from human-induced climate warming, droughts, and wildfires across California.

Genomic data reveal local endemism in Southern California Rubber Boas (Serpentes: Boidae, Charina) and the critical need for enhanced conservation actions

The mountains of southern California represent unique, isolated ecosystems that support distinct high-elevation habitats found nowhere else in the area. Analyses of several moisture-dependent species across these sky-islands indicate they exist as locally endemic lineages that occur across these fragmented mountains ranges. The Rubber Boa is a semi-fossorial snake species that is widely distributed in the cooler and more moist ecoregions regions of western North America, including isolated populations across southern California mountain ranges. We developed a genomic and ecological dataset to examine genetic diversity within Rubber Boas and to determine if the endemic Southern Rubber Boa represents a distinct lineage. We quantified current and future habitat suitability under a range of climate change scenarios, and discuss the possible environmental threats facing these unique montane isolates. Our results support four major lineages within Rubber Boas, with genetic breaks that are consistent with biogeographic boundaries observed in other co-distributed, cool-temperature, moisture adapted species. Our data support previous studies that the Southern Rubber Boa is an independent evolutionary unit and now includes multiple locally endemic sky-island populations, restricted to isolated mountain tops and ranges across southern California. Analyses of future habitat suitability indicate that many of these sky-island populations will lose most of their suitable habitat over the next 70 years given predicted increases in drought, rising temperatures, and wildfires. Collectively these data emphasize the critical conservation needs of these montane ecosystems in southern California under current and projected climate change conditions.

Phylogeography of an endemic California silkmoth genus suggests the importance of an unheralded central California province in generating regional endemic biodiversity

The California Floristic province is a biodiversity hotspot. Understanding the phylogeographic patterns that exist in this diverse region is essential to understanding its evolution and for guiding conservation efforts. Calosaturnia is a charismatic silkmoth genus endemic to large portions of the region with three described species, C. mendocino, C. walterorum, and C. albofasciata. We sampled all three species from across their ranges, sequenced 1463 bp of mitochondrial COI and 1941 bp of nuclear DNA from three genes, and reconstructed phylogenetic relationships and estimated divergence times within the lineages. All three species show pronounced evidence of isolation and, in two cases, secondary reconnection. An unexpected monophyletic mtDNA lineage was found in the Central Coast region, in a region thought to represent an intergrade between C. mendocino and C. walterorum. Our genetic data also significantly revise previous hypotheses as to which species occur in which regions, suggesting that historical ecological changes around four Ma ago isolated some lineages, and a secondary isolation event two Ma ago led to isolation of populations both in the Central Coast region and the southern Sierra Nevada. Our results add to a currently under-appreciated pattern suggesting that coastal Central California is not a transition zone between Northern and Southern California Floristic Province faunas but rather its own unique, periodically isolated, biogeographic region. They also suggest cryptic diversity may be present in many other currently unrecognized groups. Additional research should focus on this central California region because many species may be highly restricted in range and in need of conservation attention.

Genome-wide genetic variation coupled with demographic and ecological niche modeling of the dusky-footed woodrat (Neotoma fuscipes) reveal patterns of deep divergence and widespread Holocene expansion across northern California

Understanding how species have responded to past climate change may help refine projections of how species and biotic communities will respond to future change. Here, we integrate estimates of genome-wide genetic variation with demographic and niche modeling to investigate the historical biogeography of an important ecological engineer: the dusky-footed woodrat, Neotoma fuscipes. We use RADseq to generate a genome-wide dataset for 71 individuals from across the geographic distribution of the species in California. We estimate population structure using several model-based methods and infer the demographic history of regional populations using a site frequency spectrum-based approach. Additionally, we use ecological niche modeling to infer current and past (Last Glacial Maximum) environmental suitability across the species' distribution. Finally, we estimate the directionality and possible spatial origins of regional population expansions. Our analyses indicate this species is subdivided into three regionally distinct populations, with the deepest divergence occurring ~1.7 million years ago across the modern-day San Francisco-Bay Delta region; a common biogeographic barrier for the flora and fauna of California. Our models of environmental suitability through time coincide with our estimates of population expansion, with relative long-term stability in the southern portion of the range, and more recent expansion into the northern end of the range. Our study illustrates how the integration of genome-wide data with spatial and demographic modeling can reveal the timing and spatial extent of historic events that determine patterns of biotic diversity and may help predict biotic response to future change.

The influence of latitude, geographic distance, and habitat discontinuities on genetic variation in a high latitude montane species

Examining the factors that influence contemporary genetic patterns is important given the alarming rate at which natural environments are changing. In particular habitat fragmentation and climate change are expected to influence the distribution and diversity of natural populations. In this study we used both mitochondrial control region (mtDNA) and microsatellite data to answer the following questions about genetic diversity and divergence in mountain chickadees (Poecile gambeli) a resident bird species in western North America: (1) Do populations exhibit similar levels of genetic diversity across the range? (2) What is the genetic affinity of western populations in Oregon and Washington? (3) Do genetic patterns exhibit isolation by distance, or are genetic patterns more heavily influenced by habitat discontinuity? We tested the effects of isolation by distance and habitat distribution on genetic structure by analyzing 266 samples from 17 sites across western Canada and the United States. We found a near significant relationship between genetic diversity and latitude, however, our results indicate that overall, latitude is not a strong predictor of genetic diversity. Our analyses of populations in Oregon and Washington revealed a mismatch between patterns detected with mtDNA and microsatellite data. In particular, Washington clustered with the Coast Range/Cascades/Rocky Mountain mtDNA group, but with populations in southern Oregon/California based on microsatellite data. These results suggest the presence of a contact zone in Washington between the two mtDNA clades Coast Range/Cascades/Rocky Mountain and southern Oregon/California clades. Finally, our study revealed a greater effect of isolation by distance than isolation by habitat for both mtDNA and microsatellite data. Overall the isolation by distance signal was greater for mtDNA than microsatellite patterns. The greater signal of isolation by distance on mtDNA patterns likely reflects the strong effects of Pleistocene glaciations in shaping genetic patterns in western North America.

Montane refugia predict population genetic structure in the Large-blotched Ensatina salamander

Understanding the biotic consequences of Pleistocene range shifts and fragmentation remains a fundamental goal in historical biogeography and evolutionary biology. Here, we combine species distribution models (SDM) from the present and two late Quaternary time periods with multilocus genetic data (mitochondrial DNA and microsatellites) to evaluate the effect of climate-induced habitat shifts on population genetic structure in the Large-blotched Ensatina (Ensatina eschscholtzii klauberi), a plethodontid salamander endemic to middle and high-elevation conifer forest in the Transverse and Peninsular Ranges of southern California and northern Baja California. A composite SDM representing the range through time predicts two disjunct refugia, one in southern California encompassing the core of the species range and the other in the Sierra San Pedro Mártir of northern Baja California at the southern limit of the species range. Based on our spatial model, we would expect a pattern of high connectivity among populations within the northern refugium and, conversely, a pattern of isolation due to long-term persistence of the Sierra San Pedro Mártir population. Our genetic results are consistent with these predictions based on the hypothetical refugia in that (i) historical measures of population connectivity among stable areas are correlated with gene flow estimates; and (ii) there is strong geographical structure between separate refugia. These results provide evidence for the role of recent climatic change in shaping patterns of population persistence and connectivity within the Transverse and Peninsular Ranges, an evolutionary hotspot.

Mitochondrial DNA variation reveals recent evolutionary history of main Boa constrictor clades

We sequenced a 1114-bp fragment of cytochrome b gene in six subspecies (115 samples) of Boa constrictor and detected 67 haplotypes. Our analyses revealed the presence of two distinct clades, one from Central America (CA) including the neighboring part of South America west of the Andes, and the other covering the rest of South America (SA). Sequence divergence between CA and SA clades is about 5-7%, which roughly corresponds to a separation at the time of uplift of the Colombian Andes following formation of the Panama Isthmus before 3.5 Myr Sequence divergence within the SA and CA clades is only 2-3%, suggesting a fairly recent spread of these clades Into their current geographic ranges. Thus, we may not be dealing with taxa with a markedly old evolutionary history. Because juveniles of B. constrictor feed mostly on small rodents, we hypothesized that spread of this species was allowed by a new food source represented by murold rodents that appeared after closure of the Panama portal. With respect to the taxonomy, B. c. imperator may be elevated to full species rank. Within the SA clade, a haplotype of Argentinian B. c. occidentalis is markedly distinct, while the remaining haplotype groups analyzed are distributed throughout large ranges and may all belong to a single nominotypic subspecies.

TOWARD A BETTER UNDERSTANDING OF THE "TRANSVERSE RANGE BREAK": LINEAGE DIVERSIFICATION IN SOUTHERN CALIFORNIA

The Transverse Ranges in southern California have been identified as having a prominent phylogeographic role. Numerous studies have identified distinct north-south and/or east-west lineage breaks involving the Transverse Ranges. However, in evaluating their findings, most authors have regarded this complex system somewhat simplistically. In this study we more deeply investigate these breaks using two approaches: first we examine the phylogeographic history of Sepedophilus castaneus (Coleoptera: Staphylinidae) and then implement a comparative phylogeography approach applying Brooks parsimony analysis to the topologies of nine additional taxa. Phylogenetic analysis, nested clade analysis, and AMOVAs for S. castaneus agree that there is a major lineage break between the eastern and western Transverse Ranges, localized between the Sierra Pelona and the San Gabriel Mountains. The comparative phylogeographic analysis supports a generally strong concordance of area relationships with geographic proximity. It is notable, however, that the Transverse Ranges as a group do not show phylogenetic cohesion, but rather they are split into three main regions: an eastern region (San Gabriel, San Bernardino, and San Jacinto Mountains), a central region (central Transverse Ranges and Sierra Pelona) that is often grouped with the Tehachapi and Sierra Nevada populations, and a western region (northwestern Transverse Ranges and Santa Ynez Mountains) that is consistently grouped with coast range areas to the north. The lineage break between east and west Transverse Ranges is attributable to the presence of a marine embayment in what is now the Santa Clara River valley 5-2.5 million years ago.

Glacial Effects on Sequence Divergence of Mitochondrial COII of Polyura eudamippus (Lepidoptera: Nymphalidae) in China

Sequence divergence of mitochondrial COII was analyzed in 50 specimens belonging to five subspecies of Polyura eudamippus (Lepidoptera: Nymphalidae) collected from southern China. There were nine haplotypes in the 405 bp of partial COII sequence. Distribution of the five subspecies was primarily consistent with the spatial distribution of haplotypes. The K (st) statistic showed genetic differentiation among these subspecies, except between the pair of P. e. kuangtungensis and P. e. formosana, which were separated by the Taiwan Strait. This is consistent with the 10,000-year history of the Taiwan Strait, not long enough for detectable differentiation. The present distribution pattern of COII haplotypes of P. eudamippus should be shaped by the alteration of Pleistocene glaciations, and Yunnan might be the refugium of P. eudamippus in the ice age, judging from the abundant haplotypes remaining. There were two routes for P. eudamippus in the postglacial expansion, one northward to Sichuan, Chongqing, and Hubei and another eastward to the southeastern coast of mainland China and Taiwan Island. Because the haplotype of butterflies on Hainan Island (P. e. whiteheadi) was completely different from that of mainland China, it was estimated that butterflies on Hainan Island might be from the Indo-China Peninsula rather than from mainland China.

Comparative phylogeography of woodland reptiles in California: repeated patterns of cladogenesis and population expansion

The ultimate goal of comparative phylogeographical analyses is to infer processes of diversification from contemporary geographical patterns of genetic diversity. When such studies are employed across diverse groups in an array of communities, it may be difficult to discover common evolutionary and ecological processes associated with diversification. In order to identify taxa that have responded in a similar fashion to historical events, we conducted comparative phylogeographical analyses on a phylogenetically and ecologically limited set of taxa. Here, we focus on a group of squamate reptiles (snakes and lizards) that share similar ecological requirements and generally occupy the same communities in the western USA. At a gross level, deep genetic division in Contia tenuis, Diadophis punctatus, Elgaria multicarinata, the Charina bottae complex, and Lampropeltis zonata are often concordant in the Transverse Ranges, the Monterey Bay and Sacramento-San Joaquin Delta region, and the southern Sierra Nevada in California. Molecular clock estimates suggest that major phyletic breaks within many of these taxa roughly coincide temporally, and may correspond to important geological events. Furthermore, significant congruence between the phylogeographies of E. multicarinata and L. zonata suggests that the succession of vicariance and dispersal events in these species progressed in concert. Such congruence suggests that E. multicarinata and L. zonata have occupied the same communities through time. However, across our entire multi-taxon data set, the sequence of branching events rarely match between sympatric taxa, indicating the importance of subtle differences in life history features as well as random processes in creating unique genetic patterns. Lastly, coalescent and noncoalescent estimates of population expansion suggest that populations in the more southerly distributed clades of C. tenuis, D. punctatus, E. multicarinata, and L. zonata have been stable, while populations in more northerly clades appear to have recently expanded. This concerted demographic response is consistent with palaeontological data and previous phylogeographical work that suggests that woodland habitat has become more restricted in southern California, but more widespread in the North during Holocene warming. Future phylogeographical work focusing on allied and ecologically associated taxa may add insight into the ecological and evolutionary processes that yield current patterns of genetic diversity.

Shifting sands and shifty lizards: molecular phylogeny and biogeography of African flat lizards (Platysaurus)

The African flat lizard genus Platysaurus is widely distributed on rock outcrops in southern Africa and is found both east and west of the Kalahari Desert and between major river drainage systems. We assembled a molecular phylogeny for the genus in order to test several biogeographic hypotheses. Sequence data were obtained from 29 specimens representing 14 taxa of Platysaurus that span the geographic range of the genus. We targeted a fragment of the mitochondrial genome comprising the 3' half of the ND4 gene and most of the flanking tRNA-HSL cluster. The edited alignment comprised 864 characters, of which 479 (55%) were variable and 461 (96%) parsimony informative. Overall, the phylogeny was well resolved and supported by high bootstrap values. Four major clades were identified comprising two to seven species: P. mitchelli and P. maculatus maculatus from the north-eastern range of the genus; P. broadleyi and P. capensis from the western range; P. imperator, P. torquatus, and P. intermedius rhodesianus; P. i. intermedius, P. monotropis, P. minor, P. i. nigrescens, P. lebomboensis, P. i. wilhelmi, and P. o. orientalis. Platysaurus has been suggested to represent a recent adaptive radiation where rapid speciation was fuelled by population fragmentation brought on by vicariant events and possibly divergent sexual selection. The traditional explanation for the radiation of the genus is that the eastern migration of the Kalahari sands fragmented populations in the Plio-Pleistocene, resulting in conditions favorable for speciation. Our genetic data strongly suggests that many of the speciation events in Platysaurus already had occurred prior to the Plio-Pleistocene. Moreover, vicariant events associated with the formation of the major river systems played an additional role in the evolution and distribution of Platysaurus species. Our topology displays long internodes and long terminal branches, suggesting that the radiation is much older than previously believed.

Species boundaries, phylogeography and conservation genetics of the red-legged frog (Rana aurora/draytonii) complex

The red-legged frog, Rana aurora, has been recognized as both a single, polytypic species and as two distinct species since its original description 150 years ago. It is currently recognized as one species with two geographically contiguous subspecies, aurora and draytonii; the latter is protected under the US Endangered Species Act. We present the results of a survey of 50 populations of red-legged frogs from across their range plus four outgroup species for variation in a phylogenetically informative, approximately 400 base pairs (bp) fragment of the mitochondrial cytochrome b gene. Our mtDNA analysis points to several major results. (1) In accord with several other lines of independent evidence, aurora and draytonii are each diagnosably distinct, evolutionary lineages; the mtDNA data indicate that they do not constitute a monophyletic group, but rather that aurora and R. cascadae from the Pacific northwest are sister taxa; (2) the range of the draytonii mtDNA clade extends about 100 km further north in coastal California than was previously suspected, and corresponds closely with the range limits or phylogeographical breaks of several codistributed taxa; (3) a narrow zone of overlap exists in southern Mendocino County between aurora and draytonii haplotypes, rather than a broad intergradation zone; and (4) the critically endangered population of draytonii in Riverside County, CA forms a distinct clade with frogs from Baja California, Mexico. The currently available evidence favours recognition of aurora and draytonii as separate species with a narrow zone of overlap in northern California.

Phylogeography of the California Thrasher (Toxostoma Redivivum) Based on Nested-Clade Analysis of Mitochondrial-DNA Variation

Distribution of genealogical lineages within a species is likely the result of a complicated series of ecological and historical events. Nested-clade analysis is specifically designed as an objective phylogeographic approach for inferring evolutionary processes on a spatial and temporal scale for small subclades within a larger set of intraspecific relationships. Here, we use nested-clade analysis as well as other phylogeographic methods to investigate the evolutionary history of California Thrasher (Toxostoma redivivum) populations. Inferences resulting from nested clade analysis suggest a history that includes past fragmentation, range expansion, and isolation-by-distance. Along with root information, those inferences enable the construction of a biogeographic scenario for this species involving general southern ancestry, an early north–south division, northward range expansion, and a southward back-expansion into an already populated southern region. Isolation-by-distance is also identified, particularly in southern California, indicating that gene flow between localities does occur but is restricted. Many conclusions drawn from this study are concordant with geologic data as well as phylogeographic scenarios drawn for other codistributed California taxa.

The evolutionary response of predators to dangerous prey: hotspots and coldspots in the geographic mosaic of coevolution between garter snakes and newts

The "geographic mosaic" approach to understanding coevolution is predicated on the existence of variable selection across the landscape of an interaction between species. A range of ecological factors, from differences in resource availability to differences in community composition, can generate such a mosaic of selection among populations, and thereby differences in the strength of coevolution. The result is a mixture of hotspots, where reciprocal selection is strong, and coldspots, where reciprocal selection is weak or absent, throughout the ranges of species. Population subdivision further provides the opportunity for nonadaptive forces, including gene flow, drift, and metapopulation dynamics, to influence the coevolutionary interaction between species. Some predicted results of this geographic mosaic of coevolution include maladapted or mismatched phenotypes, maintenance of high levels of polymorphism, and prevention of stable equilibrium trait combinations. To evaluate the potential for the geographic mosaic to influence predator-prey coevolution, we investigated the geographic pattern of genetically determined TTX resistance in the garter snake Thamnophis sirtalis over much of the range of its ecological interaction with toxic newts of genus Taricha. We assayed TTX resistance in over 2900 garter snakes representing 333 families from 40 populations throughout western North America. Our results provide dramatic evidence that geographic structure is an important component in coevolutionary interactions between predators and prey. Resistance levels vary substantially (over three orders of magnitude) among populations and over short distances. The spatial array of variation is consistent with two areas of intense evolutionary response by predators ("hotspots") surrounded by clines of decreasing resistance. Some general predictions of the geographic mosaic process are supported, including clinal variation in phenotypes, polymorphism in some populations, and divergent outcomes of the interaction between predator and prey. Conversely, our data provide little support for one of the major predictions, mismatched values of interacting traits. Two lines of evidence suggest selection is paramount in determining population variation in resistance. First, phylogenetic information indicates that two hotspots of TTX resistance have evolved independently. Second, in the one region that TTX levels in prey have been quantified, resistance and toxicity levels match almost perfectly over a wide phenotypic and geographic range. However, these results do not preclude the role the nonadaptive forces in generating the overall geographic mosaic of TTX resistance. Much work remains to fill in the geographic pattern of variation among prey populations and, just as importantly, to explore the variation in the ecology of the interaction that occurs within populations.

Molecular phylogeography of common garter snakes (Thamnophis sirtalis) in western North America: implications for regional historical forces

Complete ND2 and partial ND4 and cytochrome b mitochondrial DNA (mtDNA) sequences were analysed to evaluate the phylogeographic patterns of common garter snakes (Thamnophis sirtalis) in western North America. This species is widely distributed throughout North America, and exhibits extensive phenotypic variation in the westernmost part of its range. The overall phylogeographic pattern based on mtDNA sequences is concordant with results from studies of other species in this region, implicating historical vicariant processes during the Pleistocene and indicating bottleneck effects of recent dispersal into postglacial habitat. Indeed, the topology is statistically consistent with the hypothesis of both southern (Great Basin and California) and northern (Haida Gwaii) refugia. Specifically, we identified genetic breaks among three major clades: Northwest Coastal populations, Intermountain populations, and all California populations. The California clade contained the only other well-supported branching patterns detected; relationships among populations within the two northern clades were indistinguishable. These molecular splits contrast sharply with all prior geographical analyses of phenotypic variation in T. sirtalis in this region. Our results suggest that the extensive phenotypic variation in western T. sirtalis has been shaped more by local evolutionary forces than by shared common ancestry. Consequently, we consider all morphologically based subspecies designations of T. sirtalis in this region invalid because they do not reflect reciprocal monophyly of the mtDNA sequences.

Molecular phylogenetics of emydine turtles: taxonomic revision and the evolution of shell kinesis

The 10 extant species of emydine turtles represent an array of morphological and ecological forms recognizable and popular among scientists and hobbyists. Nevertheless, the phylogenetic affinities of most emydines remain contentious. Here, we examine the evolutionary relationships of emydine turtles using 2092 bp of DNA encoding the mitochondrial genes cyt b, ND4, and adjacent tRNAs. These data contain 339 parsimony informative characters that we use to erect hypotheses of relationships for the Emydinae. Both maximum parsimony and maximum likelihood methods yield a monophyletic Emydinae in which all but three nodes are well resolved. Emys orbicularis, Emydoidea blandingii, and Clemmys marmorata form a monophyletic clade, as do the species of Terrapene. Clemmys muhlenbergii and Clemmys insculpta form a third monophyletic group that may be sister to all other emydines. Clemmys guttata is problematic and probably related to Terrapene. Based on this phylogeny, and previous molecular work on the group, we suggest the following taxonomic revisions: (1) Clemmys should be restricted to a single species, C. guttata. (2) Calemys should be resurrected for C. muhlenbergii and C. insculpta. (3) Emys should be expanded to include three species: E. orbicularis, E. blandingii, and E. marmorata. Furthermore, our analyses show that neither kinetic-shelled nor akinetic-shelled emydines form monophyletic groups. Therefore, shell kinesis was either independently gained in Emys and Terrapene or secondarily lost in E. marmorata and C. guttata. Parsimony, paleontological evidence, and the multiple origins of shell kinesis in related turtle lineages (especially geoemydines) support the independent origin of plastral kinesis.