Rapid communication: Thirty-eight polymorphic microsatellite markers for mapping in rainbow trout

Conservation Genetics of Mediterranean Brown Trout in Central Italy (Latium): A Multi-Marker Approach

Brown trout is considered a complex of incipient species, including several phylogenetic lineages, whose natural distribution in the Mediterranean area has been altered, since the beginning of the 1900s, by massive introductions of domestic strains of Atlantic origin to support fisheries. Introduced trout naturalize in new suitable environments and extensively hybridize with native populations. Here, we characterized putatively neutral and adaptive genetic variability and population structure of Mediterranean brown trout from six river catchments in central peninsular Italy, as revealed by both mitochondrial (Control Region) and nuclear (microsatellites, LDH-C1, major histocompatibility complex) markers. We quantified the admixture of wild populations with hatchery strains and evaluated the effects of domestic trout introductions on shaping population genetics. Our analyses indicated: (1) a composite picture of genetic variability in the area, with the presence of all native Mediterranean trout mitochondrial lineages (“Adriatic”, “Mediterranean”, “marmoratus”), various frequencies of allochthonous genotypes and different rates of introgression among sampling sites; (2) asymmetric mito-nuclear introgression; (3) increasing nuclear marker diversity with increasing levels of admixture across populations; (4) strong population structure coupled with relatively low effective population size. Data allowed the identification of five management units and we propose specific actions to support ongoing and future conservation strategies within the examined area.

Comparative genetic variability of pink salmon from different parts of their range: native Pacific, artificially introduced White Sea and naturally invasive Atlantic Scottish rivers

Trans-oceanic movement, stocking and subsequent establishment of Pacific pink salmon (Oncorhynchus gorbuscha) into the Atlantic White Sea area have resulted in their spreading further across the northern Atlantic, with spawning being reported in a number of regions within this area. Such expansions of non-native species bring potential risks to the ecosystems in question. It has not yet been established if the spawning events of pink salmon observed are the result of self-sustaining populations in these areas, or are because of repeated invasions of strayers from the White Sea stocks. In 2017 pink salmon were observed in a number of Scottish rivers in historically large numbers. This study set out to examine genetic variation in these fish and compare this to fish in Pacific founder regions and the White Sea translocated populations. A total of 286 samples from Scotland, the Atlantic White Sea, the Pacific Okhotsk region and Northern Pacific Bering Sea were screened using a 1018 bp sequenced region of the Cytochrome b mtDNA gene and 205 of these samples for 13 microsatellites. Significant bottleneck and founder effects were observed in the White Sea stocks in both mitochondrial and nuclear DNA, including loss of diversity and changes in haplotype and allele proportions. Scottish fish were indistinguishable from White Sea populations and as such it was not possible to determine if the fish were strayers from this region or returning fish from previous spawning events in Scotland. Therefore, although the fish caught in Scotland have their origins in the White Sea population, it may not be easy to determine whether self-sustaining populations have, or are becoming, established in the UK using genetic analysis and other techniques may need to be used.

Genetic structure of rainbow trout Oncorhynchus mykiss (Salmoniformes, Salmonidae) from aquaculture by DNA-markers

The rational use of valuable fish species from aquaculture is difficult to implement without knowledge of the state of the genetic structure of local stocks. Different types of DNA markers can be used to achieve the goals of selection and breeding work. The genetic structure of a local stock of rainbow trout (Oncorhynchus mykiss Walbaum, 1792) (Salmoniformes, Salmonidae) farmed in Ukraine was studied using DNA-markers: microsatellite (SSR-markers – simple-sequence repeats-markers) and intermicrosatellite (ISSR – inter-simple sequence repeat). Five fragments of trinucleotide microsatellite motifs with a single anchor nucleotide at the 3'-end were used as a primer for analysis by the ISSR-PCR method. Totally, 85 amplicons were obtained across the five loci, of which 92.9% were polymorphic. The total number of alleles ranged from 10 (marker (ACC)₆G) to 23 (marker (AGC)₆G). The following monomorphic amplicons were determined for the studied local stock of rainbow trout: according to marker (CTC)₆C – 770 and 520 bp bands, for the marker (GAG)₆C – 345, 295 and 260 bp, and for the marker (AGC)₆C – 350 bp. The average number of polymorphic bands per locus was 15.8. The selected ISSR primers had a level of polymorphic information content above the average. The most effective markers for molecular-genetic analysis of rainbow trout were (AGC)₆G and (AGC)₆C according to the percentage of polymorphic bands, marker index, effective multiplex ratio and resolving power. The selected ISSR loci allow the genetic structure of the studied local stock to be characterized using the total and the effective number of alleles per locus (Na and Ne were 1.9 and 1.4, respectively), the Shannon index (average value I was 0.4) and the unbiased expected heterozygosity (mean uHe = 0.3). Microsatellite-based analysis showed features of the genetic structure of the local stock of rainbow trout at six microsatellite loci (OMM 1032, OMM 1077, OMM 1088, Str 15, Str 60, Str 73). Allelic diversity was established and alleles with the highest frequency and most typical for the given stock were identified. The Shannon index and unbiased expected heterozygosity were determined using SSR-markers and were 1.42 and 0.79, respectively. This depicts the complexity of the population structure, a high level of genetic diversity and indicates a high level of heterozygosity of local stock. The “gene pool profile” established as a result of ISSR-PCR in the future will help to differentiate local stocks of rainbow trout in aquaculture of Ukraine. Microsatellite markers provide the ability to determine individual features of genetic variation of local populations and to conduct the management of genetic resources on fish farms.

