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, F_ST 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 F_ST = 0.038) and contemporary data sets (global F_ST = 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 F_ST 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.     

Are QST–FST comparisons for natural populations meaningful?

Comparisons between putatively neutral genetic differentiation amongst populations, F_ST, and quantitative genetic variation, Q_ST, are increasingly being used to test for natural selection. However, we find that approximately half of the comparisons that use only data from wild populations confound phenotypic and genetic variation. We urge the use of a clear distinction between narrow‐sense Q_ST, which can be meaningfully compared with F_ST, and phenotypic divergence measured between populations, P_ST, which is inadequate for comparisons in the wild. We also point out that an unbiased estimate of Q_ST can be found using the so‐called ‘animal model’ of quantitative genetics.     

Inference of Potential Genetic Risks Associated with Large-Scale Releases of Red Sea Bream in Kanagawa Prefecture, Japan Based on Nuclear and Mitochondrial DNA Analysis

Since 1978, millions of hatchery-reared red sea bream (Pagrus major) juveniles have been released in Sagami Bay and Tokyo Bay in Kanagawa Prefecture, Japan. The stock enhancement program has contributed to total catch; however, no information regarding the genetic interactions with wild counterparts is available. Here, we combined 15 microsatellite loci and mitochondrial D-loop sequencing to characterize the genetic resources of red sea bream in Sagami Bay and Tokyo Bay and to elucidate the potential harmful genetic effects associated with fish releases. Both types of markers evidenced higher levels of genetic diversity in wild samples (SB and TB) compared with offspring before stocking (H07 and H08) as well as a hatchery-released sample recaptured in Sagami Bay (HR). Microsatellite F _ST estimates and Bayesian clustering analysis found significant genetic differences among samples (F _ST?=?0.013–0.054), except for the two wild samples (F _ST?=?0.002) and HR vs. H07 (F _ST?=?0.007). On the other hand, mitochondrial-based Ф _ST suggested haplotypic similarity between SB, H07, and HR. The low effective number of females contributing to the offspring over multiple generations may be responsible for the lack of haplotypic differentiation. Moreover, the putative hatchery origin to three fish (8 %) without deformity in the inter-nostril epidermis was inferred for the first time. Our results showed the usefulness of combining nuclear and mitochondrial markers to elucidate genetic interactions between hatchery-released and wild red sea bream and warned about potential harmful genetic effects should interbreeding takes place.     

Detection of farm fox and hybrid genotypes among wild arctic foxes in Scandinavia

In Scandinavia, farmed arctic foxes frequently escape from farms, raising concern about hybridization with the endangered wild population. This study was performed to find a genetic marker to distinguish escaped farm foxes from wild Scandinavian foxes. Microsatellite and mitochondrial control region variation were analyzed in 41 farm foxes. The results were compared with mitochondrial and microsatellite data from the wild population in Scandinavia. The farm foxes were genetically distinct from the wild foxes (F _ST=0.254, P < 0.00001) and all farm foxes had a single control region haplotype different from those observed in the wild population. We developed a method based on Restriction Fragment Length Polymorphism (RFLP) on the mitochondrial control region to differentiate between farmed and wild arctic foxes. This test was subsequently successfully used on 25 samples from free-ranging foxes, of which four had a suspected farm origin. All four of the suspected foxes, and none of the others, carried the farm fox haplotype. Three of these were successfully genotyped for all eleven microsatellite loci. A population assignment test and a Bayesian Markov Chain Monte Carlo analysis indicated that two of these individuals were escaped farm foxes, and that the third possibly was a hybrid between a farmed and a wild arctic fox.     

Contrasting patterns of morphological and neutral genetic divergence among geographically proximate populations of sockeye salmon Oncorhynchus nerka in Lake Aleknagik,Alaska

Levels of differentiation in morphological traits (age at maturity, body length at age, egg mass and body depth) and spawning time were examined in sockeye salmon Oncorhynchus nerka from three geographically proximate but physically distinct creeks in Lake Aleknagik, Alaska. Happy Creek fish had significantly greater values for most measured morphological traits, and Eagle Creek fish spawned significantly later than fish in the other creeks. Phenotypic differentiation between creeks, measured as P_ST, was then compared with microsatellite marker differentiation between creeks, measured as F_ST. No correlations were apparent between P_ST and F_ST values, and P_ST values were generally significantly larger than zero (P_ST= 0·0018–0·31) whereas F_ST values were not (F_ST=?0·0004 to 0·0016). The insignificant pair‐wise F_ST values between creek samples indicated that gene flow occurs between creeks, assuming the creek populations have reached migration–drift equilibrium. However, the strong homing behaviour of sockeye salmon precludes a scenario in which fish from the three creeks constitute a single population that segregates by body size. Rather, significant phenotypic differentiation suggests that strong divergent selection occurs on the phenotypic traits despite the homogenizing effects of gene flow.     

Microgeographic structure of Anopheles gambiae in western Kenya based on mtDNA and microsatellite loci

The population genetic structure of the Anopheles gambiae in western Kenya was studied using length variation at five microsatellite loci and sequence variation in a 648-nt mtDNA fragment. Mosquitoes were collected from houses in villages spanning up to 50 km distance, The following questions were answered, (i) Are mosquitoes in a house more related genetically to each other than mosquitoes between houses? (ii) What degree of genetic differentiation occurs on these geographical scales? (iii) How consistent are the results obtained with both types of genetic markers? At the house level, no differentiation was detected by F_ST and R_ST, and the band sharing index test revealed no significant associations of alleles across loci. Likewise, indices of kinship based on mtDNA haplotypes in houses were even lower than in the pooled sample. Therefore, the hypothesis that mosquitoes in a house are more related genetically was rejected. At increasing geographical scales, microsatellite allele distributions were similar among all population samples and no subdivision of the gene pool was detected by F_ST or R_ST. Likewise, estimates of haplotype divergence of mtDNA between populations were not higher than the within population estimates, and mtDNA-based F_ST values were not significantly different from zero. That sequence variation in mtDNA provided matching results with microsatellite loci (while high genetic variation was observed in all loci), suggested that this pattern represents the whole genome. The minimum area associated with a deme of A. gambiae in western Kenya is therefore larger than 50 km in diameter.     

Population genomics of wild and laboratory zebrafish (Danio rerio)

Understanding a wider range of genotype–phenotype associations can be achieved through ecological and evolutionary studies of traditional laboratory models. Here, we conducted the first large‐scale geographic analysis of genetic variation within and among wild zebrafish (Danio rerio) populations occurring in Nepal, India, and Bangladesh, and we genetically compared wild populations to several commonly used lab strains. We examined genetic variation at 1832 polymorphic EST‐based single nucleotide polymorphisms (SNPs) and the cytb mitochondrial gene in 13 wild populations and three lab strains. Natural populations were subdivided into three major mitochondrial DNA clades with an average among‐clade sequence divergence of 5.8%. SNPs revealed five major evolutionarily and genetically distinct groups with an overall F_ST of 0.170 (95% CI 0.105–0.254). These genetic groups corresponded to discrete geographic regions and appear to reflect isolation in refugia during past climate cycles. We detected 71 significantly divergent outlier loci (3.4%) and nine loci (0.5%) with significantly low F_ST values. Valleys of reduced heterozygosity, consistent with selective sweeps, surrounded six of the 71 outliers (8.5%). The lab strains formed two additional groups that were genetically distinct from all wild populations. An additional subset of outlier loci was consistent with domestication selection within lab strains. Substantial genetic variation that exists in zebrafish as a whole is missing from lab strains that we analysed. A combination of laboratory and field studies that incorporates genetic variation from divergent wild populations along with the wealth of molecular information available for this model organism provides an opportunity to advance our understanding of genetic influences on phenotypic variation for a vertebrate species.     

Population Structure and Stock Identification of Eulachon (Thaleichthys pacificus), an Anadromous Smelt, in the Pacific Northwest

The genetic structure of eulachon (Thaleichthys pacificus) populations was examined in an analysis of variation of 14 microsatellite loci representing approximately 1900 fish from 9 sites between the Columbia River and Cook Inlet, Alaska. Significant genetic differentiation occurred among the putative populations. The mean F_ST for all loci was 0.0046, and there was a significant correlation between population genetic differentiation (F_ST) and geographic distance. Simulated mixed-stock samples comprising populations from different regions suggested that variation at microsatellite loci provided reasonably accurate estimates of stock composition for potential fishery samples. Marine sampling indicated that immature eulachons from different rivers, during the 2 to 3 years of prespawning life in offshore marine waters, do not mix thoroughly. For eulachons captured incidentally in offshore trawl fisheries, there was a clear geographic cline in relative abundance of eulachons from different geographic areas. The sample from northern British Columbia was dominated by northern and central coastal populations of British Columbia, the sample from central British Columbia was composed of eulachons from all regions, and the sample from southern British Columbia was dominated by Columbia River and Fraser River populations. These results have implications for the management of trawl fisheries and conservation of spawning populations in some rivers where abundance is at historically low levels.     

Genetic rescue of an inbred captive population of the critically endangered Puerto Rican crested toad (<Emphasis Type="Italic">Peltophryne lemur</Emphasis>) by mixing lineages

The Puerto Rican crested toad (Peltophryne lemur) is currently composed of a single wild population on the south coast of Puerto Rico and two captive populations founded by animals from the northern and southern coasts. The main factors contributing to its decline are habitat loss, inundation of breeding ponds during storms, and impacts of invasive species. Recovery efforts have been extensive, involving captive breeding and reintroductions, habitat restoration, construction of breeding ponds, and public education. To guide future conservation efforts, genetic variation and differentiation were assessed for the two captive colonies and the remaining wild population using the mitochondrial control region and six novel microsatellite loci. Only two moderately divergent mitochondrial haplotypes were found, with one fixed in each of the southern and northern lineages. Moderate genetic variation exists for microsatellite loci in all three groups. The captive southern population has not diverged substantially from the wild population at microsatellite loci (F _ST = 0.03), whereas there is little allelic overlap between the northern and southern lineages at five of six loci (F _ST > 0.3). Despite this differentiation, they are no more divergent than many populations of other amphibian species. As the northern breeding colony may not remain viable due to its small size and inbred nature, it is recommended that a third breeding colony be established in which northern and southern individuals are combined. This will preserve any northern adaptive traits that may exist, and provide animals for release in the event that the pure northern lineage becomes extirpated.     

Preliminary data on the variability of three microsatellite loci in Pacific <Emphasis Type="Italic">Gadus macrocephalus</Emphasis> and Atlantic <Emphasis Type="Italic">G. morhua</Emphasis> cod (Gadidae)

Variability of microsatellite DNA loci Gmo3, Gmo34, and Gmo35 is studied in samples of Pacific cod Gadus macrocephalus and Atlantic cod G. morhua. The results show high values of identity of the samples within the North Pacific basin (0.9766–0.9924) and within the Northeast Atlantic basin (0.9580). Based on the pairwise assessment of genetic differentiation, the F _ST values are significantly different in all variants between the samples of Pacific and Atlantic cod (F _ST = 0.5235–0.6719, p < 0.001). Within the basins, the significant differences in the frequencies of main alleles are revealed in the loci Gmo3 and Gmo34 for the samples from the Pacific and Atlantic oceans, respectively.     

Genetic diversity and population differentiation in the yellow drum Nibea albiflora along the coast of the China Sea

Nibea albiflora (yellow drum) is an important seafood fish species in East Asia. We explored the population genetic variation of N. albiflora along the coastal waters of the China Sea using microsatellite markers to facilitate a selective breeding programme that is undertaken in China. A total of 256 alleles were detected at 12 loci in four wild populations. A high level of genetic diversity was observed with the mean number of alleles and the observed and expected heterozygosity in each population ranging from 7.917 to 14.083, 0.701 to 0.764 and 0.765 to 0.841, respectively. Pairwise fixation index (F_ST) analysis indicated significant but weak genetic differentiation among populations from four localities (F_ST?=?0.030, P?<?0.01), which was also confirmed by analysis of molecular variance (AMOVA). Significant genetic differentiation was detected between Ningde and the other populations (F_ST?=?0.047–0.056, P?<?0.01). Structure analysis suggested that N. albiflora within the examined range might be composed of two stocks. The data of the present study revealed high genetic diversity and low genetic differentiation among the N. albiflora populations along the coast of the China Sea. This baseline information could be valuable for future selective breeding programmes of N. albiflora.     

Population structure and pattern of geographic variation in Muscari comosum along its range of distribution

Genetic variation for six loci in 37 populations of Muscari comosum L. (Liliaeeae) is surveyed. One locus is monomorphic and identical in all the populations. The remaining loci are polymorphic. Although the GOT-1 and GOT-3 loci show a pronounced heterozygote deficit explained by selection acting upon these loci (or on genes linked to them), the remaining loci nearly conform to Hardy-Weinberg proportions. The overall pattern shows a low level of heterozygote deficit (F_IS=0.08) explained by the mixed mating system. The organization of genetic variation shows a low level of interpopulation differentiation (F_ST or G_ST=0.04). At the same time, autocorrelation analysis shows no pattern of geographical variation. It is concluded that gene flow and selection interact to produce the overall pattern of genetic variation.     

Comparison of Four Mendelian Loci of the California Sea Hare (Aplysia californica) from Populations of the Coast of California and the Sea of Cortez

Aplysia californica is a species widely used in neurobiology, and specimens are collected from a wide range of places along its distribution range. A. californica is endemic to the coast of California and the Gulf of California. On the west coast, this is an unusual distribution range relative to other benthic species from that region. Four polymorphic nuclear Mendelian markers were identified (three single-copy nuclear DNA loci and one microsatellite) for an initial survey of genetic variation of wild populations. F _ST values not significantly different from 0 (overall F _ST= 0.0148) suggest there was no geographic genetic population subdivision in 177 individuals examined. Received December 3, 1999; accepted March 3, 2000.     

Targeted resequencing reveals geographical patterns of differentiation for loci implicated in parallel evolution

Parallel divergence and speciation provide evidence for the role of divergent selection in generating biological diversity. Recent studies indicate that parallel phenotypic divergence may not have the same genetic basis in different geographical locations – ‘outlier loci’ (loci potentially affected by divergent selection) are often not shared among parallel instances of phenotypic divergence. However, limited sharing may be due, in part, to technical issues if false‐positive outliers occur. Here, we test this idea in the marine snail Littorina saxatilis, which has evolved two partly isolated ecotypes (adapted to crab predation vs. wave action) in multiple locations independently. We argue that if the low extent of sharing observed in earlier studies in this system is due to sampling effects, we expect outliers not to show elevated F_ST when sequenced in new samples from the original locations and also not to follow predictable geographical patterns of elevated F_ST. Following a hierarchical sampling design (within vs. between country), we applied capture sequencing, targeting outliers from earlier studies and control loci. We found that outliers again showed elevated levels of F_ST in their original location, suggesting they were not generated by sampling effects. Outliers were also likely to show increased F_ST in geographically close locations, which may be explained by higher levels of gene flow or shared ancestral genetic variation compared with more distant locations. However, in contrast to earlier findings, we also found some outlier types to show elevated F_ST in geographically distant locations. We discuss possible explanations for this unexpected result.     

Background selection and FST: Consequences for detecting local adaptation

Background selection is a process whereby recurrent deleterious mutations cause a decrease in the effective population size and genetic diversity at linked loci. Several authors have suggested that variation in the intensity of background selection could cause variation in F_ST across the genome, which could confound signals of local adaptation in genome scans. We performed realistic simulations of DNA sequences, using recombination maps from humans and sticklebacks, to investigate how variation in the intensity of background selection affects F_ST and other statistics of population differentiation in sexual, outcrossing species. We show that, in populations connected by gene flow, Weir and Cockerham's (1984; Evolution, 38 , 1358) estimator of F_ST is largely insensitive to locus‐to‐locus variation in the intensity of background selection. Unlike F_ST, however, d_XY is negatively correlated with background selection. Moreover, background selection does not greatly affect the false‐positive rate in F_ST outlier studies in populations connected by gene flow. Overall, our study indicates that background selection will not greatly interfere with finding the variants responsible for local adaptation.     

Replicate divergence between and within sounds in a marine fish: the copper rockfish (Sebastes caurinus)

Understanding the factors that influence larval dispersal and connectivity among marine populations is critical to the conservation and sustainable management of marine resources. We assessed genetic subdivision among ten populations of copper rockfish (Sebastes caurinus) representing paired samples of outer coast and the heads of inlets in five replicate sounds on the west coast of Vancouver Island, British Columbia, using 17 microsatellite DNA loci. Overall, subdivision (F_ST) was low (F_ST = 0.031, P < 0.001), but consistently higher between paired coast and head of inlet sites (mean F_ST = 0.047, P < 0.001) compared to among the five coast sites (mean F_ST = ?0.001, P > 0.5) or among the five head of inlet sites (mean F_ST = 0.026, P < 0.001). Heterozygosity, allelic richness and estimates of effective population size were also lower in head of inlet sites than in coast sites. Bayesian analysis identified two genetic groups across all samples, a single genetic group among only coast samples, two genetic groups among head of inlet samples and two genetic groups within each sound analysed separately. Head of inlet copper rockfish tended to be shorter with lower condition factors and grew more slowly than coast sites fish. Reduced physical connectivity and selection against immigrants in contrasting outer coast–head of inlet environments likely contribute to the evolution of population structure of copper rockfish. Based on genetic connectivity, coast sites appear to be better served by existing marine protected areas than are head of inlet sites.     

A Beta-mixture model for assessing genetic population structure

Summary An important fraction of recently generated molecular data is dominant markers. They contain substantial information about genetic variation but dominance makes it impossible to apply standard techniques to calculate measures of genetic differentiation, such as F‐statistics. In this article, we propose a new Bayesian beta‐mixture model that more accurately describes the genetic structure from dominant markers and estimates multiple F_ST s from the sample. The model also has important application for codominant markers and single‐nucleotide polymorphism (SNP) data. The number of F_ST is assumed unknown beforehand and follows a random distribution. The reversible jump algorithm is used to estimate the unknown number of multiple F_ST s. We evaluate the performance of three split proposals and the overall performance of the proposed model based on simulated dominant marker data. The model could reliably identify and estimate a spectrum of degrees of genetic differentiation present in multiple loci. The estimates of F_ST s also incorporate uncertainty about the magnitude of within‐population inbreeding coefficient. We illustrate the method with two examples, one using dominant marker data from a rare orchid and the other using codominant marker data from human populations.     

Identifying footprints of selection in stocked brown trout populations: a spatio‐temporal approach

Studies of interactions between farmed and wild salmonid fishes have suggested reduced fitness of farmed strains in the wild, but evidence for selection at the genic level is lacking. We studied three brown trout populations in Denmark which have been significantly admixed with stocked hatchery trout (19–64%), along with two hatchery strains used for stocking. The wild populations were represented by contemporary samples (2000–2006) and two of them by historical samples (1943–1956). We analysed 61 microsatellite loci, nine of which showed putative functional relationships [expressed sequence tag (EST)‐linked or quantitative trait loci]. F_ST‐based outlier tests provided support for diversifying selection at chromosome regions marked by three loci, two anonymous and one EST‐linked. Patterns of differentiation suggested that the loci were candidates for being under diversifying hitch‐hiking selection in hatchery vs. wild environments. Analysis of hatchery strain admixture proportions showed that in one wild population, two of the loci showed significantly lower admixture proportions than the putatively neutral loci, implying contemporary selection against alleles introduced by hatchery strain trout. In the most strongly admixed population, however, there was no evidence for selection, possibly because of immigration by stocked trout overcoming selection against hatchery‐derived alleles or supportive breeding practices allowing hatchery strain trout to escape natural selection. To our knowledge, this is the first study demonstrating footprints of selection in wild salmonid populations subject to spawning intrusion by farmed fish.     

Identifying outlier loci in admixed and in continuous populations using ancestral population differentiation statistics

Finding genetic signatures of local adaptation is of great interest for many population genetic studies. Common approaches to sorting selective loci from their genomic background focus on the extreme values of the fixation index, F_ST, across loci. However, the computation of the fixation index becomes challenging when the population is genetically continuous, when predefining subpopulations is a difficult task, and in the presence of admixed individuals in the sample. In this study, we present a new method to identify loci under selection based on an extension of the F_ST statistic to samples with admixed individuals. In our approach, F_ST values are computed from the ancestry coefficients obtained with ancestry estimation programs. More specifically, we used factor models to estimate F_ST, and we compared our neutrality tests with those derived from a principal component analysis approach. The performances of the tests were illustrated using simulated data and by re‐analysing genomic data from European lines of the plant species Arabidopsis thaliana and human genomic data from the population reference sample, POPRES.     

Evaluation of Genetic Variation in the Clown Knifefish, Chitala chitala, Using Allozymes, RAPD, and Microsatellites

Twenty-seven enzyme systems, six random amplified polymorphic DNA (RAPD) primers, and two microsatellite loci were tested to determine intraspecific divergence in the natural population of the endangered Indian featherback fish, Chitala chitala, for the first time. The 262 samples of C. chitala were collected from six riverine locations in India: the Satluj, Ganga (Ghagra, Bhagirathi, and Brahmaputra), Mahanadi, and Narmada river systems. The analysis revealed population subdivisions, with an F_ST value from 0.1235 (95% confidence 0.0868–0.1621) for RAPD and a combined F_ST of 0.0344 (95% confidence 0.0340–0.0350) for microsatellite loci. An analysis of 38 allozyme loci did not reveal any polymorphism in the samples from any of the riverine localities; a possible explanation for this could be that the ancestors of Chitala could have faced a population reduction in prehistoric periods, as low allozyme variation is also reported for other species of Chitala from south Asia.