The influence of historical landscape change on genetic variation and population structure of a terrestrial salamander (Plethodon cinereus)

Genetic divergence correlates with the contemporary landscape in populations of Slimy Salamander (Plethodon glutinosus) species complex across the lower Piedmont and Coastal Plain of the southeastern United States

A primary goal of landscape genetics is to elucidate factors associated with genetic structure among populations. Among the important patterns identified have been isolation by distance (IBD), isolation by barrier (IBB), and isolation by environment (IBE). We tested hypotheses relating each of these possible patterns to genetic divergence in the Slimy Salamander (Plethodon glutinosus (Green, 1818)) species complex across the lower Piedmont and Coastal Plain of Georgia, USA, and adjacent areas of South Carolina, USA. We sequenced 2148 total bp, including three regions of the mitochondrial genome and a nuclear intron, and related genetic distance to GIS-derived surrogate variables representing possible IBD (geographic distance), IBE (principal components of 19 climate variables, watershed, and normalized difference vegetation index (NDVI)), and IBB (streams of fourth order and higher). Multiple matrix regression with randomization analysis indicated significant relationships between genetic distance and two principal components of climate, as well as NDVI. These results support roles for environment (IBE) in helping to drive genetic divergence in this group of salamanders. The absence of a significant influence of IBD and IBB was surprising. It is possible that the signal effects of geographic distance and barriers on genetic divergence may have been erased by more recent responses to the environment.

Assessment of intra and interregional genetic variation in the Eastern Red-backed Salamander, Plethodon cinereus, via analysis of novel microsatellite markers

The red-backed salamander (Plethodon cinereus) has long-served as a model system in ecology, evolution, and behavior, and studies surveying molecular variation in this species have become increasingly common over the past decade. However, difficulties are commonly encountered when extending microsatellite markers to populations that are unstudied from a genetic perspective due to high levels of genetic differentiation across this species’ range. To ameliorate this issue, we used 454 pyrosequencing to identify hundreds of microsatellite loci. We then screened 40 of our top candidate loci in populations in Virginia, Pennsylvania, and Ohio—including an isolated island population ~ 4.5 km off the shore of Lake Erie (South Bass Island). We identified 25 loci that are polymorphic in a well-studied region of Virginia and 11 of these loci were polymorphic in populations located in the genetically unstudied regions of Ohio and Pennsylvania. Use of these loci to examine patterns of variation within populations revealed that South Bass Island has low diversity in comparison to other sites. However, neither South Bass Island nor isolated populations around Cleveland are inbred. Assessment of variation between populations revealed three well defined genetic clusters corresponding to Virginia, mainland Ohio/Pennsylvania, and South Bass Island. Comparisons of our results to those of others working in various parts of the range are consistent with the idea that differentiation is lower in regions that were once glaciated. However, these comparisons also suggest that well differentiated isolated populations in the formerly glaciated portion of the range are not uncommon. This work provides novel genetic resources that will facilitate population genetic studies in a part of the red-backed salamander’s range that has not previously been studied in this manner. Moreover, this work refines our understanding of how neutral variation is distributed in this ecologically important organism.

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.

Current and historical drivers of landscape genetic structure differ in core and peripheral salamander populations

With predicted decreases in genetic diversity and greater genetic differentiation at range peripheries relative to their cores, it can be difficult to distinguish between the roles of current disturbance versus historic processes in shaping contemporary genetic patterns. To address this problem, we test for differences in historic demography and landscape genetic structure of coastal giant salamanders (Dicamptodon tenebrosus) in two core regions (Washington State, United States) versus the species' northern peripheral region (British Columbia, Canada) where the species is listed as threatened. Coalescent-based demographic simulations were consistent with a pattern of post-glacial range expansion, with both ancestral and current estimates of effective population size being much larger within the core region relative to the periphery. However, contrary to predictions of recent human-induced population decline in the less genetically diverse peripheral region, there was no genetic signature of population size change. Effects of current demographic processes on genetic structure were evident using a resistance-based landscape genetics approach. Among core populations, genetic structure was best explained by length of the growing season and isolation by resistance (i.e. a ‘flat’ landscape), but at the periphery, topography (slope and elevation) had the greatest influence on genetic structure. Although reduced genetic variation at the range periphery of D. tenebrosus appears to be largely the result of biogeographical history rather than recent impacts, our analyses suggest that inherent landscape features act to alter dispersal pathways uniquely in different parts of the species' geographic range, with implications for habitat management.

Simulating Pattern-Process Relationships to Validate Landscape Genetic Models

Landscapes may resist gene flow and thereby give rise to a pattern of genetic isolation within a population. The mechanism by which a landscape resists gene flow can be inferred by evaluating the relationship between landscape models and an observed pattern of genetic isolation. This approach risks false inferences because researchers can never feasibly test all plausible alternative hypotheses. In this paper, rather than infer the process of gene flow from an observed genetic pattern, we simulate gene flow and determine if the simulated genetic pattern is related to the observed empirical genetic pattern. This is a form of inverse modeling and can be used to independently validate a landscape genetic model. In this study, we used this approach to validate a model of landscape resistance based on elevation, landcover, and roads that was previously related to genetic isolation among mountain goats (Oreamnos americanus) inhabiting the Cascade Range, Washington (USA). The strong relationship between the empirical and simulated patterns of genetic isolation we observed provides independent validation of the resistance model and demonstrates the utility of this approach in supporting landscape genetic inferences.