Genetic drift and rapid evolution of viviparity in insular fire salamanders (Salamandra salamandra)

Continental islands offer an excellent opportunity to investigate adaptive processes and to time microevolutionary changes that precede macroevolutionary events. We performed a population genetic study of the fire salamander (Salamandra salamandra), a species that displays unique intraspecific diversity of reproductive strategies, to address the microevolutionary processes leading to phenotypic and genetic differentiation of island, coastal and interior populations. We used eight microsatellite markers to estimate genetic diversity, population structure and demographic parameters in viviparous insular populations and ovoviviparous coastal and interior populations. Our results show considerable genetic differentiation (F(ST) range: 0.06-0.27), and no clear signs of gene flow among populations, except between the large and admixed interior populations. We find no support for island colonization by rafting or intentional/accidental anthropogenic introductions, indicating that rising sea levels were responsible for isolation of the island populations approximately 9000 years ago. Our study provides evidence of rapid genetic differentiation between island and coastal populations, and rapid evolution of viviparity driven by climatic selective pressures on island populations, geographic isolation with genetic drift, or a combination of these factors. Studies of these viviparous island populations in early stages of divergence help us better understand the microevolutionary processes involved in rapid phenotypic shifts.     

Population genetics of Anopheles koliensis through Papua New Guinea: New cryptic species and landscape topography effects on genetic connectivity

New Guinea is a topographically and biogeographically complex region that supports unique endemic fauna. Studies describing the population connectivity of species through this region are scarce. We present a population and landscape genetic study on the endemic malaria‐transmitting mosquito, Anopheles koliensis (Owen). Using mitochondrial and nuclear sequence data, as well as microsatellites, we show the evidence of geographically discrete population structure within Papua New Guinea (PNG). We also confirm the existence of three rDNA ITS2 genotypes within this mosquito and assess reproductive isolation between individuals carrying different genotypes. Microsatellites reveal the clearest population structure and show four clear population units. Microsatellite markers also reveal probable reproductive isolation between sympatric populations in northern PNG with different ITS2 genotypes, suggesting that these populations may represent distinct cryptic species. Excluding individuals belonging to the newly identified putative cryptic species (ITS2 genotype 3), we modeled the genetic differences between A. koliensis populations through PNG as a function of terrain and find that dispersal is most likely along routes with low topographic relief. Overall, these results show that A. koliensis is made up of geographically and genetically discrete populations in Papua New Guinea with landscape topography being important in restricting dispersal.     

Spatial genetic structure and landscape connectivity in black bears: Investigating the significance of using different land cover datasets and classifications in landscape genetics analyses

Landscape genetic analyses allow detection of fine‐scale spatial genetic structure (SGS) and quantification of effects of landscape features on gene flow and connectivity. Typically, analyses require generation of resistance surfaces. These surfaces characteristically take the form of a grid with cells that are coded to represent the degree to which landscape or environmental features promote or inhibit animal movement. How accurately resistance surfaces predict association between the landscape and movement is determined in large part by (a) the landscape features used, (b) the resistance values assigned to features, and (c) how accurately resistance surfaces represent landscape permeability. Our objective was to evaluate the performance of resistance surfaces generated using two publicly available land cover datasets that varied in how accurately they represent the actual landscape. We genotyped 365 individuals from a large black bear population (Ursus americanus) in the Northern Lower Peninsula (NLP) of Michigan, USA at 12 microsatellite loci, and evaluated the relationship between gene flow and landscape features using two different land cover datasets. We investigated the relative importance of land cover classification and accuracy on landscape resistance model performance. We detected local spatial genetic structure in Michigan''s NLP black bears and found roads and land cover were significantly correlated with genetic distance. We observed similarities in model performance when different land cover datasets were used despite 21% dissimilarity in classification between the two land cover datasets. However, we did find the performance of land cover models to predict genetic distance was dependent on the way the land cover was defined. Models in which land cover was finely defined (i.e., eight land cover classes) outperformed models where land cover was defined more coarsely (i.e., habitat/non‐habitat or forest/non‐forest). Our results show that landscape genetic researchers should carefully consider how land cover classification changes inference in landscape genetic studies.     

植物景观遗传学研究进展

