Outbreeding reduces survival during metamorphosis in a headwater stream salamander

Assessing direct fitness effects of individual genetic diversity is challenging due to the intensive and long-term data needed to quantify survival and reproduction in the wild. But resolving these effects is necessary to determine how inbreeding and outbreeding influence eco-evolutionary processes. We used 8 years of capture-recapture data and single nucleotide polymorphism genotypes for 1906 individuals to test for effects of individual heterozygosity on stage-specific survival probabilities in the salamander Gyrinophilus porphyriticus. The life cycle of G. porphyriticus includes an aquatic larval stage followed by metamorphosis into a semi-aquatic adult stage. In our study populations, the larval stage lasts 6-10 years, metamorphosis takes several months, and lifespan can reach 20 years. Previous studies showed that metamorphosis is a sensitive life stage, leading us to predict that fitness effects of individual heterozygosity would occur during metamorphosis. Consistent with this prediction, monthly probability of survival during metamorphosis declined with multi-locus heterozygosity (MLH), from 0.38 at the lowest MLH (0.10) to 0.06 at the highest MLH (0.38), a reduction of 84%. Body condition of larvae also declined significantly with increasing MLH. These relationships were consistent in the three study streams. With evidence of localised inbreeding within streams, these results suggest that outbreeding disrupts adaptations in pre-metamorphic and metamorphic individuals to environmental gradients along streams, adding to evidence that headwater streams are hotspots of microgeographic adaptation. Our results also underscore the importance of incorporating life history in analyses of the fitness effects of individual genetic diversity and suggest that metamorphosis and similar discrete life stage transitions may be critical periods of viability selection.

Stage-Specific Demographic Effects of Hydrologic Variation in a Stream Salamander

AbstractWe lack a strong understanding of how organisms with complex life histories respond to climate variation. Many stream-associated species have multistage life histories that are likely to influence the demographic consequences of floods and droughts. However, tracking stage-specific demographic responses requires high-resolution, long-term data that are rare. We used 8 years of capture-recapture data for the headwater stream salamander Gyrinophilus porphyriticus to quantify the effects of flooding and drying magnitude on stage-specific vital rates and population growth. Drying reduced larval recruitment but increased the probability of metamorphosis (i.e., adult recruitment). Flooding reduced adult recruitment but had no effect on larval recruitment. Larval and adult survival declined with flooding but were unaffected by drying. Annual population growth rates (λ) declined with flooding and drying. Lambda also declined over the study period (2012-2021), although mean λ was 1.0 over this period. Our results indicate that G. porphyriticus populations are resilient to hydrologic variation because of compensatory effects on recruitment of larvae versus adults (i.e., reproduction vs. metamorphosis). Complex life cycles may enable this resilience to climate variation by creating opportunities for compensatory demographic responses across stages. However, more frequent and intense hydrologic variation in the latter half of this study contributed to a decline in λ over time, suggesting that increasing environmental variability poses a threat even when demographic compensation occurs.

Individual growth rates and size at metamorphosis increase with watershed area in a stream salamander

A fundamental goal of ecology is to understand how the physical environment influences intraspecific variability in life history and, consequently, fitness. In streams, discharge and associated habitat conditions change along a continuum from intermittency to permanence: Headwater streams typically have smaller watersheds and are thus more prone to drying than higher-order streams with larger watersheds and more consistent discharge. However, few empirical studies have assessed life history and associated population responses to this continuum in aquatic organisms. We tested the prediction that individual growth, rate of development, and population growth increase with watershed area in the long-lived stream salamander Gyrinophilus porphyriticus, where we use watershed area as a proxy for hydrologic intermittence. To address this hypothesis, we used 8 years of mark-recapture data from 53 reaches across 10 headwater streams in New Hampshire, USA. Individual growth rates and mean size at metamorphosis increased with watershed area for watersheds from 0.12 to 1.66 km2 . Population growth rates increased with watershed area; however, this result was not statistically significant at our sample size. Mean age of metamorphosis did not vary across watershed areas. Lower individual growth rates and smaller sizes at metamorphosis likely contributed to reduced lifetime fecundity and population growth in reaches with the smallest watershed areas and highest vulnerability to drought. These responses suggest that as droughts increase due to climate change, headwater specialists in hydrologically intermittent environments will experience a reduction in fitness due to smaller body sizes or other growth-related mechanisms.

Source-sink dynamics within a complex life history