The genetic status of the Hungarian brown trout populations: exploration of a blind spot on the European map of Salmo trutta studies

BACKGROUND

Analyses of the control region sequences of European brown trout populations' mitrochondrial DNA have revealed five main evolutionary lineages (Atlantic, Danubian, Mediterranean, Adriatic, Marble) mostly relating to the main water basins; however, the hybridization between lineages were increasingly reported. Due to the hydrogeography of Hungary, wild populations should theoretically belong to the Danubian lineage, however, this has not been verified by genetic studies.

METHODS

In our study multiple molecular marker sets (mitochondrial sequence, microsatellites, PCR-RFLP of nuclear markers and sex marker) were used to investigate the genetic composition and population genetics of the brown trout populations in two broodstocks, six wild streams in Hungary and one Serbian population.

RESULTS

The admixture of Atlantic and Danubian lineages in these populations, except the Serbian population with pure Danubian origin, was observed by control region sequences of mitochondrial DNA and PCR-RFLP markers in the nuclear genome, and one unpublished Danubian haplotype was found in Hungarian populations. A sex-specific marker revealed equal gender ratio in broodstocks and Kemence stream, whereas in other wild streams the proportion of female individuals were less than 50%. Structure and principal component analyses based on the alleles of microsatellite loci also revealed overlapping populations, however the populations were still significantly different from each other and were mostly in Hardy-Weinberg equilibrium.

DISCUSSION

Stocking and migration can have a significant genetic impact on trout populations of wild streams, however there are no guidelines or common practices for stocking of small streams in Hungary, thus the genetic background of these populations should be considered when developing conservation actions.

Developmental competence of eggs produced by rainbow trout Doubled Haploids (DHs) and generation of the clonal lines

Poor quality eggs produced by the fully homozygous doubled haploids (DHs) may impair generation of clonal lines in fish species. In the present research, gynogenetic development of rainbow trout (Oncorhynchus mykiss) was induced in eggs originated from the DH females. Eggs were activated with the UV-irradiated grayling (Thymallus thymallus) spermatozoa and subjected to the high hydrostatic pressure (HHP) shock to provide diploid clonal individuals. Only two of four DH females produced eggs that were successfully activated by the irradiated spermatozoa and subsequently developed into the gynogenetic embryos. Survival rates of rainbow trout from the clonal lines equalled 21.5% and 19.8% during embryogenesis and decreased after hatching to 18.6% and 14.9%, respectively. Some of the dead rainbow trout clones collected between hatching and swim-up stage were emaciated and exhibited spinal deformities including scoliosis. Provided results confirmed limited developmental competences of eggs produced by rainbow trout DH females. Clonal rainbow trout developing in such eggs exhibited reduced survival and increased frequency of the body abnormalities.

Low levels of hybridization between sympatric Arctic char (Salvelinus alpinus) and Dolly Varden char (Salvelinus malma) highlights their genetic distinctiveness and ecological segregation