植物景观遗传学是新兴的景观遗传学交叉学科的一个重要研究方向。目前植物景观遗传学的研究虽落后于动物,但其在生物多样性保护方面具有的巨大潜力不可忽视。从景观特征对遗传结构、环境因素对适应性遗传变异影响两个方面,系统综述了近十年来国际上植物景观遗传学的研究焦点和研究进展,比较了植物景观遗传学与动物景观遗传学研究在研究设计和研究方法上的异同,并基于将来植物景观遗传学由对空间遗传结构的描述发展为对景观遗传效应的量化分析及预测的发展框架,具体针对目前景观特征与遗传结构研究设计的系统性差、遗传结构与景观格局在时间上的误配、适应性位点与环境变量的模糊匹配、中性遗传变异与适应性遗传变异研究的分隔、景观与遗传关系分析方法的局限等五个方面提出了研究对策。     

Spatial genetic analysis reveals high connectivity of tiger (Panthera tigris) populations in the Satpura–Maikal landscape of Central India

We investigated the spatial genetic structure of the tiger meta‐population in the Satpura–Maikal landscape of central India using population‐ and individual‐based genetic clustering methods on multilocus genotypic data from 273 individuals. The Satpura–Maikal landscape is classified as a global‐priority Tiger Conservation Landscape (TCL) due to its potential for providing sufficient habitat that will allow the long‐term persistence of tigers. We found that the tiger meta‐population in the Satpura–Maikal landscape has high genetic variation and very low genetic subdivision. Individual‐based Bayesian clustering algorithms reveal two highly admixed genetic populations. We attribute this to forest connectivity and high gene flow in this landscape. However, deforestation, road widening, and mining may sever this connectivity, impede gene exchange, and further exacerbate the genetic division of tigers in central India.     

Trapped within the city: integrating demography,time since isolation and population‐specific traits to assess the genetic effects of urbanization

Urbanization is a severe form of habitat fragmentation that can cause many species to be locally extirpated and many others to become trapped and isolated within an urban matrix. The role of drift in reducing genetic diversity and increasing genetic differentiation is well recognized in urban populations. However, explicit incorporation and analysis of the demographic and temporal factors promoting drift in urban environments are poorly studied. Here, we genotyped 15 microsatellites in 320 fire salamanders from the historical city of Oviedo (Est. 8th century) to assess the effects of time since isolation, demographic history (historical effective population size; N_e) and patch size on genetic diversity, population structure and contemporary N_e. Our results indicate that urban populations of fire salamanders are highly differentiated, most likely due to the recent N_e declines, as calculated in coalescence analyses, concomitant with the urban development of Oviedo. However, urbanization only caused a small loss of genetic diversity. Regression modelling showed that patch size was positively associated with contemporary N_e, while we found only moderate support for the effects of demographic history when excluding populations with unresolved history. This highlights the interplay between different factors in determining current genetic diversity and structure. Overall, the results of our study on urban populations of fire salamanders provide some of the very first insights into the mechanisms affecting changes in genetic diversity and population differentiation via drift in urban environments, a crucial subject in a world where increasing urbanization is forecasted.     

Increased genetic structuring of isolated Salamandra salamandra populations (Caudata: Salamandridae) at the margins of the Carpathian Mountains

The traditional southern Pleistocene refugia hypothesis in Europe has lately been challenged for several animal and plant species. The Carpathian Basin, especially at the marginal regions, is one of the recently recognized biodiversity hotspots in Europe. Marginal populations are prone to have lower genetic diversity and higher genetic differentiation than central populations. Here, we examined one mitochondrial DNA fragment (D‐loop) and nine nuclear (microsatellite) loci to describe the genetic diversity and phylogeographical pattern of fire salamander (Salamandra salamandra) populations in the Carpathian Basin with focusing on the southern margins of the Western Carpathians, where isolated populations of this species are present. Analyses of microsatellites indicated reduced genetic diversity for most of the isolated populations. Based on the mitochondrial DNA, only two haplotypes were found, whereas the analyses with the nuclear markers revealed a more recent genetic split between Western (Alpine) and Eastern (Carpathian) populations, and separated the Apuseni Mountains population (part of the Western Carpathians). Using approximate Bayesian computation analyses, we identified the most probable colonization scenario for the isolated North Hungarian Carpathian Basin populations. The split between isolated salamander populations from the central populations in the Carpathian Mountains dates back to the beginning of the Late Pleistocene, while the split between most of the Hungarian populations can be associated with the Last Glacial Maximum. We found evidence for long‐time isolation between the marginal Carpathian Basin and central populations. Our results also show that S. salamandra survived glacial periods in the temperate forests of north‐east Pannonia (North Hungarian Mountains), confirming that the Carpathian Basin served as important northerly refugia during the Pleistocene climatic oscillations.     

