Anti Vasemägi

Challenges for identifying functionally important genetic variation: the promise of combining complementary research strategies

Strategies for the identification of functional genetic variation underlying phenotypic traits of ecological and evolutionary importance have received considerable attention in the literature recently. This paper aims to bring together and compare the relative strengths and limitations of various potentially useful research strategies for dissecting functionally important genetic variation in a wide range of organisms. We briefly explore the relative strengths and limitations of traditional and emerging approaches and evaluate their potential use in free‐living populations. While it is likely that much of the progress in functional genetic analyses will rely on progress in traditional model species, it is clear that with prudent choices of methods and appropriate sampling designs, much headway can be also made in a diverse range of species. We suggest that combining research approaches targeting different functional and biological levels can potentially increase understanding the genetic basis of ecological and evolutionary processes both in model and non‐model organisms.

Matrilinear phylogeography of Atlantic salmon ( Salmo salar L.) in Europe and postglacial colonization of the Baltic Sea area

Sixty‐four samples from 46 salmon populations totalling 2369 specimens were used for polymerase chain reaction–restriction fragment length polymorphism (PCR–RFLP) analysis of the mitochondrial ND1 region. The final analyses included 3095 specimens from 60 populations in Northern Europe. A subsample was analysed by RFLP of ND3/4/5/6. Representative RFLP haplotypes from different parts of the distribution area were sequenced and the phylogeny of European haplotypes and their relations to the North American lineage was described. The four common European haplotypes derive from the ancestral ND1‐BBBA (rooting the European clade to the North American) by one‐step substitutions: AAAA < AABA < BBBA > BBBB. The Swedish west‐coast populations differ from the geographically close southern Baltic, indicating absence of inward and limited outward gene flow through the Danish straits during the last 8000 years. Within the Baltic Sea, only three ND1 haplotypes were detected and there was no variation for ND3/4/5/6. In the whole southern Baltic and in lakes Vänern, Ladoga and Onega the haplotype AABA dominated. Proposed postglacial colonization routes to the Baltic Sea are discussed in relation to the haplotype distribution pattern.

SNP-array reveals genome-wide patterns of geographical and potential adaptive divergence across the natural range of Atlantic salmon (Salmo salar)

Atlantic salmon (Salmo salar) is one of the most extensively studied fish species in the world due to its significance in aquaculture, fisheries and ongoing conservation efforts to protect declining populations. Yet, limited genomic resources have hampered our understanding of genetic architecture in the species and the genetic basis of adaptation to the wide range of natural and artificial environments it occupies. In this study, we describe the development of a medium‐density Atlantic salmon single nucleotide polymorphism (SNP) array based on expressed sequence tags (ESTs) and genomic sequencing. The array was used in the most extensive assessment of population genetic structure performed to date in this species. A total of 6176 informative SNPs were successfully genotyped in 38 anadromous and freshwater wild populations distributed across the species natural range. Principal component analysis clearly differentiated European and North American populations, and within Europe, three major regional genetic groups were identified for the first time in a single analysis. We assessed the potential for the array to disentangle neutral and putative adaptive divergence of SNP allele frequencies across populations and among regional groups. In Europe, secondary contact zones were identified between major clusters where endogenous and exogenous barriers could be associated, rendering the interpretation of environmental influence on potentially adaptive divergence equivocal. A small number of markers highly divergent in allele frequencies (outliers) were observed between (multiple) freshwater and anadromous populations, between northern and southern latitudes, and when comparing Baltic populations to all others. We also discuss the potential future applications of the SNP array for conservation, management and aquaculture.

Extensive immigration from compensatory hatchery releases into wild Atlantic salmon population in the Baltic sea: spatio-temporal analysis over 18 years

Genetic homogenization has been recognized as a serious threat in an increasing number of species, including many salmonid fishes. We assessed the rate and impact of immigration from the main hatchery stocks of Atlantic salmon in the Gulf of Bothnia into one of the largest wild salmon populations in the Baltic Sea, the River Vindelälven, within a temporal framework of 18 years (from 1985–2003). We provide genetic evidence based on mtDNA and microsatellite markers, using mixed-stock analysis, that a large proportion (66%) of fin-damaged spawners (n=181) caught in the Ume/Vindelälven during 1997–2003 originated from the hatcheries in the Rivers Ångermanälven, Luleälven and Ljusnan. The maximum-likelihood estimate of immigration rate from these hatcheries into the wild Vindelälven population was 0.068 (95% CI 0.021–0.128) over the studied time period (1985–2003) and reached up to a quarter (m=0.249, 95% CI 0.106–0.419) of the total population during 1993–2000. This resulted in significant (P<0.01) genetic homogenization trend between the wild Vindelälven population and hatchery stocks of the Ångermanälven and Luleälven. Our results demonstrate extensive straying from geographically distant hatchery releases into wild salmon population and emphasize the genetic risks associated with current large-scale stocking practices in the Baltic Sea.

