SNP genotyping by the 5'-nuclease reaction: advances in high-throughput genotyping with nonmodel organisms

Hatchery supplementation provides a demographic boost but alters age composition of sockeye salmon in Auke Lake, Southeast Alaska

Evaluating salmon hatchery supplementation programs requires assessing not only program objectives but identifying potential risks to wild populations as well. Such evaluations can be hampered by difficulty in distinguishing between hatchery- and wild-born returning adults. Here, we conducted 3 years (2011-2013) of experimental hatchery supplementation of sockeye salmon in Auke Lake, Juneau, Alaska where a permanent weir allows sampling and genotyping of every returning adult (2008-2019). We identified both hatchery- and wild-born returning adults with parentage assignment, quantified the productivity (adult offspring/spawner) of hatchery spawners relative to that of wild spawners, and compared run timing, age, and size at age between hatchery- and wild-born adults. Hatchery-spawning females produced from approximately six to 50 times more returning adults than did naturally spawning females. Supplementation had no discernable effect on run timing and limited consequences for size at age, but we observed a distinct shift to younger age at maturity in the hatchery-born individuals in all three brood years. The shift appeared to be driven by hatchery-born fish being more likely to emigrate after one, rather than two, years in the lake but the cause is unknown. In cases when spawning or incubation habitat is limiting sockeye salmon production, hatchery supplementation can be effective for enhancing the number of returning adult fish but not without the risk of phenotypic change in the recipient population, which can be an undesired outcome of hatchery supplementation. This study adds to a growing body of evidence suggesting that phenotypic change within a single generation of captive spawning might be widespread in salmon hatchery programs.

Pharmacogenomics Beyond Single Common Genetic Variants: The Way Forward

Interindividual variability in genes encoding drug-metabolizing enzymes, transporters, receptors, and human leukocyte antigens has a major impact on a patient's response to drugs with regard to efficacy and safety. Enabled by both technological and conceptual advances, the field of pharmacogenomics is developing rapidly. Major progress in omics profiling methods has enabled novel genotypic and phenotypic characterization of patients and biobanks. These developments are paralleled by advances in machine learning, which have allowed us to parse the immense wealth of data and establish novel genetic markers and polygenic models for drug selection and dosing. Pharmacogenomics has recently become more widespread in clinical practice to personalize treatment and to develop new drugs tailored to specific patient populations. In this review, we provide an overview of the latest developments in the field and discuss the way forward, including how to address the missing heritability, develop novel polygenic models, and further improve the clinical implementation of pharmacogenomics.

Uniparental and biparental markers unravel invasion pathways, population admixture, and genetic structure in naturalized rainbow trout (Oncorhynchus mykiss)

The present study combined uniparental mtDNA and biparental SNPs to illuminate the invasion and colonization pathways of rainbow trout, Oncorhynchus mykiss, one of the world's most widespread invasive fishes, that has been intensively propagated in Chile, South America. The specific aims of the study were (i) to evaluate potential donor populations, which could be either from the species' native range in North America or from introduced populations in Europe, by comparing mtDNA D-loop/control region haplotypes; and (ii) to assess the factors that have shaped genetic diversity and contemporary genetic structure of rainbow trout populations introduced to Chile through SNP genotyping. The authors comprehensively sampled 24 sites in 12 basins ranging from the High Andean Plateau (Altiplano, 18° S) to northern Patagonia (41° S). Results of the mtDNA data of naturalized trout populations from rivers in the Altiplano (northern Chile) differed from those collected in central and southern Chile, suggesting an origin from North American hatcheries. Naturalized trout populations in central and southern Chile, on the contrary, shared haplotypes with specimens found in European hatcheries. The southern and central Chile populations also contained rare haplotypes, possibly indicating potential spread through aquaculture escapes. Results of the SNP analysis revealed higher allelic richness for trout sampled in sites influenced by commercial aquaculture than sites without commercial aquaculture, likely due to increased admixture between aquaculture broodstock and naturalized trout. The analysis further uncovered some complex patterns of divergent trout populations with low genetic diversity as well as increased relatedness between individuals from isolated sites, suggesting possible local populations. A comprehensive characterization of genetic diversity and structure of rainbow trout should help identify management areas that may augment socioeconomic benefits while preventing the spread and further impacts on biodiversity.

Molecular sexing of degraded DNA from elephants and mammoths: a genotyping assay relevant both to conservation biology and to paleogenetics

It is important to determine the sex of elephants from their samples-faeces from the field or seized ivory-for forensic reasons or to understand population demography and genetic structure. Molecular sexing methods developed in the last two decades have often shown limited efficiency, particularly in terms of sensitivity and specificity, due to the degradation of DNA in these samples. These limitations have also prevented their use with ancient DNA samples of elephants or mammoths. Here we propose a novel TaqMan-MGB qPCR assay to address these difficulties. We designed it specifically to allow the characterization of the genetic sex for highly degraded samples of all elephantine taxa (elephants and mammoths). In vitro experiments demonstrated a high level of sensitivity and low contamination risks. We applied this assay in two actual case studies where it consistently recovered the right genotype for specimens of known sex a priori. In the context of a modern conservation survey of African elephants, it allowed determining the sex for over 99% of fecal samples. In a paleogenetic analysis of woolly mammoths, it produced a robust hypothesis of the sex for over 65% of the specimens out of three PCR replicates. This simple, rapid, and cost-effective procedure makes it readily applicable to large sample sizes.

Discriminating Among Pacific Salmon, Rainbow Trout, and Atlantic Salmon Species Using Common Genetic Screening Methods

