Genetic adaptation despite high gene flow in a range-expanding population.

Signals of natural selection can be quickly eroded in high gene flow systems, curtailing efforts to understand how and when genetic adaptation occurs in the ocean. This long-standing, unresolved topic in ecology and evolution has renewed importance because changing environmental conditions are driving range expansions that may necessitate rapid evolutionary responses. One example occurs in Kellet's whelk (Kelletia kelletii), a common subtidal gastropod with an ~40- to 60-day pelagic larval duration that expanded their biogeographic range northwards in the 1970s by over 300 km. To test for genetic adaptation, we performed a series of experimental crosses with Kellet's whelk adults collected from their historical (HxH) and recently expanded range (ExE), and conducted RNA-Seq on offspring that we reared in a common garden environment. We identified 2770 differentially expressed genes (DEGs) between 54 offspring samples with either only historical range (HxH offspring) or expanded range (ExE offspring) ancestry. Using SNPs called directly from the DEGs, we assigned samples of known origin back to their range of origin with unprecedented accuracy for a marine species (92.6% and 94.5% for HxH and ExE offspring, respectively). The SNP with the highest predictive importance occurred on triosephosphate isomerase (TPI), an essential metabolic enzyme involved in cold stress response. TPI was significantly upregulated and contained a non-synonymous mutation in the expanded range. Our findings pave the way for accurately identifying patterns of dispersal, gene flow and population connectivity in the ocean by demonstrating that experimental transcriptomics can reveal mechanisms for how marine organisms respond to changing environmental conditions.

Detection of contractility changes in the heart from arterial blood pressure data using symmetric Projection Attractor Reconstruction.

The potential for unintended drug induced changes in cardiac contractility is a major concern in medicines development. Whilst direct left ventricular pressure (LVP) measurement is the gold standard for measuring cardiac contractility in vivo, it is resource intensive and poses a welfare burden on research animals. In contrast, arterial blood pressure (BP) measurement has fewer challenges. Symmetric Projection Attractor Reconstruction (SPAR) is a signal processing technique which transforms physiological time-series signals into a corresponding visual image ('attractor'), amplifying morphology changes within physiological waveforms. It was hypothesized that SPAR analysis of BP signals would provide a surrogate measure of cardiac contractility by specifically amplifying the maximum slope of the systolic upstroke. BP (abdominal aorta) signals obtained from beagle dogs, treated with positive and negative inotropes, were retrospectively analysed to identify signal features that correlated with the maximum upslope of the LVP signal from simultaneously acquired LVP recordings. SPAR transformation of BP signals quantified drug induced changes in the maximum slope of the systolic upstroke. We identified key SPAR metrics that provided >0.8 correlation with the LVP maximum upslope, outperforming the BP systolic upstroke alone. This was observed for all 4 different drugs, doses and time points evaluated across studies. Thus, we conclude that the SPAR measures derived from the BP signal could be used as a first pass in vivo screen to flag any risk of drug induced changes in cardiac contractility during the conduct of non-clinical medicines development, potentially reducing or replacing the need to perform direct left ventricular measurements.

Next-generation data filtering in the genomics era.

Genomic data are ubiquitous across disciplines, from agriculture to biodiversity, ecology, evolution and human health. However, these datasets often contain noise or errors and are missing information that can affect the accuracy and reliability of subsequent computational analyses and conclusions. A key step in genomic data analysis is filtering - removing sequencing bases, reads, genetic variants and/or individuals from a dataset - to improve data quality for downstream analyses. Researchers are confronted with a multitude of choices when filtering genomic data; they must choose which filters to apply and select appropriate thresholds. To help usher in the next generation of genomic data filtering, we review and suggest best practices to improve the implementation, reproducibility and reporting standards for filter types and thresholds commonly applied to genomic datasets. We focus mainly on filters for minor allele frequency, missing data per individual or per locus, linkage disequilibrium and Hardy-Weinberg deviations. Using simulated and empirical datasets, we illustrate the large effects of different filtering thresholds on common population genetics statistics, such as Tajima's D value, population differentiation (FST), nucleotide diversity (π) and effective population size (Ne).

