Use of genotyping by sequencing data to develop a high-throughput and multifunctional SNP panel for conservation applications in Pacific lamprey

Practical application of the linkage disequilibrium method for estimating contemporary effective population size: A review

The method to estimate contemporary effective population size (Ne ) based on patterns of linkage disequilibrium (LD) at unlinked loci has been widely applied to natural and managed populations. The underlying model makes many simplifying assumptions, most of which have been evaluated in numerous studies published over the last two decades. Here, these performance evaluations are reviewed and summarized, with a focus on information that facilitates practical application to real populations in nature. Potential sources of bias that are discussed include calculation of r2 (a measure of LD), adjustments for sampling error, physical linkage, age structure, migration and spatial structure, mutation and selection, mating systems, changes in abundance, rare alleles, missing data, genotyping errors, data filtering choices and methods for combining multiple Ne estimates. Factors that affect precision are reviewed, including pseudoreplication that limits the information gained from large genomics datasets, constraints imposed by small samples of individuals, and the challenges in obtaining robust estimates for large populations. Topics that merit further research include the potential to weight r2 values by allele frequency, lump samples of individuals, use genotypic likelihoods rather than called genotypes, prune large LD values and apply the method to species practising partial monogamy.

Pedigree analysis and estimates of effective breeding size characterize sea lamprey reproductive biology

The sea lamprey (Petromyzon marinus) is an invasive species in the Great Lakes and the focus of a large control and assessment program. Current assessment methods provide information on the census size of spawning adult sea lamprey in a small number of streams, but information characterizing reproductive success of spawning adults is rarely available. We used RAD-capture sequencing to genotype single nucleotide polymorphism (SNP) loci for ~1600 sea lamprey larvae collected from three streams in northern Michigan (Black Mallard, Pigeon, and Ocqueoc Rivers). Larval genotypes were used to reconstruct family pedigrees, which were combined with Gaussian mixture analyses to identify larval age classes for estimation of spawning population size. Two complementary estimates of effective breeding size (N b), as well as the extrapolated minimum number of spawners (N s), were also generated for each cohort. Reconstructed pedigrees highlighted inaccuracies of cohort assignments from traditionally used mixture analyses. However, combining genotype-based pedigree information with length-at-age assignment of cohort membership greatly improved cohort identification accuracy. Population estimates across all three streams sampled in this study indicate a small number of successfully spawning adults when barriers were in operation, implying that barriers limited adult spawning numbers but were not completely effective at blocking access to spawning habitats. Thus, the large numbers of larvae present in sampled systems were a poor indicator of spawning adult abundance. Overall, pedigree-based N b and N s estimates provide a promising and rapid assessment tool for sea lamprey and other species.

Whole genome resequencing data enables a targeted SNP panel for conservation and aquaculture of Oreochromis cichlid fishes

Cichlid fish of the genus Oreochromis form the basis of the global tilapia aquaculture and fisheries industries. Broodstocks for aquaculture are often collected from wild populations, which in Africa may be from locations containing multiple Oreochromis species. However, many species are difficult to distinguish morphologically, hampering efforts to maintain good quality farmed strains. Additionally, non-native farmed tilapia populations are known to be widely distributed across Africa and to hybridize with native Oreochromis species, which themselves are important for capture fisheries. The morphological identification of these hybrids is particularly unreliable. Here, we describe the development of a single nucleotide polymorphism (SNP) genotyping panel from whole-genome resequencing data that enables targeted species identification in Tanzania. We demonstrate that an optimized panel of 96 genome-wide SNPs based on F_ST outliers performs comparably to whole genome resequencing in distinguishing species and identifying hybrids. We also show this panel outperforms microsatellite-based and phenotype-based classification methods. Case studies indicate several locations where introduced aquaculture species have become established in the wild, threatening native Oreochromis species. The novel SNP markers identified here represent an important resource for assessing broodstock purity in hatcheries and helping to conserve unique endemic biodiversity.

Development and Validation of a Multi-Locus PCR-HRM Method for Species Identification in Mytilus Genus with Food Authenticity Purposes

