The influence of altitude and topography on genetic structure in the long-toed salamander (Ambystoma macrodactulym)

Evolutionary history of a cold-adapted limnephilid caddisfly: Effects of climate change and topography on genetic structure

The distribution of organisms is influenced by complex factors such as the phylogenetic evolutionary histories of species, the physiological and ecological characteristics of organisms, climate, and geographical and geohistorical features. In this study, we focused on a caddisfly, Asynarchus sachalinensis (Trichoptera: Limnephilidae), which has adapted to cold habitats. From phylogeographic analyses based on the mitochondrial DNA (mtDNA) cytochrome c oxidase subunit I (COI) and 16S rRNA regions and the nuclear DNA (nDNA) 18S rRNA, 28S rRNA, carbamoyl-phosphate synthetase (CAD), elongation factor-1 alpha (EF1-α), and RNA polymerase II (POLII) regions, two distinct genetic clades were detected. Clade I was shown to be widely distributed from Sakhalin to Honshu, whereas Clade II was only distributed within Honshu. The distributions of these clades overlapped in Honshu. The habitats were located at relatively lower altitudes for Clade I and higher altitudes for Clade II. The divergence time of these clades was estimated to be during the Pleistocene, indicating that repeated climatic changes facilitated distributional shifts. Haplotype network and demographic analyses based on the mtDNA COI region showed contrasting genetic structures in the two clades. It was indicated that the population sizes of Clade I had expanded rapidly in a recent period, whereas Clade II had maintained stable population sizes. The habitats of Clade II were typically isolated and scattered at high altitudes, resulting in restricted migration and dispersal because of their discontinuous "Sky Island" habitats. The habitats of Clade I were located at relatively low altitudes, and it was assumed that the populations were continuous, which resulted in a higher frequency of migration and dispersal between populations. Thus, differences in the spatial scale of the adapted habitats of each clade may have resulted in different patterns of population connectivity and fragmentation associated with repeated climatic changes during the Pleistocene. Our study provided new insight into the distributional patterns of cold-adapted aquatic insects in the Japanese Archipelago. Furthermore, the distributional shifts predicted by ecological niche modeling under future climatic change conditions were different for each clade. Therefore, different principles are required in the assessment of each clade to predict temporal changes in their distributions.

Conservation genomics of an endangered montane amphibian reveals low population structure, low genomic diversity and selection pressure from disease

Wildlife diseases are a major global threat to biodiversity. Boreal toads (Anaxyrus [Bufo] boreas) are a state-endangered species in the southern Rocky Mountains of Colorado and New Mexico, and a species of concern in Wyoming, largely due to lethal skin infections caused by the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). We performed conservation and landscape genomic analyses using single nucleotide polymorphisms from double-digest, restriction site-associated DNA sequencing in combination with the development of the first boreal toad (and first North American toad) reference genome to investigate population structure, genomic diversity, landscape connectivity and adaptive divergence. Genomic diversity (π = 0.00034-0.00040) and effective population sizes (Ne = 8.9-38.4) were low, likely due to post-Pleistocene founder effects and Bd-related population crashes over the last three decades. Population structure was also low, likely due to formerly high connectivity among a higher density of geographically proximate populations. Boreal toad gene flow was facilitated by low precipitation, cold minimum temperatures, less tree canopy, low heat load and less urbanization. We found >8X more putatively adaptive loci related to Bd intensity than to all other environmental factors combined, and evidence for genes under selection related to immune response, heart development and regulation and skin function. These data suggest boreal toads in habitats with Bd have experienced stronger selection pressure from disease than from other, broad-scale environmental variations. These findings can be used by managers to conserve and recover the species through actions including reintroduction and supplementation of populations that have declined due to Bd.

Modeling flood susceptibility zones using hybrid machine learning models of an agricultural dominant landscape of India

Flooding events are determining a significant amount of damages, in terms of economic loss and also casualties in Asia and Pacific areas. Due to complexity and ferocity of severe flooding, predicting flood-prone areas is a difficult task. Thus, creating flood susceptibility maps at local level is though challenging but an inevitable task. In order to implement a flood management plan for the Balrampur district, an agricultural dominant landscape of India, and strengthen its resilience, flood susceptibility modeling and mapping are carried out. In the present study, three hybrid machine learning (ML) models, namely, fuzzy-ANN (artificial neural network), fuzzy-RBF (radial basis function), and fuzzy-SVM (support vector machine) with 12 topographic, hydrological, and other flood influencing factors were used to determine flood-susceptible zones. To ascertain the relationship between the occurrences and flood influencing factors, correlation attribute evaluation (CAE) and multicollinearity diagnostic tests were used. The predictive power of these models was validated and compared using a variety of statistical techniques, including Wilcoxon signed-rank, t-paired tests and receiver operating characteristic (ROC) curves. Results show that fuzzy-RBF model outperformed other hybrid ML models for modeling flood susceptibility, followed by fuzzy-ANN and fuzzy-SVM. Overall, these models have shown promise in identifying flood-prone areas in the basin and other basins around the world. The outcomes of the work would benefit policymakers and government bodies to capture the flood-affected areas for necessary planning, action, and implementation.

Computer vision for assessing species color pattern variation from web-based community science images