Genomic DNA extraction optimization and validation for genome sequencing using the marine gastropod Kellet's whelk.

Next-generation sequencing technologies, such as Nanopore MinION, Illumina Hiseq and Novaseq, and PacBio Sequel II, hold immense potential for advancing genomic research on non-model organisms, including the vast majority of marine species. However, application of these technologies to marine invertebrate species is often impeded by challenges in extracting and purifying their genomic DNA due to high polysaccharide content and other secondary metabolites. In this study, we help resolve this issue by developing and testing DNA extraction protocols for Kellet's whelk (Kelletia kelletii), a subtidal gastropod with ecological and commercial importance, by comparing four DNA extraction methods commonly used in marine invertebrate studies. In our comparison of extraction methods, the Salting Out protocol was the least expensive, produced the highest DNA yields, produced consistent high DNA quality, and had low toxicity. We validated the protocol using an independent set of tissue samples, then applied it to extract high-molecular-weight (HMW) DNA from over three thousand Kellet's whelk tissue samples. The protocol demonstrated scalability and, with added clean-up, suitability for RAD-seq, GT-seq, as well as whole genome sequencing using both long read (ONT MinION) and short read (Illumina NovaSeq) sequencing platforms. Our findings offer a robust and versatile DNA extraction and clean-up protocol for supporting genomic research on non-model marine organisms, to help mediate the under-representation of invertebrates in genomic studies.

Rapid genetic adaptation to a novel ecosystem despite a large founder event.

Introduced and invasive species make excellent natural experiments for investigating rapid evolution. Here, we describe the effects of genetic drift and rapid genetic adaptation in pink salmon (Oncorhynchus gorbuscha) that were accidentally introduced to the Great Lakes via a single introduction event 31 generations ago. Using whole-genome resequencing for 134 fish spanning five sample groups across the native and introduced range, we estimate that the source population's effective population size was 146,886 at the time of introduction, whereas the founding population's effective population size was just 72-a 2040-fold decrease. As expected with a severe founder event, we show reductions in genome-wide measures of genetic diversity, specifically a 37.7% reduction in the number of SNPs and an 8.2% reduction in observed heterozygosity. Despite this decline in genetic diversity, we provide evidence for putative selection at 47 loci across multiple chromosomes in the introduced populations, including missense variants in genes associated with circadian rhythm, immunological response and maturation, which match expected or known phenotypic changes in the Great Lakes. For one of these genes, we use a species-specific agent-based model to rule out genetic drift and conclude our results support a strong response to selection occurring in a period gene (per2) that plays a predominant role in determining an organism's daily clock, matching large day length differences experienced by introduced salmon during important phenological periods. Together, these results inform how populations might evolve rapidly to new environments, even with a small pool of standing genetic variation.

Dispersive currents explain patterns of population connectivity in an ecologically and economically important fish.

How to identify the drivers of population connectivity remains a fundamental question in ecology and evolution. Answering this question can be challenging in aquatic environments where dynamic lake and ocean currents coupled with high levels of dispersal and gene flow can decrease the utility of modern population genetic tools. To address this challenge, we used RAD-Seq to genotype 959 yellow perch (Perca flavescens), a species with an ~40-day pelagic larval duration (PLD), collected from 20 sites circumscribing Lake Michigan. We also developed a novel, integrative approach that couples detailed biophysical models with eco-genetic agent-based models to generate "predictive" values of genetic differentiation. By comparing predictive and empirical values of genetic differentiation, we estimated the relative contributions for known drivers of population connectivity (e.g., currents, behavior, PLD). For the main basin populations (i.e., the largest contiguous portion of the lake), we found that high gene flow led to low overall levels of genetic differentiation among populations (F ST = 0.003). By far the best predictors of genetic differentiation were connectivity matrices that were derived from periods of time when there were strong and highly dispersive currents. Thus, these highly dispersive currents are driving the patterns of population connectivity in the main basin. We also found that populations from the northern and southern main basin are slightly divergent from one another, while those from Green Bay and the main basin are highly divergent (F ST = 0.11). By integrating biophysical and eco-genetic models with genome-wide data, we illustrate that the drivers of population connectivity can be identified in high gene flow systems.

