Ecological correlates of population genetic structure: a comparative approach using a vertebrate metacommunity

What Have Long-Term Field Studies Taught Us About Population Dynamics?

Long-term studies have been crucial to the advancement of population biology, especially our understanding of population dynamics. We argue that this progress arises from three key characteristics of long-term research. First, long-term data are necessary to observe the heterogeneity that drives most population processes. Second, long-term studies often inherently lead to novel insights. Finally, long-term field studies can serve as model systems for population biology, allowing for theory and methods to be tested under well-characterized conditions. We illustrate these ideas in three long-term field systems that have made outsized contributions to our understanding of population ecology, evolution, and conservation biology. We then highlight three emerging areas to which long-term field studies are well positioned to contribute in the future: ecological forecasting, genomics, and macrosystems ecology. Overcoming the obstacles associated with maintaining long-term studies requires continued emphasis on recognizing the benefits of such studies to ensure that long-term research continues to have a substantial impact on elucidating population biology.

Reliable effective number of breeders/adult census size ratios in seasonal-breeding species: Opportunity for integrative demographic inferences based on capture-mark-recapture data and multilocus genotypes

The ratio of the effective number of breeders (Nb) to the adult census size (Na), Nb/Na, approximates the departure from the standard capacity of a population to maintain genetic diversity in one reproductive season. This information is relevant for assessing population status, understanding evolutionary processes operating at local scales, and unraveling how life-history traits affect these processes. However, our knowledge on Nb/Na ratios in nature is limited because estimation of both parameters is challenging. The sibship frequency (SF) method is adequate for reliable Nb estimation because it is based on sibship and parentage reconstruction from genetic marker data, thereby providing demographic inferences that can be compared with field-based information. In addition, capture-mark-recapture (CMR) robust design methods are well suited for Na estimation in seasonal-breeding species. We used tadpole genotypes of three pond-breeding amphibian species (Epidalea calamita, Hyla molleri, and Pelophylax perezi, n = 73-96 single-cohort tadpoles/species genotyped at 15-17 microsatellite loci) and candidate parental genotypes (n = 94-300 adults/species) to estimate Nb by the SF method. To assess the reliability of Nb estimates, we compared sibship and parentage inferences with field-based information and checked for the convergence of results in replicated subsampled analyses. Finally, we used CMR data from a 6-year monitoring program to estimate annual Na in the three species and calculate the Nb/Na ratio. Reliable ratios were obtained for E. calamita (Nb/Na = 0.18-0.28) and P. perezi (0.5), but in H. molleri, Na could not be estimated and genetic information proved insufficient for reliable Nb estimation. Integrative demographic studies taking full advantage of SF and CMR methods can provide accurate estimates of the Nb/Na ratio in seasonal-breeding species. Importantly, the SF method provides results that can be readily evaluated for reliability. This represents a good opportunity for obtaining robust demographic inferences with wide applications for evolutionary and conservation research.

Differential influences of local subpopulations on regional diversity and differentiation for greater sage-grouse (Centrocercus urophasianus)

The distribution of spatial genetic variation across a region can shape evolutionary dynamics and impact population persistence. Local population dynamics and among-population dispersal rates are strong drivers of this spatial genetic variation, yet for many species we lack a clear understanding of how these population processes interact in space to shape within-species genetic variation. Here, we used extensive genetic and demographic data from 10 subpopulations of greater sage-grouse to parameterize a simulated approximate Bayesian computation (ABC) model and (i) test for regional differences in population density and dispersal rates for greater sage-grouse subpopulations in Wyoming, and (ii) quantify how these differences impact subpopulation regional influence on genetic variation. We found a close match between observed and simulated data under our parameterized model and strong variation in density and dispersal rates across Wyoming. Sensitivity analyses suggested that changes in dispersal (via landscape resistance) had a greater influence on regional differentiation, whereas changes in density had a greater influence on mean diversity across all subpopulations. Local subpopulations, however, varied in their regional influence on genetic variation. Decreases in the size and dispersal rates of central populations with low overall and net immigration (i.e. population sources) had the greatest negative impact on genetic variation. Overall, our results provide insight into the interactions among demography, dispersal and genetic variation and highlight the potential of ABC to disentangle the complexity of regional population dynamics and project the genetic impact of changing conditions.

