Population genomic structure of a widespread, urban-dwelling mammal: The eastern grey squirrel (Sciurus carolinensis).

Urbanization is a persistent and widespread driver of global environmental change, potentially shaping evolutionary processes due to genetic drift and reduced gene flow in cities induced by habitat fragmentation and small population sizes. We tested this prediction for the eastern grey squirrel (Sciurus carolinensis), a common and conspicuous forest-dwelling rodent, by obtaining 44K SNPs using reduced representation sequencing (ddRAD) for 403 individuals sampled across the species' native range in eastern North America. We observed moderate levels of genetic diversity, low levels of inbreeding, and only a modest signal of isolation-by-distance. Clustering and migration analyses show that estimated levels of migration and genetic connectivity were higher than expected across cities and forested areas, specifically within the eastern portion of the species' range dominated by urbanization, and genetic connectivity was less than expected within the western range where the landscape is fragmented by agriculture. Landscape genetic methods revealed greater gene flow among individual squirrels in forested regions, which likely provide abundant food and shelter for squirrels. Although gene flow appears to be higher in areas with more tree cover, only slight discontinuities in gene flow suggest eastern grey squirrels have maintained connected populations across urban areas in all but the most heavily fragmented agricultural landscapes. Our results suggest urbanization shapes biological evolution in wildlife species depending strongly on the composition and habitability of the landscape matrix surrounding urban areas.

Population genomic analyses support sympatric origins of parapatric morphs in a salamander.

In numerous clades, divergent sister species have largely non-overlapping geographic ranges. This pattern presumably arises because species diverged in allopatry or parapatry, prior to a subsequent contact. Here, we provide population-genomic evidence for the opposite scenario: previously sympatric ecotypes that have spatially separated into divergent monomorphic populations over large geographic scales (reverse sympatric scenario). We analyzed a North American salamander (Plethodon cinereus) with two color morphs that are broadly sympatric: striped (redback) and unstriped (leadback). Sympatric morphs can show considerable divergence in other traits, and many Plethodon species are fixed for a single morph. Long Island (New York) is unusual in having many pure redback and leadback populations that are spatially separated, with pure redback populations in the west and pure leadbacks in the east. Previous work showed that these pure-morph populations were genetically, morphologically, and ecologically divergent. Here, we performed a coalescent-based analysis of new data from 88,696 single-nucleotide polymorphisms to address the origins of these populations. This analysis strongly supports the monophyly of Long Island populations and their subsequent divergence into pure redback and pure leadback populations. Taken together, these results suggest that the formerly sympatric mainland morphs separated into parapatric populations on Long Island, reversing the conventional speciation scenario.

SNAPSHOT USA 2020: A second coordinated national camera trap survey of the United States during the COVID-19 pandemic.

Managing wildlife populations in the face of global change requires regular data on the abundance and distribution of wild animals, but acquiring these over appropriate spatial scales in a sustainable way has proven challenging. Here we present the data from Snapshot USA 2020, a second annual national mammal survey of the USA. This project involved 152 scientists setting camera traps in a standardized protocol at 1485 locations across 103 arrays in 43 states for a total of 52,710 trap-nights of survey effort. Most (58) of these arrays were also sampled during the same months (September and October) in 2019, providing a direct comparison of animal populations in 2 years that includes data from both during and before the COVID-19 pandemic. All data were managed by the eMammal system, with all species identifications checked by at least two reviewers. In total, we recorded 117,415 detections of 78 species of wild mammals, 9236 detections of at least 43 species of birds, 15,851 detections of six domestic animals and 23,825 detections of humans or their vehicles. Spatial differences across arrays explained more variation in the relative abundance than temporal variation across years for all 38 species modeled, although there are examples of significant site-level differences among years for many species. Temporal results show how species allocate their time and can be used to study species interactions, including between humans and wildlife. These data provide a snapshot of the mammal community of the USA for 2020 and will be useful for exploring the drivers of spatial and temporal changes in relative abundance and distribution, and the impacts of species interactions on daily activity patterns. There are no copyright restrictions, and please cite this paper when using these data, or a subset of these data, for publication.

SNAPSHOT USA 2019: a coordinated national camera trap survey of the United States.