The five most common species of Pacific salmon, Rainbow Trout (steelhead) Oncorhynchus spp., and Atlantic Salmon Salmo salar intermingle in the North Pacific Ocean and its freshwater tributaries. Efficient morphological methods for distinguishing among these species are sometimes limited by condition of the specimen (degraded or missing morphology), life history stage, or training of the observer. Researchers have successfully applied various genetic methods to distinguish among these species when morphological analyses are not possible, but they cannot easily incorporate these methods into standard fish and wildlife population monitoring analysis workflows. Here we test five 5′–3′ exonuclease (TaqMan) assays developed from mitochondrial genes and provide novel methods that take advantage of TaqMan output to distinguish among these species. We found that combinations of as few as two of the five assays were adequate to distinguish all species. TaqMan chemistry is designed to interrogate a single nucleotide locus. We also explore the basis for the variation in the observed scatter plot distributions (variation in florescent signals) and show that this variation is due to nucleotide diversity in and near the probe site. Because the SNPs underlying the assays developed here are all physically close to one another along the mitochondrial genome, the potential exists to develop a single DNA sequence-based assay to discriminate among salmon species. This single assay can be added to a genotyping-by-sequencing panel to identify and exclude nontarget species from analyses.

Genetic signals of artificial and natural dispersal linked to colonization of South America by non-native Chinook salmon (Oncorhynchus tshawytscha)

Genetics data have provided unprecedented insights into evolutionary aspects of colonization by non-native populations. Yet, our understanding of how artificial (human-mediated) and natural dispersal pathways of non-native individuals influence genetic metrics, evolution of genetic structure, and admixture remains elusive. We capitalize on the widespread colonization of Chinook salmon Oncorhynchus tshawytscha in South America, mediated by both dispersal pathways, to address these issues using data from a panel of polymorphic SNPs. First, genetic diversity and the number of effective breeders (Nb) were higher among artificial than natural populations. Contemporary gene flow was common between adjacent artificial and natural and adjacent natural populations, but uncommon between geographically distant populations. Second, genetic structure revealed four distinct clusters throughout the Chinook salmon distributional range with varying levels of genetic connectivity. Isolation by distance resulted from weak differentiation between adjacent artificial and natural and between natural populations, with strong differentiation between distant Pacific Ocean and Atlantic Ocean populations, which experienced strong genetic drift. Third, genetic mixture analyses revealed the presence of at least six donor geographic regions from North America, some of which likely hybridized as a result of multiple introductions. Relative propagule pressure or the proportion of Chinook salmon propagules introduced from various geographic regions according to government records significantly influenced genetic mixtures for two of three artificial populations. Our findings support a model of colonization in which high-diversity artificial populations established first; some of these populations exhibited significant admixture resulting from propagule pressure. Low-diversity natural populations were likely subsequently founded from a reduced number of individuals.

Contrasting genetic metrics and patterns among naturalized rainbow trout (Oncorhynchus mykiss) in two Patagonian lakes differentially impacted by trout aquaculture

Different pathways of propagation and dispersal of non‐native species into new environments may have contrasting demographic and genetic impacts on established populations. Repeated introductions of rainbow trout (Oncorhynchus mykiss) to Chile in South America, initially through stocking and later through aquaculture escapes, provide a unique setting to contrast these two pathways. Using a panel of single nucleotide polymorphisms, we found contrasting genetic metrics and patterns among naturalized trout in Lake Llanquihue, Chile's largest producer of salmonid smolts for nearly 50 years, and Lake Todos Los Santos (TLS), a reference lake where aquaculture has been prohibited by law. Trout from Lake Llanquihue showed higher genetic diversity, weaker genetic structure, and larger estimates for the effective number of breeders (N_b) than trout from Lake TLS. Trout from Lake TLS were divergent from Lake Llanquihue and showed marked genetic structure and a significant isolation‐by‐distance pattern consistent with secondary contact between documented and undocumented stocking events in opposite shores of the lake. Multiple factors, including differences in propagule pressure, origin of donor populations, lake geomorphology, habitat quality or quantity, and life history, may help explain contrasting genetic metrics and patterns for trout between lakes. We contend that high propagule pressure from aquaculture may not only increase genetic diversity and N_b via demographic effects and admixture, but also may impact the evolution of genetic structure and increase gene flow, consistent with findings from artificially propagated salmonid populations in their native and naturalized ranges.

Discovery and characterization of single nucleotide polymorphisms in two anadromous alosine fishes of conservation concern

Freshwater habitat alteration and marine fisheries can affect anadromous fish species, and populations fluctuating in size elicit conservation concern and coordinated management. We describe the development and characterization of two sets of 96 single nucleotide polymorphism (SNP) assays for two species of anadromous alosine fishes, alewife and blueback herring (collectively known as river herring), that are native to the Atlantic coast of North America. We used data from high-throughput DNA sequencing to discover SNPs and then developed molecular genetic assays for genotyping sets of 96 individual loci in each species. The two sets of assays were validated with multiple populations that encompass both the geographic range and the known regional genetic stocks of both species. The SNP panels developed herein accurately resolved the genetic stock structure for alewife and blueback herring that was previously identified using microsatellites and assigned individuals to regional stock of origin with high accuracy. These genetic markers, which generate data that are easily shared and combined, will greatly facilitate ongoing conservation and management of river herring including genetic assignment of marine caught individuals to stock of origin.

Non-parallel divergence across freshwater and marine three-spined stickleback Gasterosteus aculeatus populations