Source-sink patch dynamics occur when movement from sources stabilizes sinks by compensating for low local vital rates. The mechanisms underlying source-sink dynamics may be complicated in species that undergo transitions between discrete life stages, particularly when stages have overlapping habitat requirements and similar movement abilities. In these species, for example, the demographic effects of movement by one stage may augment or offset the effects of movement by another stage. We used a stream salamander system to investigate patch dynamics within this form of complex life history. Specifically, we tested the hypothesis that the salamander Gyrinophilus porphyriticus experiences source-sink dynamics in riffles and pools, the dominant geomorphic patch types in headwater streams. We estimated stage-specific survival probabilities in riffles and pools and stage-specific movement probabilities between the two patch types using 8 years of capture-recapture data on 4491 individuals, including premetamorphic larvae and postmetamorphic adults. We then incorporated survival and movement probabilities into a stage-structured, two-patch model to determine the demographic interactions between riffles and pools. Monthly survival probabilities of both stages were higher in pools than in riffles. Larvae were more likely to move from riffles to pools, but adults were more likely to move from pools to riffles, despite experiencing much lower survival in riffles. In simulations, eliminating interpatch movements by both stages indicated that riffles are sinks that rely on immigration from pools for stability. Allowing only larvae to move stabilized both patch types, but allowing only adults to move destabilized pools due to the demographic cost of adult emigration. These results indicated that larval movement not only stabilizes riffles, but also offsets the destabilizing effects of maladaptive adult movement. Similar patch dynamics may emerge in any structured population in which movement and local vital rates differ by age, size, or stage. Addressing these forms of internal demographic structure in patch dynamics analyses will help to refine and advance general understanding of spatial ecology.

Differences in Corticosterone Release Rates of Larval Spring Salamanders (Gyrinophilus porphyriticus) in Response to Native Fish Presence

Simple Summary

In amphibians, glucocorticoid hormones play a key role in the response to predation stress. Predators can directly affect prey via injury and death, but they can also have indirect effects due to the activity of glucocorticoids. The regulation of glucocorticoids can differ between populations that have co-evolved with predators and those that have not. We measured glucocorticoids at baseline and in response to a novel stressor in free-living larval salamanders that either live with or without fish predators naturally. We found that salamanders living with fish predators had lower measures of glucocorticoids than those without fish predators. Our study indicates that predator presence alters glucocorticoid regulation, which may allow species to better cope with native and introduced predators.

Abstract

Invasive fish predators are an important factor causing amphibian declines and may have direct and indirect effects on amphibian survival. For example, early non-lethal exposure to these stressors may reduce survival in later life stages, especially in biphasic species. In amphibians, the glucocorticoid hormone corticosterone is released by the hypothalamo–pituitary–interrenal axis (HPI), as an adaptive physiological response to environmental stressors. The corticosterone response (baseline and response to acute stressors) is highly flexible and context dependent, and this variation can allow individuals to alter their phenotype and behavior with environmental changes, ultimately increasing survival. We sampled larvae of the spring salamander (Gyrinophilus porphyriticus) from two streams that each contained predatory brook trout (Slavelinus fontinalis) in the lower reaches and no predatory brook trout in the upper reaches. We measured baseline and stress-induced corticosterone release rates of larvae from the lower and upper reaches using a non-invasive water-borne hormone assay. We hypothesized that corticosterone release rates would differ between larvae from fish-present reaches and larvae from fish-free reaches. We found that baseline and stressor-induced corticosterone release rates were downregulated in larvae from reaches with fish predators. These results indicate that individuals from reaches with predatory trout are responding to fish predators by downregulating corticosterone while maintaining an active HPI axis. This may allow larvae more time to grow before metamorphosing, while also allowing them to physiologically respond to novel stressors. However, prolonged downregulation of corticosterone release rates can impact growth in post-metamorphic individuals.

Sex-related differences in aging rate are associated with sex chromosome system in amphibians

Sex‐related differences in mortality are widespread in the animal kingdom. Although studies have shown that sex determination systems might drive lifespan evolution, sex chromosome influence on aging rates have not been investigated so far, likely due to an apparent lack of demographic data from clades including both XY (with heterogametic males) and ZW (heterogametic females) systems. Taking advantage of a unique collection of capture–recapture datasets in amphibians, a vertebrate group where XY and ZW systems have repeatedly evolved over the past 200 million years, we examined whether sex heterogamy can predict sex differences in aging rates and lifespans. We showed that the strength and direction of sex differences in aging rates (and not lifespan) differ between XY and ZW systems. Sex‐specific variation in aging rates was moderate within each system, but aging rates tended to be consistently higher in the heterogametic sex. This led to small but detectable effects of sex chromosome system on sex differences in aging rates in our models. Although preliminary, our results suggest that exposed recessive deleterious mutations on the X/Z chromosome (the “unguarded X/Z effect”) or repeat‐rich Y/W chromosome (the “toxic Y/W effect”) could accelerate aging in the heterogametic sex in some vertebrate clades.

Hybridization alters growth and migratory life-history expression of native trout

Human-mediated hybridization threatens many native species, but the effects of introgressive hybridization on life-history expression are rarely quantified, especially in vertebrates. We quantified the effects of non-native rainbow trout admixture on important life-history traits including growth and partial migration behavior in three populations of westslope cutthroat trout over five years. Rainbow trout admixture was associated with increased summer growth rates in all populations and decreased spring growth rates in two populations with cooler spring temperatures. These results indicate that non-native admixture may increase growth under warmer conditions, but cutthroat trout have higher growth rates during cooler periods. Non-native admixture consistently increased expression of migratory behavior, suggesting that there is a genomic basis for life-history differences between these species. Our results show that effects of interspecific hybridization on fitness traits can be the product of genotype-by-environment interactions even when there are minor differences in environmental optima between hybridizing species. These results also indicate that while environmentally mediated traits like growth may play a role in population-level consequences of admixture, strong genetic influences on migratory life-history differences between these species likely explains the continued spread of non-native hybridization at the landscape-level, despite selection against hybrids at the population-level.