Negative frequency dependent selection unites ecology and evolution.

From genes to communities, understanding how diversity is maintained remains a fundamental question in biology. One challenging to identify, yet potentially ubiquitous, mechanism for the maintenance of diversity is negative frequency dependent selection (NFDS), which occurs when entities (e.g., genotypes, life history strategies, species) experience a per capita reduction in fitness with increases in relative abundance. Because NFDS allows rare entities to increase in frequency while preventing abundant entities from excluding others, we posit that negative frequency dependent selection plays a central role in the maintenance of diversity. In this review, we relate NFDS to coexistence, identify mechanisms of NFDS (e.g., mutualism, predation, parasitism), review strategies for identifying NFDS, and distinguish NFDS from other mechanisms of coexistence (e.g., storage effects, fluctuating selection). We also emphasize that NFDS is a key place where ecology and evolution intersect. Specifically, there are many examples of frequency dependent processes in ecology, but fewer cases that link this process to selection. Similarly, there are many examples of selection in evolution, but fewer cases that link changes in trait values to negative frequency dependence. Bridging these two well-developed fields of ecology and evolution will allow for mechanistic insights into the maintenance of diversity at multiple levels.

Large genetic divergence underpins cryptic local adaptation across ecological and evolutionary gradients.

Environmentally covarying local adaptation is a form of cryptic local adaptation in which the covariance of the genetic and environmental effects on a phenotype obscures the divergence between locally adapted genotypes. Here, we systematically document the magnitude and drivers of the genetic effect (VG) for two forms of environmentally covarying local adaptation: counter- and cogradient variation. Using a hierarchical Bayesian meta-analysis, we calculated the overall effect size of VG as 1.05 and 2.13 for populations exhibiting countergradient or cogradient variation, respectively. These results indicate that the genetic contribution to phenotypic variation represents a 1.05 to 2.13 s.d. change in trait value between the most disparate populations depending on if populations are expressing counter- or cogradient variation. We also found that while there was substantial variance among abiotic and biotic covariates, the covariates with the largest mean effects were temperature (2.41) and gamete size (2.81). Our results demonstrate the pervasiveness and large genetic effects underlying environmentally covarying local adaptation in wild populations and highlight the importance of accounting for these effects in future studies.

Evolutionary rescue in host-pathogen systems.

Natural populations encounter a variety of threats that can increase their risk of extinction. Populations can avoid extinction through evolutionary rescue (ER), which occurs when an adaptive, genetic response to selection allows a population to recover from an environmental change that would otherwise cause extinction. While the traditional framework for ER was developed with abiotic risk factors in mind, ER may also occur in response to a biotic source of demographic change, such as the introduction of a novel pathogen. We first describe how ER in response to a pathogen differs from the traditional ER framework; density-dependent transmission, pathogen evolution, and pathogen extinction can change the strength of selection imposed by a pathogen and make host population persistence more likely. We also discuss several variables that affect traditional ER (abundance, genetic diversity, population connectivity, and community composition) that also directly affect disease risk resulting in diverse outcomes for ER in host-pathogen systems. Thus, generalizations developed in studies of traditional ER may not be relevant for ER in response to the introduction of a pathogen. Incorporating pathogens into the framework of ER will lead to a better understanding of how and when populations can avoid extinction in response to novel pathogens.

Rapid genetic adaptation to recently colonized environments is driven by genes underlying life history traits.