Node‐based measures of connectivity in genetic networks

At‐site environmental conditions can have strong influences on genetic connectivity, and in particular on the immigration and settlement phases of dispersal. However, at‐site processes are rarely explored in landscape genetic analyses. Networks can facilitate the study of at‐site processes, where network nodes are used to model site‐level effects. We used simulated genetic networks to compare and contrast the performance of 7 node‐based (as opposed to edge‐based) genetic connectivity metrics. We simulated increasing node connectivity by varying migration in two ways: we increased the number of migrants moving between a focal node and a set number of recipient nodes, and we increased the number of recipient nodes receiving a set number of migrants. We found that two metrics in particular, the average edge weight and the average inverse edge weight, varied linearly with simulated connectivity. Conversely, node degree was not a good measure of connectivity. We demonstrated the use of average inverse edge weight to describe the influence of at‐site habitat characteristics on genetic connectivity of 653 American martens (Martes americana) in Ontario, Canada. We found that highly connected nodes had high habitat quality for marten (deep snow and high proportions of coniferous and mature forest) and were farther from the range edge. We recommend the use of node‐based genetic connectivity metrics, in particular, average edge weight or average inverse edge weight, to model the influences of at‐site habitat conditions on the immigration and settlement phases of dispersal.     

Dispersal of viviparity across contact zones in Iberian populations of fire salamanders (Salamandra) inferred from discordance of genetic and morphological traits

Abstract.— We used partial sequences of the cytochrome b mitochondrial DNA (mtDNA) gene, obtained from 76 individuals representing 45 populations of Iberian Salamandra salamandra plus 15 sequences of additional species of Salamandra and related genera, to investigate contact zones. These zones, identified by earlier allozymic and morphological analyses, are between populations of viviparous (S. s. bernardezi and S. s. fastuosa ) and ovoviviparous (S. s. gallaica and S. s. terrestris ) salamanders. The distribution of mtDNA and nuclear markers is mostly concordant at one contact zone (between S. s. gallaica and S. s. bernardezi ), but at another (between S. s. fastuosa and S. s. terrestris ) the markers are offset by about 250 km. The observed geographic variation fits a model of mtDNA capture. Among the potential mechanisms responsible for such discordance we favor a combination of range shifts due to climatic fluctuations and biased genetic admixture across moving contact zones. We apply our findings to the issue of possible homoplasy in the evolution of viviparity and conclude that viviparity likely arose only once within S. salamandra .     

Directed dispersal by rotational shepherding supports landscape genetic connectivity in a calcareous grassland plant

Directed dispersal by animal vectors has been found to have large effects on the structure and dynamics of plant populations adapted to frugivory. Yet, empirical data are lacking on the potential of directed dispersal by rotational grazing of domestic animals to mediate gene flow across the landscape. Here, we investigated the potential effect of large‐flock shepherding on landscape‐scale genetic structure in the calcareous grassland plant Dianthus carthusianorum, whose seeds lack morphological adaptations to dispersal to animals or wind. We found a significant pattern of genetic structure differentiating population within grazed patches of three nonoverlapping shepherding systems and populations of ungrazed patches. Among ungrazed patches, we found a strong and significant effect of isolation by distance (r = 0.56). In contrast, genetic distance between grazed patches within the same herding system was unrelated to geographical distance but significantly related to distance along shepherding routes (r = 0.44). This latter effect of connectivity along shepherding routes suggests that gene flow is spatially restricted occurring mostly between adjacent populations. While this study used nuclear markers that integrate gene flow by pollen and seed, the significant difference in the genetic structure between ungrazed patches and patches connected by large‐flock shepherding indicates the potential of directed seed dispersal by sheep across the landscape.     

Spatial genetic analysis of two polyploid macrophytes reveals high connectivity in a modified wetland

相似文献   

Watershed characteristics shape the landscape genetics of brook stickleback (Culaea inconstans) in shallow prairie lakes