DNA-based methods using informative markers such as single nucleotide polymorphism (SNPs) are suitable for reliable species identification (SI) needed to enforce compliance with seafood labelling regulations (EU No.1379/2013). We developed a panel of 10 highly informative SNPs to be genotyped by PCR-High resolution melting (HRM) for SI in the Mytilus genus through in silico and in vitro stages. Its fitness for purpose and concordance were assessed by an internal validation process and by the transference to a second laboratory. The method was applicable to identify M. chilensis, M. edulis, M. galloprovincialis and M. trossulus mussels, fresh, frozen and canned with brine, oil and scallop sauce, but not in preserves containing acetic acid (wine vinegar) and tomato sauce. False-positive and negative rates were zero. Sensitivity, expressed as limit of detection (LOD), ranged between 5 and 8 ng/μL. The method was robust against small variations in DNA quality, annealing time and temperature, primer concentration, reaction volume and HRM kit. Reference materials and 220 samples were tested in an inter-laboratory assay obtaining an “almost perfect agreement” (κ = 0.925, p < 0.001). In conclusion, the method was suitable for the intended use and to be applied in the seafood industry.

Unlocking the potential of a validated single nucleotide polymorphism array for genomic monitoring of trade in cheetahs (Acinonyx jubatus)

Cheetahs (Acinonyx jubatus) are listed as vulnerable on the International Union for Conservation of Nature Red List of Threatened Species. Threats include loss of habitat, human-wildlife conflict and illegal wildlife trade. In South Africa, the export of wild cheetah is a restricted activity under the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES), however, limited legal trade is permitted of animals born to captive parents. To effectively monitor the legal and illegal trade in South Africa, it was thus essential to develop a validated molecular test. Here, we designed a single nucleotide polymorphism (SNP) array for cheetah from Double Digest Restriction Associated DNA sequencing data for individual identification and parentage testing. In order to validate the array, unrelated individuals and 16 family groups consisting of both parents and one to three offspring were genotyped using the Applied Biosystems™ QuantStudio™ 12K Flex Real-Time PCR System. In addition, parentage assignments were compared to microsatellite data. Cross-species amplification was tested in various felids and cheetah sub-species in order to determine the utility of the SNP array in other species. We obtained successful genotyping results for 218 SNPs in cheetah (A. j. jubatus) with an optimal DNA input concentration ranging from 10 to 30 ng/µl. The combination of SNPs had a higher resolving power for individual identification compared to microsatellites and provided high assignment accuracy in known pedigrees. Cross-species amplification in other felids was determined to be limited. However, the SNP array demonstrated a clear genetic discrimination of two cheetah subspecies tested here. We conclude that the described SNP array is suitable for accurate parentage assignment and provides an important traceability tool for forensic investigations of cheetah trade.

Population genomics for wildlife conservation and management

Biodiversity is under threat worldwide. Over the past decade, the field of population genomics has developed across nonmodel organisms, and the results of this research have begun to be applied in conservation and management of wildlife species. Genomics tools can provide precise estimates of basic features of wildlife populations, such as effective population size, inbreeding, demographic history and population structure, that are critical for conservation efforts. Moreover, population genomics studies can identify particular genetic loci and variants responsible for inbreeding depression or adaptation to changing environments, allowing for conservation efforts to estimate the capacity of populations to evolve and adapt in response to environmental change and to manage for adaptive variation. While connections from basic research to applied wildlife conservation have been slow to develop, these connections are increasingly strengthening. Here we review the primary areas in which population genomics approaches can be applied to wildlife conservation and management, highlight examples of how they have been used, and provide recommendations for building on the progress that has been made in this field.

Genomic islands of divergence infer a phenotypic landscape in Pacific lamprey

High rates of dispersal can breakdown coadapted gene complexes. However, concentrated genomic architecture (i.e., genomic islands of divergence) can suppress recombination to allow evolution of local adaptations despite high gene flow. Pacific lamprey (Entosphenus tridentatus) is a highly dispersive anadromous fish. Observed trait diversity and evidence for genetic basis of traits suggests it may be locally adapted. We addressed whether concentrated genomic architecture could influence local adaptation for Pacific lamprey. Using two new whole genome assemblies and genotypes from 7,716 single nucleotide polymorphism (SNP) loci in 518 individuals from across the species range, we identified four genomic islands of divergence (on chromosomes 01, 02, 04, and 22). We determined robust phenotype-by-genotype relationships by testing multiple traits across geographic sites. These trait associations probably explain genomic divergence across the species' range. We genotyped a subset of 302 broadly distributed SNPs in 2,145 individuals for association testing for adult body size, sexual maturity, migration distance and timing, adult swimming ability, and larval growth. Body size traits were strongly associated with SNPs on chromosomes 02 and 04. Moderate associations also implicated SNPs on chromosome 01 as being associated with variation in female maturity. Finally, we used candidate SNPs to extrapolate a heterogeneous spatiotemporal distribution of these predicted phenotypes based on independent data sets of larval and adult collections. These maturity and body size results guide future elucidation of factors driving regional optimization of these traits for fitness. Pacific lamprey is culturally important and imperiled. This research addresses biological uncertainties that challenge restoration efforts.