With the accelerating pace of global change, it is imperative that we obtain rapid inventories of the status and distribution of wildlife for ecological inferences and conservation planning. To address this challenge, we launched the SNAPSHOT USA project, a collaborative survey of terrestrial wildlife populations using camera traps across the United States. For our first annual survey, we compiled data across all 50 states during a 14-week period (17 August-24 November of 2019). We sampled wildlife at 1,509 camera trap sites from 110 camera trap arrays covering 12 different ecoregions across four development zones. This effort resulted in 166,036 unique detections of 83 species of mammals and 17 species of birds. All images were processed through the Smithsonian's eMammal camera trap data repository and included an expert review phase to ensure taxonomic accuracy of data, resulting in each picture being reviewed at least twice. The results represent a timely and standardized camera trap survey of the United States. All of the 2019 survey data are made available herein. We are currently repeating surveys in fall 2020, opening up the opportunity to other institutions and cooperators to expand coverage of all the urban-wild gradients and ecophysiographic regions of the country. Future data will be available as the database is updated at eMammal.si.edu/snapshot-usa, as will future data paper submissions. These data will be useful for local and macroecological research including the examination of community assembly, effects of environmental and anthropogenic landscape variables, effects of fragmentation and extinction debt dynamics, as well as species-specific population dynamics and conservation action plans. There are no copyright restrictions; please cite this paper when using the data for publication.

Inferring introgression using RADseq and DFOIL : Power and pitfalls revealed in a case study of spiny lizards (Sceloporus).

Introgression is now commonly reported in studies across the Tree of Life, aided by recent advancements in data collection and analysis. Nevertheless, researchers working with nonmodel species lacking reference genomes may be stymied by a mismatch between available resources and methodological demands. In this study, we demonstrate a fast and simple approach for inferring introgression using RADseq data, and apply it to a case study involving spiny lizards (Sceloporus) from northeastern México. First, we find evidence for recurrent mtDNA introgression between the two focal species based on patterns of mito-nuclear discordance. We then test for nuclear introgression by exhaustively applying the "five-taxon" D-statistic (DFOIL ) to all relevant individuals sampled for RADseq data. In our case, this exhaustive approach (dubbed "ExDFOIL ") entails testing up to ~250,000 unique four-taxon combinations of individuals across species. To facilitate use of this ExDFOIL approach, we provide scripts for many relevant tasks, including the selection of appropriate four-taxon combinations, execution of DFOIL tests in parallel and visualization of introgression results in phylogenetic and geographic space. Using ExDFOIL , we find evidence for ancient introgression between the focal species. Furthermore, we reveal geographic variation in patterns of introgression that is consistent with patterns of mito-nuclear discordance and with recurrent introgression. Overall, our study demonstrates that the combination of DFOIL and RADseq data can effectively detect introgression under a variety of sampling conditions (for individuals, populations and loci). Importantly, we also find evidence that batch-specific error and linkage in RADseq data may mislead inferences of introgression under certain conditions.

Parapatric divergence of sympatric morphs in a salamander: incipient speciation on Long Island?

Speciation is often categorized based on geographic modes (allopatric, parapatric or sympatric). Although it is widely accepted that species can arise in allopatry and then later become sympatrically or parapatrically distributed, patterns in the opposite direction are also theoretically possible (e.g. sympatric lineages or ecotypes becoming parapatric), but such patterns have not been shown at a macrogeographic scale. Here, we analyse genetic, climatic, ecological and morphological data and show that two typically sympatric colour morphs of the salamander Plethodon cinereus (redback and leadback) appear to have become parapatrically distributed on Long Island, New York, with pure-redback populations in the west and pure-leadback populations in the east (and polymorphic populations in between and on the mainland). In addition, the pure-leadback populations in eastern Long Island are genetically, ecologically and morphologically divergent from both mainland and other Long Island populations, suggesting the possibility of incipient speciation. This parapatric separation seems to be related to the different ecological preferences of the two morphs, preferences which are present on the mainland and across Long Island. These results potentially support the idea that spatial segregation of sympatric ecotypes may sometimes play an important part in parapatric speciation.

How does climate change cause extinction?