BACKGROUND: Uncovering the mechanisms underlying rapid genetic adaptation can provide insight into adaptive evolution and shed light on conservation, invasive species control, and natural resource management. However, it can be difficult to experimentally explore rapid adaptation due to the challenges associated with propagating and maintaining species in captive environments for long periods of time. By contrast, many introduced species have experienced strong selection when colonizing environments that differ substantially from their native range and thus provide a "natural experiment" for studying rapid genetic adaptation. One such example occurred when sea lamprey (Petromyzon marinus), native to the northern Atlantic, naturally migrated into Lake Champlain and expanded their range into the Great Lakes via man-made shipping canals. RESULTS: Utilizing 368,886 genome-wide single nucleotide polymorphisms (SNPs), we calculated genome-wide levels of genetic diversity (i.e., heterozygosity and π) for sea lamprey collected from native (Connecticut River), native but recently colonized (Lake Champlain), and invasive (Lake Michigan) populations, assessed genetic differentiation between all populations, and identified candidate genes that responded to selection imposed by the novel environments. We observed a 14 and 24% reduction in genetic diversity in Lake Michigan and Lake Champlain populations, respectively, compared to individuals from the Connecticut River, suggesting that sea lamprey populations underwent a genetic bottleneck during colonization. Additionally, we identified 121 and 43 outlier genes in comparisons between Lake Michigan and Connecticut River and between Lake Champlain and Connecticut River, respectively. Six outlier genes that contained synonymous SNPs in their coding regions and two genes that contained nonsynonymous SNPs may underlie the rapid evolution of growth (i.e., GHR), reproduction (i.e., PGR, TTC25, STARD10), and bioenergetics (i.e., OXCT1, PYGL, DIN4, SLC25A15). CONCLUSIONS: By identifying the genomic basis of rapid adaptation to novel environments, we demonstrate that populations of invasive species can be a useful study system for understanding adaptive evolution. Furthermore, the reduction in genome-wide levels of genetic diversity associated with colonization coupled with the identification of outlier genes underlying key life history traits known to have changed in invasive sea lamprey populations (e.g., growth, reproduction) illustrate the utility in applying genomic approaches for the successful management of introduced species.

Drug reformulation for a neglected disease. The NANOHAT project to develop a safer more effective sleeping sickness drug.

BACKGROUND: Human African trypanosomiasis (HAT or sleeping sickness) is caused by the parasite Trypanosoma brucei sspp. The disease has two stages, a haemolymphatic stage after the bite of an infected tsetse fly, followed by a central nervous system stage where the parasite penetrates the brain, causing death if untreated. Treatment is stage-specific, due to the blood-brain barrier, with less toxic drugs such as pentamidine used to treat stage 1. The objective of our research programme was to develop an intravenous formulation of pentamidine which increases CNS exposure by some 10-100 fold, leading to efficacy against a model of stage 2 HAT. This target candidate profile is in line with drugs for neglected diseases inititative recommendations. METHODOLOGY: To do this, we evaluated the physicochemical and structural characteristics of formulations of pentamidine with Pluronic micelles (triblock-copolymers of polyethylene-oxide and polypropylene oxide), selected candidates for efficacy and toxicity evaluation in vitro, quantified pentamidine CNS delivery of a sub-set of formulations in vitro and in vivo, and progressed one pentamidine-Pluronic formulation for further evaluation using an in vivo single dose brain penetration study. PRINCIPAL FINDINGS: Screening pentamidine against 40 CNS targets did not reveal any major neurotoxicity concerns, however, pentamidine had a high affinity for the imidazoline2 receptor. The reduction in insulin secretion in MIN6 ß-cells by pentamidine may be secondary to pentamidine-mediated activation of ß-cell imidazoline receptors and impairment of cell viability. Pluronic F68 (0.01%w/v)-pentamidine formulation had a similar inhibitory effect on insulin secretion as pentamidine alone and an additive trypanocidal effect in vitro. However, all Pluronics tested (P85, P105 and F68) did not significantly enhance brain exposure of pentamidine. SIGNIFICANCE: These results are relevant to further developing block-copolymers as nanocarriers, improving BBB drug penetration and understanding the side effects of pentamidine.

