María Quintela

Three Decades of Farmed Escapees in the Wild: A Spatio-Temporal Analysis of Atlantic Salmon Population Genetic Structure throughout Norway

Each year, hundreds of thousands of domesticated farmed Atlantic salmon escape into the wild. In Norway, which is the world’s largest commercial producer, many native Atlantic salmon populations have experienced large numbers of escapees on the spawning grounds for the past 15–30 years. In order to study the potential genetic impact, we conducted a spatio-temporal analysis of 3049 fish from 21 populations throughout Norway, sampled in the period 1970–2010. Based upon the analysis of 22 microsatellites, individual admixture, FST and increased allelic richness revealed temporal genetic changes in six of the populations. These changes were highly significant in four of them. For example, 76% and 100% of the fish comprising the contemporary samples for the rivers Vosso and Opo were excluded from their respective historical samples at P = 0.001. Based upon several genetic parameters, including simulations, genetic drift was excluded as the primary cause of the observed genetic changes. In the remaining 15 populations, some of which had also been exposed to high numbers of escapees, clear genetic changes were not detected. Significant population genetic structuring was observed among the 21 populations in the historical (global FST = 0.038) and contemporary data sets (global FST = 0.030), although significantly reduced with time (P = 0.008). This reduction was especially distinct when looking at the six populations displaying temporal changes (global FST dropped from 0.058 to 0.039, P = 0.006). We draw two main conclusions: 1. The majority of the historical population genetic structure throughout Norway still appears to be retained, suggesting a low to modest overall success of farmed escapees in the wild; 2. Genetic introgression of farmed escapees in native salmon populations has been strongly population-dependent, and it appears to be linked with the density of the native population.

Can balancing selection on MHC loci counteract genetic drift in small fragmented populations of black grouse

The ability of natural populations to adapt to new environmental conditions is crucial for their survival and partly determined by the standing genetic variation in each population. Populations with higher genetic diversity are more likely to contain individuals that are better adapted to new circumstances than populations with lower genetic diversity. Here, we use both neutral and major histocompatibility complex (MHC) markers to test whether small and highly fragmented populations hold lower genetic diversity than large ones. We use black grouse as it is distributed across Europe and found in populations with varying degrees of isolation and size. We sampled 11 different populations; five continuous, three isolated, and three small and isolated. We tested patterns of genetic variation in these populations using three different types of genetic markers: nine microsatellites and 21 single nucleotide polymorphisms (SNPs) which both were found to be neutral, and two functional MHC genes that are presumably under selection. The small isolated populations displayed significantly lower neutral genetic diversity compared to continuous populations. A similar trend, but not as pronounced, was found for genotypes at MHC class II loci. Populations were less divergent at MHC genes compared to neutral markers. Measures of genetic diversity and population genetic structure were positively correlated among microsatellites and SNPs, but none of them were correlated to MHC when comparing all populations. Our results suggest that balancing selection at MHC loci does not counteract the power of genetic drift when populations get small and fragmented.

Decreased TBT pollution and changing bioaccumulation pattern in gastropods imply butyltin desorption from sediments.

Two monitoring surveys were repeated in Galicia (NW Spain): one in 2003 concerned the rocky shore gastropod Nucella lapillus (19 populations, the reference campaign in 1996), the other dealt in 2005 with the infaunal snail Nassarius reticulatus (25 sites, the previous one in 2000). Samples were subject to a standard protocol to determine the concentrations of butyltins (BTs) in tissues. Results show that pollution in most populations has considerably decreased over the last decade: for N. lapillus the mean descent ranged from 37% (for tributyltin -TBT-) to 66% (for monobutyltin -MBT-), and TBT concentrations were on average halved in N. reticulatus. However, derivatives generally increased in this latter species, to the extent that dibutyltin -DBT- in several 2005 samples exceeded the aggregate of all three BTs in 2000. As a consequence, a major change in the bioaccumulation patterns becomes evident, particularly when computing the butyltin degradation index [BDI: (DBT+MBT)/TBT]. This shift is most marked at sites where pollution has always been lowest, and it shows significant negative correlation between both gastropods. Since sources other than antifouling paints are not important in the area, it is proposed that observations are due to BT desorption from sediments through some interplay involving the different characteristics of the chemicals and the contrasting biology of the animals.