This work investigated whether multiple freshwater populations of three-spined stickleback Gasterosteus aculeatus in different freshwater catchments in the Jutland Peninsula, Denmark, derived from the same marine populations show repeated adaptive responses. A total of 327 G. aculeatus collected at 13 sampling locations were screened for genetic variation using a combination of 70 genes putatively under selection and 26 neutral genes along with a marker linked to the ectodysplasin gene (eda), which is strongly correlated with plate armour morphs in the species. A highly significant genetic differentiation was found that was higher among different freshwater samples than between marine-freshwater samples. Tests for selection between marine and freshwater populations showed a very low degree of parallelism and no single nucleotide polymorphism was detected as outlier in all freshwater-marine pairwise comparisons, including the eda. This suggests that G. aculeatus is not necessarily the prime example of parallel local adaptation suggested in much of the literature and that important exceptions exist (i.e. the Jutland Peninsula). While marine populations in the results described here showed a high phenotype-genotype correlation at eda, a low association was found for most of the freshwater populations. The most extreme case was found in the freshwater Lake Hald where all low-plated phenotypes were either homozygotes for the allele supposed to be associated with completely plated morphs or heterozygotes, but none were homozygotes for the putative low-plated allele. Re-examination of data from seven G. aculeatus studies agrees in showing a high but partial association between phenotype-genotype at eda in G. aculeatus freshwater populations and that mismatches occur everywhere in the European regions studied (higher in some areas, i.e. Denmark). This is independent of the eda marker used.

Historical DNA documents long-distance natal homing in marine fish

The occurrence of natal homing in marine fish remains a fundamental question in fish ecology as its unequivocal demonstration requires tracking of individuals from fertilization to reproduction. Here, we provide evidence of long-distance natal homing (>1000 km) over more than 60 years in Atlantic cod (Gadus morhua), through genetic analysis of archived samples from marked and recaptured individuals. Using a high differentiation single-nucleotide polymorphism assay, we demonstrate that the vast majority of cod tagged in West Greenland and recaptured on Icelandic spawning grounds belonged to the Iceland offshore population, strongly supporting a hypothesis of homing. The high degree of natal fidelity observed provides the evolutionary settings for development of locally adapted populations in marine fish and emphasize the need to consider portfolio effects in marine fisheries management strategies.

Temporal Genetic Variance and Propagule-Driven Genetic Structure Characterize Naturalized Rainbow Trout (Oncorhynchus mykiss) from a Patagonian Lake Impacted by Trout Farming

Knowledge about the genetic underpinnings of invasions—a theme addressed by invasion genetics as a discipline—is still scarce amid well documented ecological impacts of non-native species on ecosystems of Patagonia in South America. One of the most invasive species in Patagonia’s freshwater systems and elsewhere is rainbow trout (Oncorhynchus mykiss). This species was introduced to Chile during the early twentieth century for stocking and promoting recreational fishing; during the late twentieth century was reintroduced for farming purposes and is now naturalized. We used population- and individual-based inference from single nucleotide polymorphisms (SNPs) to illuminate three objectives related to the establishment and naturalization of Rainbow Trout in Lake Llanquihue. This lake has been intensively used for trout farming during the last three decades. Our results emanate from samples collected from five inlet streams over two seasons, winter and spring. First, we found that significant intra- population (temporal) genetic variance was greater than inter-population (spatial) genetic variance, downplaying the importance of spatial divergence during the process of naturalization. Allele frequency differences between cohorts, consistent with variation in fish length between spring and winter collections, might explain temporal genetic differences. Second, individual-based Bayesian clustering suggested that genetic structure within Lake Llanquihue was largely driven by putative farm propagules found at one single stream during spring, but not in winter. This suggests that farm broodstock might migrate upstream to breed during spring at that particular stream. It is unclear whether interbreeding has occurred between “pure” naturalized and farm trout in this and other streams. Third, estimates of the annual number of breeders (N_b) were below 73 in half of the collections, suggestive of genetically small and recently founded populations that might experience substantial genetic drift. Our results reinforce the notion that naturalized trout originated recently from a small yet genetically diverse source and that farm propagules might have played a significant role in the invasion of Rainbow Trout within a single lake with intensive trout farming. Our results also argue for proficient mitigation measures that include management of escapes and strategies to minimize unintentional releases from farm facilities.

Archived DNA reveals fisheries and climate induced collapse of a major fishery

Fishing and climate change impact the demography of marine fishes, but it is generally ignored that many species are made up of genetically distinct locally adapted populations that may show idiosyncratic responses to environmental and anthropogenic pressures. Here, we track 80 years of Atlantic cod (Gadus morhua) population dynamics in West Greenland using DNA from archived otoliths in combination with fish population and niche based modeling. We document how the interacting effects of climate change and high fishing pressure lead to dramatic spatiotemporal changes in the proportions and abundance of different genetic populations, and eventually drove the cod fishery to a collapse in the early 1970s. Our results highlight the relevance of fisheries management at the level of genetic populations under future scenarios of climate change.

Genetic evidence for ecological divergence in kokanee salmon

The evolution of locally adapted phenotypes among populations that experience divergent selective pressures is a central mechanism for generating and maintaining biodiversity. Recently, the advent of high-throughput DNA sequencing technology has provided tools for investigating the genetic basis of this process in natural populations of nonmodel organisms. Kokanee, the freshwater form of sockeye salmon (Oncorhynchus nerka), occurs as two reproductive ecotypes, which differ in spawning habitat (tributaries vs. shorelines); however, outside of the spawning season the two ecotypes co-occur in many lakes and lack diagnostic morphological characteristics. We used restriction site-associated DNA (RAD) sequencing to identify 6145 SNPs and genotype kokanee from multiple spawning sites in Okanagan Lake (British Columbia, Canada). Outlier tests revealed 18 loci putatively under divergent selection between ecotypes, all of which exhibited temporally stable allele frequencies within ecotypes. Six outliers were annotated to sequences in the NCBI database, two of which matched genes associated with early development. There was no evidence for neutral genetic differentiation; however, outlier loci demonstrated significant structure with respect to ecotype and had high assignment accuracy in mixed composition simulations. The absence of neutral structure combined with a small number of highly divergent outlier loci is consistent with theoretical predictions for the early stages of ecological divergence. These outlier loci were then applied to a realistic fisheries scenario in which additional RAD sequencing was used to genotype kokanee collected by trawl in Okanagan Lake, providing preliminary evidence that this approach may be an effective tool for conservation and management.