Openly available community science digital vouchers provide a wealth of data to study phenotypic change across space and time. However, extracting phenotypic data from these resources requires significant human effort. Here, we demonstrate a workflow and computer vision model for automatically categorizing species color pattern from community science images. Our work is focused on documenting the striped/unstriped color polymorphism in the Eastern Red-backed Salamander (Plethodon cinereus). We used an ensemble convolutional neural network model to analyze this polymorphism in 20,318 iNaturalist images. Our model was highly accurate (∼98%) despite image heterogeneity. We used the resulting annotations to document extensive niche overlap between morphs, but wider niche breadth for striped morphs at the range-wide scale. Our work showcases key design principles for using machine learning with heterogeneous community science image data to address questions at an unprecedented scale.

•

We built a deep learning model to group color morphs from community science images

•

Our model achieved 98% accuracy for classifying striped and unstriped salamanders

•

We used our model to classify >20,000 images and built morph-specific niche models

•

We then determined if Red-backed salamanders niche partition at a range-wide scale

Topographic Complexity Facilitates Persistence Compared to Signals of Contraction and Expansion in the Adjacent Subdued Landscape

Topographically heterogeneous areas are likely to act as refugia for species because they facilitate survival during regional climatic stress due to availability of a range of microenvironments. The Stirling Ranges are a topographically complex area in the generally subdued and ancient landscape of south-western Australia. We investigated the influence of these landscape features on the evolutionary history of the rare woody shrub, Banksia brownii through a combined approach using phylogeographic analysis of sequence data from three chloroplast sequences, the trnV–ndhC, trnQ–rps16, and rpl32–ndhF intergenic spacer regions, and species distribution modeling. The Stirling Ranges showed high genetic diversity and differentiation among populations consistent with localized persistence and maintenance of large populations in an area that species distribution modeling identified as providing habitat stability at the Last Glacial Maximum as well as under warmer conditions. In contrast, populations in the adjacent subdued lowlands showed signals of low diversity, suggesting contraction, and subsequent expansion from localized refugia in the west. Cool summers are an important climatic variable for the species and species distribution modeling showed suitable habitat identified at the LGM suggesting expansion at this time following likely contraction during earlier warmer climatic oscillations. The isolated, coastal population at Vancouver Peninsula showed low diversity but no differentiation and it may have been established in more recent historical times, possibly through Aboriginal movement of seed. Our analysis of B. brownii highlights the complex evolutionary history of the species and the influence of topographic complexity and habitat heterogeneity in this global biodiversity hotspot.

Processes that drive the population structuring of Jenynsia lineata (Cyprinidontiformes, Anablepidae) in the La Plata Basin

The distribution of genetic diversity across a species distribution range is rarely homogeneous, as the genetic structure among populations is related to the degree of isolation among them, such as isolation by distance, isolation by barrier, and isolation by environment. Jenynsia lineata is a small viviparous fish that inhabits a wide range of habitats in South America. To decipher the isolation processes that drive population structuring in J. lineata, we analyzed 221 sequences of the mitochondrial cytochrome c oxidase I gene (COI), from 19 localities. Then, we examined the influence of the three most common types of isolation in order to explain the genetic variation found in this species.Our results revealed a marked structuration, with three groups: (a) La Plata/Desaguadero Rivers (sampling sites across Argentina, Uruguay, and Southern Brazil), (b) Central Argentina, and (c) Northern Argentina. A distance-based redundancy analysis, including the explanatory variables geographical distances, altitude, latitude, and basin, was able to explain up to 65% of the genetic structure. A variance partitioning analysis showed that the two most important variables underlying the structuration in J. lineata were altitude (isolation by environment) and type of basin (isolation by barrier).Our results show that in this species, the processes of population diversification are complex and are not limited to a single mechanism. The processes that play a prominent role in this study could explain the high rate of diversity that characterizes freshwater fish species. And these processes in turn are the basis for possible speciation events.

Modeling habitat connectivity in support of multiobjective species movement: An application to amphibian habitat systems

Reasoning about the factors underlying habitat connectivity and the inter-habitat movement of species is essential to many areas of biological inquiry. In order to better describe and understand the ways in which the landscape may support species movement, an increasing amount of research has focused on identification of paths or corridors that may be important in providing connectivity among habitat. The least-cost path problem has proven to be an instrumental analytical tool in this sense. A complicating aspect of such path identification methods is how to best reconcile and integrate the array of criteria or objectives that species may consider in traversal of a landscape. In cases where habitat connectivity is thought to be influenced or guided by multiple objectives, numerous solutions to least-cost path problems can exist, representing tradeoffs between the objectives. In practice though, identification of these solutions can be very challenging and as such, only a small proportion of them are typically examined leading to a weak characterization of habitat connectivity. To address this computational challenge, a multiobjective optimization framework is proposed. A generalizable multiobjective least-cost path model is first detailed. A non-inferior set estimation (MONISE) algorithm for identifying supported efficient solutions to the multiobjective least-cost path model is then described. However, it is well known that unsupported efficient solutions (which are equally important) can also exist, but are typically ignored given that they are more difficult to identify. Thus, to enable the identification of the full set of efficient solutions (supported and unsupported) to the multiobjective model, a multi-criteria labeling algorithm is then proposed. The developed framework is applied to assess different conceptualizations of habitat connectivity supporting amphibian movement in a wetland system. The results highlight the range of tradeoffs in characterizations of connectivity that can exist when multiple objectives are thought to contribute to movement decisions and that the number of unsupported efficient solutions (which are typically ignored) can vastly outweigh that of the supported efficient solutions.