Incipient resistance to an effective pesticide results from genetic adaptation and the canalization of gene expression.

The resistance of pest species to chemical controls has vast ecological, economic, and societal costs. In most cases, resistance is only detected after spreading throughout an entire population. Detecting resistance in its incipient stages, by comparison, provides time to implement preventative strategies. Incipient resistance can be detected by coupling standard toxicology assays with large-scale gene expression experiments. We apply this approach to a system where an invasive parasite, sea lamprey (Petromyzon marinus), has been treated with the highly effective pesticide 3-trifluoromethyl-4-nitrophenol (TFM) for 60 years. Toxicological experiments revealed that lamprey from treated populations did not have higher survival to TFM exposure than lamprey from untreated populations, demonstrating that full-fledged resistance has not yet evolved. In contrast, we find hundreds of genes differentially expressed in response to TFM in the population with the longest history of exposure, many of which relate to TFM's primary mode of action, the uncoupling of oxidative phosphorylation, and subsequent depletion of ATP. Three genes critical to oxidative phosphorylation, ATP5PB, PLCB1, and NDUFA9, were nearly fixed for alternative alleles in comparisons of SNPs between treated and untreated populations (FST  > 5 SD from the mean). ATP5PB encodes subunit b of ATP synthase and an additional subunit, ATP5F1B, was canalized for high expression in treated populations, but remained plastic in response to TFM treatment in individuals from the untreated population. These combined genomic and transcriptomic results demonstrate that an adaptive, genetic response to TFM is likely driving incipient resistance in a damaging pest species.

Characterization of a Y-specific duplication/insertion of the anti-Mullerian hormone type II receptor gene based on a chromosome-scale genome assembly of yellow perch, Perca flavescens.

Yellow perch, Perca flavescens, is an ecologically and economically important species native to a large portion of the northern United States and southern Canada and is also a promising candidate species for aquaculture. However, no yellow perch reference genome has been available to facilitate improvements in both fisheries and aquaculture management practices. By combining Oxford Nanopore Technologies long-reads, 10X Genomics Illumina short linked reads and a chromosome contact map produced with Hi-C, we generated a high-continuity chromosome-scale yellow perch genome assembly of 877.4 Mb. It contains, in agreement with the known diploid chromosome yellow perch count, 24 chromosome-size scaffolds covering 98.8% of the complete assembly (N50 = 37.4 Mb, L50 = 11). We also provide a first characterization of the yellow perch sex determination locus that contains a male-specific duplicate of the anti-Mullerian hormone type II receptor gene (amhr2by) inserted at the proximal end of the Y chromosome (chromosome 9). Using this sex-specific information, we developed a simple PCR genotyping assay which accurately differentiates XY genetic males (amhr2by+ ) from XX genetic females (amhr2by- ). Our high-quality genome assembly is an important genomic resource for future studies on yellow perch ecology, toxicology, fisheries and aquaculture research. In addition, characterization of the amhr2by gene as a candidate sex-determining gene in yellow perch provides a new example of the recurrent implication of the transforming growth factor beta pathway in fish sex determination, and highlights gene duplication as an important genomic mechanism for the emergence of new master sex determination genes.

Among-family variation in survival and gene expression uncovers adaptive genetic variation in a threatened fish.