Fast and highly specific DNA-based multiplex detection on a solid support

Highly specific and fast multiplex detection methods are essential to conduct reasonable DNA-based diagnostics and are especially important to characterise infectious diseases. More than 1000 genetic targets such as antibiotic resistance genes, virulence factors and phylogenetic markers have to be identified as fast as possible to facilitate the correct treatment of a patient. In the present work, we developed a novel ligation-based DNA probe concept that was combined with the microarray technology and used it for the detection of bacterial pathogens. The novel linear chain (LNC) probes identified all tested species correctly within 1 h based on their 16S rRNA gene in a 25-multiplex reaction. Genomic DNA was used directly as template in the ligation reaction identifying as little as 10(7) cells without any pre-amplification. The high specificity was further demonstrated characterising a single nucleotide polymorphism leading to no false positive fluorescence signals of the untargeted single nucleotide polymorphism (SNP) variants. In comparison to conventional microarray probes, the sensitivity of the novel LNC3 probes was higher by a factor of 10 or more. In summary, we present a fast, simple, highly specific and sensitive multiplex detection method adaptable for a wide range of applications.

A low-density SNP array for analyzing differential selection in freshwater and marine populations of threespine stickleback (Gasterosteus aculeatus)

Background

The threespine stickleback (Gasterosteus aculeatus) has become an important model species for studying both contemporary and parallel evolution. In particular, differential adaptation to freshwater and marine environments has led to high differentiation between freshwater and marine stickleback populations at the phenotypic trait of lateral plate morphology and the underlying candidate gene Ectodysplacin (EDA). Many studies have focused on this trait and candidate gene, although other genes involved in marine-freshwater adaptation may be equally important. In order to develop a resource for rapid and cost efficient analysis of genetic divergence between freshwater and marine sticklebacks, we generated a low-density SNP (Single Nucleotide Polymorphism) array encompassing markers of chromosome regions under putative directional selection, along with neutral markers for background.

Results

RAD (Restriction site Associated DNA) sequencing of sixty individuals representing two freshwater and one marine population led to the identification of 33,993 SNP markers. Ninety-six of these were chosen for the low-density SNP array, among which 70 represented SNPs under putatively directional selection in freshwater vs. marine environments, whereas 26 SNPs were assumed to be neutral. Annotation of these regions revealed several genes that are candidates for affecting stickleback phenotypic variation, some of which have been observed in previous studies whereas others are new.

Conclusions

We have developed a cost-efficient low-density SNP array that allows for rapid screening of polymorphisms in threespine stickleback. The array provides a valuable tool for analyzing adaptive divergence between freshwater and marine stickleback populations beyond the well-established candidate gene Ectodysplacin (EDA).

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-867) contains supplementary material, which is available to authorized users.

Temporally isolated lineages of pink salmon reveal unique signatures of selection on distinct pools of standing genetic variation

A species' genetic diversity bears the marks of evolutionary processes that have occurred throughout its history. However, robust detection of selection in wild populations is difficult and often impeded by lack of replicate tests. Here, we investigate selection in pink salmon (Oncorhynchus gorbuscha) using genome scans coupled with inference from a haploid-assisted linkage map. Pink salmon have a strict 2-year semelparous life history which has resulted in temporally isolated (allochronic) lineages that remain sympatric through sharing of spawning habitats in alternate years. The lineages differ in a range of adaptive traits, suggesting different genetic backgrounds. We used genotyping by sequencing of haploids to generate a high-density linkage map with 7035 loci and screened an existing panel of 8036 loci for signatures of selection. The linkage map enabled identification of novel genomic regions displaying signatures of parallel selection shared between lineages. Furthermore, 24 loci demonstrated divergent selection and differences in genetic diversity between lineages, suggesting that adaptation in the 2 lineages has arisen from different pools of standing genetic variation. Findings have implications for understanding asynchronous population abundances as well as predicting future ecosystem impacts from lineage-specific responses to climate change.

Detection and mapping of QTL for temperature tolerance and body size in Chinook salmon (Oncorhynchus tshawytscha) using genotyping by sequencing

Understanding how organisms interact with their environments is increasingly important for conservation efforts in many species, especially in light of highly anticipated climate changes. One method for understanding this relationship is to use genetic maps and QTL mapping to detect genomic regions linked to phenotypic traits of importance for adaptation. We used high-throughput genotyping by sequencing (GBS) to both detect and map thousands of SNPs in haploid Chinook salmon (Oncorhynchus tshawytscha). We next applied this map to detect QTL related to temperature tolerance and body size in families of diploid Chinook salmon. Using these techniques, we mapped 3534 SNPs in 34 linkage groups which is consistent with the haploid chromosome number for Chinook salmon. We successfully detected three QTL for temperature tolerance and one QTL for body size at the experiment-wide level, as well as additional QTL significant at the chromosome-wide level. The use of haploids coupled with GBS provides a robust pathway to rapidly develop genomic resources in nonmodel organisms; these QTL represent preliminary progress toward linking traits of conservation interest to regions in the Chinook salmon genome.

Assessing patterns of hybridization between North Atlantic eels using diagnostic single-nucleotide polymorphisms