Stepping into the past to conserve the future: Archived skin swabs from extant and extirpated populations inform genetic management of an endangered amphibian

Moving animals on a landscape through translocations and reintroductions is an important management tool used in the recovery of endangered species, particularly for the maintenance of population genetic diversity and structure. Management of imperiled amphibian species rely heavily on translocations and reintroductions, especially for species that have been brought to the brink of extinction by habitat loss, introduced species, and disease. One striking example of amphibian declines and associated management efforts is in California's Sequoia and Kings Canyon National Parks with the mountain yellow-legged frog species complex (Rana sierrae/muscosa). Mountain yellow-legged frogs have been extirpated from more than 93% of their historic range, and limited knowledge of their population genetics has made long-term conservation planning difficult. To address this, we used 598 archived skin swabs from both extant and extirpated populations across 48 lake basins to generate a robust Illumina-based nuclear amplicon data set. We found that samples grouped into three main genetic clusters, concordant with watershed boundaries. We also found evidence for historical gene flow across watershed boundaries with a north-to-south axis of migration. Finally, our results indicate that genetic diversity is not significantly different between populations with different disease histories. Our study offers specific management recommendations for imperiled mountain yellow-legged frogs and, more broadly, provides a population genetic framework for leveraging minimally invasive samples for the conservation of threatened species.

Genetic diversity in natural populations of noble crayfish (Astacus astacus L.) in north-western Poland on the basis of combined SSR and AFLP data

Background

Conservation of noble crayfish (Astacus astacus) populations is becoming particularly important since the number of individuals is rapidly declining across the distribution range of the species in Europe. Five crayfish populations in northwestern Poland have been constantly monitored for two decades. However, the genetic structure of these populations has not been analysed, although this information is important to devise effective conservation strategies.

Methods

Noble crayfish were collected in the autumn of 2014 by scuba diving in Lakes Graniczne, Babinki, Biwakowe, Sęki and Kwisno, all of which are situated in the Bytów Lakeland of northwestern Poland. Genetic diversity of the five populations was assessed based on allele variability in nine SSR regions and six AFLP primer combinations.

Results

Microsatellite results analysed with AMOVA showed that the diversity between populations corresponds to 18% of total variability, which was confirmed by similar results obtained using AFLP. Additionally, significant genetic diversity was revealed by high average F_ST values. All of the studied crayfish populations significantly deviated from the expected Hardy–Weinberg genetic equilibrium and were characterised by negative values of inbreeding coefficient (F_IS).

Discussion

The invariably negative inbreeding coefficients (F_IS) suggest a low number of mating individuals, a possible consequence of the phenomenon known as genetic bottleneck. However, additional comprehensive analyses are needed to assess the genetic structure, origin and vulnerability of the remaining populations of noble crayfish in the Bytów Lakeland of northwestern Poland, which have high conservation value and are particularly important as a live genetic bank for breeding and restitution purposes.

Genetic diversity increases with depth in red gorgonian populations of the Mediterranean Sea and the Atlantic Ocean

Background

In the ocean, the variability of environmental conditions found along depth gradients exposes populations to contrasting levels of perturbation, which can be reflected in the overall patterns of species genetic diversity. At shallow sites, resource availability may structure large, persistent and well-connected populations with higher levels of diversity. In contrast, the more extreme conditions, such as thermal stress during heat waves, can lead to population bottlenecks and genetic erosion, inverting the natural expectation. Here we examine how genetic diversity varies along depth for a long-lived, important ecosystem-structuring species, the red gorgonian, Paramuricea clavata.

Methods

We used five polymorphic microsatellite markers to infer differences in genetic diversity and differentiation, and to detect bottleneck signs between shallow and deeper populations across the Atlantic Ocean and the Mediterranean Sea. We further explored the potential relationship between depth and environmental gradients (temperature, ocean currents, productivity and slope) on the observed patterns of diversity by means of generalized linear mixed models.

Results

An overall pattern of higher genetic diversity was found in the deeper sites of the Atlantic Ocean and the Mediterranean Sea. This pattern was largely explained by bottom temperatures, with a linear pattern of decreasing genetic diversity with increasing thermal stress. Genetic differentiation patterns showed higher gene flow within sites (i.e., shallow vs. deeper populations) than between sites. Recent genetic bottlenecks were found in two populations of shallow depths.

Discussion

Our results highlight the role of deep refugial populations safeguarding higher and unique genetic diversity for marine structuring species. Theoretical regression modelling demonstrated how thermal stress alone may reduce population sizes and diversity levels of shallow water populations. In fact, the examination of time series on a daily basis showed the upper water masses repeatedly reaching lethal temperatures for P. clavata. Differentiation patterns showed that the deep richer populations are isolated. Gene flow was also inferred across different depths; however, not in sufficient levels to offset the detrimental effects of surface environmental conditions on genetic diversity. The identification of deep isolated areas with high conservation value for the red gorgonian represents an important step in the face of ongoing and future climate changes.

