Mikhail Ozerov

Microsatellite-based genetic diversity and population structure of domestic sheep in northern Eurasia

BackgroundIdentification of global livestock diversity hotspots and their importance in diversity maintenance is essential for making global conservation efforts. We screened 52 sheep breeds from the Eurasian subcontinent with 20 microsatellite markers. By estimating and weighting differently within- and between-breed genetic variation our aims were to identify genetic diversity hotspots and prioritize the importance of each breed for conservation, respectively. In addition we estimated how important within-species diversity hotspots are in livestock conservation.ResultsBayesian clustering analysis revealed three genetic clusters, termed Nordic, Composite and Fat-tailed. Southern breeds from close to the region of sheep domestication were more variable, but less genetically differentiated compared with more northern populations. Decreasing weight for within-breed diversity component led to very high representation of genetic clusters or regions containing more diverged breeds, but did not increase phenotypic diversity among the high ranked breeds. Sampling populations throughout 14 regional groups was suggested for maximized total genetic diversity.ConclusionsDuring initial steps of establishing a livestock conservation program populations from the diversity hot-spot area are the most important ones, but for the full design our results suggested that approximately equal population presentation across environments should be considered. Even in this case, higher per population emphasis in areas of high diversity is appropriate. The analysis was based on neutral data, but we have no reason to think the general trend is limited to this type of data. However, a comprehensive valuation of populations should balance production systems, phenotypic traits and available genetic information, and include consideration of probability of success.

High Gyrodactylus salaris infection rate in triploid Atlantic salmon Salmo salar.

We describe an unusually high infection rate of Gyrodactylus salaris Malmberg in juvenile Atlantic salmon Salmo salar L. of Baltic Sea origin, which are generally believed to be more resistant to G. salaris than East Atlantic salmon populations. Based on analyses of mitochondrial (complete cytochrome oxidase 1 [CO1] gene, 1548 bp) and nuclear (ADNAM1, 435 bp; internal transcribed spacer [ITS] rDNA region, 1232 bp) DNA fragments, the closest relatives of the characterized Estonian G. salaris strain were parasites found off the Swedish west coast and in Raasakka hatchery, Iijoki (Baltic Sea, Finland). Analyses of 14 microsatellite loci of the host S. salarrevealed that approximately 40% of studied fish were triploids. We subsequently identified triploid Atlantic salmon of Baltic origin as more susceptible to G. salaris infection than their diploid counterparts, possibly due to compromised complement-dependent immune pathways in triploid salmon. This is in accordance with earlier studies that have shown elevated susceptibility of triploids to various viral or bacterial pathogens, and represents one of the first reports of increased susceptibility of triploid salmonid fish to an ectoparasite. However, further experimental work is needed to determine whether triploid Atlantic salmon is generally more susceptible to G. salaris compared to their diploid counterparts, irrespective of the particular triploidization method and population of origin.

Cost-effective genome-wide estimation of allele frequencies from pooled DNA in Atlantic salmon (Salmo salar L.)

BackgroundNew sequencing technologies have tremendously increased the number of known molecular markers (single nucleotide polymorphisms; SNPs) in a variety of species. Concurrently, improvements to genotyping technology have now made it possible to efficiently genotype large numbers of genome-wide distributed SNPs enabling genome wide association studies (GWAS). However, genotyping significant numbers of individuals with large number of SNPs remains prohibitively expensive for many research groups. A possible solution to this problem is to determine allele frequencies from pooled DNA samples, such ‘allelotyping’ has been presented as a cost-effective alternative to individual genotyping and has become popular in human GWAS. In this article we have tested the effectiveness of DNA pooling to obtain accurate allele frequency estimates for Atlantic salmon (Salmo salar L.) populations using an Illumina SNP-chip.ResultsIn total, 56 Atlantic salmon DNA pools from 14 populations were analyzed on an Atlantic salmon SNP-chip containing probes for 5568 SNP markers, 3928 of which were bi-allelic. We developed an efficient quality control filter which enables exclusion of loci showing high error rate and minor allele frequency (MAF) close to zero. After applying multiple quality control filters we obtained allele frequency estimates for 3631 bi-allelic loci. We observed high concordance (r > 0.99) between allele frequency estimates derived from individual genotyping and DNA pools. Our results also indicate that even relatively small DNA pools (35 individuals) can provide accurate allele frequency estimates for a given sample.ConclusionsDespite of higher level of variation associated with array replicates compared to pool construction, we suggest that both sources of variation should be taken into account. This study demonstrates that DNA pooling allows fast and high-throughput determination of allele frequencies in Atlantic salmon enabling cost-efficient identification of informative markers for discrimination of populations at various geographical scales, as well as identification of loci controlling ecologically and economically important traits.

Genetic structure of freshwater Atlantic salmon (Salmo salar L.) populations from the lakes Onega and Ladoga of northwest Russia and implications for conservation