Navigating the pitfalls and promise of landscape genetics

The field of landscape genetics has been evolving rapidly since its emergence in the early 2000s. New applications, techniques and criticisms of techniques appear like clockwork with each new journal issue. The developments are an encouraging, and at times bewildering, sign of progress in an exciting new field of study. However, we suggest that the rapid expansion of landscape genetics has belied important flaws in the development of the field, and we add an air of caution to this breakneck pace of expansion. Specifically, landscape genetic studies often lose sight of the fundamental principles and complex consequences of gene flow, instead favouring simplistic interpretations and broad inferences not necessarily warranted by the data. Here, we describe common pitfalls that characterize such studies, and provide practical guidance to improve landscape genetic investigation, with careful consideration of inferential limits, scale, replication, and the ecological and evolutionary context of spatial genetic patterns. Ultimately, the utility of landscape genetics will depend on translating the relationship between gene flow and landscape features into an understanding of long-term population outcomes. We hope the perspective presented here will steer landscape genetics down a more scientifically sound and productive path, garnering a field that is as informative in the future as it is popular now.

A Comparative Analysis of Genetic Diversity and Structure in Jaguars (Panthera onca), Pumas (Puma concolor), and Ocelots (Leopardus pardalis) in Fragmented Landscapes of a Critical Mesoamerican Linkage Zone

With increasing anthropogenic impact and landscape change, terrestrial carnivore populations are becoming more fragmented. Thus, it is crucial to genetically monitor wild carnivores and quantify changes in genetic diversity and gene flow in response to these threats. This study combined the use of scat detector dogs and molecular scatology to conduct the first genetic study on wild populations of multiple Neotropical felids coexisting across a fragmented landscape in Belize, Central America. We analyzed data from 14 polymorphic microsatellite loci in 1053 scat samples collected from wild jaguars (Panthera onca), pumas (Puma concolor), and ocelots (Leopardus pardalis). We assessed levels of genetic diversity, defined potential genetic clusters, and examined gene flow for the three target species on a countrywide scale using a combination of individual- and population-based analyses. Wild felids in Belize showed moderate levels of genetic variation, with jaguars having the lowest diversity estimates (HE = 0.57 ± 0.02; AR = 3.36 ± 0.09), followed by pumas (HE = 0.57 ± 0.08; AR = 4.20 ± 0.16), and ocelots (HE = 0.63 ± 0.03; AR = 4.16 ± 0.08). We observed low to moderate levels of genetic differentiation for all three target species, with jaguars showing the lowest degree of genetic subdivision across the country, followed by ocelots and pumas. Although levels of genetic diversity and gene flow were still fairly high, we detected evidence of fine-scale genetic subdivision, indicating that levels of genetic connectivity for wild felids in Belize are likely to decrease if habitat loss and fragmentation continue at the current rate. Our study demonstrates the value of understanding fine-scale patterns of gene flow in multiple co-occurring felid species of conservation concern, which is vital for wildlife movement corridor planning and prioritizing future conservation and management efforts within human-impacted landscapes.

How Do Species Interactions Affect Evolutionary Dynamics Across Whole Communities?

Theories of how species evolve in changing environments mostly consider single species in isolation or pairs of interacting species. Yet all organisms live in diverse communities containing many hundreds of species. This review discusses how species interactions influence the evolution of constituent species across whole communities. When species interactions are weak or inconsistent, evolutionary dynamics should be predictable by factors identified by single-species theory. Stronger species interactions, however, can alter evolutionary outcomes and either dampen or promote evolution of constituent species depending on the number of species and the distribution of interaction strengths across the interaction network. Genetic interactions, such as horizontal gene transfer, might also affect evolutionary outcomes. These evolutionary mechanisms in turn affect whole-community properties, such as the level of ecosystem functioning. Successful management of both ecosystems and focal species requires new understanding of evolutionary interactions across whole communities.