Speciation with gene flow: Evidence from a complex of alpine butterflies (Coenonympha, Satyridae)

Until complete reproductive isolation is achieved, the extent of differentiation between two diverging lineages is the result of a dynamic equilibrium between genetic isolation and mixing. This is especially true for hybrid taxa, for which the degree of isolation in regard to their parental species is decisive in their capacity to rise as a new and stable entity. In this work, we explored the past and current patterns of hybridization and divergence within a complex of closely related butterflies in the genus Coenonympha in which two alpine species, C. darwiniana and C. macromma, have been shown to result from hybridization between the also alpine C. gardetta and the lowland C. arcania. By testing alternative scenarios of divergence among species, we show that gene flow has been uninterrupted throughout the speciation process, although leading to different degrees of current genetic isolation between species in contact zones depending on the pair considered. Nonetheless, at broader geographic scale, analyses reveal a clear genetic differentiation between hybrid lineages and their parental species, pointing out to an advanced stage of the hybrid speciation process. Finally, the positive correlation observed between ecological divergence and genetic isolation among these butterflies suggests a potential role for ecological drivers during their speciation processes.

North-facing slopes and elevation shape asymmetric genetic structure in the range-restricted salamander Plethodon shenandoah

Species with narrow environmental tolerances are often distributed within fragmented patches of suitable habitat, and dispersal among these subpopulations can be difficult to directly observe. Genetic data can help quantify gene flow between localities, which is especially important for vulnerable species with a disjunct range. The Shenandoah salamander (Plethodon shenandoah) is a federally endangered species known only from three mountaintops in Virginia, USA. To reconstruct the evolutionary history and population connectivity of this species, we generated both mitochondrial and nuclear data using sequence capture from individuals collected across all three mountaintops. Applying population and landscape genetic methods, we found strong population structure that was independent of geographic distance. Both the nuclear markers and mitochondrial genomes indicated a deep split between the most southern population and the genetically similar central and northern populations. Although there was some mitochondrial haplotype-splitting between the central and northern populations, there was admixture in nuclear markers. This is indicative of either a recent split or current male-biased dispersal among mountain isolates. Models of landscape resistance found that dispersal across north-facing slopes at mid-elevation levels best explain the observed genetic structure among populations. These unexpected results highlight the importance of incorporating landscape features in understanding and predicting the movement and fragmentation of this range-restricted salamander species across space.

Genetic diversity and demography of the critically endangered Roberts' false brook salamander (Pseudoeurycea robertsi) in Central Mexico

Land use changes are threatening the maintenance of biodiversity. Genetic diversity is one of the main indicators of biological diversity and is highly important as it shapes the capability of populations to respond to environmental changes. We studied eleven populations of Pseudoeurycea robertsi, a micro-endemic and critically endangered species from the Nevado de Toluca Volcano, a mountain that is part of the Trans-Mexican Volcanic Belt, Mexico. We sequenced the mitochondrial cytochrome b gene from 71 individuals and genotyped 9 microsatellites from 150 individuals. Our results based on the cytochrome b showed two divergent lineages, with moderate levels of genetic diversity and a recently historical demographic expansion. Microsatellite-based results indicated low levels of heterozygosity for all populations and few alleles per locus, as compared with other mole salamander species. We identified two genetically differentiated subpopulations with a significant level of genetic structure. These results provide fundamental data for the development of management plans and conservation efforts for this critically endangered species.

Jailed in the mountains: Genetic diversity and structure of an endemic newt species across the Pyrenees

The Pyrenees represent a natural laboratory for biogeographic, evolutionary and ecological research of mountain fauna as a result of the high variety of habitats and the profound effect of the glacial and interglacial periods. There is a paucity of studies providing a detailed insight into genetic processes and better knowledge on the patterns of genetic diversity and how they are maintained under high altitude conditions. This is of particular interest when considering the course of past climate conditions and glaciations in a species which is considered site tenacious, with long generation times. Here we analyzed the genetic patterns of diversity and structure of the endemic Pyrenean brook newt (Calotriton asper) along its distribution range, with special emphasis on the distinct habitat types (caves, streams, and lakes), and the altitudinal and geographical ranges, using a total set of 900 individuals from 44 different localities across the Pyrenean mountain range genotyped for 19 microsatellite loci. We found evidence for a negative longitudinal and positive altitudinal gradient of genetic diversity in C. asper populations. The fact that genetic diversity was markedly higher westwards is in accordance with other Pyrenean species. However, the impact of altitudinal gradient on the genetic diversity seems to differ from other species, and mostly from other amphibians. We found that lower altitudes can act as a barrier probably because the lowlands do not provide a suitable habitat for C. asper. Regarding the distinct habitat types, caves had significantly lower values of genetic diversity compared to streams or lakes. The mean F_ST value was relatively high (0.304) with maximum values as high as 0.771, suggesting a highly structured total population. Indeed, populations were grouped into five subclusters, the eastern populations (cluster 1) remained grouped into two subclusters and the central-western Pyrenees (cluster 2) into three subclusters. The increase of isolation with geographical distance is consistent with the population structure detected. In conclusion, C. asper seems to be adapted to high altitude mountain habitats, and its genetic diversity is higher in the western Pyrenees. In terms of conservation priority, we consider more relevant the populations that represent a reservoir of genetic diversity.