The two North Atlantic eel species, the European eel (Anguilla anguilla) and the American eel (Anguilla rostrata), spawn in partial sympatry in the Sargasso Sea, providing ample opportunity to interbreed. In this study, we used a RAD (Restriction site Associated DNA) sequencing approach to identify species-specific diagnostic single-nucleotide polymorphisms (SNPs) and design a low-density array that combined with screening of a diagnostic mitochondrial DNA marker. Eels from Iceland (N=159) and from the neighboring Faroe Islands (N=29) were genotyped, along with 94 larvae (49 European and 45 American eel) collected in the Sargasso Sea. Our SNP survey showed that the majority of Icelandic eels are pure European eels but there is also an important contribution of individuals of admixed ancestry (10.7%). Although most of the hybrids were identified as F1 hybrids from European eel female × American eel male crosses, backcrosses were also detected, including a first-generation backcross (F1 hybrid × pure European eel) and three individuals identified as second-generation backcrosses originating from American eel × F1 hybrid backcrosses interbreeding with pure European eels. In comparison, no hybrids were observed in the Faroe Islands, the closest bodies of land to Iceland. It is possible that hybrids show an intermediate migratory behaviour between the two parental species that ultimately brings hybrid larvae to the shores of Iceland, situated roughly halfway between the Sargasso Sea and Europe. Only two hybrids were observed among Sargasso Sea larvae, both backcrosses, but no F1 hybrids, that points to temporal variation in the occurrence of hybridization.

Cryptic species in plant-parasitic nematodes

This paper summarises the current knowledge concerning cryptic species of plant-parasitic nematode and briefly reviews the different methods available for their detection and characterisation. Cryptic species represent an important component of biodiversity, such speciation being common among plant-parasitic nematodes and occurring in diverse groups with different life history traits, including the spiral, virus vector, root-lesion and false root-knot nematodes. Cryptic species are important for a number of reasons, including food security, quarantine, non-chemical management technologies and species conservation, and should not be ignored. The magnitude of the phenomenon is largely unknown, but the available data on plant-parasitic nematodes demonstrate that reliance on morphology alone for species delimitation seriously underestimates the total number of taxa. Future research should focus on appropriately designed case studies using combined approaches, including large-scale, whole sample analyses by next-generation sequencing or proteomics in order to be able to answer the many questions that still remain.

Infinium assay for large-scale SNP genotyping applications

Genotyping variants in the human genome has proven to be an efficient method to identify genetic associations with phenotypes. The distribution of variants within families or populations can facilitate identification of the genetic factors of disease. Illumina's panel of genotyping BeadChips allows investigators to genotype thousands or millions of single nucleotide polymorphisms (SNPs) or to analyze other genomic variants, such as copy number, across a large number of DNA samples. These SNPs can be spread throughout the genome or targeted in specific regions in order to maximize potential discovery. The Infinium assay has been optimized to yield high-quality, accurate results quickly. With proper setup, a single technician can process from a few hundred to over a thousand DNA samples per week, depending on the type of array. This assay guides users through every step, starting with genomic DNA and ending with the scanning of the array. Using propriety reagents, samples are amplified, fragmented, precipitated, resuspended, hybridized to the chip, extended by a single base, stained, and scanned on either an iScan or Hi Scan high-resolution optical imaging system. One overnight step is required to amplify the DNA. The DNA is denatured and isothermally amplified by whole-genome amplification; therefore, no PCR is required. Samples are hybridized to the arrays during a second overnight step. By the third day, the samples are ready to be scanned and analyzed. Amplified DNA may be stockpiled in large quantities, allowing bead arrays to be processed every day of the week, thereby maximizing throughput.

Secondary contact and changes in coastal habitat availability influence the nonequilibrium population structure of a salmonid (Oncorhynchus keta)

Numerous empirical studies have reported lack of migration–drift equilibrium in wild populations. Determining the causes of nonequilibrium population structure is challenging because different evolutionary processes acting at a variety of spatiotemporal scales can produce similar patterns. Studies of contemporary populations in northern latitudes suggest that nonequilibrium population structure is probably caused by recent colonization of the region after the last Pleistocene ice age ended ∼13 000 years ago. The chum salmon's (Oncorhynchus keta) range was fragmented by dramatic environmental changes during the Pleistocene. We investigated the population structure of chum salmon on the North Alaska Peninsula (NAP) and, using both empirical data and simulations, evaluated the effects of colonization timing and founder population heterogeneity on patterns of genetic differentiation. We screened 161 single nucleotide polymorphisms and found evidence of nonequilibrium population structure when the slope of the isolation-by-distance relationship was examined at incremental spatial scales. In addition, simulations suggested that this pattern closely matched models of recent colonization of the NAP by secondary contact. Our results agree with geological and archaeological data indicating that the NAP was a dynamic landscape that may have been more recently colonized than during the last deglaciation because of dramatic changes in coastal hydrology over the last several thousand years.

How stock of origin affects performance of individuals across a meta-ecosystem: an example from sockeye salmon

Connectivity among diverse habitats can buffer populations from adverse environmental conditions, influence the functioning of meta-ecosystems, and ultimately affect the reliability of ecosystem services. This stabilizing effect on populations is proposed to derive from complementarity in growth and survival conditions experienced by individuals in the different habitats that comprise meta-ecosystems. Here we use the fine scale differentiation of salmon populations between diverse lake habitats to assess how rearing habitat and stock of origin affect the body condition of juvenile sockeye salmon. We use genetic markers (single nucleotide polymorphisms) to assign individuals of unknown origin to stock group and in turn characterize ecologically relevant attributes across habitats and stocks. Our analyses show that the body condition of juvenile salmon is related to the productivity of alternative habitats across the watershed, irrespective of their stock of origin. Emigrants and residents with genetic origins in the high productivity lake were also differentiated by their body condition, poor and high respectively. These emigrants represented a substantial proportion of juvenile sockeye salmon rearing in the lower productivity lake habitat. Despite emigrants originating from the more productive lake, they did not differ in body condition from the individuals spawned in the lower productivity, recipient habitat. Genetic tools allowed us to assess the performance of different stocks groups across the diverse habitats comprising their meta-ecosystem. The ability to characterize the ecological consequences of meta-ecosystem connectivity can help develop strategies to protect and restore ecosystems and the services they provide to humans.