Investigating the consequences of landscape features on population genetic patterns is increasingly important to elucidate the ecological factors governing connectivity between populations and predicting the evolutionary consequences of landscapes. Small prairie lakes in Alberta, Canada, and the brook stickleback (Culaea inconstans) that inhabit them, provide a unique aquatic system whereby populations are highly isolated from one another. These heterogeneous and extreme environments are prone to winterkills, an event whereby most of the fish die and frequent bottlenecks occur. In this study, we characterized the genetic population structure of brook stickleback among several lakes, finding that the species is hierarchically influenced by within‐lake characteristics in small‐scale watersheds. Landscape genetic analyses of the role of spatial features found support for basin characteristics associated with genetic diversity and bottlenecks in 20% of the sampled lakes. These results suggest that brook stickleback population genetic patterns may be driven, at least in part, by ecological processes that accelerate genetic drift and landscape patterns associated with reduced dispersal. Collectively, these results reinforce the potential importance of connectivity in the maintenance of genetic diversity, especially in fragmented landscapes.     

The influence of breeding phenology on the genetic structure of four pond‐breeding salamanders

Understanding metapopulation dynamics requires knowledge about local population dynamics and movement in both space and time. Most genetic metapopulation studies use one or two study species across the same landscape to infer population dynamics; however, using multiple co‐occurring species allows for testing of hypotheses related to different life history strategies. We used genetic data to study dispersal, as measured by gene flow, in three ambystomatid salamanders (Ambystoma annulatum , A. maculatum , and A. opacum ) and the Central Newt (Notophthalmus viridescens louisianensis ) on the same landscape in Missouri, USA . While all four salamander species are forest dependent organisms that require fishless ponds to reproduce, they differ in breeding phenology and spatial distribution on the landscape. We use these differences in life history and distribution to address the following questions: (1) Are there species‐level differences in the observed patterns of genetic diversity and genetic structure? and (2) Is dispersal influenced by landscape resistance? We detected two genetic clusters in A. annulatum and A. opacum on our landscape; both species breed in the fall and larvae overwinter in ponds. In contrast, no structure was evident in A. maculatum and N. v. louisianensis , species that breed during the spring. Tests for isolation by distance were significant for the three ambystomatids but not for N. v. louisianensis . Landscape resistance also contributed to genetic differentiation for all four species. Our results suggest species‐level differences in dispersal ability and breeding phenology are driving observed patterns of genetic differentiation. From an evolutionary standpoint, the observed differences in dispersal distances and genetic structure between fall breeding and spring breeding species may be a result of the trade‐off between larval period length and size at metamorphosis which in turn may influence the long‐term viability of the metapopulation. Thus, it is important to consider life history differences among closely related and ecologically similar species when making management decisions.     

Range‐wide analysis of genetic structure in a widespread,highly mobile species (Odocoileus hemionus) reveals the importance of historical biogeography

Highly mobile species that thrive in a wide range of habitats are expected to show little genetic differentiation across their range. A limited but growing number of studies have revealed that patterns of broad‐scale genetic differentiation can and do emerge in vagile, continuously distributed species. However, these patterns are complex and often shaped by both historical and ecological factors. Comprehensive surveys of genetic variation at a broad scale and at high resolution are useful for detecting cryptic spatial genetic structure and for investigating the relative roles of historical and ecological processes in structuring widespread, highly mobile species. In this study, we analysed 10 microsatellite loci from over 1900 samples collected across the full range of mule deer (Odocoileus hemionus), one of the most widely distributed and abundant of all large mammal species in North America. Through both individual‐ and population‐based analyses, we found evidence for three main genetic lineages, one corresponding to the ‘mule deer’ morphological type and two to the ‘black‐tailed deer’ type. Historical biogeographic events likely are the primary drivers of genetic divergence in this species; boundaries of the three lineages correspond well with predictions based on Pleistocene glacial cycles, and substructure within each lineage demonstrates island vicariance. However, across large geographic areas, including the entire mule deer lineage, we found that genetic variation fit an isolation‐by‐distance pattern rather than discrete clusters. A lack of genetic structure across wide geographic areas of the continental west indicates that ecological processes have not resulted in restrictions to gene flow sufficient for spatial genetic structure to emerge. Our results have important implications for our understanding of evolutionary mechanisms of divergence, as well as for taxonomy, conservation and management.     