Variation in among-family transcriptional responses to different environmental conditions can help to identify adaptive genetic variation, even prior to a selective event. Coupling differential gene expression with formal survival analyses allows for the disentanglement of treatment effects, required for understanding how individuals plastically respond to environmental stressors, from the adaptive genetic variation responsible for differential survival. We combined these two approaches to investigate responses to an emerging conservation issue, thiamine (vitamin B1 ) deficiency, in a threatened population of Atlantic salmon (Salmo salar). Thiamine is an essential vitamin that is increasingly limited in many ecosystems. In Lake Champlain, Atlantic salmon cannot acquire thiamine in sufficient quantities to support natural reproduction; fertilized eggs must be reared in hatcheries and treated with supplemental thiamine. We evaluated transcriptional responses (via RNA sequencing) to thiamine treatment across families and found 3,616 genes differentially expressed between control (no supplemental thiamine) and treatment individuals. Fewer genes changed expression equally across families (i.e., additively) than exhibited genotype × environment interactions in response to thiamine. Differentially expressed genes were related to known physiological effects of thiamine deficiency, including oxidative stress, cardiovascular irregularities and neurological abnormalities. We also identified 1,446 putatively adaptive genes that were strongly associated with among-family survival in the absence of thiamine treatment, many of which related to neurogenesis and visual perception. Our results highlight the utility of coupling RNA sequencing with formal survival analyses to identify candidate genes that underlie the among-family variation in survival required for an adaptive response to natural selection.

Rapid resistance to pesticide control is predicted to evolve in an invasive fish.

Xenobiotic resistance is commonly found in species with short generation times such as bacteria, annual plants, and insects. Nevertheless, the fundamental evolutionary principles that govern the spread of resistance alleles hold true for species with longer generation times. One such example could occur with sea lamprey (Petromyzon marinus), a parasitic invasive species in the Laurentian Great Lakes that decimated native fish populations prior to its control with the pesticide 3-trifluoromethyl-4-nitrophenol (TFM). Since the 1950s, tributaries have been treated annually with TFM, where treatments effectively remove most, but not all, larval sea lamprey. We developed an eco-genetic model of sea lamprey to examine factors affecting the evolution of resistance and found that resistance alleles rapidly rise to fixation after 40-80 years of treatment, despite the species' relatively long generation time (4-7 years). The absence of natal homing allows resistant individuals to spread quickly throughout the entire system, but also makes the early detection of resistance challenging. High costs of resistance and density independent reproduction can delay, but not prevent, the onset of resistance. These results illustrate that sea lamprey have the potential to evolve resistance to their primary control agent in the near future, highlighting the urgent need for alternative controls.

Inbreeding load and inbreeding depression estimated from lifetime reproductive success in a small, dispersal-limited population.

The fitness consequences of inbreeding and the individual behaviors that prevent its detrimental effects can be challenging to document in wild populations. Here, we use field and molecular data from a 17-year study of banner-tailed kangaroo rats (Dipodomys spectabilis) to quantify the relationship between inbreeding, mate kinship, and lifetime reproductive success. Using a pedigree that was reconstructed using genetic and field data within a Bayesian framework (median probability of parental assignment = 0.92, mean pedigree depth = 6 generations), we estimated both inbreeding coefficients and kinship between individuals that produced offspring (mean inbreeding coefficient = 0.07, mean mate kinship = 0.08). We also used the pedigree, in combination with census data, to generate a series of fitness estimates, ranging from survival to reproductive maturity to lifetime reproductive success. We found that the population's inbreeding load was low to moderate (0.98-4.66 haploid lethal equivalents) and increased with the time frame over which fitness was estimated (lowest for survival to maturity, highest for adult-to-adult reproductive success). Fitness decreased with increasing inbreeding coefficients. For example, lifetime reproductive success was reduced by 24% for individuals with inbreeding coefficients greater than twice the population mean. Within full sibling pairs, the sibling with less-related mates produced an average of 30% more offspring over its lifetime. These data further illustrate that inbreeding can have a negative effect on lifetime reproductive success.

Long-term demographic and genetic effects of releasing captive-born individuals into the wild.