The Tanggula Mountains enhance population divergence in Carex moorcroftii: a dominant sedge on the Qinghai-Tibetan Plateau

High-altitude mountains are often geographic barriers to gene flow and play important roles in shaping population divergence. The central Qinghai-Tibetan Plateau (QTP) stands the location of the Tanggula Mountains (TM). We use the TM as a case, using Carex moorcroftii, a dominant species on the QTP to test the effects of geographic barriers on plant population divergence. We sampled 18 C. moorcroftii populations along a north-south transect crossing the TM to investigate the correlations of genetic variation and morphological traits with climate variables. The results showed this species holds high genetic diversity (H_e = 0.58) and the surveyed populations can be genetically clustered into two groups: populations from the north face of TM, and the other from the south. Gene flow between populations within groups is higher than those between groups. The traits, number and mass of seeds, mass of root and infructescence significantly varied among populations. Mantel-tests detected a weak but significantly positive correlation between genetic and geographic (R² = 0.107, p = 0.032) and climatic distance (R² = 0.162, p = 0.005), indicating both isolation by distance and isolation by environment. These findings together suggest high-altitude mountains of TM interrupt habitat continuity, result in distinct climatic conditions on both sides, increasing population divergence of plant species.

Genetic variability and structure of an isolated population of Ambystoma altamirani, a mole salamander that lives in the mountains of one of the largest urban areas in the world

Amphibians are globally threatened by habitat loss and fragmentation; species within the order Ambystoma are not the exception, as there are 18 species of mole salamanders in México, of which 16 are endemic and all species are under some national or international status of protection. The mole salamander, Ambystoma altamirani is a microendemic species, which is distributed in central México, within the trans-Mexican volcanic belt, and is one of the most threatened species due to habitat destruction and the introduction of exotic species. Nine microsatellite markers were used to determine the genetic structure, genetic variability, effective population size, presence of bottlenecks and inbreeding coefficient of one population of A. altamirani to generate information which might help to protect and conserve this threatened species. We found two genetic subpopulations with significant level of genetic structure (F_ST = 0.005) and high levels of genetic variability (H_o = 0.883; H_e = 0.621); we also found a small population size (N_e = 8.8), the presence of historical (M = 0.486) and recent bottlenecks under IAM and TPM models, with a low, but significant coefficient of inbreeding (F_IS = -0.451). This information will help us to raise conservation strategies of this microendemic mole salamander species.

Distribution and space use of seed-dispersing rodents in central Pyrenees: implications for genetic diversity, conservation and plant recruitment

The function and conservation of many forest ecosystems depend on the distribution and diversity of the community of rodents that consume and disperse seeds. The habitat preferences and interactions are especially relevant in alpine systems where such granivorous rodents reach the southernmost limit of their distribution and are especially sensitive to global warming. We analyzed the community of granivorous rodents in the Pyrenees, one of the southernmost mountain ranges of Europe. Rodent species were identified by DNA with particular attention to the Apodemus species, which are prominent seed-dispersing rodents in Europe. We confirmed for the first time the presence of the yellow-necked mouse, Apodemus flavicollis, in central Pyrenees, a typical Eurosiberian species that reaches its southernmost distribution limit in this area. We also found the wood mouse, Apodemus sylvaticus, a related species more tolerant to Mediterranean environments. Both rodents were spatially segregated by altitude. A. sylvaticus was rare at high altitudes, which might cause the genetic differentiation between populations of the different valleys reported here. We also found other seed consumers like dormice, Elyomis quercinus, and voles, Myodes glareolus, with marked habitat preferences. We suggest that population isolation among valleys may increase the genetic diversity of rodents, like A. sylvaticus. We also highlight the potential threat that global warming may represent for species linked to high-altitude refuges at the southern edge of its distribution, like Apodemus flavicollis. Finally, we discuss how this threat may have a dimension in the conservation of alpine forests dispersed by these rodent populations.

Comparative landscape genetics of pond-breeding amphibians in Mediterranean temporal wetlands: The positive role of structural heterogeneity in promoting gene flow

Comparative landscape genetics studies can provide key information to implement cost-effective conservation measures favouring a broad set of taxa. These studies are scarce, particularly in Mediterranean areas, which include diverse but threatened biological communities. Here, we focus on Mediterranean wetlands in central Iberia and perform a multi-level, comparative study of two endemic pond-breeding amphibians, a salamander (Pleurodeles waltl) and a toad (Pelobates cultripes). We genotyped 411 salamanders from 20 populations and 306 toads from 16 populations at 18 and 16 microsatellite loci, respectively, and identified major factors associated with population connectivity through the analysis of three sets of variables potentially affecting gene flow at increasingly finer levels of spatial resolution. Topographic, land use/cover, and remotely sensed vegetation/moisture indices were used to derive optimized resistance surfaces for the two species. We found contrasting patterns of genetic structure, with stronger, finer scale genetic differentiation in Pleurodeles waltl, and notable differences in the role of fine-scale patterns of heterogeneity in vegetation cover and water content in shaping patterns of regional genetic structure in the two species. Overall, our results suggest a positive role of structural heterogeneity in population connectivity in pond-breeding amphibians, with habitat patches of Mediterranean scrubland and open oak woodlands ("dehesas") facilitating gene flow. Our study highlights the usefulness of remotely sensed continuous variables of land cover, vegetation and water content (e.g., NDVI, NDMI) in conservation-oriented studies aimed at identifying major drivers of population connectivity.

