Effects of slope and riparian habitat connectivity on gene flow in an endangered Panamanian frog,Atelopus varius

Contrasting environmental drivers of genetic and phenotypic divergence in an Andean poison frog (Epipedobates anthonyi)

Phenotypic and genetic divergence are shaped by the homogenizing effects of gene flow and the differentiating processes of genetic drift and local adaptation. Herein, we examined the mechanisms that underlie phenotypic (size and color) and genetic divergence in 35 populations (535 individuals) of the poison frog Epipedobates anthonyi along four elevational gradients (0-1800 m asl) in the Ecuadorian Andes. We found phenotypic divergence in size and color despite relatively low genetic divergence at neutral microsatellite loci. Genetic and phenotypic divergence were both explained by landscape resistance between sites (isolation-by-resistance, IBR), likely due to a cold and dry mountain ridge between the northern and southern elevational transects that limits dispersal and separates two color morphs. Moreover, environmental differences among sites also explained genetic and phenotypic divergence, suggesting isolation-by-environment (IBE). When northern and southern transects were analyzed separately, genetic divergence was predicted either by distance (isolation-by-distance, IBD; northern) or environmental resistance between sites (IBR; southern). In contrast, phenotypic divergence was primarily explained by environmental differences among sites, supporting the IBE hypothesis. These results indicate that although distance and geographic barriers are important drivers of population divergence, environmental variation has a two-fold effect on population divergence. On the one hand, landscape resistance between sites reduces gene flow (IBR), while on the other hand, environmental differences among sites exert divergent selective pressures on phenotypic traits (IBE). Our work highlights the importance of studying both genetic and phenotypic divergence to better understand the processes of population divergence and speciation along ecological gradients.

Divergent regional evolutionary histories of a devastating global amphibian pathogen

Emerging infectious diseases are a pressing threat to global biological diversity. Increased incidence and severity of novel pathogens underscores the need for methodological advances to understand pathogen emergence and spread. Here, we use genetic epidemiology to test, and challenge, key hypotheses about a devastating zoonotic disease impacting amphibians globally. Using an amplicon-based sequencing method and non-invasive samples we retrospectively explore the history of the fungal pathogen Batrachochytrium dendrobatidis (Bd) in two emblematic amphibian systems: the Sierra Nevada of California and Central Panama. The hypothesis in both regions is the hypervirulent Global Panzootic Lineage of Bd (BdGPL) was recently introduced and spread rapidly in a wave-like pattern. Our data challenge this hypothesis by demonstrating similar epizootic signatures can have radically different underlying evolutionary histories. In Central Panama, our genetic data confirm a recent and rapid pathogen spread. However, BdGPL in the Sierra Nevada has remarkable spatial structuring, high genetic diversity and a relatively older history inferred from time-dated phylogenies. Thus, this deadly pathogen lineage may have a longer history in some regions than assumed, providing insights into its origin and spread. Overall, our results highlight the importance of integrating observed wildlife die-offs with genetic data to more accurately reconstruct pathogen outbreaks.

Comparative analysis of spatial genetic structures in sympatric populations of two riparian plants, Saxifraga acerifolia and Saxifraga fortunei

PREMISE

The genetic structure between plant populations is facilitated by the spatial population arrangement and limited dispersal of seed and pollen. Saxifraga acerifolia, a local endemic species in Japan, is a habitat specialist that is confined to waterfalls in riparian environments. Its sister species, Saxifraga fortunei, is a generalist that is widely distributed along riverbanks. Here, we examined sympatric populations of the two Saxifraga species to test whether the differences in habitat preference and colonization process influenced regional and local genetic structures.

METHODS

To reveal genetic structures, we examined chloroplast microsatellite variations and genome-wide nucleotide polymorphisms obtained by genotyping by sequencing. We also estimated the gene flow among and within populations and performed landscape genetic analyses to evaluate seed and pollen movement and the extent of genetic isolation related to geographic distance and/or habitat differences.

RESULTS

We found strong genetic structure in the specialist S. acerifolia, even on a small spatial scale (<1 km part); each population on a different waterfall in one river system had a completely different predominant haplotype. By contrast, the generalist S. fortunei showed no clear genetic differentiation.

CONCLUSIONS

Our findings suggest that the level of genetic isolation was increased in S. acerifolia by the spatially fragmented habitat and limited seed and pollen dispersal over waterfalls. Habitat differentiation between the sister taxa could have contributed to the different patterns of gene flow and then shaped the contrasting genetic structures.

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.

