Discovery and characterization of novel genetic markers for use in the management of Lahontan cutthroat trout (Oncorhynchus clarkii henshawi)

The Lahontan cutthroat trout (Oncorhynchus clarkii henshawi) is threatened by habitat destruction, over‐harvest and hybridization with nonnative trout. Currently, three Geographic Management Units (GMUs) are recognized within the taxon. Here, we describe a suite of 68 single‐nucleotide polymorphism (SNP) genetic markers for use in the study and management of Lahontan cutthroat trout and a closely related subspecies, the Paiute cutthroat trout (O. c. seleneris). These include markers variable within the two subspecies (n = 35), diagnostic for the two subspecies (n = 23) and diagnostic for Yellowstone cutthroat trout (O. c. bouvieri) and other closely related subspecies (n = 10). Sixty‐three markers were discovered by Sanger sequencing of 171 EST loci in an ascertainment panel including Lahontan cutthroat trout from four populations representing all GMUs. Five markers were identified in a secondary sequencing effort with a single population of Lahontan cutthroat trout. TaqMan assays were validated on six Lahontan cutthroat trout populations and a diverse panel of other trout. Over 90% of the markers variable in Lahontan cutthroat trout were polymorphic in at least two populations, and 66% were variable within all three GMUs. All Lahontan diagnostic markers were also fixed for the Lahontan allele in Paiute cutthroat trout. Most of the Yellowstone diagnostic markers can also be used for this purpose in other cutthroat trout subspecies. This is the first set of SNP markers to be developed for Lahontan cutthroat trout, and will be an important tool for conservation and management.     

A suite of twelve single nucleotide polymorphism markers for detecting introgression between cutthroat and rainbow trout

A suite of 12 subspecies and species-specific single nucleotide polymorphism (species-specific SNP) markers was developed to distinguish rainbow trout (RT) Oncorhynchus mykiss from the four major subspecies of cutthroat trout: westslope cutthroat trout (WCT) Oncorhynchus clarki lewisi, Yellowstone cutthroat trout (YCT) Oncorhynchus clarki bouvieri, coastal cutthroat trout (CCT) Oncorhynchus clarki clarki, Lahontan cutthroat trout (LCT) Oncorhynchus clarki henshawi, and their hybrids. Several of the markers were linked to help strengthen hybrid determinations, and sex-specific species-specific SNP assays were also developed.     

Ten species specific microsatellite loci for Lahontan cutthroat trout,Oncorhynchus clarki henshawi

Ten polymorphic microsatellite loci (containing tri and tetra‐nucleotide repeats) were developed for the Lahontan cutthroat trout (Oncorhynchus clarki henshawi), a subspecies of cutthroat trout listed as threatened under the United States Endangered Species Act. Polymorphism was assessed for 445 individuals from 12 populations representing eight watersheds spread throughout the three Distinct Population Segments defined for this subspecies. All loci were polymorphic (X? = 19, range 7–26 alleles). All loci were in Hardy–Weinberg equilibrium (HWE) except for one locus (OCH 9) in a single population (P < 0.00014 after Bonferroni correction for multiple tests).     

Assessing thermal adaptation using family‐based association and FST outlier tests in a threatened trout species

Discovering genetic markers associated with phenotypic or ecological characteristics can improve our understanding of adaptation and guide conservation of key evolutionary traits. The Lahontan cutthroat trout (Oncorhynchus clarkii henshawi) of the northern Great Basin Desert, USA, demonstrated exceptional tolerance to high temperatures in the desert lakes where it resided historically. This trait is central to a conservation hatchery effort to protect the genetic legacy of the nearly extinct lake ecotype. We genotyped full‐sibling families from this conservation broodstock and samples from the only two remaining, thermally distinct, native lake populations at 4,644 new single nucleotide polymorphisms (SNPs). Family‐based genome‐wide association testing of the broodstock identified nine and 26 SNPs associated with thermal tolerance (p < 0.05 and p < 0.1), measured in a previous thermal challenge experiment. Genes near the associated SNPs had complex functions related to immunity, growth, metabolism and ion homeostasis. Principal component analysis using the thermotolerance‐related SNPs showed unexpected divergence between the conservation broodstock and the native lake populations at these loci. F_ST outlier tests on the native lake populations identified 18 loci shared between two or more of the tests, with two SNPs identified by all three tests (p < 0.01); none overlapped with loci identified by association testing in the broodstock. A recent history of isolation and the complex genetic and demographic backgrounds of Lahontan cutthroat trout probably limited our ability to find shared thermal tolerance loci. Our study extends the still relatively rare application of genomic tools testing for markers associated with important phenotypic or environmental characteristics in species of conservation concern.     