Ecology has contrasting effects on genetic variation within species versus rates of molecular evolution across species in water beetles

Comparative analysis is a potentially powerful approach to study the effects of ecological traits on genetic variation and rate of evolution across species. However, the lack of suitable datasets means that comparative studies of correlates of genetic traits across an entire clade have been rare. Here, we use a large DNA-barcode dataset (5062 sequences) of water beetles to test the effects of species ecology and geographical distribution on genetic variation within species and rates of molecular evolution across species. We investigated species traits predicted to influence their genetic characteristics, such as surrogate measures of species population size, latitudinal distribution and habitat types, taking phylogeny into account. Genetic variation of cytochrome oxidase I in water beetles was positively correlated with occupancy (numbers of sites of species presence) and negatively with latitude, whereas substitution rates across species depended mainly on habitat types, and running water specialists had the highest rate. These results are consistent with theoretical predictions from nearly-neutral theories of evolution, and suggest that the comparative analysis using large databases can give insights into correlates of genetic variation and molecular evolution.

Temporal population genetic instability in range-edge western toads, Anaxyrus boreas

In this article, we address the temporal stability of population genetic structure in a range-edge population that is undergoing continual, short-distance colonization events. We sampled western toad, Anaxyrus boreas, breeding populations over 2 seasons near their northern range limit in southeast Alaska. We sampled 20 ponds each during the summers of 2008 and 2009, with 14 ponds sampled in both summers. We found considerable turnover in the population genetic relationships among ponds in those 2 seasons, as well as biologically meaningful genetic differentiation between years within some ponds. We found relatively consistent relationships between major population centers, whereas the relationships between the central ponds and smaller, outlying populations differed year to year. This finding indicates that multiple years of genetic sampling may be important for understanding the genetic landscape of some populations.

Physiological vagility and its relationship to dispersal and neutral genetic heterogeneity in vertebrates

Vagility is the inherent power of movement by individuals. Vagility and the available duration of movement determine the dispersal distance individuals can move to interbreed, which affects the fine-scale genetic structure of vertebrate populations. Vagility and variation in population genetic structure are normally explained by geographic variation and not by the inherent power of movement by individuals. We present a new, quantitative definition for physiological vagility that incorporates aerobic capacity, body size, body temperature and the metabolic cost of transport, variables that are independent of the physical environment. Physiological vagility is the speed at which an animal can move sustainably based on these parameters. This meta-analysis tests whether this definition of physiological vagility correlates with empirical data for maximal dispersal distances and measured microsatellite genetic differentiation with distance {[F(ST)/[1-F(ST))]/ln distance} for amphibians, reptiles, birds and mammals utilizing three locomotor modes (running, flying, swimming). Maximal dispersal distance and physiological vagility increased with body mass for amphibians, reptiles and mammals utilizing terrestrial movement. The relative slopes of these relationships indicate that larger individuals require longer movement durations to achieve maximal dispersal distances. Both physiological vagility and maximal dispersal distance were independent of body mass for flying vertebrates. Genetic differentiation with distance was greatest for terrestrial locomotion, with amphibians showing the greatest mean and variance in differentiation. Flying birds, flying mammals and swimming marine mammals showed the least differentiation. Mean physiological vagility of different groups (class and locomotor mode) accounted for 98% of the mean variation in genetic differentiation with distance in each group. Genetic differentiation with distance was not related to body mass. The physiological capacity for movement (physiological vagility) quantitatively predicts genetic isolation by distance in the vertebrates examined.

Contemporary effective population and metapopulation size (N e and meta-N e): comparison among three salmonids inhabiting a fragmented system and differing in gene flow and its asymmetries