Predicting larval metamorphosis of Pacific lamprey Entosphenus tridentatus through measurements of length and mass

The utility of length and mass measurements to predict the larval metamorphosis of Pacific lamprey Entosphenus tridentatus was evaluated. During 2004-2008, larval E. tridentatus were collected from Cedar Creek (Washington, USA) in either the spring or autumn, measured for total length and total mass, reared in captivity and monitored for metamorphosis. The minimum total length, total mass and condition factor of larvae that were observed to go through metamorphosis were 102 mm, 2.0 g and 1.52, respectively. Logistic models indicated that total length and condition factor in both spring and autumn were the most significant variables for predicting metamorphosis of Pacific lamprey during the subsequent summer. Mass in the autumn also appeared important to predict whether metamorphosis occurred in the subsequent summer. Collectively, all models using specific minimums of total length, total mass or condition factor of larvae as criteria for them to metamorphose were sometimes (5 of 14 cases) able to predict the percentage of larvae that would metamorphose but rarely (1 of 12 cases) able to predict which individual larvae would metamorphose. Similar to other anadromous species of lampreys, the size and condition of larval E. tridentatus have utility for predicting metamorphic fate.

One panel to rule them all: DArTcap genotyping for population structure, historical demography, and kinship analyses, and its application to a threatened shark

With recent advances in sequencing technology, genomic data are changing how important conservation management decisions are made. Applications such as Close-Kin Mark-Recapture demand large amounts of data to estimate population size and structure, and their full potential can only be realised through ongoing improvements in genotyping strategies. Here we introduce DArTcap, a cost-efficient method that combines DArTseq and sequence capture, and illustrate its use in a high resolution population analysis of Glyphis garricki, a rare, poorly known and threatened euryhaline shark. Clustering analyses and spatial distribution of kin pairs from four different regions across northern Australia and one in Papua New Guinea, representing its entire known range, revealed that each region hosts at least one distinct population. Further structuring is likely within Van Diemen Gulf, the region that included the most rivers sampled, suggesting additional population structuring would be found if other rivers were sampled. Coalescent analyses and spatially explicit modelling suggest that G. garricki experienced a recent range expansion during the opening of the Gulf of Carpentaria following the conclusion of the Last Glacial Maximum. The low migration rates between neighbouring populations of a species that is found only in restricted coastal and riverine habitats show the importance of managing each population separately, including careful monitoring of local and remote anthropogenic activities that may affect their environments. Overall we demonstrated how a carefully chosen SNP panel combined with DArTcap can provide highly accurate kinship inference and also support population structure and historical demography analyses, therefore maximising cost-effectiveness.

Helping decision making for reliable and cost-effective 2b-RAD sequencing and genotyping analyses in non-model species

High-throughput sequencing has revolutionized population and conservation genetics. RAD sequencing methods, such as 2b-RAD, can be used on species lacking a reference genome. However, transferring protocols across taxa can potentially lead to poor results. We tested two different IIB enzymes (AlfI and CspCI) on two species with different genome sizes (the loggerhead turtle Caretta caretta and the sharpsnout seabream Diplodus puntazzo) to build a set of guidelines to improve 2b-RAD protocols on non-model organisms while optimising costs. Good results were obtained even with degraded samples, showing the value of 2b-RAD in studies with poor DNA quality. However, library quality was found to be a critical parameter on the number of reads and loci obtained for genotyping. Resampling analyses with different number of reads per individual showed a trade-off between number of loci and number of reads per sample. The resulting accumulation curves can be used as a tool to calculate the number of sequences per individual needed to reach a mean depth ≥20 reads to acquire good genotyping results. Finally, we demonstrated that selective-base ligation does not affect genomic differentiation between individuals, indicating that this technique can be used in species with large genome sizes to adjust the number of loci to the study scope, to reduce sequencing costs and to maintain suitable sequencing depth for a reliable genotyping without compromising the results. Here, we provide a set of guidelines to improve 2b-RAD protocols on non-model organisms with different genome sizes, helping decision-making for a reliable and cost-effective genotyping.

RAPTURE (RAD capture) panel facilitates analyses characterizing sea lamprey reproductive ecology and movement dynamics