Six diagnostic single nucleotide polymorphism markers for detecting introgression between cutthroat and rainbow trouts

Ten primer pairs were screened to develop single nucleotide polymorphism (SNP) TaqMan assays that will distinguish California golden trout and some rainbow trouts (Oncorhynchus mykiss sspp., O. m. aguabonita) from the Paiute and Lahontan cutthroat trouts (Oncorhynchus clarkii seleniris, O. c. henshawi). From these 10 primer pairs, one mitochondrial and five nuclear fixed SNP differences were discovered and developed into TaqMan assays. These six assays will be useful for characterizing and monitoring hybridization between these groups. Additional Oncorhynchus clarkii sspp. and Oncorhynchus mykiss sspp. were assayed to determine if these assays are useful in closely related species.     

Landscape topography and the distribution of Lahontan cutthroat trout (Oncorhynchus clarki henshawi) in a high desert stream

Lahontan cutthroat trout, Oncorhynchus clarki henshawi, are currently limited in their distribution to a patchwork of small isolated populations, the result of habitat degradation and natural variation in landscape and in-stream conditions. The objectives of this study were to determine if landscape topography influences trout distribution, and if water temperatures control this response. The work was carried out in a sub-basin of the Quinn River system, McDermitt Creek, which drains the sagebrush desert of southeastern Oregon and northern Nevada. Headwater tributaries of this creek consist of alternating canyon-confined and valley bounded reaches. Trout within these systems are challenged by low discharge and high temperatures during the summer, and anchor ice during the winter. Contiguous whole stream surveys were used to look at trout distribution during the summer of 2003 and spring and fall of 2004. Our results suggested that topography can affect trout distribution. Trout numbers were highest in areas with greater numbers of nick-points (the transition zones between less confined and more confined valley segments) and greater valley confinement. Additionally, in the downstream portion of our headwater reaches, more trout were found in nick-points than expected based on the availability of this habitat type. Our data suggest that hyporheic inputs may be high in such areas, thus providing trout with shelter from warm water in the summer, anchor ice in the winter, and shallow stream depths during all seasons. Spatial occurrence of these areas of refugia can be taken into consideration when planning land use activities and restoration efforts. Further research is required to confirm that topography can affect the distribution of Lahontan cutthroat trout in other systems, and to better understand the mechanisms behind these patterns.     

Next-generation RAD sequencing identifies thousands of SNPs for assessing hybridization between rainbow and westslope cutthroat trout

The increased numbers of genetic markers produced by genomic techniques have the potential to both identify hybrid individuals and localize chromosomal regions responding to selection and contributing to introgression. We used restriction-site-associated DNA sequencing to identify a dense set of candidate SNP loci with fixed allelic differences between introduced rainbow trout (Oncorhynchus mykiss) and native westslope cutthroat trout (Oncorhynchus clarkii lewisi). We distinguished candidate SNPs from homeologs (paralogs resulting from whole-genome duplication) by detecting excessively high observed heterozygosity and deviations from Hardy-Weinberg proportions. We identified 2923 candidate species-specific SNPs from a single Illumina sequencing lane containing 24 barcode-labelled individuals. Published sequence data and ongoing genome sequencing of rainbow trout will allow physical mapping of SNP loci for genome-wide scans and will also provide flanking sequence for design of qPCR-based TaqMan(?) assays for high-throughput, low-cost hybrid identification using a subset of 50-100 loci. This study demonstrates that it is now feasible to identify thousands of informative SNPs in nonmodel species quickly and at reasonable cost, even if no prior genomic information is available.     

Characterization of tetranucleotide microsatellites for Rio Grande cutthroat trout and rainbow trout,and their cross‐amplification in other cutthroat trout subspecies

We describe the isolation and characterization of 12 tetranucleotide microsatellites for Rio Grande cutthroat trout (Oncorhynchus clarkii virginalis) and rainbow trout (Oncorhynchus mykiss), and subsequently investigate their performance in Colorado River cutthroat trout (Oncorhynchus clarkii pleuriticus), greenback cutthroat trout (Oncorhynchus clarkii stomias) and Yellowstone cutthroat trout (Oncorhynchus clarki bouvieri). All 12 loci are polymorphic in all subspecies of O. clarkii examined.     

Development and characterization of novel tetra‐, tri‐, and dinucleotide microsatellite markers in rainbow trout (Oncorhynchus mykiss)