We estimated local and metapopulation effective sizes ([Formula: see text] and meta-[Formula: see text]) for three coexisting salmonid species (Salmo salar, Salvelinus fontinalis, Salvelinus alpinus) inhabiting a freshwater system comprising seven interconnected lakes. First, we hypothesized that [Formula: see text] might be inversely related to within-species population divergence as reported in an earlier study (i.e., FST: S. salar> S. fontinalis> S. alpinus). Using the approximate Bayesian computation method implemented in ONeSAMP, we found significant differences in [Formula: see text] ([Formula: see text]) between species, consistent with a hierarchy of adult population sizes ([Formula: see text]). Using another method based on a measure of linkage disequilibrium (LDNE: [Formula: see text]), we found more finite [Formula: see text] values for S. salar than for the other two salmonids, in line with the results above that indicate that S. salar exhibits the lowest [Formula: see text] among the three species. Considering subpopulations as open to migration (i.e., removing putative immigrants) led to only marginal and non-significant changes in [Formula: see text], suggesting that migration may be at equilibrium between genetically similar sources. Second, we hypothesized that meta-[Formula: see text] might be significantly smaller than the sum of local [Formula: see text]s (null model) if gene flow is asymmetric, varies among subpopulations, and is driven by common landscape features such as waterfalls. One 'bottom-up' or numerical approach that explicitly incorporates variable and asymmetric migration rates showed this very pattern, while a number of analytical models provided meta-[Formula: see text] estimates that were not significantly different from the null model or from each other. Our study of three species inhabiting a shared environment highlights the importance and utility of differentiating species-specific and landscape effects, not only on dispersal but also in the demography of wild populations as assessed through local [Formula: see text]s and meta-[Formula: see text]s and their relevance in ecology, evolution and conservation.

Divergent landscape effects on population connectivity in two co-occurring amphibian species

The physical and environmental attributes of landscapes often shape patterns of population connectivity by influencing dispersal and gene flow. Landscape effects on movement are typically evaluated for single species. However, inferences from multiple species are required for multi-species management strategies increasingly being applied in conservation. In this study, I compared the spatial genetic patterns of two amphibian species across the northeastern United States and estimated the influence of specific landscape features on the observed genetic structure. The spotted salamander (Ambystoma maculatum) and wood frog (Rana sylvatica) share many ecological attributes related to habitat use, phenology and site fidelity. However, I hypothesized that important differences in their movement patterns and life history would create distinct genetic patterns for each species. Using 14 microsatellite loci, I tested for differences in the level of genetic differentiation between the two species across 22 breeding ponds. The effects of eight landscape features were also estimated by evaluating 32 landscape resistance models. Spotted salamanders exhibited significantly higher genetic differentiation than wood frogs. Different landscape features were also identified as potential drivers of the genetic patterns in each species, with little overlap in model support between species. Collectively, these results provide strong evidence that these two amphibian species interact with the landscape in measurably different ways. The distinct genetic patterns observed are consistent with key differences in movement ability and life history between A. maculatum and R. sylvatica. These results highlight the importance of considering more than one species when assessing the impacts of the landscape matrix on population connectivity, even for ecologically similar species within the same habitats.

Machine learning identifies specific habitats associated with genetic connectivity in Hyla squirella

The goal of this study was to identify and differentiate the influence of multiple habitat types that span a spectrum of suitability for Hyla squirella, a widespread frog species that occurs in a broad range of habitat types. We collected microsatellite data from 675 samples representing 20 localities from the southeastern USA and used machine-learning methodologies to identify significant habitat features associated with genetic structure. In simulation, we confirm that our machine-learning algorithm can successfully identify landscape features responsible for generating between-population genetic differentiation, suggesting that it can be a useful hypothesis-generating tool for landscape genetics. In our study system, we found that H. squirella were spatially structured and models including specific habitat types (i.e. upland oak forest and urbanization) consistently explained more variation in genetic distance (median pR(2)  = 47.78) than spatial distance alone (median pR(2)  = 23.81). Moreover, we estimate the relative importance that spatial distance, upland oak and urbanized habitat have in explaining genetic structure of H. squirella. We discuss how these habitat types may mechanistically facilitate dispersal in H. squirella. This study provides empirical support for the hypothesis that habitat-use can be an informative correlate of genetic differentiation, even for species that occur in a wide range of habitats.

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.