Because of continued habitat destruction and species extirpations, the need to use captive breeding for conservation purposes has been increasing steadily. However, the long-term demographic and genetic effects associated with releasing captive-born individuals with varied life histories into the wild remain largely unknown. To address this question, we developed forward-time, agent-based models for 4 species with long-running captive-breeding and release programs: coho salmon (Oncorhynchus kisutch), golden lion tamarin (Leontopithecus rosalia), western toad (Anaxyrus boreas), and Whooping Crane (Grus americana). We measured the effects of supplementation by comparing population size and neutral genetic diversity in supplemented populations to the same characteristics in unaltered populations 100 years after supplementation ended. Releasing even slightly less fit captive-born individuals to supplement wild populations typically resulted in reductions in population sizes and genetic diversity over the long term when the fitness reductions were heritable (i.e., due to genetic adaptation to captivity) and populations continued to be regulated by density-dependent mechanisms over time. Negative effects for species with longer life spans and lower rates of population replacement were smaller than for species with shorter life spans and higher rates of population replacement. Programs that released captive-born individuals over fewer years or that avoided breeding individuals with captive ancestry had smaller reductions in population size and genetic diversity over the long term. Relying on selection in the wild to remove individuals with reduced fitness mitigated some negative demographic effects, but at a substantial cost to neutral genetic diversity. Our results suggest that conservation-focused captive-breeding programs should take measures to prevent even small amounts of genetic adaptation to captivity, quantitatively determine the minimum number of captive-born individuals to release each year, and fully account for the interactions among genetic adaptation to captivity, population regulation, and life-history variation.

Efectos Genéticos y Demográficos a Largo Plazo de la Liberación de Individuos Nacidos en Cautiverio Resumen Debido a la continua destrucción de hábitats y a la extracción de especies, la necesidad de utilizar la reproducción en cautiverio con motivos de conservación ha aumentado constantemente. Sin embargo, los efectos genéticos y demográficos a largo plazo asociados con la liberación de individuos nacidos en cautiverio con historias de vida variadas permanecen en gran parte desconocidos. Para responder a la pregunta anterior desarrollamos modelos de tiempo futuro basados en agentes para cuatro especies con programas de reproducción en cautiverio y liberación con una larga duración: el salmón plateado (Oncorhynchus kisutch), la marmoseta dorada (Leontopithecus rosalia), el sapo occidental (Anaxyrus boreas), y la grulla trompetera (Grus americana). Medimos los efectos de la suplementación al comparar el tamaño poblacional y la diversidad genética neutra en las poblaciones suplementadas con las mismas características en poblaciones sin alteraciones 100 años después de la suplementación. La liberación de individuos criados en cautiverio con una mínima aptitud física como suplemento para las poblaciones silvestres resultó típicamente en la reducción del tamaño poblacional y de la diversidad genética a largo plazo cuando la reducción en la aptitud fue heredable (es decir, debido a la adaptación genética al cautiverio) y las poblaciones siguieron reguladas a lo largo del tiempo por los mecanismos dependientes de la densidad. Los efectos negativos para las especies con ciclos de vida mayores y tasas más bajas de reemplazo poblacional fueron menores que para aquellas especies con ciclos de vida más cortos y tasas más altas de reemplazo poblacional. Los programas que liberaron individuos criados en cautiverio durante menos años o que evitaron reproducir individuos con ascendencia en cautiverio tuvieron reducciones menores en el tamaño poblacional y en la diversidad genética a largo plazo. La dependencia de la selección silvestre para extirpar a los individuos con aptitud física reducida mitigó algunos efectos demográficos negativos, pero a un precio sustancial para la diversidad genética neutra. Nuestros resultados sugieren que los programas de reproducción en cautiverio enfocados en la conservación deberían tomar medidas para prevenir las más mínimas cantidades de adaptación genética al cautiverio, determinar cuantitativamente el número mínimo de individuos nacidos en cautiverio para liberar cada año, y compensar totalmente las interacciones entre la adaptación genética al cautiverio, la regulación poblacional, y la variación en historias de vida.

Integrating larval connectivity with local demography reveals regional dynamics of a marine metapopulation.