We discuss the development and characterization of 40 polymorphic rainbow trout (Oncorhynchus mykiss) microsatellite loci. We used enriched libraries to isolate 14 dinucleotide, seven trinucleodide, eight compound di/tetranucleotide, and 11 tetranucleotide loci. These markers will be useful for selective breeding via marker‐assisted selection, population genetics studies, parentage analysis, and have already been used for genome mapping.     

Isolation of 15 single nucleotide polymorphisms from coastal steelhead, Oncorhynchus mykiss (Salmonidae)

We describe 15 single nucleotide polymorphisms (SNPs) isolated in coastal California populations of steelhead (Oncorhynchus mykiss). SNP loci were developed using a 'gene-targeted' approach, which involved the development of primers from functional genes in O. mykiss that were deposited in GenBank or in the published literature. These markers show a wide range of variability in three coastal steelhead populations, and will be useful in population genetic studies and in pedigree reconstruction. Potential applications include evaluation of population structure, introgression between native and hatchery trout, and evaluating reproductive success.     

Development and evaluation of 200 novel SNP assays for population genetic studies of westslope cutthroat trout and genetic identification of related taxa

DNA sequence data were collected and screened for single nucleotide polymorphisms (SNPs) in westslope cutthroat trout (Oncorhynchus clarki lewisi) and also for substitutions that could be used to genetically discriminate rainbow trout (O. mykiss) and cutthroat trout, as well as several cutthroat trout subspecies. In total, 260 expressed sequence tag‐derived loci were sequenced and allelic discrimination genotyping assays developed from 217 of the variable sites. Another 50 putative SNPs in westslope cutthroat trout were identified by restriction‐site‐associated DNA sequencing, and seven of these were developed into assays. Twelve O. mykiss SNP assays that were variable within westslope cutthroat trout and 12 previously published SNP assays were also included in downstream testing. A total of 241 assays were tested on six westslope cutthroat trout populations (N = 32 per population), as well as collections of four other cutthroat trout subspecies and a population of rainbow trout. All assays were evaluated for reliability and deviation from Hardy–Weinberg and linkage equilibria. Poorly performing and duplicate assays were removed from the data set, and the remaining 200 assays were used in tests of population differentiation. The remaining markers easily distinguished the various subspecies tested, as evidenced by mean G_ST of 0.74. A smaller subset of the markers (N = 86; average G_ST = 0.40) was useful for distinguishing the six populations of westslope cutthroat trout. This study increases by an order of magnitude the number of genetic markers available for the study of westslope cutthroat trout and closely related taxa and includes many markers in genes (developed from ESTs).     

Characterization and comparison of microsatellites derived from repeat-enriched libraries and expressed sequence tags

The construction of high-density linkage maps for use in identifying loci underlying important traits requires the development of large numbers of polymorphic genetic markers spanning the entire genome at regularly spaced intervals. As part of our efforts to develop markers for rainbow trout (Oncorhynchus mykiss), we performed a comparison of allelic variation between microsatellite markers developed from expressed sequence tag (EST) data and anonymous markers identified from repeat-enriched libraries constructed from genomic DNA. A subset of 70 markers (37 from EST databases and 33 from repeat enriched libraries) was characterized with respect to polymorphism information content (PIC), number of alleles, repeat number, locus duplication within the genome and ability to amplify in other salmonid species. Higher PIC was detected in dinucleotide microsatellites derived from ESTs than anonymous markers (72.7% vs. 54.0%). In contrast, dinucleotide repeat numbers were higher for anonymous microsatellites than for EST derived microsatellites (27.4 vs.18.1). A higher rate of cross-species amplification was observed for EST microsatellites. Approximately half of each marker type was duplicated within the genome. Unlike single-copy markers, amplification of duplicated microsatellites in other salmonids was not correlated to phylogenetic distance. Genomic microsatellites proved more useful than EST derived microsatellites in discriminating among the salmonids. In total, 428 microsatellite markers were developed in this study for mapping and population genetic studies in rainbow trout.     

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.     

Characterization of twenty‐four microsatellite markers for rainbow trout (Oncorhynchus mykiss)

Twenty‐four new microsatellite markers were developed for genome mapping and population genetics studies in rainbow trout (Oncorhynchus mykiss). The amount of polymorphism, percentage heterozygosity and ability of each marker to amplify genomic DNA from other salmonids were recorded. Seven markers were observed to be duplicated in the rainbow trout genome by containing more than one allele in homozygous (clonal) fish.     

Quantitative trait loci (QTLs) associated with resistance/susceptibility to infectious pancreatic necrosis virus (IPNV) in rainbow trout (Oncorhynchus mykiss)