Freshwater Atlantic salmon (Salmo salar L.) populations are a living example of adaptation to the changing conditions caused by glacial cycles. The uniqueness of these populations is emphasized by almost complete resistance to the dangerous parasite Gyrodactylus salaris. In Europe, freshwater salmon populations occur primarily in north-western Russia in the republic of Karelia. These systems include Lakes Ladoga and Onega, the two largest lakes in Europe, each of which harbours a number of freshwater salmon spawning rivers. We used microsatellite markers to study the genetic structure and temporal stability of 11 freshwater salmon populations in Russian Karelia. Populations clustered according their region of origin. Although temporal variation in allele frequencies was observed in the majority of temporal comparisons, various lines of evidence demonstrated that this influence was relatively minor compared to spatial variation that explained eight times more of the variability than temporal variation. Temporal stability tended to occur in populations from rivers with a higher linear lake coefficient. The high level of genetic structuring observed in both lake systems and the apparent low level of migration between populations suggests that treating each river as a separate management unit is recommended. In addition, as the number of populations is large, the best strategy for such fine scale management would be to ensure that the level of natural reproduction in each river is sufficient to sustain the population. A prioritization strategy for population conservation based on estimating the relative roles of different evolutionary forces shaping the gene pools highlighted a number of populations where further monitoring is warranted and also identified populations which could be prioritized for conservation as living gene banks in the event that conservation resources are limited. This prioritization agreed well with the occurrence of temporal (in)stability.

Finding Markers That Make a Difference: DNA Pooling and SNP-Arrays Identify Population Informative Markers for Genetic Stock Identification

Genetic stock identification (GSI) using molecular markers is an important tool for management of migratory species. Here, we tested a cost-effective alternative to individual genotyping, known as allelotyping, for identification of highly informative SNPs for accurate genetic stock identification. We estimated allele frequencies of 2880 SNPs from DNA pools of 23 Atlantic salmon populations using Illumina SNP-chip. We evaluated the performance of four common strategies (global F ST, pairwise F ST, Delta and outlier approach) for selection of the most informative set of SNPs and tested their effectiveness for GSI compared to random sets of SNP and microsatellite markers. For the majority of cases, SNPs selected using the outlier approach performed best followed by pairwise F ST and Delta methods. Overall, the selection procedure reduced the number of SNPs required for accurate GSI by up to 53% compared with randomly chosen SNPs. However, GSI accuracy was more affected by populations in the ascertainment group rather than the ranking method itself. We demonstrated for the first time the compatibility of different large-scale SNP datasets by compiling the largest population genetic dataset for Atlantic salmon to date. Finally, we showed an excellent performance of our top SNPs on an independent set of populations covering the main European distribution range of Atlantic salmon. Taken together, we demonstrate how combination of DNA pooling and SNP arrays can be applied for conservation and management of salmonids as well as other species.

Genomewide introgressive hybridization patterns in wild Atlantic salmon influenced by inadvertent gene flow from hatchery releases

Many salmonid fish populations are threatened by genetic homogenization, primarily due to introgressive hybridization with hatchery‐reared conspecifics. By applying genomewide analysis using two molecular marker types (1986 SNPs and 17 microsatellites), we assessed the genetic impacts of inadvertent gene flow via straying from hatchery releases on wild populations of Atlantic salmon in the Gulf of Finland, Baltic Sea, over 16 years (1996–2012). Both microsatellites and SNPs revealed congruent population genetic structuring, indicating that introgression changed the genetic make‐up of wild populations by increasing genetic diversity and reducing genetic divergence. However, the degree of genetic introgression varied among studied populations, being higher in the eastern part and lower in the western part of Estonia, which most likely reflects the history of past stocking activities. Using kernel smoothing and permutation testing, we detected considerable heterogeneity in introgression patterns across the genome, with a large number of regions exhibiting nonrandom introgression widely dispersed across the genome. We also observed substantial variation in nonrandom introgression patterns within populations, as the majority of genomic regions showing elevated or reduced introgression were not consistently detected among temporal samples. This suggests that recombination, selection and stochastic processes may contribute to complex nonrandom introgression patterns. Our results suggest that (i) some genomic regions in Atlantic salmon are more vulnerable to introgressive hybridization, while others show greater resistance to unidirectional gene flow; and (ii) the hybridization of previously separated populations leads to complex and dynamic nonrandom introgression patterns that most likely have functional consequences for indigenous populations.

Generation of a neutral FST baseline for testing local adaptation on gill raker number within and between European whitefish ecotypes in the Baltic Sea basin