Understanding the extent of interspecific hybridization and how ecological segregation may influence hybridization requires comprehensively sampling different habitats over a range of life history stages. Arctic char (Salvelinus alpinus) and Dolly Varden (S. malma) are recently diverged salmonid fishes that come into contact in several areas of the North Pacific where they occasionally hybridize. To better quantify the degree of hybridization and ecological segregation between these taxa, we sampled over 700 fish from multiple lake (littoral and profundal) and stream sites in two large, interconnected southwestern Alaskan lakes. Individuals were genotyped at 12 microsatellite markers, and genetic admixture (Q) values generated through Bayesian-based clustering revealed hybridization levels generally lower than reported in a previous study (<0.6% to 5% of samples classified as late-generation hybrids). Dolly Varden and Arctic char tended to make different use of stream habitats with the latter apparently abandoning streams for lake habitats after 2-3 years of age. Our results support the distinct biological species status of Dolly Varden and Arctic char and suggest that ecological segregation may be an important factor limiting opportunities for hybridization and/or the ecological performance of hybrid char.

The evolutionary jigsaw puzzle of the surviving trout (Salmo trutta L. complex) diversity in the Italian region. A multilocus Bayesian approach

Mediterranean trout populations display a diversity of phenotypes, representing a valuable model for the study of adaptation and a puzzling dilemma for taxonomists and biogeographers, which is further entangled by the widespread introgression of allochthonous genes. In this paper we analysed DNA polymorphism at multiple loci (sequence variation of the mitochondrial control region and eight nuclear fragments and length variation at eleven nuclear microsatellite loci) in representative samples of the autochthonous taxonomic diversity described in Italian trout populations (Salmo marmoratus, S. carpio, S. cenerinus, S. cettii and S. fibreni) and in samples from hatchery-originated strains of Atlantic S. trutta. We employed model-based clustering and Approximate Bayesian Computation in order to: (i) describe the phylogeographic structure of Italian autochthonous trout populations; (ii) evaluate a set of evolutionary/biogeographic models. The inclusion of hatchery-originated strains allowed to account for man-mediated allochthonous introgression in Italian populations. Our results (i) showed that the analysed sample consists of two main autochthonous evolutionary lineages, including the marble trout populations on one side ('marble' lineage) and the three peninsular populations of S. cettii, S. cenerinus and S. fibreni on the other side ('peninsular' lineage); (ii) indicated that S. carpio originated from a 'peninsular' population, with a possible, limited contribution from the 'marble' lineage; (iii) pointed out that the 'marble' lineage started diverging before the separation of the 'peninsular' lineage from Atlantic S. trutta; (iv) suggested that a model of divergence involving gene flow from the 'peninsular' population into the ancestral gene pool of 'marble' trout is most consistent with the genetic data; (v) provided evidence that the autochthonous trout gene pools in the Tyrrhenian and Adriatic basins of the Italian peninsula started diverging very recently (most likely after the last glacial maximum).

Genetic diversity and conservation of the Prespa trout in the Balkans

The Balkans are known to have a high level of biodiversity and endemism. No less than 15 taxa have been recorded in salmonids of the Salmo genus. Among them, the Prespa trout is found in only four river systems flowing into Lake Macro Prespa, three in the Former Yugoslav Republic of Macedonia and one in Greece. This is the first comprehensive survey of all streams located within the Macro Prespa Basin, encompassing the whole taxon range. A large genetic sample of 536 Prespa trout was collected mainly between 2005 and 2007. The sampling included 59 individuals from the Golema river system, 93 from the Kranska, 260 from the Brajcinska, 119 from the Agios Germanos, and five individuals from the lake itself. These specimens were analyzed with six microsatellite markers and by sequencing the mitochondrial control region. Nuclear data were examined through multidimensional analysis and assignment tests. Five clusters were detected by assignment: Golema, Kranska, Brajcinska upstream, Rzanska Brajcinska tributary and Brajcinska downstream. Most of these river systems thus hosted differentiated Prespa trout populations (with past gene flows likely dating before the construction of dams), except Agios Germanos, which was found to be composed of 5% to 32% of each cluster. Among the five trout individuals from the lake, four originated from Kranska River and one was admixed. Supported parsimonious hypotheses are proposed to explain these specificities. Conservation of this endemic taxon should take these results into account. No translocation should be performed between different tributaries of the lake and preservation of the Brajcinska populations should address the upstream-downstream differentiation described.

Winning the invasion roulette: escapes from fish farms increase admixture and facilitate establishment of non-native rainbow trout