A garter snake transcriptome: pyrosequencing, de novo assembly, and sex-specific differences

Background

The reptiles, characterized by both diversity and unique evolutionary adaptations, provide a comprehensive system for comparative studies of metabolism, physiology, and development. However, molecular resources for ectothermic reptiles are severely limited, hampering our ability to study the genetic basis for many evolutionarily important traits such as metabolic plasticity, extreme longevity, limblessness, venom, and freeze tolerance. Here we use massively parallel sequencing (454 GS-FLX Titanium) to generate a transcriptome of the western terrestrial garter snake (Thamnophis elegans) with two goals in mind. First, we develop a molecular resource for an ectothermic reptile; and second, we use these sex-specific transcriptomes to identify differences in the presence of expressed transcripts and potential genes of evolutionary interest.

Results

Using sex-specific pools of RNA (one pool for females, one pool for males) representing 7 tissue types and 35 diverse individuals, we produced 1.24 million sequence reads, which averaged 366 bp in length after cleaning. Assembly of the cleaned reads from both sexes with NEWBLER and MIRA resulted in 96,379 contigs containing 87% of the cleaned reads. Over 34% of these contigs and 13% of the singletons were annotated based on homology to previously identified proteins. From these homology assignments, additional clustering, and ORF predictions, we estimate that this transcriptome contains ~13,000 unique genes that were previously identified in other species and over 66,000 transcripts from unidentified protein-coding genes. Furthermore, we use a graph-clustering method to identify contigs linked by NEWBLER-split reads that represent divergent alleles, gene duplications, and alternatively spliced transcripts. Beyond gene identification, we identified 95,295 SNPs and 31,651 INDELs. From these sex-specific transcriptomes, we identified 190 genes that were only present in the mRNA sequenced from one of the sexes (84 female-specific, 106 male-specific), and many highly variable genes of evolutionary interest.

Conclusions

This is the first large-scale, multi-organ transcriptome for an ectothermic reptile. This resource provides the most comprehensive set of EST sequences available for an individual ectothermic reptile species, increasing the number of snake ESTs 50-fold. We have identified genes that appear to be under evolutionary selection and those that are sex-specific. This resource will assist studies on gene expression and comparative genomics, and will facilitate the study of evolutionarily important traits at the molecular level.

Landscape genetics: where are we now?

Landscape genetics has seen rapid growth in number of publications since the term was coined in 2003. An extensive literature search from 1998 to 2008 using keywords associated with landscape genetics yielded 655 articles encompassing a vast array of study organisms, study designs and methodology. These publications were screened to identify 174 studies that explicitly incorporated at least one landscape variable with genetic data. We systematically reviewed this set of papers to assess taxonomic and temporal trends in: (i) geographic regions studied; (ii) types of questions addressed; (iii) molecular markers used; (iv) statistical analyses used; and (v) types and nature of spatial data used. Overall, studies have occurred in geographic regions proximal to developed countries and more commonly in terrestrial vs. aquatic habitats. Questions most often focused on effects of barriers and/or landscape variables on gene flow. The most commonly used molecular markers were microsatellites and amplified fragment length polymorphism (AFLPs), with AFLPs used more frequently in plants than animals. Analysis methods were dominated by Mantel and assignment tests. We also assessed differences among journals to evaluate the uniformity of reporting and publication standards. Few studies presented an explicit study design or explicit descriptions of spatial extent. While some landscape variables such as topographic relief affected most species studied, effects were not universal, and some species appeared unaffected by the landscape. Effects of habitat fragmentation were mixed, with some species altering movement paths and others unaffected. Taken together, although some generalities emerged regarding effects of specific landscape variables, results varied, thereby reinforcing the need for species-specific work. We conclude by: highlighting gaps in knowledge and methodology, providing guidelines to authors and reviewers of landscape genetics studies, and suggesting promising future directions of inquiry.