Anthropogenic climate change is predicted to be a major cause of species extinctions in the next 100 years. But what will actually cause these extinctions? For example, will it be limited physiological tolerance to high temperatures, changing biotic interactions or other factors? Here, we systematically review the proximate causes of climate-change related extinctions and their empirical support. We find 136 case studies of climatic impacts that are potentially relevant to this topic. However, only seven identified proximate causes of demonstrated local extinctions due to anthropogenic climate change. Among these seven studies, the proximate causes vary widely. Surprisingly, none show a straightforward relationship between local extinction and limited tolerances to high temperature. Instead, many studies implicate species interactions as an important proximate cause, especially decreases in food availability. We find very similar patterns in studies showing decreases in abundance associated with climate change, and in those studies showing impacts of climatic oscillations. Collectively, these results highlight our disturbingly limited knowledge of this crucial issue but also support the idea that changing species interactions are an important cause of documented population declines and extinctions related to climate change. Finally, we briefly outline general research strategies for identifying these proximate causes in future studies.

How is the rate of climatic-niche evolution related to climatic-niche breadth?

The rate of climatic-niche evolution is important to many research areas in ecology, evolution, and conservation biology, including responses of species to global climate change, spread of invasive species, speciation, biogeography, and patterns of species richness. Previous studies have implied that clades with higher rates of climatic-niche evolution among species should have species with narrower niche breadths, but there is also evidence suggesting the opposite pattern. However, the relationships between rate and breadth have not been explicitly analyzed. Here, we examine the relationships between the rate of climatic-niche evolution and climatic-niche breadth using phylogenetic and climatic data for 250 species in the salamander family Plethodontidae, a group showing considerable variation in both rates of climatic-niche evolution and climatic-niche breadths. Contrary to some expectations, we find no general relationship between climatic-niche breadth and the rate of climatic-niche evolution. Climatic-niche breadths for some ecologically important climatic variables considered separately (temperature seasonality and annual precipitation) do show significant relationships with the rate of climatic-niche evolution, but rates are faster in clades in which species have broader (not narrower) niche breadths. In summary, our results show that narrower niche breadths are not necessarily associated with faster rates of niche evolution.

What are the consequences of combining nuclear and mitochondrial data for phylogenetic analysis? Lessons from Plethodon salamanders and 13 other vertebrate clades.

BACKGROUND: The use of mitochondrial DNA data in phylogenetics is controversial, yet studies that combine mitochondrial and nuclear DNA data (mtDNA and nucDNA) to estimate phylogeny are common, especially in vertebrates. Surprisingly, the consequences of combining these data types are largely unexplored, and many fundamental questions remain unaddressed in the literature. For example, how much do trees from mtDNA and nucDNA differ? How are topological conflicts between these data types typically resolved in the combined-data tree? What determines whether a node will be resolved in favor of mtDNA or nucDNA, and are there any generalities that can be made regarding resolution of mtDNA-nucDNA conflicts in combined-data trees? Here, we address these and related questions using new and published nucDNA and mtDNA data for Plethodon salamanders and published data from 13 other vertebrate clades (including fish, frogs, lizards, birds, turtles, and mammals). RESULTS: We find widespread discordance between trees from mtDNA and nucDNA (30-70% of nodes disagree per clade), but this discordance is typically not strongly supported. Despite often having larger numbers of variable characters, mtDNA data do not typically dominate combined-data analyses, and combined-data trees often share more nodes with trees from nucDNA alone. There is no relationship between the proportion of nodes shared between combined-data and mtDNA trees and relative numbers of variable characters or levels of homoplasy in the mtDNA and nucDNA data sets. Congruence between trees from mtDNA and nucDNA is higher on branches that are longer and deeper in the combined-data tree, but whether a conflicting node will be resolved in favor mtDNA or nucDNA is unrelated to branch length. Conflicts that are resolved in favor of nucDNA tend to occur at deeper nodes in the combined-data tree. In contrast to these overall trends, we find that Plethodon have an unusually large number of strongly supported conflicts between data types, which are generally resolved in favor of mtDNA in the combined-data tree (despite the large number of nuclear loci sampled). CONCLUSIONS: Overall, our results from 14 vertebrate clades show that combined-data analyses are not necessarily dominated by the more variable mtDNA data sets. However, given cases like Plethodon, there is also the need for routine checking of incongruence between mtDNA and nucDNA data and its impacts on combined-data analyses.