Metamorphosis in an Era of Increasing Climate Variability

Most animals have complex life cycles including metamorphosis or other discrete life stage transitions, during which individuals may be particularly vulnerable to environmental stressors. With climate change, individuals will be exposed to increasing thermal and hydrologic variability during metamorphosis, which may affect survival and performance through physiological, behavioral, and ecological mechanisms. Furthermore, because metamorphosis entails changes in traits and vital rates, it is likely to play an important role in how populations respond to increasing climate variability. To identify mechanisms underlying population responses and associated trait and life history evolution, we need new approaches to estimating changes in individual traits and performance throughout metamorphosis, and we need to integrate metamorphosis as an explicit life stage in analytical models.

Long-term survival probability, not current habitat quality, predicts dispersal distance in a stream salamander

Dispersal evolves as an adaptive mechanism to optimize individual fitness across the landscape. Specifically, dispersal represents a mechanism to escape fitness costs resulting from changes in environmental conditions. Decades of empirical work suggest that individuals use local habitat cues to make movement decisions, but theory predicts that dispersal can also evolve as a fixed trait, independent of local conditions, in environments characterized by a history of stochastic spatiotemporal variation. Until now, however, both conditional and fixed models of dispersal evolution have primarily been evaluated using emigration data (stay vs. leave), and not dispersal distances: a more comprehensive measure of dispersal. Our goal was to test whether conditional or fixed models of dispersal evolution predict variation in dispersal distance in the stream salamander Gyrinophilus porphyriticus. We quantified variation in habitat conditions using measures of salamander performance from 4 yr of spatially explicit, capture-mark-recapture (CMR) data across three headwater streams in the Hubbard Brook Experimental Forest in central New Hampshire, USA. We used body condition as an index of local habitat quality that individuals may use to make dispersal decisions, and survival probability estimated from multistate CMR models as an index of mortality risk resulting from the long-term history of environmental variation. We found that dispersal distances increased with declining survival probability, indicating that salamanders disperse further in risky environments. Dispersal distances were unrelated to spatial variation in body condition, suggesting that salamanders do not base dispersal distance decisions on local habitat quality. Our study provides the first empirical support for fixed models of dispersal evolution, which predict that dispersal evolves in response to a history of spatiotemporal environmental variation, rather than as a conditional response to current habitat conditions. More broadly, this study underscores the value of assessing alternative scales of environmental variation to gain a more complete and balanced understanding of dispersal evolution.

The Landscape Ecology of Rivers: from Patch-Based to Spatial Network Analyses

Purpose of Review

We synthesize recent methodological and conceptual advances in the field of riverscape ecology, emphasizing areas of synergy with current research in landscape ecology.

Recent Findings

Recent advances in riverscape ecology highlight the need for spatially explicit examinations of how network structure influences ecological pattern and process, instead of the simple linear (upstream-downstream) view. Developments in GIS, remote sensing, and computer technologies already offer powerful tools for the application of patch- and gradient-based models for characterizing abiotic and biotic heterogeneity across a range of spatial and temporal scales. Along with graph-based analyses and spatial statistical stream network models (i.e., geostatistical modelling), these approaches offer improved capabilities for quantifying spatial and temporal heterogeneity and connectivity relationships, thereby allowing for rigorous and high-resolution analyses of pattern, process, and scale relationships.

Summary

Spatially explicit network approaches are able to quantify and predict biogeochemical, hydromorphological, and ecological patterns and processes more precisely than models based on longitudinal or lateral riverine gradients alone. Currently, local habitat characteristics appear to be more important than spatial effects in determining population and community dynamics, but this conclusion may change with direct quantification of the movement of materials, energy, and organisms along channels and across ecosystem boundaries—a key to improving riverscape ecology. Coupling spatially explicit riverscape models with optimization approaches will improve land protection and water management efforts, and help to resolve the land sharing vs. land sparing debate.

A distance-performance trade-off in the phenotypic basis of dispersal

Across taxa, individuals vary in how far they disperse, with most individuals staying close to their origin and fewer dispersing long distances. Costs associated with dispersal (e.g., energy, risk) are widely believed to trade off with benefits (e.g., reduced competition, increased reproductive success) to influence dispersal propensity. However, this framework has not been applied to understand variation in dispersal distance, which is instead generally attributed to extrinsic environmental factors. We alternatively hypothesized that variation in dispersal distances results from trade-offs associated with other aspects of locomotor performance. We tested this hypothesis in the stream salamander Gyrinophilus porphyriticus and found that salamanders that dispersed farther in the field had longer forelimbs but swam at slower velocities under experimental conditions. The reduced swimming performance of long-distance dispersers likely results from drag imposed by longer forelimbs. Longer forelimbs may facilitate moving longer distances, but the proximate costs associated with reduced swimming performance may help to explain the rarity of long-distance dispersal. The historical focus on environmental drivers of dispersal distances misses the importance of individual traits and associated trade-offs among traits affecting locomotion.

Spatial capture-recapture reveals age- and sex-specific survival and movement in stream amphibians