The use of Nucella lapillus (L.) transplanted in cages to monitor tributyltin (TBT) pollution

Specimens of the gastropod Nucella lapillus were transplanted from their natural habitat to sites affected by tributyltin (TBT) pollution to diverse degrees; individuals were held in cages that also included mussels as food. Subsamples were taken from the cages and the natural site at different time intervals (1, 2 and 5 months) to compare growth, imposex (a set of sexual anomalies caused by TBT) and butyltin concentrations accumulated in female tissues. Individuals in cages experienced negligible mortality and considerable growth. In addition, while the Relative Penis Size Index (RPSI, an index of imposex) remained consistent in the natural site through the experiment, it increased markedly at the implant sites; the same applies to the butyltin residues. These results are discussed to conclude that cage transplanting of N. lapillus is a suitable technique that can quickly (i.e. 1 month) detect differences in the level of TBT contamination among sites; however, extrapolation of experimental results to the field seems to require a longer exposure time (at least 5 months). Finally, the use of such an approach is deemed helpful to interpret the relationship between N. lapillus population dynamics and TBT pollution.

Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.

Fine‐grained adaptive divergence in an amphibian: genetic basis of phenotypic divergence and the role of nonrandom gene flow in restricting effective migration among wetlands

Adaptive ecological differentiation among sympatric populations is promoted by environmental heterogeneity, strong local selection and restricted gene flow. High gene flow, on the other hand, is expected to homogenize genetic variation among populations and therefore prevent local adaptation. Understanding how local adaptation can persist at the spatial scale at which gene flow occurs has remained an elusive goal, especially for wild vertebrate populations. Here, we explore the roles of natural selection and nonrandom gene flow (isolation by breeding time and habitat choice) in restricting effective migration among local populations and promoting generalized genetic barriers to neutral gene flow. We examined these processes in a network of 17 breeding ponds of the moor frog Rana arvalis, by combining environmental field data, a common garden experiment and data on variation in neutral microsatellite loci and in a thyroid hormone receptor (TRβ) gene putatively under selection. We illustrate the connection between genotype, phenotype and habitat variation and demonstrate that the strong differences in larval life history traits observed in the common garden experiment can result from adaptation to local pond characteristics. Remarkably, we found that haplotype variation in the TRβ gene contributes to variation in larval development time and growth rate, indicating that polymorphism in the TRβ gene is linked with the phenotypic variation among the environments. Genetic distance in neutral markers was correlated with differences in breeding time and environmental differences among the ponds, but not with geographical distance. These results demonstrate that while our study area did not exceed the scale of gene flow, ecological barriers constrained gene flow among contrasting habitats. Our results highlight the roles of strong selection and nonrandom gene flow created by phenological variation and, possibly, habitat preferences, which together maintain genetic and phenotypic divergence at a fine‐grained spatial scale.

Local divergence of thermal reaction norms among amphibian populations is affected by pond temperature variation

Although temperature variation is known to cause large‐scale adaptive divergence, its potential role as a selective factor over microgeographic scales is less well‐understood. Here, we investigated how variation in breeding pond temperature affects divergence in multiple physiological (thermal performance curve and critical thermal maximum [CTmax]) and life‐history (thermal developmental reaction norms) traits in a network of Rana arvalis populations. The results supported adaptive responses to face two main constraints limiting the evolution of thermal adaptation. First, we found support for the faster–slower model, indicating an adaptive response to compensate for the thermodynamic constraint of low temperatures in colder environments. Second, we found evidence for the generalist–specialist trade‐off with populations from colder and less thermally variable environments exhibiting a specialist phenotype performing at higher rates but over a narrower range of temperatures. By contrast, the local optimal temperature for locomotor performance and CTmax did not match either mean or maximum pond temperatures. These results highlight the complexity of the adaptive multiple‐trait thermal responses in natural populations, and the role of local thermal variation as a selective force driving diversity in life‐history and physiological traits in the presence of gene flow.