Functional connectivity and home range inferred at a microgeographic landscape genetics scale in a desert‐dwelling rodent

Gene flow in animals is limited or facilitated by different features within the landscape matrix they inhabit. The landscape representation in landscape genetics (LG) is traditionally modeled as resistance surfaces (RS), where novel optimization approaches are needed for assigning resistance values that adequately avoid subjectivity. Also, desert ecosystems and mammals are scarcely represented in LG studies. We addressed these issues by evaluating, at a microgeographic scale, the effect of landscape features on functional connectivity of the desert‐dwelling Dipodomys merriami. We characterized genetic diversity and structure with microsatellites loci, estimated home ranges and movement of individuals using telemetry—one of the first with rodents, generated a set of individual and composite environmental surfaces based on hypotheses of variables influencing movement, and assessed how these variables relate to individual‐based gene flow. Genetic diversity and structure results evidenced a family‐induced pattern driven by first‐order‐related individuals, notably determining landscape genetic inferences. The vegetation cover and soil resistance optimized surface (NDVI) were the best‐supported model and a significant predictor of individual genetic distance, followed by humidity and NDVI+humidity. Based on an accurate definition of thematic resolution, we also showed that vegetation is better represented as continuously (vs. categorically) distributed. Hence, with a nonsubjective optimization framework for RS and telemetry, we were able to describe that vegetation cover, soil texture, and climatic variables influence D. merriami's functional connectivity at a microgeographic scale, patterns we could further explain based on the home range, habitat use, and activity observed between sexes. We describe the relationship between environmental features and some aspects of D. merriami‘s behavior and physiology.     

Comparative landscape genetics reveals differential effects of environment on host and pathogen genetic structure in Tasmanian devils (Sarcophilus harrisii) and their transmissible tumour

Genetic structure in host species is often used to predict disease spread. However, host and pathogen genetic variation may be incongruent. Understanding landscape factors that have either concordant or divergent influence on host and pathogen genetic structure is crucial for wildlife disease management. Devil facial tumour disease (DFTD) was first observed in 1996 and has spread throughout almost the entire Tasmanian devil geographic range, causing dramatic population declines. Whereas DFTD is predominantly spread via biting among adults, devils typically disperse as juveniles, which experience low DFTD prevalence. Thus, we predicted little association between devil and tumour population structure and that environmental factors influencing gene flow differ between devils and tumours. We employed a comparative landscape genetics framework to test the influence of environmental factors on patterns of isolation by resistance (IBR) and isolation by environment (IBE) in devils and DFTD. Although we found evidence for broad‐scale costructuring between devils and tumours, we found no relationship between host and tumour individual genetic distances. Further, the factors driving the spatial distribution of genetic variation differed for each. Devils exhibited a strong IBR pattern driven by major roads, with no evidence of IBE. By contrast, tumours showed little evidence for IBR and a weak IBE pattern with respect to elevation in one of two tumour clusters we identify herein. Our results warrant caution when inferring pathogen spread using host population genetic structure and suggest that reliance on environmental barriers to host connectivity may be ineffective for managing the spread of wildlife diseases. Our findings demonstrate the utility of comparative landscape genetics for identifying differential factors driving host dispersal and pathogen transmission.     

Dealing with uncertainty in landscape genetic resistance models: a case of three co‐occurring marsupials

Landscape genetics lacks explicit methods for dealing with the uncertainty in landscape resistance estimation, which is particularly problematic when sample sizes of individuals are small. Unless uncertainty can be quantified, valuable but small data sets may be rendered unusable for conservation purposes. We offer a method to quantify uncertainty in landscape resistance estimates using multimodel inference as an improvement over single model‐based inference. We illustrate the approach empirically using co‐occurring, woodland‐preferring Australian marsupials within a common study area: two arboreal gliders (Petaurus breviceps, and Petaurus norfolcensis) and one ground‐dwelling antechinus (Antechinus flavipes). First, we use maximum‐likelihood and a bootstrap procedure to identify the best‐supported isolation‐by‐resistance model out of 56 models defined by linear and non‐linear resistance functions. We then quantify uncertainty in resistance estimates by examining parameter selection probabilities from the bootstrapped data. The selection probabilities provide estimates of uncertainty in the parameters that drive the relationships between landscape features and resistance. We then validate our method for quantifying uncertainty using simulated genetic and landscape data showing that for most parameter combinations it provides sensible estimates of uncertainty. We conclude that small data sets can be informative in landscape genetic analyses provided uncertainty can be explicitly quantified. Being explicit about uncertainty in landscape genetic models will make results more interpretable and useful for conservation decision‐making, where dealing with uncertainty is critical.     