Altitudinal gradients, biogeographic history and microhabitat adaptation affect fine-scale spatial genetic structure in African and Neotropical populations of an ancient tropical tree species

The analysis of fine-scale spatial genetic structure (FSGS) within populations can provide insights into eco-evolutionary processes. Restricted dispersal and locally occurring genetic drift are the primary causes for FSGS at equilibrium, as described in the isolation by distance (IBD) model. Beyond IBD expectations, spatial, environmental or historical factors can affect FSGS. We examined FSGS in seven African and Neotropical populations of the late-successional rain forest tree Symphonia globulifera L. f. (Clusiaceae) to discriminate the influence of drift-dispersal vs. landscape/ecological features and historical processes on FSGS. We used spatial principal component analysis and Bayesian clustering to assess spatial genetic heterogeneity at SSRs and examined its association with plastid DNA and habitat features. African populations (from Cameroon and São Tomé) displayed a stronger FSGS than Neotropical populations at both marker types (mean Sp = 0.025 vs. Sp = 0.008 at SSRs) and had a stronger spatial genetic heterogeneity. All three African populations occurred in pronounced altitudinal gradients, possibly restricting animal-mediated seed dispersal. Cyto-nuclear disequilibria in Cameroonian populations also suggested a legacy of biogeographic history to explain these genetic patterns. Conversely, Neotropical populations exhibited a weaker FSGS, which may reflect more efficient wide-ranging seed dispersal by Neotropical bats and other dispersers. The population from French Guiana displayed an association of plastid haplotypes with two morphotypes characterized by differential habitat preferences. Our results highlight the importance of the microenvironment for eco-evolutionary processes within persistent tropical tree populations.

Catchments catch all in South African coastal lowlands: topography and palaeoclimate restricted gene flow in Nymania capensis (Meliaceae)-a multilocus phylogeographic and distribution modelling approach

BACKGROUND

This study investigates orbitally-forced range dynamics at a regional scale by exploring the evolutionary history of Nymania capensis (Meliaceae) across the deeply incised landscapes of the subescarpment coastal lowlands of South Africa; a region that is home to three biodiversity hotspots (Succulent Karoo, Fynbos, and Maputaland-Pondoland-Albany hotspots).

METHODS

A range of methods are used including: multilocus phylogeography (chloroplast and high- and low-copy nuclear DNA), molecular dating and species distribution modelling (SDM).

RESULTS

The results support an 'evolutionarily distinct catchment' hypothesis where: (1) different catchments contain genetically distinct lineages, (2) limited genetic structuring was detected within basins whilst high structuring was detected between basins, and (3) within primary catchment populations display a high degree of genealogical lineage sorting. In addition, the results support a glacial refugia hypothesis as: (a) the timing of chloroplast lineage diversification is restricted to the Pleistocene in a landscape that has been relatively unchanged since the late Pliocene, and (b) the projected LGM distribution of suitable climate for N. capensis suggest fragmentation into refugia that correspond to the current phylogeographic populations.

DISCUSSION

This study highlights the interaction of topography and subtle Pleistocene climate variations as drivers limiting both seed and pollen flow along these lowlands. This lends support to the region's large-scale conservation planning efforts, which used catchments as foundational units for conservation as these are likely to be evolutionarily significant units.

Population Structure and Evolution after Speciation of the Hokkaido Salamander (Hynobius retardatus)

The Hokkaido salamander (Hynobius retardatus) is endemic to Hokkaido Island, Japan, and shows intriguing flexible phenotypic plasticity and regional morphological diversity. However, to date, allozymes and partial mitochondria DNA sequences have provided only an outline of its demographic histories and the pattern of its genetic diversification. To understand the finer details of the population structure of this species and its evolution since speciation, we genotyped five regional populations by using 12 recently developed microsatellite polymorphic markers. We found a clear population structure with low gene flow among the five populations, but a close genetic relationship between the Teshio and Kitami populations. Our demographic analysis suggested that Teshio and Erimo had the largest effective population sizes among the five populations. These findings regarding the population structure and demography of H. retardatus improve our understanding of the faunal phylogeography on Hokkaido Island and also provide fundamental genetic information that will be useful for future studies.

Contrasting evolutionary histories of the legless lizards slow worms (Anguis) shaped by the topography of the Balkan Peninsula

Background