Comparative landscape genetics of two pond-breeding amphibian species in a highly modified agricultural landscape

Evaluating fine-scale population structure of multiple species in the same landscape increases our ability to identify common patterns as well as discern ecological differences among species' landscape genetic relationships. In the Palouse bioregion of northern Idaho, USA, 99% of the native prairie has been converted to nonirrigated agriculture and exotic grasslands. Columbia spotted frogs (Rana luteiventris) and long-toed salamanders (Ambystoma macrodactylum) in this area breed almost entirely in artificial ponds on private land. We used genetic distances (F(ST) and D(c)) derived from eight microsatellite loci in 783 samples to evaluate the relationships among sympatric breeding populations (N = 20 and 26) of these species in a 213-km(2) landscape. Both species showed a pattern of isolation by distance that was not improved when distance was measured along drainages instead of topographically corrected straight lines (P < 0.01). After testing for autocorrelation among genetic distances, we used an information theoretic approach to model landscape resistance based on slope, soil type, solar insolation, and land cover, and multi-model inference to rank the resistance of landscape surfaces to dispersal (represented by genetic distance). For both species, urban and rural developed land cover provided the highest landscape resistances. Resistance values for long-toed salamanders followed a moisture gradient where forest provided the least resistance, while agriculture and shrub/clearcut provided the least resistance for Columbia spotted frogs. Comparative landscape genetics can be a powerful tool for detecting similarities and differences between codistributed species, and resulting models can be used to predict species-specific responses to landscape change.

Landscape genetics of high mountain frog metapopulations

Explaining functional connectivity among occupied habitats is crucial for understanding metapopulation dynamics and species ecology. Landscape genetics has primarily focused on elucidating how ecological features between observations influence gene flow. Functional connectivity, however, may be the result of both these between-site (landscape resistance) landscape characteristics and at-site (patch quality) landscape processes that can be captured using network based models. We test hypotheses of functional connectivity that include both between-site and at-site landscape processes in metapopulations of Columbia spotted frogs (Rana luteiventris) by employing a novel justification of gravity models for landscape genetics (eight microsatellite loci, 37 sites, n = 441). Primarily used in transportation and economic geography, gravity models are a unique approach as flow (e.g. gene flow) is explained as a function of three basic components: distance between sites, production/attraction (e.g. at-site landscape process) and resistance (e.g. between-site landscape process). The study system contains a network of nutrient poor high mountain lakes where we hypothesized a short growing season and complex topography between sites limit R. luteiventris gene flow. In addition, we hypothesized production of offspring is limited by breeding site characteristics such as the introduction of predatory fish and inherent site productivity. We found that R. luteiventris connectivity was negatively correlated with distance between sites, presence of predatory fish (at-site) and topographic complexity (between-site). Conversely, site productivity (as measured by heat load index, at-site) and growing season (as measured by frost-free period between-sites) were positively correlated with gene flow. The negative effect of predation and positive effect of site productivity, in concert with bottleneck tests, support the presence of source-sink dynamics. In conclusion, gravity models provide a powerful new modelling approach for examining a wide range of both basic and applied questions in landscape genetics.

Quantifying Bufo boreas connectivity in Yellowstone National Park with landscape genetics