Landscape genetics of northern crested newt Triturus cristatus populations in a contrasting natural and human-impacted boreal forest

Among vertebrates, amphibians currently have the highest proportion of threatened species worldwide, mainly through loss of habitat, leading to increased population isolation. Smaller amphibian populations may lose more genetic diversity, and become more dependent on immigration for survival. Investigations of landscape factors and patterns mediating migration and population genetic differentiation are fundamental for knowledge-based conservation. The pond-breeding northern crested newt (Triturus cristatus) populations are decreasing throughout Europe, and are a conservation concern. Using microsatellites, we studied the genetic structure of the northern crested newt in a boreal forest ecosystem containing two contrasting landscapes, one subject to recent change and habitat loss by clear-cutting and roadbuilding, and one with little anthropogenic disturbance. Newts from 12 breeding ponds were analyzed for 13 microsatellites and 7 landscape and spatial variables. With a Maximum-likelihood population-effects model we investigated important landscape factors potentially explaining genetic patterns. Results indicate that intervening landscape factors between breeding ponds, explain the genetic differentiation in addition to an isolation-by-distance effect. Geographic distance, gravel roads, and south/south-west facing slopes reduced landscape permeability and increased genetic differentiation for these newts. The effect was opposite for streams, presumably being more favorable for newt dispersal. Populations within or bordering on old growth forest had a higher allelic richness than populations in managed forest outside these areas. Old growth forest areas may be important source habitats in the conservation of northern crested newt populations.

Amphibian fungal panzootic causes catastrophic and ongoing loss of biodiversity

Anthropogenic trade and development have broken down dispersal barriers, facilitating the spread of diseases that threaten Earth's biodiversity. We present a global, quantitative assessment of the amphibian chytridiomycosis panzootic, one of the most impactful examples of disease spread, and demonstrate its role in the decline of at least 501 amphibian species over the past half-century, including 90 presumed extinctions. The effects of chytridiomycosis have been greatest in large-bodied, range-restricted anurans in wet climates in the Americas and Australia. Declines peaked in the 1980s, and only 12% of declined species show signs of recovery, whereas 39% are experiencing ongoing decline. There is risk of further chytridiomycosis outbreaks in new areas. The chytridiomycosis panzootic represents the greatest recorded loss of biodiversity attributable to a disease.

A landscape genetic analysis of swamp rabbits (Sylvilagus aquaticus) suggests forest canopy cover enhances gene flow in an agricultural matrix

Habitat loss and fragmentation pose a continued and immediate threat to wildlife and create a persistent need for ecological information at the landscape scale to guide conservation efforts. Landscape features influence population connectivity for many species and genetic analyses can be employed to determine which of these features are most important. Because population connectivity through dispersal is important to the persistence of swamp rabbits (Sylvilagus aquaticus (Bachman, 1837)) at the northern edge of their range, we used a landscape genetic approach to relate gene flow to landscape features that may impact dispersal success. We tested resistance values for attributes of land cover, watercourse corridors, canopy cover, and roads and used causal modeling and redundancy analysis to relate these representations of landscapes to genetic distance for swamp rabbits in southern Illinois, USA. Models that included canopy cover had the strongest correlations with genetic distance and were supported by our methods whereas other models were not. We concluded that high tree canopy cover enhances gene flow and landscape connectivity for swamp rabbits in southern Illinois. Our study provides important empirical evidence that landscape variables may impact the habitat connectivity of swamp rabbits. Preserving dispersal routes for swamp rabbits should focus on improving canopy cover, in both bottomland and upland, to connect suitable habitat.

Inferring spatial patterns and drivers of population divergence of Neolitsea sericea (Lauraceae), based on molecular phylogeography and landscape genomics

The relative roles of geography, climate and ecology in driving population divergence and (incipient) speciation has so far been largely neglected in studies addressing the evolution of East Asia's island flora. Here, we employed chloroplast and ribosomal DNA sequences and restriction site-associated DNA sequencing (RADseq) loci to investigate the phylogeography and drivers of population divergence of Neolitsea sericea. These data sets support the subdivision of N. sericea populations into the Southern and Northern lineages across the 'Tokara gap'. Two distinct sublineages were further identified for the Northern lineage of N. sericea from the RADseq data. RADseq was also used along with approximate Bayesian computation to show that the current distribution and differentiation of N. sericea populations resulted from a combination of relatively ancient migration and successive vicariant events that likely occurred during the mid to late Pleistocene. Landscape genomic analyses showed that, apart from geographic barriers, barrier, potentially local adaptation to different climatic conditions appears to be one of the major drivers for lineage diversification of N. sericea.