Genomic tools are lacking for invasive and native populations of sea lamprey (Petromyzon marinus). Our objective was to discover single nucleotide polymorphism (SNP) loci to conduct pedigree analyses to quantify reproductive contributions of adult sea lampreys and dispersion of sibling larval sea lampreys of different ages in Great Lakes tributaries. Additional applications of data were explored using additional geographically expansive samples. We used restriction site-associated DNA sequencing (RAD-Seq) to discover genetic variation in Duffins Creek (DC), Ontario, Canada, and the St. Clair River (SCR), Michigan, USA. We subsequently developed RAD capture baits to genotype 3,446 RAD loci that contained 11,970 SNPs. Based on RAD capture assays, estimates of variance in SNP allele frequency among five Great Lakes tributary populations (mean F ST 0.008; range 0.00-0.018) were concordant with previous microsatellite-based studies; however, outlier loci were identified that contributed substantially to spatial population genetic structure. At finer scales within streams, simulations indicated that accuracy in genetic pedigree reconstruction was high when 200 or 500 independent loci were used, even in situations of high spawner abundance (e.g., 1,000 adults). Based on empirical collections of larval sea lamprey genotypes, we found that age-1 and age-2 families of full and half-siblings were widely but nonrandomly distributed within stream reaches sampled. Using the genomic scale set of SNP loci developed in this study, biologists can rapidly genotype sea lamprey in non-native and native ranges to investigate questions pertaining to population structuring and reproductive ecology at previously unattainable scales.

Genotyping-in-Thousands by sequencing (GT-seq) panel development and application to minimally invasive DNA samples to support studies in molecular ecology

Minimally invasive sampling (MIS) is widespread in wildlife studies; however, its utility for massively parallel DNA sequencing (MPS) is limited. Poor sample quality and contamination by exogenous DNA can make MIS challenging to use with modern genotyping-by-sequencing approaches, which have been traditionally developed for high-quality DNA sources. Given that MIS is often more appropriate in many contexts, there is a need to make such samples practical for harnessing MPS. Here, we test the ability for Genotyping-in-Thousands by sequencing (GT-seq), a multiplex amplicon sequencing approach, to effectively genotype minimally invasive cloacal DNA samples collected from the Western Rattlesnake (Crotalus oreganus), a threatened species in British Columbia, Canada. As there was no previous genetic information for this species, an optimized panel of 362 SNPs was selected for use with GT-seq from a de novo restriction site-associated DNA sequencing (RADseq) assembly. Comparisons of genotypes generated within and among RADseq and GT-seq for the same individuals found low rates of genotyping error (GT-seq: 0.50%; RADseq: 0.80%) and discordance (2.57%), the latter likely due to the different genotype calling models employed. GT-seq mean genotype discordance between blood and cloacal swab samples collected from the same individuals was also minimal (1.37%). Estimates of population diversity parameters were similar across GT-seq and RADseq data sets, as were inferred patterns of population structure. Overall, GT-seq can be effectively applied to low-quality DNA samples, minimizing the inefficiencies presented by exogenous DNA typically found in minimally invasive samples and continuing the expansion of molecular ecology and conservation genetics in the genomics era.

Long-term evaluation of fitness and demographic effects of a Chinook Salmon supplementation program

While the goal of supplementation programs is to provide positive, population-level effects for species of conservation concern, these programs can also present an inherent fitness risk when captive-born individuals are fully integrated into the natural population. In order to evaluate the long-term effects of a supplementation program and estimate the demographic and phenotypic factors influencing the fitness of a threatened population of Chinook Salmon (Oncorhynchus tshawytscha), we genotyped tissue samples spanning a 19-year period (1998-2016) to generate pedigrees from adult fish returning to Johnson Creek, Idaho, USA. We expanded upon previous estimates of relative reproductive success (RRS) to include grandparentage analyses and used generalized linear models to determine whether origin (hatchery or natural) or phenotypic traits (timing of arrival to spawning grounds, body length, and age) significantly predicted reproductive success (RS) across multiple years. Our results provide evidence that this supplementation program with 100% natural-origin broodstock provided a long-term demographic boost to the population (mean of 4.56 times in the first generation and mean of 2.52 times in the second generation). Overall, when spawning in nature, hatchery-origin fish demonstrated a trend toward lower RS compared to natural-origin fish (p < 0.05). However, when hatchery-origin fish successfully spawned with natural-origin fish, they had similar RS compared to natural by natural crosses (first-generation mean hatchery by natural cross RRS = 1.11 females, 1.13 males; second-generation mean hatchery by natural cross RRS = 1.03 females, 1.08 males). While origin, return year, and body length were significant predictors of fitness for both males and females (p < 0.05), return day was significant for males but not females (p > 0.05). These results indicate that supplementation programs that reduce the potential for genetic adaptation to captivity can be effective at increasing population abundance while limiting long-term fitness effects on wild populations.