A major objective of ecology is to understand how ecological processes limit population connectivity and species' distributions. By spatially quantifying ecological components driving functional connectivity, we can understand why some locally suitable habitats are unoccupied, resulting in observed discontinuities in distribution. However, estimating connectivity may be difficult due to population stochasticity and violations of assumptions of parametric statistics. To address these issues, we present a novel application of Random Forests to landscape genetic data. We address the effects of three key ecological components on Bufo boreas connectivity in Yellowstone National Park: ecological process, scale, and hierarchical organization. Habitat permeability, topographic morphology, and temperature-moisture regime are all significant ecological processes associated with B. boreas connectivity. Connectivity was influenced by growing-season precipitation, 1988 Yellowstone fires, cover, temperature, impervious surfaces (roads and development), and topographic complexity (56% variation explained). We found that habitat permeability generally operates on fine scales, while topographic morphology and temperature-moisture regime operate across multiple scales, thus demonstrating the importance of cross-scale analysis for ecological interpretation. In a hierarchical analysis, we were able to explain more variation within genetic clusters as identified using Structure (a Bayesian algorithm) (74%; dispersal cover, growing-season precipitation, impervious surfaces) as opposed to between genetic clusters (45%; ridgelines, hot, dry slopes, length of hot season, and annual precipitation). Finally, the analytical methods we developed are powerful and can be applied to any species or system with appropriate landscape and genetic data.

Habitat continuity and geographic distance predict population genetic differentiation in giant kelp

Isolation by distance (IBD) models are widely used to predict levels of genetic connectivity as a function of Euclidean distance, and although recent studies have used GIS-landscape ecological approaches to improve the predictability of spatial genetic structure, few if any have addressed the effect of habitat continuity on gene flow. Landscape effects on genetic connectivity are even less understood in marine populations, where habitat mapping is particularly challenging. In this study, we model spatial genetic structure of a habitat-structuring species, the giant kelp Macrocystis pyrifera, using highly variable microsatellite markers. GIS mapping was used to characterize habitat continuity and distance between sampling sites along the mainland coast of the Santa Barbara Channel, and their roles as predictors of genetic differentiation were evaluated. Mean dispersal distance (sigma) and effective population size (Ne) were estimated by comparing our IBD slope with those from simulations incorporating habitat continuity and spore dispersal characteristics of the study area. We found an allelic richness of 7-50 alleles/locus, which to our knowledge is the highest reported for macroalgae. The best regression model relating genetic distance to habitat variables included both geographic distance and habitat continuity, which were respectively, positively and negatively related to genetic distance. Our results provide strong support for a dependence of gene flow on both distance and habitat continuity and elucidate the combination of Ne and a that explained genetic differentiation.

Interaction of landscape and life history attributes on genetic diversity, neutral divergence and gene flow in a pristine community of salmonids

Landscape genetics holds promise for the forecasting of spatial patterns of genetic diversity based on key environmental features. Yet, the degree to which inferences based on single species can be extended to whole communities is not fully understood. We used a pristine and spatially structured community of three landlocked salmonids (Salvelinus fontinalis, Salmo salar, and Salvelinus alpinus) from Gros Morne National Park (Newfoundland, Canada) to test several predictions on the interacting effects of landscape and life history variation on genetic diversity, neutral divergence, and gene flow (m, migration rate). Landscape factors consistently influenced multispecies genetic patterns: (i) waterfalls created strong dichotomies in genetic diversity and divergence between populations above and below them in all three salmonids; (ii) contemporary m decreased with waterway distance in all three species, while neutral genetic divergence (theta) increased with waterway distance, albeit in only two taxa; (iii) river flow generally produced downstream-biased m between populations when waterfalls separated these, but not otherwise. In contrast, we expected differential life history to result in a hierarchy of neutral divergence (S. salar > S. fontinalis > S. alpinus) based on disparities in dispersal abilities and population size from previous mark-recapture studies. Such hierarchy additionally matched varying degrees of spatial genetic structure among species revealed through individual-based analyses. We conclude that, whereas key landscape attributes hold power to predict multispecies genetic patterns in equivalent communities, they are likely to interact with species-specific life history attributes such as dispersal, demography, and ecology, which will in turn affect holistic conservation strategies.