Permanent Genetic Resources added to the Molecular Ecology Resources Database 1 February 2010-31 March 2010.

This article documents the addition of 142 microsatellite marker loci to the Molecular Ecology Resources database. Loci were developed for the following species: Agriophyllum squarrosum, Amazilia cyanocephala, Batillaria attramentaria, Fungal strain CTeY1 (Ascomycota), Gadopsis marmoratus, Juniperus phoenicea subsp. turbinata, Liriomyza sativae, Lupinus polyphyllus, Metschnikowia reukaufii, Puccinia striiformis and Xylocopa grisescens. These loci were cross‐tested on the following species: Amazilia beryllina, Amazilia candida, Amazilia rutila, Amazilia tzacatl, Amazilia violiceps, Amazilia yucatanensis, Campylopterus curvipennis, Cynanthus sordidus, Hylocharis leucotis, Juniperus brevifolia, Juniperus cedrus, Juniperus osteosperma, Juniperus oxycedrus, Juniperus thurifera, Liriomyza bryoniae, Liriomyza chinensis, Liriomyza huidobrensis and Liriomyza trifolii.

Molecular evolutionary and population genomic analysis of the nine-spined stickleback using a modified restriction-site-associated DNA tag approach

In recent years, the explosion of affordable next generation sequencing technology has provided an unprecedented opportunity to conduct genome‐wide studies of adaptive evolution in organisms previously lacking extensive genomic resources. Here, we characterize genome‐wide patterns of variability and differentiation using pooled DNA from eight populations of the nine‐spined stickleback (Pungitius pungitius L.) from marine, lake and pond environments. We developed a novel genome complexity reduction protocol, defined as paired‐end double restriction‐site‐associated DNA (PE dRAD), to maximize read coverage at sequenced locations. This allowed us to identify over 114 000 short consensus sequences and 15 000 SNPs throughout the genome. A total of 6834 SNPs mapped to a single position on the related three‐spined stickleback genome, allowing the detection of genomic regions affected by divergent and balancing selection, both between species and between freshwater and marine populations of the nine‐spined stickleback. Gene ontology analysis revealed 15 genomic regions with elevated diversity, enriched for genes involved in functions including immunity, chemical stimulus response, lipid metabolism and signalling pathways. Comparisons of marine and freshwater populations identified nine regions with elevated differentiation related to kidney development, immunity and MAP kinase pathways. In addition, our analysis revealed that a large proportion of the identified SNPs mapping to LG XII is likely to represent alternative alleles from divergent X and Y chromosomes, rather than true autosomal markers following Mendelian segregation. Our work demonstrates how population‐wide sequencing and combining inter‐ and intra‐specific RAD analysis can uncover genome‐wide patterns of differentiation and adaptations in a non‐model species.

Beyond MHC: signals of elevated selection pressure on Atlantic salmon (Salmo salar) immune-relevant loci.

Using Atlantic salmon (Salmo salar) as a model system, we investigated whether 18 microsatellites tightly linked to immune-relevant genes have experienced different selection pressures than 76 loci with no obvious association with immune function. Immune-relevant loci were identified as outliers by two outlier tests significantly more often than nonimmune linked loci (22% vs. 1.6%). In addition, the allele frequencies of immune relevant markers were more often correlated with latitude and temperature. Combined, these results support the hypothesis that immune-relevant loci more frequently exhibit footprints of selection than other loci. They also indicate that the correlation between immune-relevant loci and latitude may be due to temperature-induced differences in pathogen-driven selection or some other environmental factor correlated with latitude.

A proteomics approach reveals divergent molecular responses to salinity in populations of European whitefish (Coregonus lavaretus)