Comparative Phylogeography of Ethiopian anurans: impact of the Great Rift Valley and Pleistocene climate change

Background

The Ethiopian highlands are a biodiversity hotspot, split by the Great Rift Valley into two distinct systems of plateaus and mountains. The Rift Valley is currently hot and dry and acts as a barrier to gene flow for highland-adapted species. It is however unlikely that the conditions in the Rift were inhospitable to highland species during the entire Pleistocene. To assess the significance of the Ethiopian Rift as a biogeographic barrier as well as the impact Pleistocene climatic changes have had on the evolution of Ethiopian organisms, we performed phylogeographic analyses and developed present and past niche models on seven anuran species with different elevational and ecological preferences.

Results

We found that highland species on the east and the west sides of the Rift are genetically differentiated and have not experienced any detectable gene flow for at least 0.4 my. In contrast, species found at elevations lower than 2500 m do not show any population structure. We also determined that highland species have lower effective population sizes than lowland species, which have experienced a large, yet gradual, demographic expansion, starting approximately half a million year ago.

Conclusions

The pattern we report here is consistent with the increasingly warmer and drier conditions of the Pleistocene in East Africa, which resulted in the expansion of savanna, the fragmentation of forests and the shrinking of highland habitats. Climatic niche models indicated that the Rift is currently non suitable for most of the studied species, but it could have been a more permeable barrier during the Last Glacial Maximum. However, considering the strong genetic structure of highland species, we hypothesize that the barrier mechanisms at the Rift are not only climatic but also topographical.

Electronic supplementary material

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

Defining landscape resistance values in least-cost connectivity models for the invasive grey squirrel: a comparison of approaches using expert-opinion and habitat suitability modelling

Least-cost models are widely used to study the functional connectivity of habitat within a varied landscape matrix. A critical step in the process is identifying resistance values for each land cover based upon the facilitating or impeding impact on species movement. Ideally resistance values would be parameterised with empirical data, but due to a shortage of such information, expert-opinion is often used. However, the use of expert-opinion is seen as subjective, human-centric and unreliable. This study derived resistance values from grey squirrel habitat suitability models (HSM) in order to compare the utility and validity of this approach with more traditional, expert-led methods. Models were built and tested with MaxEnt, using squirrel presence records and a categorical land cover map for Cumbria, UK. Predictions on the likelihood of squirrel occurrence within each land cover type were inverted, providing resistance values which were used to parameterise a least-cost model. The resulting habitat networks were measured and compared to those derived from a least-cost model built with previously collated information from experts. The expert-derived and HSM-inferred least-cost networks differ in precision. The HSM-informed networks were smaller and more fragmented because of the higher resistance values attributed to most habitats. These results are discussed in relation to the applicability of both approaches for conservation and management objectives, providing guidance to researchers and practitioners attempting to apply and interpret a least-cost approach to mapping ecological networks.

Reproductive isolation related to mimetic divergence in the poison frog Ranitomeya imitator

In a mimetic radiation--when a single species evolves to resemble different model species--mimicry can drive within-species morphological diversification, and, potentially, speciation. While mimetic radiations have occurred in a variety of taxa, their role in speciation remains poorly understood. We study the Peruvian poison frog Ranitomeya imitator, a species that has undergone a mimetic radiation into four distinct morphs. Using a combination of colour-pattern analysis, landscape genetics and mate-choice experiments, we show that a mimetic shift in R. imitator is associated with a narrow phenotypic transition zone, neutral genetic divergence and assortative mating, suggesting that divergent selection to resemble different model species has led to a breakdown in gene flow between these two populations. These results extend the effects of mimicry on speciation into a vertebrate system and characterize an early stage of speciation where reproductive isolation between mimetic morphs is incomplete but evident.

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.

Population structure and landscape genetics of two endangered frog species of genus Odorrana: different scenarios on two islands