Temporal dynamics of genetic clines of invasive European green crab (Carcinus maenas) in eastern North America

Two genetically distinct lineages of European green crabs (Carcinus maenas) were independently introduced to eastern North America, the first in the early 19th century and the second in the late 20th century. These lineages first came into secondary contact in southeastern Nova Scotia, Canada (NS), where they hybridized, producing latitudinal genetic clines. Previous studies have documented a persistent southward shift in the clines of different marker types, consistent with existing dispersal and recruitment pathways. We evaluated current clinal structure by quantifying the distribution of lineages and fine-scale hybridization patterns across the eastern North American range (25 locations, ~39 to 49°N) using informative single nucleotide polymorphisms (SNPs; n = 96). In addition, temporal changes in the genetic clines were evaluated using mitochondrial DNA and microsatellite loci (n = 9-11) over a 15-year period (2000-2015). Clinal structure was consistent with prior work demonstrating the existence of both northern and southern lineages with a hybrid zone occurring between southern New Brunswick (NB) and southern NS. Extensive later generation hybrids were detected in this region and in southeastern Newfoundland. Temporal genetic analysis confirmed the southward progression of clines over time; however, the rate of this progression was slower than predicted by forecasting models, and current clines for all marker types deviated significantly from these predictions. Our results suggest that neutral and selective processes contribute to cline dynamics, and ultimately, highlight how selection, hybridization, and dispersal can collectively influence invasion success.

Guidelines for planning genomic assessment and monitoring of locally adaptive variation to inform species conservation

Identifying and monitoring locally adaptive genetic variation can have direct utility for conserving species at risk, especially when management may include actions such as translocations for restoration, genetic rescue, or assisted gene flow. However, genomic studies of local adaptation require careful planning to be successful, and in some cases may not be a worthwhile use of resources. Here, we offer an adaptive management framework to help conservation biologists and managers decide when genomics is likely to be effective in detecting local adaptation, and how to plan assessment and monitoring of adaptive variation to address conservation objectives. Studies of adaptive variation using genomic tools will inform conservation actions in many cases, including applications such as assisted gene flow and identifying conservation units. In others, assessing genetic diversity, inbreeding, and demographics using selectively neutral genetic markers may be most useful. And in some cases, local adaptation may be assessed more efficiently using alternative approaches such as common garden experiments. Here, we identify key considerations of genomics studies of locally adaptive variation, provide a road map for successful collaborations with genomics experts including key issues for study design and data analysis, and offer guidelines for interpreting and using results from genomic assessments to inform monitoring programs and conservation actions.

The conservation physiology toolbox: status and opportunities

For over a century, physiological tools and techniques have been allowing researchers to characterize how organisms respond to changes in their natural environment and how they interact with human activities or infrastructure. Over time, many of these techniques have become part of the conservation physiology toolbox, which is used to monitor, predict, conserve, and restore plant and animal populations under threat. Here, we provide a summary of the tools that currently comprise the conservation physiology toolbox. By assessing patterns in articles that have been published in 'Conservation Physiology' over the past 5 years that focus on introducing, refining and validating tools, we provide an overview of where researchers are placing emphasis in terms of taxa and physiological sub-disciplines. Although there is certainly diversity across the toolbox, metrics of stress physiology (particularly glucocorticoids) and studies focusing on mammals have garnered the greatest attention, with both comprising the majority of publications (>45%). We also summarize the types of validations that are actively being completed, including those related to logistics (sample collection, storage and processing), interpretation of variation in physiological traits and relevance for conservation science. Finally, we provide recommendations for future tool refinement, with suggestions for: (i) improving our understanding of the applicability of glucocorticoid physiology; (ii) linking multiple physiological and non-physiological tools; (iii) establishing a framework for plant conservation physiology; (iv) assessing links between environmental disturbance, physiology and fitness; (v) appreciating opportunities for validations in under-represented taxa; and (vi) emphasizing tool validation as a core component of research programmes. Overall, we are confident that conservation physiology will continue to increase its applicability to more taxa, develop more non-invasive techniques, delineate where limitations exist, and identify the contexts necessary for interpretation in captivity and the wild.

Fine-scale temperature-associated genetic structure between inshore and offshore populations of sea scallop (Placopecten magellanicus)