Osmoregulation is a vital physiological function for fish, as it helps maintain a stable intracellular concentration of ions in environments of variable salinities. We focused on a primarily freshwater species, the European whitefish (Coregonus lavaretus), to investigate the molecular mechanisms underlying salinity tolerance and examine whether these mechanisms differ between genetically similar populations that spawn in freshwater vs. brackishwater environments. A common garden experiment involving 27 families in two populations and five salinity treatments together with a large‐scale, high‐resolution mass spectrometry experiment that quantified 1500 proteins was conducted to assess phenotypic and proteomic responses during early development, from fertilization until hatching, in the studied populations. The populations displayed drastically different phenotypic and proteomic responses to salinity. Freshwater‐spawning whitefish showed a significantly higher mortality rate in higher salinity treatments. Calcium, an ion involved in osmotic stress sensing, had a central role in the observed proteomic responses. Brackishwater‐spawning fish were capable of viable osmoregulation, which was modulated by cortisol, an important seawater‐adaptation hormone in teleost fish. Several proteins were identified to play key roles in osmoregulation, most importantly a highly conserved cytokine, tumour necrosis factor, whereas calcium receptor activities were associated with salinity adaptation. These results imply that individuals from these populations are most likely adapted to their local environments, even though the baseline level of genetic divergence between them is low (FSTu2003=u20030.049). They also provide clues for choosing candidate loci for studying the molecular basis of salinity adaptation in other species. Further, our approach provides an example of how proteomic methods can be successfully used to obtain novel insights into the molecular mechanisms behind adaptation in non‐model organism.

Historical and recent genetic bottlenecks in European grayling, Thymallus thymallus

Sharp declines in population size, known as genetic bottlenecks, increase the level of inbreeding and reduce genetic diversity threatening population sustainability in both short- and long-term. We evaluated the presence, severity and approximate time of bottlenecks in 34 European grayling (Thymallus thymallus) populations covering the majority of the species distribution using microsatellite markers. We identified footprints of population decline in all grayling populations using the M ratio test. In contrast to earlier simulation studies assuming isolated populations, forward simulations allowing low levels of migration demonstrated that bottleneck footprints measured using the M ratio can persist within small populations much longer (up to thousands of generations) than previously anticipated. Using a coalescence approach, the beginning of population reduction was dated back to 1,000–10,000xa0years ago which suggests that the extremely low M ratio in European grayling is most likely caused by the last glaciation and subsequent post-glacial recolonization processes. In contrast to the M ratio, two alternative methods for bottleneck detection identified more recent bottlenecks in six populations and thus, from a conservation perspective, these populations warrant future monitoring. Based on a single time-point analysis using approximate Bayesian computation methodology, all grayling populations exhibited very small effective population sizes with the majority of Ne estimates below 50. Taken together, our results demonstrate the predominate role of genetic drift in European grayling populations in the short term but also emphasize the importance of gene flow counteracting the effects of genetic drift and loss of variation over longer evolutionary timescales.

Discovery and application of insertion-deletion (INDEL) polymorphisms for QTL mapping of early life-history traits in Atlantic salmon

BackgroundFor decades, linkage mapping has been one of the most powerful and widely used approaches for elucidating the genetic architecture of phenotypic traits of medical, agricultural and evolutionary importance. However, successful mapping of Mendelian and quantitative phenotypic traits depends critically on the availability of fast and preferably high-throughput genotyping platforms. Several array-based single nucleotide polymorphism (SNP) genotyping platforms have been developed for genetic model organisms during recent years but most of these methods become prohibitively expensive for screening large numbers of individuals. Therefore, inexpensive, simple and flexible genotyping solutions that enable rapid screening of intermediate numbers of loci (~75-300) in hundreds to thousands of individuals are still needed for QTL mapping applications in a broad range of organisms.ResultsHere we describe the discovery of and application of insertion-deletion (INDEL) polymorphisms for cost-efficient medium throughput genotyping that enables analysis of >75 loci in a single automated sequencer electrophoresis column with standard laboratory equipment. Genotyping of INDELs requires low start-up costs, includes few standard sample handling steps and is applicable to a broad range of species for which expressed sequence tag (EST) collections are available. As a proof of principle, we generated a partial INDEL linkage map in Atlantic salmon (Salmo salar) and rapidly identified a number of quantitative trait loci (QTLs) affecting early life-history traits that are expected to have important fitness consequences in the natural environment.ConclusionsThe INDEL genotyping enabled fast coarse-mapping of chromosomal regions containing QTL, thus providing an efficient means for characterization of genetic architecture in multiple crosses and large pedigrees. This enables not only the discovery of larger number of QTLs with relatively smaller phenotypic effect but also provides a cost-effective means for evaluation of the frequency of segregating QTLs in outbred populations which is important for further understanding how genetic variation underlying phenotypic traits is maintained in the wild.