Contrasting patterns of spatial genetic structure in sympatric rock-dwelling cichlid fishes

The cichlid fish of Lake Tanganyika in Eastern Africa are a celebrated example of both ecological and species diversification. Because population subdivision is likely to play an important role in the speciation process, understanding how habitat features interact with species' demographic, behavioral and ecological attributes to influence gene flow and population divergence may help explain the causes of high species richness in this and other systems. Here, we test the roles of isolation-by-habitat and isolation-by-distance in generating fine-scale population genetic structure in three sympatric species of habitat-restricted cichlids in Lake Tanganyika. Using multilocus microsatellite genotypes, we contrast patterns of population differentiation in these habitat specialists along a mosaic coastline of both favorable and unfavorable habitat. Despite their close phylogenetic relationship and shared habitat affinity, these species show striking differences in their pattern of genetic subdivision within the same geographical region, suggesting substantially different patterns of gene flow. In particular, two trophically specialized species exhibit much more restricted gene flow over sandy habitat than a trophically opportunistic species. This result suggests that ecological and behavioral traits have a strong influence on the scale and degree of population subdivision, a finding that has potentially important implications for understanding differential propensities for diversification among lineages and phylogenetic patterns of diversity.

The effects of maternal corticosterone levels on offspring behavior in fast- and slow-growth garter snakes (Thamnophis elegans)

During embryonic development, viviparous offspring are exposed to maternally circulating hormones. Maternal stress increases offspring exposure to corticosterone and this hormonal exposure has the potential to influence developmental, morphological and behavioral traits of the resulting offspring. We treated pregnant female garter snakes (Thamnophis elegans) with low levels of corticosterone after determining both natural corticosterone levels in the field and pre-treatment levels upon arrival in the lab. Additional measurements of plasma corticosterone were taken at days 1, 5, and 10 during the 10-day exposure, which occurred during the last third of gestation (of 4-month gestation). These pregnant snakes were from replicate populations of fast- and slow-growth ecotypes occurring in Northern California, with concomitant short and long lifespans. Field corticosterone levels of pregnant females of the slow-growth ecotype were an order of magnitude higher than fast-growth dams. In the laboratory, corticosterone levels increased over the 10 days of corticosterone manipulation for animals of both ecotypes, and reached similar plateaus for both control and treated dams. Despite similar plasma corticosterone levels in treated and control mothers, corticosterone-treated dams produced more stillborn offspring and exhibited higher total reproductive failure than control dams. At one month of age, offspring from fast-growth females had higher plasma corticosterone levels than offspring from slow-growth females, which is opposite the maternal pattern. Offspring from corticosterone-treated mothers, although unaffected in their slither speed, exhibited changes in escape behaviors and morphology that were dependent upon maternal ecotype. Offspring from corticosterone-treated fast-growth females exhibited less anti-predator reversal behavior; offspring from corticosterone-treated slow-growth females exhibited less anti-predator tail lashing behavior.

Adaptive divergence within and between ecotypes of the terrestrial garter snake, Thamnophis elegans, assessed with F(ST)-Q(ST) comparisons

Populations of the terrestrial garter snake (Thamnophis elegans) around Eagle Lake in California exhibit dramatic ecotypic differentiation in life history, colouration and morphology across distances as small as a few kilometres. We assayed the role of selection in ecotypic differentiation in T. elegans using F(ST)-Q(ST) analysis and identified selective agents using direct and indirect observations. We extended the conventional implementation of the F(ST)-Q(ST) approach by using three-level analyses of genetic and phenotypic variance to assess the role of selection in differentiating populations both within and between ecotypes. These results suggest that selection has driven differentiation between as well as within ecotypes, and in the presence of moderate to high gene flow. Our findings are discussed in the context of previous correlational selection analyses which revealed stabilizing and correlational selection for some of the traits examined.