In the northwest Atlantic Ocean, sea scallop (Placopecten magellanicus) has been characterized by a latitudinal genetic cline with a breakpoint between northern and southern genetic clusters occurring at ~45°N along eastern Nova Scotia, Canada. Using 96 diagnostic single-nucleotide polymorphisms (SNPs) capable of discriminating between northern and southern clusters, we examined fine-scale genetic structure of scallops among 27 sample locations, spanning the largest geographic range evaluated in this species to date (~37-51°N). Here, we confirmed previous observations of northern and southern groups, but we show that the boundary between northern and southern clusters is not a discrete latitudinal break. Instead, at latitudes near the previously described boundary, we found unexpected patterns of fine-scale genetic structure occurring between inshore and offshore sites. Scallops from offshore sites, including St. Pierre Bank and the eastern Scotian Shelf, clustered with southern stocks, whereas inshore sites at similar latitudes clustered with northern stocks. Our analyses revealed significant genetic divergence across small spatial scales (i.e., 129-221 km distances), and that spatial structure over large and fine scales was strongly associated with temperature during seasonal periods of thermal minima. Clear temperature differences between inshore and offshore locations may explain the fine-scale structuring observed, such as why southern lineages of scallop occur at higher latitudes in deeper, warmer offshore waters. Our study supports growing evidence that fine-scale population structure in marine species is common, often environmentally associated, and that consideration of environmental and genomic data can significantly enhance the identification of marine diversity and management units.

A SNP-Enabled Assessment of Genetic Diversity, Evolutionary Relationships and the Identification of Candidate Genes in Chrysanthemum

Varieties of the economically important ornamental species chrysanthemum have been bred to fit a number of market niches, but the genetic basis and evolutionary relationships among various cultivated types are poorly understood. Here, a DNA marker-based analysis of 199 chrysanthemum entries representing each of the five cultivated types is presented. A set of >90,000 single nucleotide polymorphisms (SNPs) associated with a minor allele frequency of at least 5% was defined, and used to perform a phylogenetic analysis which corresponded well with the phenotypic classification. The analysis revealed that the small-flowered types, spray cut chrysanthemum (SCC) and potted and ground chrysanthemum (PGC), are more closely related to the wild progenitor species (WC) than are the large-flowered ones, disbud cut chrysanthemum (DCC) and traditional chrysanthemum (TC); and the PGC type was closest. Some 550 genetic regions appeared to have experienced selection in the separation of potted and ground-cover types from disbud cut types, and that between potted and ground-cover types from traditional types. A genome-wide association analysis revealed that seven SNPs lying within six genes were predictive of three important traits (ray floret type, cultivated type and flower shape), but no association with flower color was detected. The study has provided a number of novel insights into evolutionary relationships, the population structure and the genetic basis of some key ornamental traits.

Growth and annual survival estimates to examine the ecology of larval lamprey and the implications of ageing error in fitting models

This study used existing western brook lamprey Lampetra richardsoni age information to fit three different growth models (i.e. von Bertalanffy, Gompertz and logistic) with and without error in age estimates. Among these growth models, there was greater support for the logistic and Gompertz models than the von Bertalanffy model, regardless of ageing error assumptions. The von Bertalanffy model, however, appeared to fit the data well enough to permit survival estimates; using length-based estimators, annual survival varied between 0·64 (95% credibility interval: 0·44-0·79) and 0·81 (0·79-0·83) depending on ageing and growth process error structure. These estimates are applicable to conservation and management of L. richardsoni and other western lampreys (e.g. Pacific lamprey Entosphenus tridentatus) and can potentially be used in the development of life-cycle models for these species. These results also suggest that estimators derived from von Bertalanffy growth models should be interpreted with caution if there is high uncertainty in age estimates.

Construction of a High-Density Genetic Map and Quantitative Trait Locus Mapping in the Manila clam Ruditapes philippinarum

Genetic linkage maps are indispensable tools in a wide range of genetic and genomic research. With the advancement of genotyping-by-sequencing (GBS) methods, the construction of a high-density linkage maps has become achievable in marine organisms lacking sufficient genomic resources, such as mollusks. In this study, high-density linkage map was constructed for an ecologically and commercially important clam species, Ruditapes philippinarum. For the consensus linkage map, a total of 9658 markers spanning 1926.98 cM were mapped to 18 sex-averaged linkage groups, with an average marker distance of 0.42 cM. Based on the high-density linkage map, ten QTLs for growth-related traits and shell color were detected. The coverage and density of the current map are sufficient for us to effectively detect QTL for segregating traits, and two QTL positions were all coincident with the closest markers. This high-density genetic linkage map reveals basic genomic architecture and will be useful for comparative genomics research, genome assembly and genetic improvement of R. philippinarum and other bivalve molluscan species.

Accurate recapture identification for genetic mark-recapture studies with error-tolerant likelihood-based match calling and sample clustering