Many ocean species exist within what are called marine metapopulations: networks of otherwise isolated local populations connected by the exchange of larval offspring. In order to manage these species as effectively as possible (e.g., by designing and implementing effective networks of marine protected areas), we must know how many offspring are produced within each local population (i.e., local demography), and where those offspring disperse (i.e., larval connectivity). Although there is much interest in estimating connectivity in the relatively simple sense of identifying the locations of spawning parents and their settling offspring, true measures of demographic connectivity that account for among-site variation in offspring production have been lacking. We combined detailed studies of local reproductive output and larval dispersal of a coral reef fish to quantify demographic connectivity within a regional metapopulation that included four widely spaced islands in the Bahamas. We present a new method for estimating demographic connectivity when the levels of dispersal among populations are inferred by the collection of genetically "tagged" offspring. We estimated that 13.3% of recruits returned to natal islands, on average (95% CI = 1.1-50.3%), that local retention was high on one of the islands (41%, 95% CI = 6.0-97.0%), and that larval connectivity was appreciable, even between islands 129 km apart (mean = 1.6%, 95% CI = 0.20-8.8%). Our results emphasize the importance of properly integrating measurements of production with measurements of connectivity. Had we not accounted for among-site variation in offspring production, our estimates of connectivity would have been inaccurate by a factor as much as 6.5. At a generational timescale, lifetime offspring production varied substantially (a fivefold difference among islands) and the importance of each island to long-term metapopulation growth was dictated by both larval production and connectivity. At the scale of our study (local populations inhabiting 5-ha reefs), the regional metapopulation could not grow without external input. However, an exploratory analysis simulating a network of four marine protected areas suggested that reserves of >65 ha each would ensure persistence of this network. Thus, integrating studies of larval connectivity and local demography hold promise for both managing and conserving marine metapopulations effectively.

Evolutionary rescue in a host-pathogen system results in coexistence not clearance.

The evolutionary rescue of host populations may prevent extinction from novel pathogens. However, the conditions that facilitate rapid evolution of hosts, in particular the population variation in host susceptibility, and the effects of host evolution in response to pathogens on population outcomes remain largely unknown. We constructed an individual-based model to determine the relationships between genetic variation in host susceptibility and population persistence in an amphibian-fungal pathogen (Batrachochytrium dendrobatidis) system. We found that host populations can rapidly evolve reduced susceptibility to a novel pathogen and that this rapid evolution led to a 71-fold increase in the likelihood of host-pathogen coexistence. However, the increased rates of coexistence came at a cost to host populations; fewer populations cleared infection, population sizes were depressed, and neutral genetic diversity was lost. Larger adult host population sizes and greater adaptive genetic variation prior to the onset of pathogen introduction led to substantially reduced rates of extinction, suggesting that populations with these characteristics should be prioritized for conservation when species are threatened by novel infectious diseases.

Rapid genetic adaptation to a novel environment despite a genome-wide reduction in genetic diversity.

Introduced species often colonize regions that have vastly different ecological and environmental conditions than those found in their native range. As such, species introductions can provide a deeper understanding into the process of adaptive evolution. In the 1880s, steelhead trout (Oncorhynchus mykiss) from California were introduced into Lake Michigan (Laurentian Great Lakes, North America) where they established naturally reproducing populations. In their native range, steelhead hatch in rivers, migrate to the ocean and return to freshwater to spawn. Steelhead in Lake Michigan continue to swim up rivers to spawn, but now treat the freshwater environment of the Great Lakes as a surrogate ocean. To examine the effects of this introduction, we sequenced the genomes of 264 fish. By comparing steelhead from Lake Michigan to steelhead from their ancestral range, we determined that the introduction led to consistent reductions in genetic diversity across all 29 chromosomes. Despite this reduction in genetic diversity, three chromosomal regions were associated with rapid genetic adaptation to the novel environment. The first region contained functional changes to ceramide kinase, which likely altered metabolic and wound-healing rates in Lake Michigan steelhead. The second and third regions encoded carbonic anhydrases and a solute carrier protein, both of which are critical for osmoregulation, and demonstrate how steelhead physiologically adapted to freshwater. Furthermore, the contemporary release of diverse hatchery strains into the lake increased genetic diversity but reduced the signature of genetic adaptation. This study illustrates that species can rapidly adapt to novel environments despite genome-wide reductions in genetic diversity.