Genomic patterns of introgression in rainbow and westslope cutthroat trout illuminated by overlapping paired-end RAD sequencing

Rapid and inexpensive methods for genomewide single nucleotide polymorphism (SNP) discovery and genotyping are urgently needed for population management and conservation. In hybridized populations, genomic techniques that can identify and genotype thousands of species-diagnostic markers would allow precise estimates of population- and individual-level admixture as well as identification of 'super invasive' alleles, which show elevated rates of introgression above the genomewide background (likely due to natural selection). Techniques like restriction-site-associated DNA (RAD) sequencing can discover and genotype large numbers of SNPs, but they have been limited by the length of continuous sequence data they produce with Illumina short-read sequencing. We present a novel approach, overlapping paired-end RAD sequencing, to generate RAD contigs of >300-400 bp. These contigs provide sufficient flanking sequence for design of high-throughput SNP genotyping arrays and strict filtering to identify duplicate paralogous loci. We applied this approach in five populations of native westslope cutthroat trout that previously showed varying (low) levels of admixture from introduced rainbow trout (RBT). We produced 77 141 RAD contigs and used these data to filter and genotype 3180 previously identified species-diagnostic SNP loci. Our population-level and individual-level estimates of admixture were generally consistent with previous microsatellite-based estimates from the same individuals. However, we observed slightly lower admixture estimates from genomewide markers, which might result from natural selection against certain genome regions, different genomic locations for microsatellites vs. RAD-derived SNPs and/or sampling error from the small number of microsatellite loci (n = 7). We also identified candidate adaptive super invasive alleles from RBT that had excessively high admixture proportions in hybridized cutthroat trout populations.

High potential for using DNA from ancient herring bones to inform modern fisheries management and conservation

Pacific herring (Clupea pallasi) are an abundant and important component of the coastal ecosystems for the west coast of North America. Current Canadian federal herring management assumes five regional herring populations in British Columbia with a high degree of exchange between units, and few distinct local populations within them. Indigenous traditional knowledge and historic sources, however, suggest that locally adapted, distinct regional herring populations may have been more prevalent in the past. Within the last century, the combined effects of commercial fishing and other anthropogenic factors have resulted in severe declines of herring populations, with contemporary populations potentially reflecting only the remnants of a previously more abundant and genetically diverse metapopulation. Through the analysis of 85 archaeological herring bones, this study attempted to reconstruct the genetic diversity and population structure of ancient herring populations using three different marker systems (mitochondrial DNA (mtDNA), microsatellites and SNPs). A high success rate (91%) of DNA recovery was obtained from the extremely small herring bone samples (often <10 mg). The ancient herring mtDNA revealed high haplotype diversity comparable to modern populations, although population discrimination was not possible due to the limited power of the mtDNA marker. Ancient microsatellite diversity was also similar to modern samples, but the data quality was compromised by large allele drop-out and stuttering. In contrast, SNPs were found to have low error rates with no evidence for deviations from Hardy-Weinberg equilibrium, and simulations indicated high power to detect genetic differentiation if loci under selection are used. This study demonstrates that SNPs may be the most effective and feasible approach to survey genetic population structure in ancient remains, and further efforts should be made to screen for high differentiation markers.This study provides the much needed foundation for wider scale studies on temporal genetic variation in herring, with important implications for herring fisheries management, Aboriginal title rights and herring conservation.

Rank and order: evaluating the performance of SNPs for individual assignment in a non-model organism

Single nucleotide polymorphisms (SNPs) are valuable tools for ecological and evolutionary studies. In non-model species, the use of SNPs has been limited by the number of markers available. However, new technologies and decreasing technology costs have facilitated the discovery of a constantly increasing number of SNPs. With hundreds or thousands of SNPs potentially available, there is interest in comparing and developing methods for evaluating SNPs to create panels of high-throughput assays that are customized for performance, research questions, and resources. Here we use five different methods to rank 43 new SNPs and 71 previously published SNPs for sockeye salmon: F(ST), informativeness (I(n)), average contribution to principal components (LC), and the locus-ranking programs BELS and WHICHLOCI. We then tested the performance of these different ranking methods by creating 48- and 96-SNP panels of the top-ranked loci for each method and used empirical and simulated data to obtain the probability of assigning individuals to the correct population using each panel. All 96-SNP panels performed similarly and better than the 48-SNP panels except for the 96-SNP BELS panel. Among the 48-SNP panels, panels created from F(ST), I(n), and LC ranks performed better than panels formed using the top-ranked loci from the programs BELS and WHICHLOCI. The application of ranking methods to optimize panel performance will become more important as more high-throughput assays become available.

Meiotic maps of sockeye salmon derived from massively parallel DNA sequencing

Background

Meiotic maps are a key tool for comparative genomics and association mapping studies. Next-generation sequencing and genotyping by sequencing are speeding the processes of SNP discovery and the development of new genetic tools, including meiotic maps for numerous species. Currently there are limited genetic resources for sockeye salmon, Oncorhynchus nerka. We develop the first dense meiotic map for sockeye salmon using a combination of novel SNPs found in restriction site associated DNA (RAD tags) and SNPs available from existing expressed sequence tag (EST) based assays.

Results

We discovered and genotyped putative SNPs in 3,430 RAD tags. We removed paralogous sequence variants leaving 1,672 SNPs; these were combined with 53 EST-based SNP genotypes for linkage mapping. The map contained 29 male and female linkage groups, consistent with the haploid chromosome number expected for sockeye salmon. The female map contains 1,057 loci spanning 4,896 cM, and the male map contains 1,118 loci spanning 4,220 cM. Regions of conservation with rainbow trout and synteny between the RAD based rainbow trout map and the sockeye salmon map were established.

Conclusions

Using RAD sequencing and EST-based SNP assays we successfully generated the first high density linkage map for sockeye salmon.