Contrasting effects of landscape features on genetic structure in different geographic regions in the ornate dragon lizard,Ctenophorus ornatus

Habitat fragmentation can have profound effects on the distribution of genetic variation within and between populations. Previously, we showed that in the ornate dragon lizard, Ctenophorus ornatus, lizards residing on outcrops that are separated by cleared agricultural land are significantly more isolated and hold less genetic variation than lizards residing on neighbouring outcrops connected by undisturbed native vegetation. Here, we extend the fine‐scale study to examine the pattern of genetic variation and population structure across the species' range. Using a landscape genetics approach, we test whether land clearing for agricultural purposes has affected the population structure of the ornate dragon lizard. We found significant genetic differentiation between outcrop populations (F_ST = 0.12), as well as isolation by distance within each geographic region. In support of our previous study, land clearing was associated with higher genetic divergences between outcrops and lower genetic variation within outcrops, but only in the region that had been exposed to intense agriculture for the longest period of time. No other landscape features influenced population structure in any geographic region. These results show that the effects of landscape features can vary across species' ranges and suggest there may be a temporal lag in response to contemporary changes in land use. These findings therefore highlight the need for caution when assessing the impact of contemporary land use practices on genetic variation and population structure.     

A comparison of individual‐based genetic distance metrics for landscape genetics

A major aim of landscape genetics is to understand how landscapes resist gene flow and thereby influence population genetic structure. An empirical understanding of this process provides a wealth of information that can be used to guide conservation and management of species in fragmented landscapes and also to predict how landscape change may affect population viability. Statistical approaches to infer the true model among competing alternatives are based on the strength of the relationship between pairwise genetic distances and landscape distances among sampled individuals in a population. A variety of methods have been devised to quantify individual genetic distances, but no study has yet compared their relative performance when used for model selection in landscape genetics. In this study, we used population genetic simulations to assess the accuracy of 16 individual‐based genetic distance metrics under varying sample sizes and degree of population genetic structure. We found most metrics performed well when sample size and genetic structure was high. However, it was much more challenging to infer the true model when sample size and genetic structure was low. Under these conditions, we found genetic distance metrics based on principal components analysis were the most accurate (although several other metrics performed similarly), but only when they were derived from multiple principal components axes (the optimal number varied depending on the degree of population genetic structure). Our results provide guidance for which genetic distance metrics maximize model selection accuracy and thereby better inform conservation and management decisions based upon landscape genetic analysis.     

Detection of environmental and morphological adaptation despite high landscape genetic connectivity in a pest grasshopper (Phaulacridium vittatum)

Widespread species that exhibit both high gene flow and the capacity to occupy heterogeneous environments make excellent models for examining local selection processes along environmental gradients. Here we evaluate the influence of temperature and landscape variables on genetic connectivity and signatures of local adaptation in Phaulacridium vittatum, a widespread agricultural pest grasshopper, endemic to Australia. With sampling across a 900‐km latitudinal gradient, we genotyped 185 P. vittatum from 19 sites at 11,408 single nucleotide polymorphisms (SNPs) using ddRAD sequencing. Despite high gene flow across sites (pairwise F_ST = 0.0003–0.08), landscape genetic resistance modelling identified a positive nonlinear effect of mean annual temperature on genetic connectivity. Urban areas and water bodies had a greater influence on genetic distance among sites than pasture, agricultural areas and forest. Together, F_ST outlier tests and environmental association analysis (EAA) detected 242 unique SNPs under putative selection, with the highest numbers associated with latitude, mean annual temperature and body size. A combination of landscape genetic connectivity analysis together with EAA identified mean annual temperature as a key driver of both neutral gene flow and environmental selection processes. Gene annotation of putatively adaptive SNPs matched with gene functions for olfaction, metabolic detoxification and ultraviolet light shielding. Our results imply that this widespread agricultural pest has the potential to spread and adapt under shifting temperature regimes and land cover change.