Infectious pancreatic necrosis (IPN) is a well-known acute viral disease of salmonid species. We have identified quantitative trait loci (QTLs) associated with resistance to this disease in rainbow trout. We searched for linkage among 51 microsatellite markers used to construct a framework linkage map in backcross families of rainbow trout (Oncorhynchus mykiss), produced by crossing IPN-resistant (YN-RT201) and -susceptible (YK-RT101) strains. Two putative QTLs affecting disease resistance were detected on chromosomes A (IPN R S-1) and C (IPN R/S-2), respectively, suggesting that this is a polygenic trait in rainbow trout. These markers have great potential for use in marker-assisted selection (MAS) for IPN resistance and provide the basis for cloning of IPN resistance genes. Clarification of the genetic bases of complex traits has broad implications for fundamental research, but will also be of practical benefit to fish breeding.     

Characterization of 38 polymorphic microsatellite markers for rainbow trout (Oncorhynchus mykiss)

Thirty‐eight new microsatellite markers were developed for genome mapping and population genetics studies in rainbow trout (Oncorhynchus mykiss). The amount of polymorphism, percentage of heterozygosity and ability of each marker to amplify genomic DNA from other salmonids were recorded. Five markers were observed to be duplicated in the rainbow trout genome by containing more than one allele in homozygous (clonal) fish.     

Distribution of temperature tolerance quantitative trait loci in Arctic charr (Salvelinus alpinus) and inferred homologies in rainbow trout (Oncorhynchus mykiss)

We searched for quantitative trait loci (QTL) affecting upper temperature tolerance (UTT) in crosses between the Nauyuk Lake and Fraser River strains of Arctic charr (Salvelinus alpinus) using survival analysis. Two QTL were detected by using two microsatellite markers after correcting for experiment-wide error. A comparative mapping approach localized these two QTL to homologous linkage groups containing UTT QTL in rainbow trout (Oncorhynchus mykiss). Additional marginal associations were detected in several families in regions homologous to those with QTL in rainbow trout. Thus, the genes underlying UTT QTL may antedate the divergence of these two species, which occurred by approximately 16 MYA. The data also indicate that one pair of homeologs (ancestrally duplicated chromosomal segments) have contained QTL in Arctic charr since the evolution of salmonids from a tetraploid ancestor 25-100 MYA. This study represents one of the first examples of comparative QTL mapping in an animal polyploid group and illustrates the fate of QTL after genome duplication and reorganization.     

HETEROZYGOSITY‐FITNESS CORRELATIONS IN RAINBOW TROUT: EFFECTS OF ALLOZYME LOCI OR ASSOCIATIVE OVERDOMINANCE?

Previous studies with rainbow trout (Oncorhynchus mykiss) have shown that increased heterozygosity at allozyme loci is correlated with several phenotypic traits associated with fitness. We expected to find a similar effect of heterozygosity at other nuclear loci if these associations are due to loci in linkage disequilibrium with the allozyme loci (i.e., associative overdominance), rather than the allozymes themselves. We examined the association between multiple locus heterozygosity and condition factor at 10 allozyme and 10 microsatellite loci. Individuals that were more heterozygous at allozyme loci had significantly greater condition factor in two hatchery cohorts of rainbow trout (1996 P = 0.006; 1997 P < 0.001). In contrast, there was no evidence at microsatellite loci that increased heterozygosity was associated with greater condition factor. Our results suggest that the observed relationship between heterozygosity at allozyme loci and condition factor in rainbow trout appears to be due to the allozyme loci themselves, rather than associative overdominance. We cannot, however, rule out that differences in the mutation process between allozymes and microsatellites may be responsible for these observations. Regardless of the underlying mechanism, these results support the view that allozymes and microsatellites are differentially affected by natural selection.     

Dinucleotide-repeat polymorphism in DNA of rainbow trout and its application in fisheries science

Dinucleotide–repeat DNA polymorphisms for rainbow trout, Oncorhynchus mykiss , are described. The potential applications of these markers in fisheries science are discussed.     

Novel molecular markers differentiate Oncorhynchus mykiss (rainbow trout and steelhead) and the O. clarki (cutthroat trout) subspecies

A suite of 26 PCR‐based markers was developed that differentiates rainbow (Oncorhynchus mykiss) and coastal cutthroat trout (O. clarki clarki). The markers also differentiated rainbow from other cutthroat trout subspecies (O. clarki), and several of the markers differentiated between cutthroat trout subspecies. This system has numerous positive attributes, including: nonlethal sampling, high species‐specificity and products that are easily identified and scored using agarose gel electrophoresis. The methodology described for developing the markers can be applied to virtually any system in which numerous markers are desired for identifying or differentiating species or subspecies.