Signatures of natural selection among lineages and habitats in Oncorhynchus mykiss

Recent advances in molecular interrogation techniques now allow unprecedented genomic inference about the role of adaptive genetic divergence in wild populations. We used high-throughput genotyping to screen a genome-wide panel of 276 single nucleotide polymorphisms (SNPs) for the economically and culturally important salmonid Oncorhynchus mykiss. Samples included 805 individuals from 11 anadromous and resident populations from the northwestern United States and British Columbia, and represented two major lineages including paired populations of each life history within single drainages of each lineage. Overall patterns of variation affirmed clear distinctions between lineages and in most instances, isolation by distance within them. Evidence for divergent selection at eight candidate loci included significant landscape correlations, particularly with temperature. High diversity of two nonsynonymous mutations within the peptide-binding region of the major histocompatibility complex (MHC) class II (DAB) gene provided signatures of balancing selection. Weak signals for potential selection between sympatric resident and anadromous populations were revealed from genome scans and allele frequency comparisons. Our results suggest an important adaptive role for immune-related functions and present a large genomic resource for future studies

Circulating micro-RNA expression profiles in early stage nonsmall cell lung cancer

Circulating micro-RNA (miR) profiles have been proposed as promising diagnostic and prognostic biomarkers for cancer, including lung cancer. We have developed methods to accurately and reproducibly measure micro-RNA levels in serum and plasma. Here, we study paired serum and plasma samples from 220 patients with early stage nonsmall cell lung cancer (NSCLC) and 220 matched controls. We use qRT-PCR to measure the circulating levels of 30 different miRs that have previously been reported to be differently expressed in lung cancer tissue. Duplicate RNA extractions were performed for 10% of all samples, and micro-RNA measurements were highly correlated among those duplicates. This demonstrates high reproducibility of our assay. The expressions of miR-146b, miR-221, let-7a, miR-155, miR-17-5p, miR-27a and miR-106a were significantly reduced in the serum of NSCLC cases, while miR-29c was significantly increased. No significant differences were observed in plasma of patients compared with controls. Overall, expression levels in serum did not correlate well with levels in plasma. In secondary analyses, reduced plasma expression of let-7b was modestly associated with worse cancer-specific mortality in all patients, and reduced serum expression of miR-223 was modestly associated with cancer-specific mortality in stage IA/B patients. MiR profiles also showed considerable differences comparing African American and European Americans. In summary, we found significant differences in miR expression when comparing cases and controls and find evidence that expression of let-7b is associated with prognosis in NSCLC.

Genetic differentiation of Alaska Chinook salmon: the missing link for migratory studies

Most information about Chinook salmon genetic diversity and life history originates from studies from the West Coast USA, western Canada and southeast Alaska; less is known about Chinook salmon from western and southcentral Alaska drainages. Populations in this large area are genetically distinct from populations to the south and represent an evolutionary legacy of unique genetic, phenotypic and life history diversity. More genetic information is necessary to advance mixed stock analysis applications for studies involving these populations. We assembled a comprehensive, open-access baseline of 45 single nucleotide polymorphisms (SNPs) from 172 populations ranging from Russia to California. We compare SNP data from representative populations throughout the range with particular emphasis on western and southcentral Alaska. We grouped populations into major lineages based upon genetic and geographic characteristics, evaluated the resolution for identifying the composition of admixtures and performed mixed stock analysis on Chinook salmon caught incidentally in the walleye pollock fishery in the Bering Sea. SNP data reveal complex genetic structure within Alaska and can be used in applications to address not only regional issues, but also migration pathways, bycatch studies on the high seas, and potential changes in the range of the species in response to climate change.

Discovery and characterization of single nucleotide polymorphisms in Chinook salmon, Oncorhynchus tshawytscha

Molecular population genetics of non-model organisms has been dominated by the use of microsatellite loci over the last two decades. The availability of extensive genomic resources for many species is contributing to a transition to the use of single nucleotide polymorphisms (SNPs) for the study of many natural populations. Here we describe the discovery of a large number of SNPs in Chinook salmon, one of the world's most important fishery species, through large-scale Sanger sequencing of expressed sequence tag (EST) regions. More than 3 Mb of sequence was collected in a survey of variation in almost 132 kb of unique genic regions, from 225 separate ESTs, in a diverse ascertainment panel of 24 salmon. This survey yielded 117 TaqMan (5' nuclease) assays, almost all from separate ESTs, which were validated in population samples from five major stocks of salmon from the three largest basins on the Pacific coast of the contiguous United States: the Sacramento, Klamath and Columbia Rivers. The proportion of these loci that was variable in each of these stocks ranged from 86.3% to 90.6% and the mean minor allele frequency ranged from 0.194 to 0.236. There was substantial differentiation between populations with these markers, with a mean F(ST) estimate of 0.107, and values for individual loci ranging from 0 to 0.592. This substantial polymorphism and population-specific differentiation indicates that these markers will be broadly useful, including for both pedigree reconstruction and genetic stock identification applications.

Multiplex preamplification PCR and microsatellite validation enables accurate single nucleotide polymorphism genotyping of historical fish scales

Incorporating historical tissues into the study of ecological, conservation and management questions can broaden the scope of population genetic research by enhancing our understanding of evolutionary processes and anthropogenic influences on natural populations. Genotyping historical and low-quality samples has been plagued by challenges associated with low amounts of template DNA and the potential for pre-existing DNA contamination among samples. We describe a two-step process designed to (i) accurately genotype large numbers of historical low-quality scale samples in a high-throughput format and (ii) screen samples for pre-existing DNA contamination. First, we describe how an efficient multiplex preamplification PCR of 45 single nucleotide polymorphisms (SNPs) can generate highly accurate genotypes with low failure and error rates in subsequent SNP genotyping reactions of individual historical scales from sockeye salmon (Oncorhynchus nerka). Second, we demonstrate how the method can be modified for the amplification of microsatellite loci to detect pre-existing DNA contamination. A total of 760 individual historical scale and 182 contemporary fin clip samples were genotyped and screened for contamination. Genotyping failure and error rates were exceedingly low and similar for both historical and contemporary samples. Pre-existing contamination in 21% of the historical samples was successfully identified by screening the amplified microsatellite loci. The advantages of automation, low failure and error rates, and ability to multiplex both the preamplification and subsequent genotyping reactions combine to make the protocol ideally suited for efficiently genotyping large numbers of potentially contaminated low-quality sources of DNA.