Life-history information sets the foundation for our understanding of ecology and conservation requirements. For many species, this information is lacking even for basic demographic rates such as survival and movement. When survival and movement estimates are available, they are often derived from mixed demographic groups and do not consider differences among life stages or sexes, which is critical, because life stages and sexes often contribute differentially to population dynamics. We used hierarchical models informed with spatial capture-mark-recapture data of Ascaphus montanus (Rocky Mountain tailed frog) in five streams and A. truei (coastal tailed frog) in one stream to estimate variation in survival and movement by sex and age, represented by size. By incorporating survival and movement into a single model, we were able to estimate both parameters with limited bias. Annual survival was similar between sexes of A. montanus [females = 0.885 (95% CI 0.614-1), males = 0.901 (0.657-1)], but was slightly higher for female A. truei [0.836 (0.560-0.993)] than for males [0.664 (0.354-0.962)]. Survival of A. montanus peaked at mid-age, suggesting that lower survival of young and actuarial senescence may influence population demographics. Our models suggest that younger A. montanus moved farther than older individuals, and that females moved farther than males in both species. Our results provide uncommon insight into age- and sex-specific rates of survival and movement that are crucial elements of life-history strategies and are important for modeling population growth and prescribing conservation actions.

Mercury bioaccumulation in temperate forest food webs associated with headwater streams

The soils and food webs associated with mid to high elevation, forested, headwater streams in northeastern North America are potential hotspots for mercury (Hg) methylation and bioaccumulation, but are not well studied. Our goals were to quantify total Hg (THg) and methyl Hg (MeHg) concentrations in soils and terrestrial food webs associated with headwater streams of northern hardwood forests to identify predictors of small-scale spatial variation in Hg bioaccumulation. We sampled soil characteristics that promote Hg methylation including pH, sulfur and calcium content, and organic matter. To assess spatial variation, we sampled at high (~700 m asl) and mid elevations (~500 m asl), both adjacent to (<1 m) and away from (>75 m) three replicate headwater streams in each of two watersheds of the White Mountains region, New Hampshire, USA. Soils of these forested watersheds differed significantly in pH and the content of calcium, sulfur, organic matter and THg. Conditions for methylation were more favorable in the upland forest sites compared to streamside sites. Significant bioaccumulation of THg occurred in all measured components of the food web, including insects, spiders, salamanders, and birds. Trophic position, as determined by δ¹⁵N, was the best predictor of both THg and MeHg bioaccumulation across the sampled taxa and was also a better predictor than spatial location. However, the degree of bioaccumulation at which MeHg significantly affects animal behavior, reproduction or survival is unknown for most taxa in terrestrial habitats, particularly for invertebrates. These findings show that Hg methylation and bioaccumulation is not limited to areas traditionally classified as wetlands or to areas with exceptionally high THg inputs, but that it is a widespread and important phenomenon in the moist deciduous forests of eastern North America.

Vive la résistance: genome-wide selection against introduced alleles in invasive hybrid zones

Evolutionary and ecological consequences of hybridization between native and invasive species are notoriously complicated because patterns of selection acting on non-native alleles can vary throughout the genome and across environments. Rapid advances in genomics now make it feasible to assess locus-specific and genome-wide patterns of natural selection acting on invasive introgression within and among natural populations occupying diverse environments. We quantified genome-wide patterns of admixture across multiple independent hybrid zones of native westslope cutthroat trout and invasive rainbow trout, the world's most widely introduced fish, by genotyping 339 individuals from 21 populations using 9380 species-diagnostic loci. A significantly greater proportion of the genome appeared to be under selection favouring native cutthroat trout (rather than rainbow trout), and this pattern was pervasive across the genome (detected on most chromosomes). Furthermore, selection against invasive alleles was consistent across populations and environments, even in those where rainbow trout were predicted to have a selective advantage (warm environments). These data corroborate field studies showing that hybrids between these species have lower fitness than the native taxa, and show that these fitness differences are due to selection favouring many native genes distributed widely throughout the genome.

Legacy introductions and climatic variation explain spatiotemporal patterns of invasive hybridization in a native trout

Hybridization between invasive and native species, a significant threat to worldwide biodiversity, is predicted to increase due to climate-induced expansions of invasive species. Long-term research and monitoring are crucial for understanding the ecological and evolutionary processes that modulate the effects of invasive species. Using a large, multidecade genetics dataset (N = 582 sites, 12,878 individuals) with high-resolution climate predictions and extensive stocking records, we evaluate the spatiotemporal dynamics of hybridization between native cutthroat trout and invasive rainbow trout, the world's most widely introduced invasive fish, across the Northern Rocky Mountains of the United States. Historical effects of stocking and contemporary patterns of climatic variation were strongly related to the spread of hybridization across space and time. The probability of occurrence, extent of, and temporal changes in hybridization increased at sites in close proximity to historical stocking locations with greater rainbow trout propagule pressure, warmer water temperatures, and lower spring precipitation. Although locations with warmer water temperatures were more prone to hybridization, cold sites were not protected from invasion; 58% of hybridized sites had cold mean summer water temperatures (<11°C). Despite cessation of stocking over 40 years ago, hybridization increased over time at half (50%) of the locations with long-term data, the vast majority of which (74%) were initially nonhybridized, emphasizing the chronic, negative impacts of human-mediated hybridization. These results show that effects of climate change on biodiversity must be analyzed in the context of historical human impacts that set ecological and evolutionary trajectories.