Isolation by distance and landscape connectivity are fundamental factors underlying speciation and evolution. To understand how landscapes affect gene flow and shape population structures, island species provide intrinsic study objects. We investigated the effects of landscapes on the population structure of the endangered frog species, Odorrana ishikawae and O. splendida, which each inhabit an island in southwest Japan. This was done by examining population structure, gene flow and demographic history of each species by analyzing 12 microsatellite loci and exploring causal environmental factors through ecological niche modeling (ENM) and the cost-distance approach. Our results revealed that the limited gene flow and multiple-population structure in O. splendida and the single-population structure in O. ishikawae were maintained after divergence of the species through ancient vicariance between islands. We found that genetic distance correlated with geographic distance between populations of both species. Our landscape genetic analysis revealed that the connectivity of suitable habitats influences gene flow and leads to the formation of specific population structures. In particular, different degrees of topographical complexity between islands are the major determining factor for shaping contrasting population structures of two species. In conclusion, our results illustrate the diversification mechanism of organisms through the interaction with space and environment. Our results also present an ENM approach for identifying the key factors affecting demographic history and population structures of target species, especially endangered species.

Differential effects of landscape-level environmental features on genetic structure in three codistributed tree species in Central America

Landscape genetic studies use spatially explicit population genetic information to determine the physical and environmental causes of population genetic structure on regional scales. Comparative studies that identify common barriers to gene flow across multiple species within a community are important to both understand the evolutionary trajectories of populations and prioritize habitat conservation. Here, we use a comparative landscape genetic approach to ask whether gradients in temperature or precipitation seasonality structure genetic variation across three codistributed tree species in Central America, or whether a simpler (geographic distance) or more complex, species-specific environmental niche model is necessary to individually explain population genetic structure. Using descriptive statistics and causal modelling, we find that different factors best explain genetic distance in each of the three species: environmental niche distance in Bursera simaruba, geographic distance in Ficus insipida and historical barriers to gene flow or cryptic reproductive barriers for Brosimum alicastrum. This study confirms suggestions from previous studies of Central American tree species that imply that population genetic structure of trees in this region is determined by complex interactions of both historical and current barriers to gene flow.

Influence of environmental heterogeneity on genetic diversity and structure in an endemic southern Californian oak

Understanding how specific environmental factors shape gene flow while disentangling their importance relative to the effects of geographical isolation is a major question in evolutionary biology and a specific goal of landscape genetics. Here, we combine information from nuclear microsatellite markers and ecological niche modelling to study the association between climate and spatial genetic structure and variability in Engelmann oak (Quercus engelmannii), a wind-pollinated species with high potential for gene flow. We first test whether genetic diversity is associated with climatic niche suitability and stability since the Last Glacial Maximum (LGM). Second, we use causal modelling to analyse the potential influence of climatic factors (current and LGM niche suitability) and altitude in the observed patterns of genetic structure. We found that genetic diversity is negatively associated with local climatic stability since the LGM, which may be due to higher immigration rates in unstable patches during favourable climatic periods and/or temporally varying selection. Analyses of spatial genetic structure revealed the presence of three main genetic clusters, a pattern that is mainly driven by two highly differentiated populations located in the northern edge of the species distribution range. After controlling for geographic distance, causal modelling analyses showed that genetic relatedness decreases with the environmental divergence among sampling sites estimated as altitude and current and LGM niche suitability. Natural selection against nonlocal genotypes and/or asynchrony in reproductive phenology may explain this pattern. Overall, this study suggests that local environmental conditions can shape patterns of genetic structure and variability even in species with high potential for gene flow and relatively small distribution ranges.

Rapid selection against inbreeding in a wild population of a rare frog

Populations that are small and isolated can be threatened through loss of fitness due to inbreeding. Nevertheless, an increased frequency of recessive homozygotes could increase the efficiency of selection against deleterious mutants, thus reducing inbreeding depression. In wild populations, observations of evolutionary changes determined by selection against inbreeding are few. We used microsatellite DNA markers to compare the genetic features of tadpoles immediately after hatch with those of metamorphosing froglets belonging to the same cohort in a small, isolated population of the threatened frog Rana latastei. Within a generation, the inbreeding coefficient (F_IS) decreased: at hatch, F_IS was significantly >0, whereas F_IS was <0 after metamorphosis. Furthermore, heterozygosity increased and allelic frequencies changed over time, resulting in the loss of genotypes at metamorphosis that were present in hatchlings. One microsatellite locus exhibited atypically large F_ST values, suggesting it might be linked to a locus under selection. These results support the hypothesis that strong selection against the most inbred genotypes occurred among early life-history stages in our population. Selective forces can promote changes that can affect population dynamics and should be considered in conservation planning.