Aquaculture is a major source of invasive aquatic species, despite the fact that cultured organisms often have low genetic diversity and tend to be maladapted to survive in the wild. Yet, to what extent aquaculture escapees become established by means of high propagule pressure and multiple origins is not clear. We analysed the genetic diversity of 15 established populations and four farmed stocks of non-native rainbow trout in Chile, a species first introduced for recreational fishing around 1900, but which has in recent decades escaped in large numbers from fish farms and become widespread. Aquaculture propagule pressure was a good predictor of the incidence of farm escapees, which represented 16% of all free-ranging rainbow trout and were present in 80% of the study rivers. Hybrids between farm escapes and established trout were present in all rivers at frequencies ranging between 7 and 69%, and population admixture was positively correlated with genetic diversity. We suggest that non-native salmonids introduced into the Southern Hemisphere could benefit from admixture because local adaptations may not have yet developed, and there may be initially little fitness loss resulting from outbreeding depression.

Estimates of linkage disequilibrium and effective population size in rainbow trout

Background

The use of molecular genetic technologies for broodstock management and selective breeding of aquaculture species is becoming increasingly more common with the continued development of genome tools and reagents. Several laboratories have produced genetic maps for rainbow trout to aid in the identification of loci affecting phenotypes of interest. These maps have resulted in the identification of many quantitative/qualitative trait loci affecting phenotypic variation in traits associated with albinism, disease resistance, temperature tolerance, sex determination, embryonic development rate, spawning date, condition factor and growth. Unfortunately, the elucidation of the precise allelic variation and/or genes underlying phenotypic diversity has yet to be achieved in this species having low marker densities and lacking a whole genome reference sequence. Experimental designs which integrate segregation analyses with linkage disequilibrium (LD) approaches facilitate the discovery of genes affecting important traits. To date the extent of LD has been characterized for humans and several agriculturally important livestock species but not for rainbow trout.

Results

We observed that the level of LD between syntenic loci decayed rapidly at distances greater than 2 cM which is similar to observations of LD in other agriculturally important species including cattle, sheep, pigs and chickens. However, in some cases significant LD was also observed up to 50 cM. Our estimate of effective population size based on genome wide estimates of LD for the NCCCWA broodstock population was 145, indicating that this population will respond well to high selection intensity. However, the range of effective population size based on individual chromosomes was 75.51 - 203.35, possibly indicating that suites of genes on each chromosome are disproportionately under selection pressures.

Conclusions

Our results indicate that large numbers of markers, more than are currently available for this species, will be required to enable the use of genome-wide integrated mapping approaches aimed at identifying genes of interest in rainbow trout.

Genetic characterisation of farmed rainbow trout in Norway: intra- and inter-strain variation reveals potential for identification of escapees

Background

The rainbow trout (Oncorhynchus mykiss) is one of the most important aquaculture species in the world, and Norway is one of the largest producers. The present study was initiated in response to a request from the Norwegian police authority to identify the farm of origin for 35 escaped rainbow trout captured in a fjord. Eleven samples, each consisting of approximately 47 fish, were collected from the three farms operating in the fjord where the escapees were captured. In order to gain a better general understanding of the genetic structure of rainbow trout strains used in Norwegian aquaculture, seven samples (47 fish per sample) were collected from six farms located outside the region where the escapees were captured. All samples, including the escapees, were genotyped with 12 microsatellite loci.

Results

All samples displayed considerable genetic variability at all loci (mean number of alleles per locus per sample ranged from 5.4–8.6). Variable degrees of genetic differentiation were observed among the samples, with pair-wise F_ST values ranging from 0–0.127. Self-assignment tests conducted among the samples collected from farms outside the fjord where the escapees were observed gave an overall correct assignment of 82.5%, demonstrating potential for genetic identification of escapees. In the "real life" assignment of the 35 captured escapees, all were excluded from two of the samples included as controls in the analysis, and 26 were excluded from the third control sample. In contrast, only 1 of the escapees was excluded from the 11 pooled samples collected on the 3 farms operating in the fjord.

Conclusion

Considerable genetic variation exists within and among rainbow trout strains farmed in Norway. Together with modern statistical methods, this will provide commercial operators with a tool to monitor breeding and fish movements, and management authorities with the ability to identify the source of escapees. The data generated in this study were used by the Norwegian police to initiate an investigation of the company operating the three farms in the fjord where escapees were observed.

A second generation genetic map for rainbow trout (Oncorhynchus mykiss)

Background

Genetic maps characterizing the inheritance patterns of traits and markers have been developed for a wide range of species and used to study questions in biomedicine, agriculture, ecology and evolutionary biology. The status of rainbow trout genetic maps has progressed significantly over the last decade due to interest in this species in aquaculture and sport fisheries, and as a model research organism for studies related to carcinogenesis, toxicology, comparative immunology, disease ecology, physiology and nutrition. We constructed a second generation genetic map for rainbow trout using microsatellite markers to facilitate the identification of quantitative trait loci for traits affecting aquaculture production efficiency and the extraction of comparative information from the genome sequences of model fish species.