Gravel-bed river floodplains are the ecological nexus of glaciated mountain landscapes

Gravel-bed river floodplains in mountain landscapes disproportionately concentrate diverse habitats, nutrient cycling, productivity of biota, and species interactions. Although stream ecologists know that river channel and floodplain habitats used by aquatic organisms are maintained by hydrologic regimes that mobilize gravel-bed sediments, terrestrial ecologists have largely been unaware of the importance of floodplain structures and processes to the life requirements of a wide variety of species. We provide insight into gravel-bed rivers as the ecological nexus of glaciated mountain landscapes. We show why gravel-bed river floodplains are the primary arena where interactions take place among aquatic, avian, and terrestrial species from microbes to grizzly bears and provide essential connectivity as corridors for movement for both aquatic and terrestrial species. Paradoxically, gravel-bed river floodplains are also disproportionately unprotected where human developments are concentrated. Structural modifications to floodplains such as roads, railways, and housing and hydrologic-altering hydroelectric or water storage dams have severe impacts to floodplain habitat diversity and productivity, restrict local and regional connectivity, and reduce the resilience of both aquatic and terrestrial species, including adaptation to climate change. To be effective, conservation efforts in glaciated mountain landscapes intended to benefit the widest variety of organisms need a paradigm shift that has gravel-bed rivers and their floodplains as the central focus and that prioritizes the maintenance or restoration of the intact structure and processes of these critically important systems throughout their length and breadth.

Human effects on ecological connectivity in aquatic ecosystems: Integrating scientific approaches to support management and mitigation

Understanding the drivers and implications of anthropogenic disturbance of ecological connectivity is a key concern for the conservation of biodiversity and ecosystem processes. Here, we review human activities that affect the movements and dispersal of aquatic organisms, including damming of rivers, river regulation, habitat loss and alteration, human-assisted dispersal of organisms and climate change. Using a series of case studies, we show that the insight needed to understand the nature and implications of connectivity, and to underpin conservation and management, is best achieved via data synthesis from multiple analytical approaches. We identify four key knowledge requirements for progressing our understanding of the effects of anthropogenic impacts on ecological connectivity: autecology; population structure; movement characteristics; and environmental tolerance/phenotypic plasticity. Structuring empirical research around these four broad data requirements, and using this information to parameterise appropriate models and develop management approaches, will allow for mitigation of the effects of anthropogenic disturbance on ecological connectivity in aquatic ecosystems.

Dispersal and selection mediate hybridization between a native and invasive species

Hybridization between native and non-native species has serious biological consequences, but our understanding of how dispersal and selection interact to influence invasive hybridization is limited. Here, we document the spread of genetic introgression between a native (Oncorhynchus clarkii) and invasive (Oncorhynchus mykiss) trout, and identify the mechanisms influencing genetic admixture. In two populations inhabiting contrasting environments, non-native admixture increased rapidly from 1984 to 2007 and was driven by surprisingly consistent processes. Individual admixture was related to two phenotypic traits associated with fitness: size at spawning and age of juvenile emigration. Fish with higher non-native admixture were larger and tended to emigrate at a younger age--relationships that are expected to confer fitness advantages to hybrid individuals. However, strong selection against non-native admixture was evident across streams and cohorts (mean selection coefficient against genotypes with non-native alleles (s) = 0.60; s.e. = 0.10). Nevertheless, hybridization was promoted in both streams by the continuous immigration of individuals with high levels of non-native admixture from other hybrid source populations. Thus, antagonistic relationships between dispersal and selection are mediating invasive hybridization between these fish, emphasizing that data on dispersal and natural selection are needed to fully understand the dynamics of introgression between native and non-native species.

Genomics and introgression: Discovery and mapping of thousands of species-diagnostic SNPs using RAD sequencing

Invasive hybridization and introgression pose a serious threat to the persistence of many native species. Understanding the effects of hybridization on native populations (e.g., fitness consequences) requires numerous species-diagnostic loci distributed genome-wide. Here we used RAD sequencing to discover thousands of single-nucleotide polymorphisms (SNPs) that are diagnostic between rainbow trout (RBT, Oncorhynchus mykiss), the world’s most widely introduced fish, and native westslope cutthroat trout (WCT, O. clarkii lewisi) in the northern Rocky Mountains, USA. We advanced previous work that identified 4,914 species-diagnostic loci by using longer sequence reads (100 bp vs. 60 bp) and a larger set of individuals (n = 84). We sequenced RAD libraries for individuals from diverse sampling sources, including native populations of WCT and hatchery broodstocks of WCT and RBT. We also took advantage of a newly released reference genome assembly for RBT to align our RAD loci. In total, we discovered 16,788 putatively diagnostic SNPs, 10,267 of which we mapped to anchored chromosome locations on the RBT genome. A small portion of previously discovered putative diagnostic loci (325 of 4,914) were no longer diagnostic (i.e., fixed between species) based on our wider survey of non-hybridized RBT and WCT individuals. Our study suggests that RAD loci mapped to a draft genome assembly could provide the marker density required to identify genes and chromosomal regions influencing selection in admixed populations of conservation concern and evolutionary interest.