Single nucleotide polymorphisms unravel hierarchical divergence and signatures of selection among Alaskan sockeye salmon (Oncorhynchus nerka) populations

Background

Disentangling the roles of geography and ecology driving population divergence and distinguishing adaptive from neutral evolution at the molecular level have been common goals among evolutionary and conservation biologists. Using single nucleotide polymorphism (SNP) multilocus genotypes for 31 sockeye salmon (Oncorhynchus nerka) populations from the Kvichak River, Alaska, we assessed the relative roles of geography (discrete boundaries or continuous distance) and ecology (spawning habitat and timing) driving genetic divergence in this species at varying spatial scales within the drainage. We also evaluated two outlier detection methods to characterize candidate SNPs responding to environmental selection, emphasizing which mechanism(s) may maintain the genetic variation of outlier loci.

Results

For the entire drainage, Mantel tests suggested a greater role of geographic distance on population divergence than differences in spawn timing when each variable was correlated with pairwise genetic distances. Clustering and hierarchical analyses of molecular variance indicated that the largest genetic differentiation occurred between populations from distinct lakes or subdrainages. Within one population-rich lake, however, Mantel tests suggested a greater role of spawn timing than geographic distance on population divergence when each variable was correlated with pairwise genetic distances. Variable spawn timing among populations was linked to specific spawning habitats as revealed by principal coordinate analyses. We additionally identified two outlier SNPs located in the major histocompatibility complex (MHC) class II that appeared robust to violations of demographic assumptions from an initial pool of eight candidates for selection.

Conclusions

First, our results suggest that geography and ecology have influenced genetic divergence between Alaskan sockeye salmon populations in a hierarchical manner depending on the spatial scale. Second, we found consistent evidence for diversifying selection in two loci located in the MHC class II by means of outlier detection methods; yet, alternative scenarios for the evolution of these loci were also evaluated. Both conclusions argue that historical contingency and contemporary adaptation have likely driven differentiation between Kvichak River sockeye salmon populations, as revealed by a suite of SNPs. Our findings highlight the need for conservation of complex population structure, because it provides resilience in the face of environmental change, both natural and anthropogenic.

Transcriptome sequencing and high-resolution melt analysis advance single nucleotide polymorphism discovery in duplicated salmonids

Until recently, single nucleotide polymorphism (SNP) discovery in nonmodel organisms faced many challenges, often depending upon a targeted-gene approach and Sanger sequencing of many individuals. The advent of next-generation sequencing technologies has dramatically improved discovery, but validating and testing SNPs for use in population studies remain labour intensive. Here, we detail a SNP discovery and validation pipeline that incorporates 454 pyrosequencing, high-resolution melt analysis (HRMA) and 5' nuclease genotyping. We generated 4.59×10(8) bp of redundant sequence from transcriptomes of two individual chum salmon, a highly valued species across the Pacific Rim. Nearly 26000 putative SNPs were identified--some as heterozygotes and some as homozygous for different nucleotides in the two individuals. For validation, we selected 202 templates containing single putative SNPs and conducted HRMA on 10 individuals from each of 19 populations from across the species range. Finally, 5' nuclease genotyping validated 37 SNPs that conformed to Hardy-Weinberg equilibrium expectations. Putative SNPs expressed as heterozygotes in an ascertainment individual had more than twice the validation rate of those homozygous for different alleles in the two fish, suggesting that many of the latter may have been paralogous sequence variants. Overall, this validation rate of 37/202 suggests that we have found more than 4500 templates containing SNPs for use in this population set. We anticipate using this pipeline to significantly expand the number of SNPs available for the studies of population structure and mixture analyses as well as for the studies of adaptive genetic variation in nonmodel organisms.

High-throughput single nucleotide polymorphism genotyping using nanofluidic Dynamic Arrays

BACKGROUND

Single nucleotide polymorphisms (SNPs) have emerged as the genetic marker of choice for mapping disease loci and candidate gene association studies, because of their high density and relatively even distribution in the human genomes. There is a need for systems allowing medium multiplexing (ten to hundreds of SNPs) with high throughput, which can efficiently and cost-effectively generate genotypes for a very large sample set (thousands of individuals). Methods that are flexible, fast, accurate and cost-effective are urgently needed. This is also important for those who work on high throughput genotyping in non-model systems where off-the-shelf assays are not available and a flexible platform is needed.

RESULTS

We demonstrate the use of a nanofluidic Integrated Fluidic Circuit (IFC) - based genotyping system for medium-throughput multiplexing known as the Dynamic Array, by genotyping 994 individual human DNA samples on 47 different SNP assays, using nanoliter volumes of reagents. Call rates of greater than 99.5% and call accuracies of greater than 99.8% were achieved from our study, which demonstrates that this is a formidable genotyping platform. The experimental set up is very simple, with a time-to-result for each sample of about 3 hours.

CONCLUSION

Our results demonstrate that the Dynamic Array is an excellent genotyping system for medium-throughput multiplexing (30-300 SNPs), which is simple to use and combines rapid throughput with excellent call rates, high concordance and low cost. The exceptional call rates and call accuracy obtained may be of particular interest to those working on validation and replication of genome- wide- association (GWA) studies.