Rapid vertebrate speciation via isolation, bottlenecks, and drift.

Speciation is often driven by selective processes like those associated with viability, mate choice, or local adaptation, and "speciation genes" have been identified in many eukaryotic lineages. In contrast, neutral processes are rarely considered as the primary drivers of speciation, especially over short evolutionary timeframes. Here, we describe a rapid vertebrate speciation event driven primarily by genetic drift. The White Sands pupfish (Cyprinodon tularosa) is endemic to New Mexico's Tularosa Basin where the species is currently managed as two Evolutionarily significant units (ESUs) and is of international conservation concern (Endangered). Whole-genome resequencing data from each ESU showed remarkably high and uniform levels of differentiation across the entire genome (global FST ≈ 0.40). Despite inhabiting ecologically dissimilar springs and streams, our whole-genome analysis revealed no discrete islands of divergence indicative of strong selection, even when we focused on an array of candidate genes. Demographic modeling of the joint allele frequency spectrum indicates the two ESUs split only ~4 to 5 kya and that both ESUs have undergone major bottlenecks within the last 2.5 millennia. Our results indicate the genome-wide disparities between the two ESUs are not driven by divergent selection but by neutral drift due to small population sizes, geographic isolation, and repeated bottlenecks. While rapid speciation is often driven by natural or sexual selection, here we show that isolation and drift have led to speciation within a few thousand generations. We discuss these evolutionary insights in light of the conservation management challenges they pose.

The endangered White Sands pupfish (Cyprinodon tularosa) genome reveals low diversity and heterogenous patterns of differentiation.

The White Sands pupfish (Cyprinodon tularosa), endemic to New Mexico in Southwestern North America, is of conservation concern due in part to invasive species, chemical pollution, and groundwater withdrawal. Herein, we developed a draft reference genome and use it to provide biological insights into the evolution and conservation of C. tularosa. We used our assembly to localize microsatellite markers previously used to demarcate evolutionary significant units (ESU), quantified genomic divergence and transposable element profiles between species, and compared C. tularosa genomic diversity related species. Our de novo assembly of PacBio Sequel II error-corrected reads resulted in a 1.08 Gb draft genome with a contig N50 of 1.4 Mb and 25,260 annotated protein coding genes, including 95% of the expected Actinopterygii conserved complete single-copy orthologues. Many of the C. tularosa microsatellite markers used for conservation assessments fell within, or near, genes and exhibited a pattern of increased heterozygosity near genic areas compared to those in intergenic regions. Nuclear alignments between these two species revealed 193 genes contained in rapidly diverging tracts; transposable element profiles were largely concordant and suggest a shared, rapid expansion of LINE and Gypsy elements. Genome-wide heterozygosity was markedly lower in C. tularosa compared to estimates from other related species, probably because of smaller long-term effective population sizes constrained by their isolated and limited habitat. Overall, these inferences provide new insights into C. tularosa that should help inform future management efforts.

Ten-year projection of white-nose syndrome disease dynamics at the southern leading-edge of infection in North America.

Predicting the emergence and spread of infectious diseases is critical for the effective conservation of biodiversity. White-nose syndrome (WNS), an emerging infectious disease of bats, has resulted in high mortality in eastern North America. Because the fungal causative agent Pseudogymnoascus destructans is constrained by temperature and humidity, spread dynamics may vary by geography. Environmental conditions in the southern part of the continent are different than the northeast, where disease dynamics are typically studied, making it difficult to predict how the disease will manifest. Herein, we modelled WNS pathogen spread in Texas based on cave densities and average dispersal distances of hosts, projecting these results out to 10 years. We parameterized a predictive model of WNS epidemiology and its effects on bat populations with observed cave environmental data. Our model suggests that bat populations in northern Texas will be more affected by WNS mortality than southern Texas. As such, we recommend prioritizing the preservation of large overwintering colonies of bats in north Texas through management actions. Our model illustrates that infectious disease spread and infectious disease severity can become uncoupled over a gradient of environmental variation and highlight the importance of understanding host, pathogen and environmental conditions across a breadth of environments.

Correction: Winter habitats of bats in Texas.

[This corrects the article DOI: 10.1371/journal.pone.0220839.].

Winter habitats of bats in Texas.