Distance, flow and PCR inhibition: eDNA dynamics in two headwater streams

Environmental DNA (eDNA) detection has emerged as a powerful tool for monitoring aquatic organisms, but much remains unknown about the dynamics of aquatic eDNA over a range of environmental conditions. DNA concentrations in streams and rivers will depend not only on the equilibrium between DNA entering the water and DNA leaving the system through degradation, but also on downstream transport. To improve understanding of the dynamics of eDNA concentration in lotic systems, we introduced caged trout into two fishless headwater streams and took eDNA samples at evenly spaced downstream intervals. This was repeated 18 times from mid-summer through autumn, over flows ranging from approximately 1-96 L/s. We used quantitative PCR to relate DNA copy number to distance from source. We found that regardless of flow, there were detectable levels of DNA at 239.5 m. The main effect of flow on eDNA counts was in opposite directions in the two streams. At the lowest flows, eDNA counts were highest close to the source and quickly trailed off over distance. At the highest flows, DNA counts were relatively low both near and far from the source. Biomass was positively related to eDNA copy number in both streams. A combination of cell settling, turbulence and dilution effects is probably responsible for our observations. Additionally, during high leaf deposition periods, the presence of inhibitors resulted in no amplification for high copy number samples in the absence of an inhibition-releasing strategy, demonstrating the necessity to carefully consider inhibition in eDNA analysis.

Combining demographic and genetic factors to assess population vulnerability in stream species

Accelerating climate change and other cumulative stressors create an urgent need to understand the influence of environmental variation and landscape features on the connectivity and vulnerability of freshwater species. Here, we introduce a novel modeling framework for aquatic systems that integrates spatially explicit, individual-based, demographic and genetic (demogenetic) assessments with environmental variables. To show its potential utility, we simulated a hypothetical network of 19 migratory riverine populations (e.g., salmonids) using a riverscape connectivity and demogenetic model (CDFISH). We assessed how stream resistance to movement (a function of water temperature, fluvial distance, and physical barriers) might influence demogenetic connectivity, and hence, population vulnerability. We present demographic metrics (abundance, immigration, and change in abundance) and genetic metrics (diversity, differentiation, and change in differentiation), and combine them into a single vulnerability index for identifying populations at risk of extirpation. We considered four realistic scenarios that illustrate the relative sensitivity of these metrics for early detection of reduced connectivity: (1) maximum resistance due to high water temperatures throughout the network, (2) minimum resistance due to low water temperatures throughout the network, (3) increased resistance at a tributary junction caused by a partial barrier, and (4) complete isolation of a tributary, leaving resident individuals only. We then applied this demogenetic framework using empirical data for a bull trout (Salvelinus confluentus) metapopulation in the upper Flathead River system, Canada and USA, to assess how current and predicted future stream warming may influence population vulnerability. Results suggest that warmer water temperatures and associated barriers to movement (e.g., low flows, dewatering) are predicted to fragment suitable habitat for migratory salmonids, resulting in the loss of genetic diversity and reduced numbers in certain vulnerable populations. This demogenetic simulation framework, which is illustrated in a web-based interactive mapping prototype, should be useful for evaluating population vulnerability in a wide variety of dendritic and fragmented riverscapes, helping to guide conservation and management efforts for freshwater species.

Robust detection of rare species using environmental DNA: the importance of primer specificity

Environmental DNA (eDNA) is being rapidly adopted as a tool to detect rare animals. Quantitative PCR (qPCR) using probe-based chemistries may represent a particularly powerful tool because of the method’s sensitivity, specificity, and potential to quantify target DNA. However, there has been little work understanding the performance of these assays in the presence of closely related, sympatric taxa. If related species cause any cross-amplification or interference, false positives and negatives may be generated. These errors can be disastrous if false positives lead to overestimate the abundance of an endangered species or if false negatives prevent detection of an invasive species. In this study we test factors that influence the specificity and sensitivity of TaqMan MGB assays using co-occurring, closely related brook trout (Salvelinus fontinalis) and bull trout (S. confluentus) as a case study. We found qPCR to be substantially more sensitive than traditional PCR, with a high probability of detection at concentrations as low as 0.5 target copies/µl. We also found that number and placement of base pair mismatches between the Taqman MGB assay and non-target templates was important to target specificity, and that specificity was most influenced by base pair mismatches in the primers, rather than in the probe. We found that insufficient specificity can result in both false positive and false negative results, particularly in the presence of abundant related species. Our results highlight the utility of qPCR as a highly sensitive eDNA tool, and underscore the importance of careful assay design.

Interactive effects of wildfire, forest management, and isolation on amphibian and parasite abundance