Genetic architecture of a species is a result of historical changes in population size and extent of distribution related to climatic and environmental factors and contemporary processes of dispersal and gene flow. Population-size and range contractions, expansions and shifts have a substantial effect on genetic diversity and intraspecific divergence, which is further shaped by gene-flow limiting barriers. The Balkans, as one of the most important sources of European biodiversity, is a region where many temperate species persisted during the Pleistocene glaciations and where high topographic heterogeneity offers suitable conditions for local adaptations of populations. In this study, we investigated the phylogeographical patterns and demographic histories of four species of semifossorial slow-worm lizards (genus Anguis) present in the Balkan Peninsula, and tested the relationship between genetic diversity and topographic heterogeneity of the inhabited ranges.

Results

We inferred phylogenetic relationships, compared genetic structure and historical demography of slow worms using nucleotide sequence variation of mitochondrial DNA. Four Anguis species with mostly parapatric distributions occur in the Balkan Peninsula. They show different levels of genetic diversity. A signature of population growth was detected in all four species but with various courses in particular populations. We found a strong correlation between genetic diversity of slow-worm populations and topographic ruggedness of the ranges (mountain systems) they inhabit. Areas with more rugged terrain harbour higher genetic diversity.

Conclusions

Phylogeographical pattern of the genus Anguis in the Balkans is concordant with the refugia-within-refugia model previously proposed for both several other taxa in the region and other main European Peninsulas. While slow-worm populations from the southern refugia mostly have restricted distributions and have not dispersed much from their refugial areas, populations from the extra-Mediterranean refugia in northern parts of the Balkans have colonized vast areas of eastern, central, and western Europe. Besides climatic historical events, the heterogeneous topography of the Balkans has also played an important role in shaping genetic diversity of slow worms.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0669-1) contains supplementary material, which is available to authorized users.

Amphibian molecular ecology and how it has informed conservation

Molecular ecology has become one of the key tools in the modern conservationist's kit. Here we review three areas where molecular ecology has been applied to amphibian conservation: genes on landscapes, within-population processes, and genes that matter. We summarize relevant analytical methods, recent important studies from the amphibian literature, and conservation implications for each section. Finally, we include five in-depth examples of how molecular ecology has been successfully applied to specific amphibian systems.

How spatio-temporal habitat connectivity affects amphibian genetic structure

Heterogeneous landscapes and fluctuating environmental conditions can affect species dispersal, population genetics, and genetic structure, yet understanding how biotic and abiotic factors affect population dynamics in a fluctuating environment is critical for species management. We evaluated how spatio-temporal habitat connectivity influences dispersal and genetic structure in a population of boreal chorus frogs (Pseudacris maculata) using a landscape genetics approach. We developed gravity models to assess the contribution of various factors to the observed genetic distance as a measure of functional connectivity. We selected (a) wetland (within-site) and (b) landscape matrix (between-site) characteristics; and (c) wetland connectivity metrics using a unique methodology. Specifically, we developed three networks that quantify wetland connectivity based on: (i) P. maculata dispersal ability, (ii) temporal variation in wetland quality, and (iii) contribution of wetland stepping-stones to frog dispersal. We examined 18 wetlands in Colorado, and quantified 12 microsatellite loci from 322 individual frogs. We found that genetic connectivity was related to topographic complexity, within- and between-wetland differences in moisture, and wetland functional connectivity as contributed by stepping-stone wetlands. Our results highlight the role that dynamic environmental factors have on dispersal-limited species and illustrate how complex asynchronous interactions contribute to the structure of spatially-explicit metapopulations.

Altitudinal Barrier to the Spread of an Invasive Species: Could the Pyrenean Chain Slow the Natural Spread of the Pinewood Nematode?

Mountain ranges may delimit the distribution of native species as well as constitute potential barriers to the spread of invasive species. The invasive pinewood nematode, Bursaphelenchus xylophilus, is a severe forest pest inducing pine wilt disease. It is vectored in Europe by a native long-horned beetle, Monochamus galloprovincialis. This study explored the potential of the Pyrenean chain to slow or prevent the natural spread of nematode-infested beetles from the Iberian Peninsula, where the nematode is established and is expanding its range, towards France and the rest of Europe. An analysis of the genetic structure and migration patterns of the beetle populations throughout the Pyrenean mountain range was combined with a spread model simulating the potential movements of nematode-infested beetles across it. The central part of the Pyrenees, which corresponds to the highest elevation zone, was shown to prevent gene flow between the French and Spanish populations of M. galloprovincialis on each side of the mountains. Conversely, strong admixture was detected between populations located on both sides of low elevation hills, and especially at the east and west extremities of the mountain range. Simulations of the spread of nematode-infested beetles under various thresholds of beetle survival and pine wilt disease expression gave results consistent with the variation in genetic make-up, suggesting that western and eastern hillsides may represent corridors favoring natural spread of the nematode from the Iberian Peninsula to France. Simulations also showed that temperature rise due to climate change may significantly reduce the extent of the barrier formed by highest elevations. Our results support the hypothesis that the Pyrenean chain represents a partial barrier to the natural spread of nematode-infested beetles. These results, which have to be considered together with potential human-assisted long-distance spread of the nematode, highlight priority zones for future pest monitoring and management programs. More generally, such an integrated approach could be used to assess the role of mountain chains in the potential spread of other invasive pests.

Conflicting effects of microhabitats on Long-toed Salamander (Ambystoma macrodactylum) movement: implications for landscape connectivity