Results

A genetic map ordering 1124 microsatellite loci spanning a sex-averaged distance of 2927.10 cM (Kosambi) and having 2.6 cM resolution was constructed by genotyping 10 parents and 150 offspring from the National Center for Cool and Cold Water Aquaculture (NCCCWA) reference family mapping panel. Microsatellite markers, representing pairs of loci resulting from an evolutionarily recent whole genome duplication event, identified 180 duplicated regions within the rainbow trout genome. Microsatellites associated with genes through expressed sequence tags or bacterial artificial chromosomes produced comparative assignments with tetraodon, zebrafish, fugu, and medaka resulting in assignments of homology for 199 loci.

Conclusion

The second generation NCCCWA genetic map provides an increased microsatellite marker density and quantifies differences in recombination rate between the sexes in outbred populations. It has the potential to integrate with cytogenetic and other physical maps, identifying paralogous regions of the rainbow trout genome arising from the evolutionarily recent genome duplication event, and anchoring a comparative map with the zebrafish, medaka, tetraodon, and fugu genomes. This resource will facilitate the identification of genes affecting traits of interest through fine mapping and positional cloning of candidate genes.

Amplification of rainbow trout microsatellites in Brachymystax lenok

Brachymystax lenok is a cold freshwater fish accustomed to inhabit relatively high concentration of dissolved oxygen and clean upper streams. Here we present 13 polymorphic microsatellite primer pairs from rainbow trout to amplify in 32 B. lenok individuals from Ussuri River of China. The number of alleles ranged from two to seven with an average of 3.9 per locus. Observed heterozygosity ranged from 0.0625 to 0.9677. One locus showed significant deviation from Hardy-Weinberg equilibrium. These 13 loci will provide a good basis for investigation of B. lenok population structure and genetic diversity in different distribution region.

Characterization of a chromosomal rearrangement responsible for producing "apparent" XY-female fall-run Chinook salmon in California

Fluorescence in situ hybridization (FISH) was used to identify the X and Y chromosomes of offspring produced by normal and "apparent" XY-female fall-run Chinook salmon (Oncorhynchus tshawytscha) from California. FISH experiments were performed using probes to 2 sex-linked loci, growth hormone pseudogene (GH-Psi), and OtY1, as well as a probe to a sex-linked microsatellite (Omy7INRA). Comparison of FISH staining patterns between the offspring produced by normal and apparent XY-females revealed that the apparent XY-female examined transmitted a "Y-like" chromosome with an attenuated OtY1 and GH-Psi signal to half of its offspring. Segregation analysis of microsatellites derived from rainbow trout (Oncorhynchus mykiss) with respect to phenotypic sex was carried out for 2 normal and 2 apparent XY-female Chinook salmon families. Inheritance patterns of Omy7INRA were consistent with this locus being closely linked to GH-Psi in males and in apparent XY-females carrying the Y-like chromosome. Another microsatellite locus (Omm1077) was closely linked to the primary sex-determining locus (SEX) in males but not to GH-Psi/OtY1 in apparent XY-females. The FISH analyses suggest that apparent XY-female fall-run Chinook salmon in California are not the product of a Y chromosome to autosome translocation. Despite the combined FISH and inheritance analyses, we were unable to differentiate between 2 alternative explanations for apparent XY-females, namely, recombination of markers between the sex chromosomes, or a Y chromosome with a dysfunctional or missing sex-determining region.

A consolidated linkage map for rainbow trout (Oncorhynchus mykiss)

Androgenetic doubled haploid progeny produced from a cross between the Oregon State University and Arlee clonal rainbow trout (Oncorhynchus mykiss) lines, used for a previous published rainbow trout map, were used to update the map with the addition of more amplified fragment length polymorphic (AFLP) markers, microsatellites, type I and allozyme markers. We have added more than 900 markers, bringing the total number to 1359 genetic markers and the sex phenotype including 799 EcoRI AFLPs, 174 PstI AFLPs, 226 microsatellites, 72 VNTR, 38 SINE markers, 29 known genes, 12 minisatellites, five RAPDs, and four allozymes. Thirty major linkage groups were identified. Synteny of linkage groups in our map with the outcrossed microsatellite map has been established for all except one linkage group in this doubled haploid cross. Putative homeologous relationships among linkage groups, resulting from the autotetraploid nature of the salmonid genome, have been revealed based on the placement of duplicated microsatellites and type I loci.