Few studies have described winter microclimate selection by bats in the southern United States. This is of particular importance as the cold-adapted fungus, Pseudogymnoascus destructans, which causes the fatal bat disease white-nose syndrome (WNS), continues to spread into southern United States. To better understand the suitability of winter bat habitats for the growth of P. destructans in this region, we collected roost temperature and vapor pressure deficit from 97 hibernacula in six ecoregions in Texas during winter 2016-17 and 2017-18. We also measured skin temperature of Rafinesque's big-eared bats (Corynorhinus townsendii), Townsend's big-eared bats (C. townsendii), big-brown bats (Eptesicus fuscus), southeastern myotis (Myotis austroriparius), cave myotis (M. velifer), tri-colored bats (Perimyotis subflavus), and Mexican free-tailed bats (Tadarida brasiliensis) during hibernation to study their use of torpor in these habitats. We found that temperatures within hibernacula were strongly correlated with external air temperatures and were often within the optimal range of temperatures for P. destructans growth. Hibernacula and skin temperatures differed among species, with Rafinesque's big-eared bats, southeastern myotis, and Mexican free-tailed bats occupying warmer microclimates and having higher torpid skin temperatures. For species that were broadly distributed throughout Texas, hibernacula and skin temperatures differed within species by ecoregion; Tri-colored bats and cave myotis in colder, northern regions occupied colder microclimates within hibernacula and exhibited colder skin temperatures, than individuals of the same species in warmer, southern regions. These data illustrate the variability in microclimates used as hibernacula by bats in Texas and suggest similar variation in susceptibility to WNS in the state. Thus, monitoring microclimates at winter roosts may help predict where WNS may develop, and where management efforts would be most effective.

External temperature and distance from nearest entrance influence microclimates of cave and culvert-roosting tri-colored bats (Perimyotis subflavus).

Many North American bat species hibernate in both natural and artificial roosts. Although hibernacula can have high internal climate stability, they still retain spatial variability in their thermal regimes, resulting in various "microclimates" throughout the roost that differ in their characteristics (e.g., temperature and air moisture). These microclimate components can be influenced by factors such as the number of entrances, the depth of the roost, and distance to the nearest entrance of the roost. Tri-colored bats are commonly found roosting in caves in winter, but they can also be found roosting in large numbers in culverts, providing the unique opportunity to investigate factors influencing microclimates of bats in both natural and artificial roost sites. As tri-colored bats are currently under consideration for federal listing, information of this type could be useful in aiding in the conservation and management of this species through a better understanding of what factors affect the microclimate near roosting bats. We collected data on microclimate temperature and microclimate actual water vapor pressure (AWVP) from a total of 760 overwintering tri-colored bats at 18 caves and 44 culverts. Using linear mixed models analysis, we found that variation in bat microclimate temperatures was best explained by external temperature and distance from nearest entrance in both caves and culverts. External temperature had a greater influence on microclimate temperatures in culverts than caves. We found that variation in microclimate AWVP was best explained by external temperature, distance from nearest entrance, and proportion from entrance (proportion of the total length of the roost from the nearest entrance) in culvert-roosting bats. Variation in microclimate AWVP was best explained by external temperature and proportion from entrance in cave-roosting bats. Our results suggest that bat microclimate temperature and AWVP are influenced by similar factors in both artificial and natural roosts, although the relative contribution of these factors differs between roost types.

Inadequate thermal refuge constrains landscape habitability for a grassland bird species.

Ecologists have long recognized the influence that environmental conditions have on abundance and range extent of animal species. We used the northern bobwhite (Colinus virginianus; hereafter bobwhite) as a model species for studying how microclimates serve as refuge against severe weather conditions. This species serves as an indicator or umbrella species for other sensitive ground-nesting, grassland obligate species. We conducted a mensurative field experiment in the rolling plains of Texas, USA, a semi-arid ecosystem on the southwestern periphery of bobwhite range, to determine whether native bunch grasses, apparently suitable for bobwhite nesting, could reduce ambient temperature below levels harmful for eggs. During the nesting season, we compared temperature and relative humidity readings at daily heat maxima (i.e., the 3 h during each day with highest temperatures) during the nesting season over the course of two years at 63 suitable nest sites paired with 63 random locations (n = 126) using two sensors at ∼10 and ∼60 cm above ground level. Mean temperature at nest height was 2.3% cooler at nest sites (35.99 °C ± 0.07 SE) compared to random locations (36.81 °C ± 0.07 SE); at ambient height, nest sites were slightly cooler (32.78 °C ± 0.06 SE) than random location (32.99 °C ± 0.06 SE). Mean relative humidity at nest sites was greater at nest height (34.53% ± 0.112 SE) and ambient height (36.22% ± 0.10 SE) compared to random locations at nest (33.35% ± 0.12 SE) and ambient height (35.75% ± 0.10 SE). Based on these results, cover at sites that appear visually suitable for nesting by bobwhites and other ground nesting birds provided adequate thermal refuge in the rolling plains by maintaining cooler, moister microclimates than surrounding non-nesting locations. Post-hoc analyses of data revealed that habitat conditions surrounding suitable nest sites strongly influenced thermal suitability of the substrate. Given that eggs of bobwhites and probably other species would experience lethal temperatures without these thermal refuges in the context of proper habitat condition, nesting vegetation is a critical component of niche space for bobwhites and other ground nesting birds in semi-arid regions. Many contemporary land uses, however, degrade or destroy bunch grasses and grassland systems, and thus decrease landscape inhabitability. Conservationists working with obligate grassland species that require bunch grasses in semi-arid regions should develop land management strategies that maximize the availability of these thermal refuges across space and time.