Projected increases in wildfire and other climate-driven disturbances will affect populations and communities worldwide, including host-parasite relationships. Research in temperate forests has shown that wildfire can negatively affect amphibians, but this research has occurred primarily outside of managed landscapes where interactions with human disturbances could result in additive or synergistic effects. Furthermore, parasites represent a large component of biodiversity and can affect host fitness and population dynamics, yet they are rarely included in studies of how vertebrate hosts respond to disturbance. To determine how wildfire affects amphibians and their parasites, and whether effects differ between protected and managed landscapes, we compared abundance of two amphibians and two nematodes relative to wildfire extent and severity around wetlands in neighboring protected and managed forests (Montana, USA). Population sizes of adult, male long-toed salamanders (Ambystoma macrodactylum) decreased with increased burn severity, with stronger negative effects on isolated populations and in managed forests. In contrast, breeding population sizes of Columbia spotted frogs (Rana luteiventris) increased with burn extent in both protected and managed protected forests. Path analysis showed that the effects of wildfire on the two species of nematodes were consistent with differences in their life history and transmission strategies and the responses of their hosts. Burn severity indirectly reduced abundance of soil-transmitted Cosmocercoides variabilis through reductions in salamander abundance. Burn severity also directly reduced C. variabilis abundance, possibly though changes in soil conditions. For the aquatically transmitted nematode Gyrinicola batrachiensis, the positive effect of burn extent on density of Columbia spotted frog larvae indirectly increased parasite abundance. Our results show that effects of wildfire on amphibians depend upon burn extent and severity, isolation, and prior land use. Through subsequent effects on the parasites, our results also reveal how changes in disturbance regimes can affect communities across trophic levels.

Rapid increases and time-lagged declines in amphibian occupancy after wildfire

Climate change is expected to increase the frequency and severity of drought and wildfire. Aquatic and moisture-sensitive species, such as amphibians, may be particularly vulnerable to these modified disturbance regimes because large wildfires often occur during extended droughts and thus may compound environmental threats. However, understanding of the effects of wildfires on amphibians in forests with long fire-return intervals is limited. Numerous stand-replacing wildfires have occurred since 1988 in Glacier National Park (Montana, U.S.A.), where we have conducted long-term monitoring of amphibians. We measured responses of 3 amphibian species to fires of different sizes, severity, and age in a small geographic area with uniform management. We used data from wetlands associated with 6 wildfires that burned between 1988 and 2003 to evaluate whether burn extent and severity and interactions between wildfire and wetland isolation affected the distribution of breeding populations. We measured responses with models that accounted for imperfect detection to estimate occupancy during prefire (0-4 years) and different postfire recovery periods. For the long-toed salamander (Ambystoma macrodactylum) and Columbia spotted frog (Rana luteiventris), occupancy was not affected for 6 years after wildfire. But 7-21 years after wildfire, occupancy for both species decreased ≥ 25% in areas where >50% of the forest within 500 m of wetlands burned. In contrast, occupancy of the boreal toad (Anaxyrus boreas) tripled in the 3 years after low-elevation forests burned. This increase in occupancy was followed by a gradual decline. Our results show that accounting for magnitude of change and time lags is critical to understanding population dynamics of amphibians after large disturbances. Our results also inform understanding of the potential threat of increases in wildfire frequency or severity to amphibians in the region.

Linking extinction-colonization dynamics to genetic structure in a salamander metapopulation

Theory predicts that founder effects have a primary role in determining metapopulation genetic structure. However, ecological factors that affect extinction-colonization dynamics may also create spatial variation in the strength of genetic drift and migration. We tested the hypothesis that ecological factors underlying extinction-colonization dynamics influenced the genetic structure of a tiger salamander (Ambystoma tigrinum) metapopulation. We used empirical data on metapopulation dynamics to make a priori predictions about the effects of population age and ecological factors on genetic diversity and divergence among 41 populations. Metapopulation dynamics of A. tigrinum depended on wetland area, connectivity and presence of predatory fish. We found that newly colonized populations were more genetically differentiated than established populations, suggesting that founder effects influenced genetic structure. However, ecological drivers of metapopulation dynamics were more important than age in predicting genetic structure. Consistent with demographic predictions from metapopulation theory, genetic diversity and divergence depended on wetland area and connectivity. Divergence was greatest in small, isolated wetlands where genetic diversity was low. Our results show that ecological factors underlying metapopulation dynamics can be key determinants of spatial genetic structure, and that habitat area and isolation may mediate the contributions of drift and migration to divergence and evolution in local populations.

Coexistence in streams: do source-sink dynamics allow salamanders to persist with fish predators?

Theory suggests that source-sink dynamics can allow coexistence of intraguild predators and prey, but empirical evidence for this coexistence mechanism is limited. We used capture-mark-recapture, genetic methods, and stable isotopes to test whether source-sink dynamics promote coexistence between stream fishes, the intraguild predator, and stream salamanders (Dicamptodon aterrimus), the intraguild prey. Salamander populations from upstream reaches without fish were predicted to maintain or supplement sink populations in downstream reaches with fish. We found instead that downstream reaches with fish were not sinks even though fish consumed salamander larvae-apparent survival, recruitment, and population growth rate did not differ between upstream and downstream reaches. There was also no difference between upstream and downstream reaches in net emigration. We did find that D. aterrimus moved frequently along streams, but believe that this is a response to seasonal habitat changes rather than intraguild predation. Our study provides empirical evidence that local-scale mechanisms are more important than dispersal dynamics to coexistence of streams salamanders and fish. More broadly, it shows the value of empirical data on dispersal and gene flow for distinguishing between local and spatial mechanisms of coexistence.

Explaining long-distance dispersal: effects of dispersal distance on survival and growth in a stream salamander