Error-tolerant likelihood-based match calling presents a promising technique to accurately identify recapture events in genetic mark-recapture studies by combining probabilities of latent genotypes and probabilities of observed genotypes, which may contain genotyping errors. Combined with clustering algorithms to group samples into sets of recaptures based upon pairwise match calls, these tools can be used to reconstruct accurate capture histories for mark-recapture modelling. Here, we assess the performance of a recently introduced error-tolerant likelihood-based match-calling model and sample clustering algorithm for genetic mark-recapture studies. We assessed both biallelic (i.e. single nucleotide polymorphisms; SNP) and multiallelic (i.e. microsatellite; MSAT) markers using a combination of simulation analyses and case study data on Pacific walrus (Odobenus rosmarus divergens) and fishers (Pekania pennanti). A novel two-stage clustering approach is demonstrated for genetic mark-recapture applications. First, repeat captures within a sampling occasion are identified. Subsequently, recaptures across sampling occasions are identified. The likelihood-based matching protocol performed well in simulation trials, demonstrating utility for use in a wide range of genetic mark-recapture studies. Moderately sized SNP (64+) and MSAT (10-15) panels produced accurate match calls for recaptures and accurate non-match calls for samples from closely related individuals in the face of low to moderate genotyping error. Furthermore, matching performance remained stable or increased as the number of genetic markers increased, genotyping error notwithstanding.

SNP discovery and genotyping using Genotyping-by-Sequencing in Pekin ducks

Genomic selection and genome-wide association studies need thousands to millions of SNPs. However, many non-model species do not have reference chips for detecting variation. Our goal was to develop and validate an inexpensive but effective method for detecting SNP variation. Genotyping by sequencing (GBS) can be a highly efficient strategy for genome-wide SNP detection, as an alternative to microarray chips. Here, we developed a GBS protocol for ducks and tested it to genotype 49 Pekin ducks. A total of 169,209 SNPs were identified from all animals, with a mean of 55,920 SNPs per individual. The average SNP density reached 1156 SNPs/MB. In this study, the first application of GBS to ducks, we demonstrate the power and simplicity of this method. GBS can be used for genetic studies in to provide an effective method for genome-wide SNP discovery.

Sequencing improves our ability to study threatened migratory species: Genetic population assignment in California's Central Valley Chinook salmon

Effective conservation and management of migratory species requires accurate identification of unique populations, even as they mix along their migratory corridors. While telemetry has historically been used to study migratory animal movement and habitat use patterns, genomic tools are emerging as a superior alternative in many ways, allowing large‐scale application at reduced costs. Here, we demonstrate the usefulness of genomic resources for identifying single‐nucleotide polymorphisms (SNPs) that allow fast and accurate identification of the imperiled Chinook salmon in the Great Central Valley of California. We show that 80 well‐chosen loci, drawn from a pool of over 11,500 SNPs developed from restriction site‐associated DNA sequencing, can accurately identify Chinook salmon runs and select populations within run. No other SNP panel for Central Valley Chinook salmon has been able to achieve the high accuracy of assignment we show here. This panel will greatly improve our ability to study and manage this ecologically, economically, and socially important species and demonstrates the great utility of using genomics to study migratory species.

The history of Old World camelids in the light of molecular genetics

Old World camels have come into the focus as sustainable livestock species, unique in their morphological and physiological characteristics and capable of providing vital products even under extreme environmental conditions. The evolutionary history of dromedary and Bactrian camels traces back to the middle Eocene (around 40 million years ago, mya), when the ancestors of Camelus emerged on the North American continent. While the genetic status of the two domestic species has long been established, the wild two-humped camel has only recently been recognized as a separate species, Camelus ferus, based on molecular genetic data. The demographic history established from genome drafts of Old World camels shows the independent development of the three species over the last 100,000 years with severe bottlenecks occurring during the last glacial period and in the recent past. Ongoing studies involve the immune system, relevant production traits, and the global population structure and domestication of Old World camels. Based on the now available whole genome drafts, specific metabolic pathways have been described shedding new light on the camels' ability to adapt to desert environments. These new data will also be at the origin for genome-wide association studies to link economically relevant phenotypes to genotypes and to conserve the diverse genetic resources in Old World camelids.

Development of a universal double-digest RAD sequencing approach for a group of nonmodel, ecologically and economically important insect and fish taxa