Understanding dispersal requires multiple lines of investigation, from the study of broad patterns of population connectivity to the identification of factors impacting movement at local scales. To determine the potential effects of different microhabitats on dispersal in the Long-toed Salamander (Ambystoma macrodactylum Baird, 1850), we experimentally evaluated mobility, moisture loss, and habitat choice in response to five common substrates (deciduous and coniferous leaf litter, grass, moss, and sand). Specifically, we examined differences in the efficiency with which salamanders moved across substrates when motivated to move. We then quantified moisture loss in each substrate and evaluated habitat preference. Our results point to a trade-off between substrates that are easily traversed and those that offer high protection against desiccation. Habitat choice appeared to balance these two aspects of performance, with salamanders favouring a substrate that offered both low resistance to movement and high protection against desiccation. This result was context-dependent, as preferences shifted towards wetter but less easily traversed substrates when supplemental cover objects were made available. Overall, our study highlights the potential for individuals to respond to a given substrate in ways that can both facilitate and limit dispersal and thus underscores the need to consider different aspects of individual performance and behaviour when studying population connectivity.

AFLPs and mitochondrial haplotypes reveal local adaptation to extreme thermal environments in a freshwater gastropod

The way environmental variation shapes neutral and adaptive genetic variation in natural populations is a key issue in evolutionary biology. Genome scans allow the identification of the genetic basis of local adaptation without previous knowledge of genetic variation or traits under selection. Candidate loci for divergent adaptation are expected to show higher F_ST than neutral loci influenced solely by random genetic drift, migration and mutation. The comparison of spatial patterns of neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection among populations living in contrasting environments. Using the gastropod Radix balthica as a system, we analyzed 376 AFLP markers and 25 mtDNA COI haplotypes for candidate loci and associations with local adaptation among contrasting thermal environments in Lake Mývatn, a volcanic lake in northern Iceland. We found that 2% of the analysed AFLP markers were under directional selection and 12% of the mitochondrial haplotypes correlated with differing thermal habitats. The genetic networks were concordant for AFLP markers and mitochondrial haplotypes, depicting distinct topologies at neutral and candidate loci. Neutral topologies were characterized by intense gene flow revealed by dense nets with edges connecting contrasting thermal habitats, whereas the connections at candidate loci were mostly restricted to populations within each thermal habitat and the number of edges decreased with temperature. Our results suggest microgeographic adaptation within Lake Mývatn and highlight the utility of genome scans in detecting adaptive divergence.

Landscape influences on dispersal behaviour: a theoretical model and empirical test using the fire salamander, Salamandra infraimmaculata

When populations reside within a heterogeneous landscape, isolation by distance may not be a good predictor of genetic divergence if dispersal behaviour and therefore gene flow depend on landscape features. Commonly used approaches linking landscape features to gene flow include the least cost path (LCP), random walk (RW), and isolation by resistance (IBR) models. However, none of these models is likely to be the most appropriate for all species and in all environments. We compared the performance of LCP, RW and IBR models of dispersal with the aid of simulations conducted on artificially generated landscapes. We also applied each model to empirical data on the landscape genetics of the endangered fire salamander, Salamandra infraimmaculata, in northern Israel, where conservation planning requires an understanding of the dispersal corridors. Our simulations demonstrate that wide dispersal corridors of the low-cost environment facilitate dispersal in the IBR model, but inhibit dispersal in the RW model. In our empirical study, IBR explained the genetic divergence better than the LCP and RW models (partial Mantel correlation 0.413 for IBR, compared to 0.212 for LCP, and 0.340 for RW). Overall dispersal cost in salamanders was also well predicted by landscape feature slope steepness (76%), and elevation (24%). We conclude that fire salamander dispersal is well characterised by IBR predictions. Together with our simulation findings, these results indicate that wide dispersal corridors facilitate, rather than hinder, salamander dispersal. Comparison of genetic data to dispersal model outputs can be a useful technique in inferring dispersal behaviour from population genetic data.

Geographic determinants of gene flow in two sister species of tropical Andean frogs

Complex interactions between topographic heterogeneity, climatic and environmental gradients, and thermal niche conservatism are commonly assumed to indicate the degree of biotic diversification in montane regions. Our aim was to investigate factors that disrupt gene flow between populations and to determine if there is evidence of downslope asymmetric migration in highland frogs with wide elevational ranges and thermal niches. We determined the role of putative impediments to gene flow (as measured by least-cost path (LCP) distances, topographic complexity, and elevational range) in promoting genetic divergence between populations of 2 tropical Andean frog sister species (Dendropsophus luddeckei, N = 114; Dendropsophus labialis, N = 74) using causal modeling and multiple matrix regression. Although the effect of geographic features was species specific, elevational range and LCP distances had the strongest effect on gene flow, with mean effect sizes (Mantel r and regression coefficients β), between 5 and 10 times greater than topographic complexity. Even though causal modeling and multiple matrix regression produced congruent results, the latter provided more information on the contribution of each geographic variable. We found moderate support for downslope migration. We conclude that the climatic heterogeneity of the landscape, the elevational distance between populations, and the inability to colonize suboptimal habitats due to thermal niche conservatism influence the magnitude of gene flow. Asymmetric migration, however, seems to be influenced by life history traits.