Divergent selection at ecologically important traits is thought to be a major factor driving phenotypic differentiation between populations. To elucidate the role of different evolutionary processes shaping the variation in gill raker number of European whitefish (Coregonus lavaretus sensu lato) in the Baltic Sea basin, we assessed the relationships between genetic and phenotypic variation among and within three whitefish ecotypes (sea spawners, river spawners and lake spawners). To generate expected neutral distribution of FST and to evaluate whether highly variable microsatellite loci resulted in deflated FST estimates compared to less variable markers, we performed population genetic simulations under finite island and hierarchical island models. The genetic divergence observed among (FCT = 0.010) and within (FST = 0.014–0.041) ecotypes was rather low. The divergence in gill raker number, however, was substantially higher between sea and river spawners compared to observed microsatellite data and simulated neutral baseline (PCT > FCT). This suggests that the differences in gill raker number between sea and river spawners are likely driven by divergent natural selection. We also found strong support for divergent selection on gill raker number among different populations of sea spawners (PST > FST), most likely caused by highly variable habitat use and diverse diet. The putative role of divergent selection within lake spawners initially inferred from empirical microsatellite data was not supported by simulated FST distributions. This work provides a first formal test of divergent selection on gill raker number in Baltic whitefish, and demonstrates the usefulness of population genetic simulations to generate informative neutral baselines for PST–FST analyses helping to disentangle the effects of stochastic evolutionary processes from natural selection.

Temporal variation of genetic composition in Atlantic salmon populations from the Western White Sea Basin: influence of anthropogenic factors?

BackgroundStudies of the temporal patterns of population genetic structure assist in evaluating the consequences of demographic and environmental changes on population stability and persistence. In this study, we evaluated the level of temporal genetic variation in 16 anadromous and 2 freshwater salmon populations from the Western White Sea Basin (Russia) using samples collected between 1995 and 2008. To assess whether the genetic stability was affected by human activity, we also evaluated the effect of fishing pressure on the temporal genetic variation in this region.ResultsWe found that the genetic structure of salmon populations in this region was relatively stable over a period of 1.5 to 2.5 generations. However, the level of temporal variation varied among geographical regions: anadromous salmon of the Kola Peninsula exhibited a higher stability compared to that of the anadromous and freshwater salmon from the Karelian White Sea coast. This discrepancy was most likely attributed to the higher census, and therefore effective, population sizes of the populations inhabiting the rivers of the Kola Peninsula compared to salmon of the Karelian White Sea coast. Importantly, changes in the genetic diversity observed in a few anadromous populations were best explained by the increased level of fishing pressure in these populations rather than environmental variation or the negative effects of hatchery escapees. The observed population genetic patterns of isolation by distance remained consistent among earlier and more recent samples, which support the stability of the genetic structure over the period studied.ConclusionsGiven the increasing level of fishing pressure in the Western White Sea Basin and the higher level of temporal variation in populations exhibiting small census and effective population sizes, further genetic monitoring in this region is recommended, particularly on populations from the Karelian rivers.

Use of sibling relationship reconstruction to complement traditional monitoring in fisheries management and conservation of brown trout.

Declining trends in the abundance of many fish urgently call for more efficient and informative monitoring methods that would provide necessary demographic data for the evaluation of existing conservation, restoration, and management actions. We investigated how genetic sibship reconstruction from young-of-the-year brown trout (Salmo trutta L.) juveniles provides valuable, complementary demographic information that allowed us to disentangle the effects of habitat quality and number of breeders on juvenile density. We studied restored (n = 15) and control (n = 15) spawning and nursery habitats in 16 brown trout rivers and streams over 2 consecutive years to evaluate the effectiveness of habitat restoration activities. Similar juvenile densities both in restored and control spawning and nursery grounds were observed. Similarly, no differences in the effective number of breeders, Nb(SA) , were detected between habitats, indicating that brown trout readily used recently restored spawning grounds. Only a weak relationship between the Nb(SA) and juvenile density was observed, suggesting that multiple factors affect juvenile abundance. In some areas, very low estimates of Nb(SA) were found at sites with high juvenile density, indicating that a small number of breeders can produce a high number of progeny in favorable conditions. In other sites, high Nb(SA) estimates were associated with low juvenile density, suggesting low habitat quality or lack of suitable spawning substrate in relation to available breeders. Based on these results, we recommend the incorporation of genetic sibship reconstruction to ongoing and future fish evaluation and monitoring programs to gain novel insights into local demographic and evolutionary processes relevant for fisheries management, habitat restoration, and conservation.

Genetic differentiation between the Old and New types of Serbian Tsigai sheep

Two Tsigai sheep populations exist in Serbia: the Old type, called Čokan, and the New type. It is assumed that the New type results from upgrading Tsigai sheep with exotic genetic material. We investigated genetic diversity and differentiation of these types by analysing 23 autosomal microsatellites. Tests for Hardy-Weinberg proportions, linkage equilibrium between genotypes across loci and the calculation of inbreeding coefficients were performed and the deficiency in the number of alleles within the Tsigai types was examined using a Wilcoxon sign-rank test. The New type displayed a higher level of genetic variability than the Čokan in terms of allele numbers, but the New Tsigai showed a pattern of heterozygosity deficiency. The positive f value for the Čokan suggests the occurrence of inbreeding in this type. The proportion of linkage disequilibrium was below that expected by chance. Exclusion of two loci in Hardy-Weinberg disequilibrium did not alter our conclusions based on the entire data set i.e. the two Tsigai types are clearly differentiated and the New Tsigai type has been influenced by crossbreeding. Therefore, the Čokan Tsigai should be considered as a distinct endangered breed in the FAO classification.