A multilocus assay reveals high nucleotide diversity and limited differentiation among Scandinavian willow grouse (Lagopus lagopus).

BackgroundThere is so far very little data on autosomal nucleotide diversity in birds, except for data from the domesticated chicken and some passerines species. Estimates of nucleotide diversity reported so far in birds have been high (~10-3) and a likely explanation for this is the generally higher effective population sizes compared to mammals. In this study, the level of nucleotide diversity has been examined in the willow grouse, a non-domesticated bird species from the order Galliformes, which also holds the chicken. The willow grouse (Lagopus lagopus) has an almost circumpolar distribution but is absent from Greenland and the north Atlantic islands. It primarily inhabits tundra, forest edge habitats and sub-alpine vegetation. Willow grouse are hunted throughout its range, and regionally it is a game bird of great cultural and economical importance.ResultsWe sequenced 18 autosomal protein coding loci from approximately 15–18 individuals per population. We found a total of 127 SNPs, which corresponds to 1 SNP every 51 bp. 26 SNPs were amino acid replacement substitutions. Total nucleotide diversity (πt) was between 1.30 × 10-4 and 7.66 × 10-3 (average πt= 2.72 × 10-3 ± 2.06 × 10-3) and silent nucleotide diversity varied between 4.20 × 10-4and 2.76 × 10-2 (average πS= 9.22 × 10-3 ± 7.43 × 10-4). The synonymous diversity is approximately 20 times higher than in humans and two times higher than in chicken. Non-synonymous diversity was on average 18 times lower than the synonymous diversity and varied between 0 and 4.90 × 10-3 (average πa= 5.08 × 10-4 ± 7.43 × 103), which suggest that purifying selection is strong in these genes. FST values based on synonymous SNPs varied between -5.60 × 10-4 and 0.20 among loci and revealed low levels of differentiation among the four localities, with an overall value of FST = 0.03 (95% CI: 0.006 – 0.057) over 60 unlinked loci. Non-synonymous SNPs gave similar results. Low levels of linkage disequilibrium were observed within genes, with an average r2 = 0.084 ± 0.110, which is expected for a large outbred population with no population differentiation. The mean per site per generation recombination parameter (ρ) was comparably high (0.028 ± 0.018), indicating high recombination rates in these genes.ConclusionWe found unusually high levels of nucleotide diversity in the Scandinavian willow grouse as well as very little population structure among localities with up to 1647 km distance. There are also low levels of linkage disequilibrium within the genes and the population recombination rate is high, which is indicative of an old panmictic population, where recombination has had time to break up any haplotype blocks. The non-synonymous nucleotide diversity is low compared with the silent, which is in agreement with effective purifying selection, possibly due to the large effective population size.

Genetic diversity and differentiation among Lagopus lagopus populations in Scandinavia and Scotland: evolutionary significant units confirmed by SNP markers

Single Nucleotide Polymorphism in four Scandinavian populations of willow grouse (Lagopus lagopus) and two Scottish populations of red grouse (Lagopus lagopus scoticus) were assessed at 13 protein‐coding loci. We found high levels of diversity, with one substitution every 55 bp as an average and a total of 76 unlinked parsimony informative SNPs. Different estimators of genetic diversity such as: number of synonymous and non‐synonymous sites, average number of alleles, number and percentage of polymorphic loci, mean nucleotide diversity (πs, πa) and gene diversity at synonymous and non‐synonymous sites showed higher diversity in the northern populations compared to southern ones. Strong levels of purifying selection found in all the populations together with neutrality tests conforming to neutral expectations agree with large effective population sizes. Assignment tests reported a clear distinction between Scandinavian and Scottish grouse suggesting the existence of two different evolutionary significant units. The divergence time between willow and red grouse ranging between 12 500 and 125 000 years, in conjunction with the presence of ‘specific’ markers for each subspecies prompt a reassessment of the taxonomical status of the Scottish red grouse.

AFLPs and Mitochondrial Haplotypes Reveal Local Adaptation to Extreme Thermal Environments in a Freshwater Gastropod

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