Long-distance dispersal (LDD) may contribute disproportionately to range expansions, the creation of new evolutionary lineages, and species persistence in human-dominated landscapes. However, because data on the individual consequences of dispersal distance are extremely limited, we have little insight on how LDD is maintained in natural populations. I used six years of spatially explicit capture-mark-recapture (CMR) data to test the prediction that individual performance increases with dispersal distance in the stream salamander Gyrinophilus porphyriticus. Dispersal distance was total distance moved along the 1-km study stream, ranging from 0 to 565 m. To quantify individual performance, I used CMR estimates of survival and individual growth rates based on change in body length. Survival and growth rates increased significantly with dispersal distance. These relationships were not confounded by pre-dispersal body condition or by ecological gradients along the stream. Individual benefits of LDD were likely caused by an increase in the upper limit of settlement site quality with dispersal distance. My results do not support the view that the fitness consequences of LDD are unpredictable and instead suggest that consistent evolutionary mechanisms may explain the prevalence of LDD in nature. They also highlight the value of direct CMR data for understanding the individual consequences of variation in dispersal distance and how that variation is maintained in natural populations.

Use of multiple dispersal pathways facilitates amphibian persistence in stream networks

Although populations of amphibians are declining worldwide, there is no evidence that salamanders occupying small streams are experiencing enigmatic declines, and populations of these species seem stable. Theory predicts that dispersal through multiple pathways can stabilize populations, preventing extinction in habitat networks. However, empirical data to support this prediction are absent for most species, especially those at risk of decline. Our mark-recapture study of stream salamanders reveals both a strong upstream bias in dispersal and a surprisingly high rate of overland dispersal to adjacent headwater streams. This evidence of route-dependent variation in dispersal rates suggests a spatial mechanism for population stability in headwater-stream salamanders. Our results link the movement behavior of stream salamanders to network topology, and they underscore the importance of identifying and protecting critical dispersal pathways when addressing region-wide population declines.

Living in the branches: population dynamics and ecological processes in dendritic networks

Spatial structure regulates and modifies processes at several levels of ecological organization (e.g. individual/genetic, population and community) and is thus a key component of complex systems, where knowledge at a small scale can be insufficient for understanding system behaviour at a larger scale. Recent syntheses outline potential applications of network theory to ecological systems, but do not address the implications of physical structure for network dynamics. There is a specific need to examine how dendritic habitat structure, such as that found in stream, hedgerow and cave networks, influences ecological processes. Although dendritic networks are one type of ecological network, they are distinguished by two fundamental characteristics: (1) both the branches and the nodes serve as habitat, and (2) the specific spatial arrangement and hierarchical organization of these elements interacts with a species' movement behaviour to alter patterns of population distribution and abundance, and community interactions. Here, we summarize existing theory relating to ecological dynamics in dendritic networks, review empirical studies examining the population- and community-level consequences of these networks, and suggest future research integrating spatial pattern and processes in dendritic systems.

Linking direct and indirect data on dispersal: isolation by slope in a headwater stream salamander

There is growing recognition of the need to incorporate information on movement behavior in landscape-scale studies of dispersal. One way to do this is by using indirect indices of dispersal (e.g., genetic differentiation) to test predictions derived from direct data on movement behavior. Mark-recapture studies documented upstream-biased movement in the salamander Gyrinophilus porphyriticus (Plethodontidae). Based on this information, we hypothesized that gene flow in G. porphyriticus is affected by the slope of the stream. Specifically, because the energy required for upstream dispersal is positively related to slope, we predicted gene flow to be negatively related to change in elevation between sampling sites. Using amplified DNA fragment length polymorphisms among tissue samples from paired sites in nine streams in the Hubbard Brook Watershed, New Hampshire, USA, we found that genetic distances between downstream and upstream sites were positively related to change in elevation over standardized 1-km distances. This pattern of isolation by slope elucidates controls on population connectivity in stream networks and underscores the potential for specific behaviors to affect the genetic structure of species at the landscape scale. More broadly, our results show the value of combining direct data on movement behavior and indirect indices to assess patterns and consequences of dispersal in spatially complex ecosystems.

Landscape-scale spatial population dynamics in human-impacted stream systems

The movement of individuals among populations can be critical in preventing local and landscape-scale species extinctions in systems exposed to human perturbation. Current understanding of spatial population dynamics in streams is largely limited to the reach scale and is therefore inadequate to address species response to spatially extensive perturbation. Using model simulations, I examined species response to perturbation in a drainage composed of multiple, hierarchically arranged stream-patches connected by in-stream and overland pathways of dispersal. Patch extinction probability, the proportion of initially occupied patches extinct after 25 years, was highly sensitive to the extent of species occupancy and perturbation within the drainage, longitudinal species distribution, perturbation decay rate and the covariance pattern of stochastic effects on colonization and extinction probabilities. Results of these simulations underscore the importance of identifying and preserving source populations and dispersal routes for stream species in human-impacted landscapes. They also highlight the vulnerability of headwater specialist taxa to anthropogenic perturbation, and the strong positive effect on species resilience of habitat rehabilitation when recolonization is possible. Efforts to conserve and manage stream species may be greatly improved by accounting for landscape-scale spatial population dynamics.