The generation of genome-scale data is critical for a wide range of questions in basic biology using model organisms, but also in questions of applied biology in nonmodel organisms (agriculture, natural resources, conservation and public health biology). Using a genome-scale approach on a diverse group of nonmodel organisms and with the goal of lowering costs of the method, we modified a multiplexed, high-throughput genomic scan technique utilizing two restriction enzymes. We analysed several pairs of restriction enzymes and completed double-digestion RAD sequencing libraries for nine different species and five genera of insects and fish. We found one particular enzyme pair produced consistently higher number of sequence-able fragments across all nine species. Building libraries off this enzyme pair, we found a range of usable SNPs between 4000 and 37 000 SNPS per species and we found a greater number of usable SNPs using reference genomes than de novo pipelines in STACKS. We also found fewer reads in the Read 2 fragments from the paired-end Illumina Hiseq run. Overall, the results of this study provide empirical evidence of the utility of this method for producing consistent data for diverse nonmodel species and suggest specific considerations for sequencing analysis strategies.

Draft genome sequence of the silver pomfret fish, Pampus argenteus

Silver pomfret, Pampus argenteus, is a fish species from coastal waters. Despite its high commercial value, this edible fish has not been sequenced. Hence, its genetic and genomic studies have been limited. We report the first draft genome sequence of the silver pomfret obtained using a Next Generation Sequencing (NGS) technology. We assembled 38.7 Gb of nucleotides into scaffolds of 350 Mb with N50 of about 1.5 kb, using high quality paired end reads. These scaffolds represent 63.7% of the estimated silver pomfret genome length. The newly sequenced and assembled genome has 11.06% repetitive DNA regions, and this percentage is comparable to that of the tilapia genome. The genome analysis predicted 16 322 genes. About 91% of these genes showed homology with known proteins. Many gene clusters were annotated to protein and fatty-acid metabolism pathways that may be important in the context of the meat texture and immune system developmental processes. The reference genome can pave the way for the identification of many other genomic features that could improve breeding and population-management strategies, and it can also help characterize the genetic diversity of P. argenteus.

A high-density linkage map for Astyanax mexicanus using genotyping-by-sequencing technology

The Mexican tetra, Astyanax mexicanus, is a unique model system consisting of cave-adapted and surface-dwelling morphotypes that diverged >1 million years (My) ago. This remarkable natural experiment has enabled powerful genetic analyses of cave adaptation. Here, we describe the application of next-generation sequencing technology to the creation of a high-density linkage map. Our map comprises more than 2200 markers populating 25 linkage groups constructed from genotypic data generated from a single genotyping-by-sequencing project. We leveraged emergent genomic and transcriptomic resources to anchor hundreds of anonymous Astyanax markers to the genome of the zebrafish (Danio rerio), the most closely related model organism to our study species. This facilitated the identification of 784 distinct connections between our linkage map and the Danio rerio genome, highlighting several regions of conserved genomic architecture between the two species despite ~150 My of divergence. Using a Mendelian cave-associated trait as a proof-of-principle, we successfully recovered the genomic position of the albinism locus near the gene Oca2. Further, our map successfully informed the positions of unplaced Astyanax genomic scaffolds within particular linkage groups. This ability to identify the relative location, orientation, and linear order of unaligned genomic scaffolds will facilitate ongoing efforts to improve on the current early draft and assemble future versions of the Astyanax physical genome. Moreover, this improved linkage map will enable higher-resolution genetic analyses and catalyze the discovery of the genetic basis for cave-associated phenotypes.

Genes predict long distance migration and large body size in a migratory fish, Pacific lamprey

Elucidation of genetic mechanisms underpinning migratory behavior could help predict how changes in genetic diversity may affect future spatiotemporal distribution of a migratory species. This ability would benefit conservation of one such declining species, anadromous Pacific lamprey (Entosphenus tridentatus). Nonphilopatric migration of adult Pacific lamprey has homogenized population-level neutral variation but has maintained adaptive variation that differentiates groups based on geography, run-timing and adult body form. To investigate causes for this adaptive divergence, we examined 647 adult lamprey sampled at a fixed location on the Columbia River and radiotracked during their subsequent upstream migration. We tested whether genetic variation [94 neutral and adaptive single nucleotide polymorphisms (SNPs) previously identified from a genomewide association study] was associated with phenotypes of migration distance, migration timing, or morphology. Three adaptive markers were strongly associated with morphology, and one marker also correlated with upstream migration distance and timing. Genes physically linked with these markers plausibly influence differences in body size, which is also consistently associated with migration distance in Pacific lamprey. Pacific lamprey conservation implications include the potential to predict an individual's upstream destination based on its genotype. More broadly, the results suggest a genetic basis for intrapopulation variation in migration distance in migratory species.

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.