Identification of single nucleotide polymorphisms in candidate genes for growth and reproduction in a nonmodel organism; the Atlantic cod,Gadus morhua

Investigating the role of genetic variation in vgll3 and six6 in the domestication of gilthead seabream (Sparus aurata Linnaeus) and European seabass (Dicentrarchus labrax Linnaeus)

Gene function conservation is crucial in molecular ecology, especially for key traits like growth and maturation in teleost fish. The vgll3 and six6 genes are known to influence age-at-maturity in Atlantic salmon, but their impact on other fish species is poorly understood. Here, we investigated the association of vgll3 and six6 in the domestication of gilthead seabream and European seabass, both undergoing selective breeding for growth-related traits in the Mediterranean. We analysed two different sets of samples using two different genotyping approaches. The first dataset comprised farmed and wild populations from Greece, genotyped for SNPs within the two genes ('gene-level genotyping'). The second dataset examined 300-600 k SNPs located in the chromosomes of the two genes, derived from a meta-analysis of a Pool-Seq experiment involving farmed and wild populations distributed widely across the Mediterranean ('chromosome-level genotyping'). The gene-level analysis revealed a statistically significant allele frequency differences between farmed and wild populations on both genes in each species. This finding was partially supported by the chromosome-level analysis, identifying highly differentiated regions may be involved in the domestication process at varying distances from the candidate genes. Noteworthy genomic features were found, such as a CpG island in gilthead seabream and novel candidate genes in European seabass, warranting further investigation. These findings support a putative role of vgll3 and six6 in the maturation and growth of gilthead seabream and European seabass, emphasizing the need for further research on their conserved function.

Genetic assignment predicts depth of benthic settlement for 0-group Atlantic cod

Atlantic cod is a keystone species that remains among the most economically important demersal fish in the North Atlantic. Throughout its distribution range, Atlantic cod is composed of populations with varying environmental preferences and migratory propensities. This life-history variation is likely to have contributed to the niche width and large population sizes of Atlantic cod, and its relative resilience to environmental change and exploitation. The Icelandic cod stock is currently managed as a single unit, but early research indicates population variation by depth and temperature and distinct offshore and inshore spawning components. Pelagic 0-group juveniles from different spawning grounds coexist in nursery areas around Iceland, but their genetic composition or habitat partitioning had not been examined post benthic settlement. In the current study we examine the genetic composition of Atlantic cod juvenile aggregations at nearshore nursery grounds in NW-Iceland and report distinct segregation by the depth of offshore and inshore juvenile cod. The physiological mechanism of this segregation is not known, but the pattern demonstrates the need to consider population structure at nursery grounds in the application of marine spatial planning and other area-based conservation tools.

Detecting population declines via monitoring the effective number of breeders (Nb )

Estimating the effective population size and effective number of breeders per year (N_b ) can facilitate early detection of population declines. We used computer simulations to quantify bias and precision of the one-sample LDNe estimator of N_b in age-structured populations using a range of published species life history types, sample sizes, and DNA markers. N_b estimates were biased by ~5%-10% when using SNPs or microsatellites in species ranging from fishes to mosquitoes, frogs, and seaweed. The bias (high or low) was similar for different life history types within a species suggesting that life history variation in populations will not influence N_b estimation. Precision was higher for 100 SNPs (H ≈ 0.30) than for 15 microsatellites (H ≈ 0.70). Confidence intervals (CIs) were occasionally too narrow, and biased high when N_b was small (N_b  < 50); however, the magnitude of bias would unlikely influence management decisions. The CIs (from LDNe) were sufficiently narrow to achieve high statistical power (≥0.80) to reject the null hypothesis that N_b  = 50 when the true N_b  = 30 and when sampling 50 individuals and 200 SNPs. Similarly, CIs were sufficiently narrow to reject N_b  = 500 when the true N_b  = 400 and when sampling 200 individuals and 5,000 loci. Finally, we present a linear regression method that provides high power to detect a decline in N_b when sampling at least five consecutive cohorts. This study provides guidelines and tools to simulate and estimate N_b for age structured populations (https://github.com/popgengui/agestrucnb/), which should help biologists develop sensitive monitoring programmes for early detection of changes in N_b and population declines.

Differential gene expression and SNP association between fast- and slow-growing turbot (Scophthalmus maximus)

Growth is among the most important traits for animal breeding. Understanding the mechanisms underlying growth differences between individuals can contribute to improving growth rates through more efficient breeding schemes. Here, we report a transcriptomic study in muscle and brain of fast- and slow-growing turbot (Scophthalmus maximus), a relevant flatfish in European and Asian aquaculture. Gene expression and allelic association between the two groups were explored. Up-regulation of the anaerobic glycolytic pathway in the muscle of fast-growing fish was observed, indicating a higher metabolic rate of white muscle. Brain expression differences were smaller and not associated with major growth-related genes, but with regulation of feeding-related sensory pathways. Further, SNP variants showing frequency differences between fast- and slow-growing fish pointed to genomic regions likely involved in growth regulation, and three of them were individually validated through SNP typing. Although different mechanisms appear to explain growth differences among families, general mechanisms seem also to be involved, and thus, results provide a set of useful candidate genes and markers to be evaluated for more efficient growth breeding programs and to perform comparative genomic studies of growth in fish and vertebrates.

Integrative Transcriptome, Genome and Quantitative Trait Loci Resources Identify Single Nucleotide Polymorphisms in Candidate Genes for Growth Traits in Turbot

Growth traits represent a main goal in aquaculture breeding programs and may be related to adaptive variation in wild fisheries. Integrating quantitative trait loci (QTL) mapping and next generation sequencing can greatly help to identify variation in candidate genes, which can result in marker-assisted selection and better genetic structure information. Turbot is a commercially important flatfish in Europe and China, with available genomic information on QTLs and genome mapping. Muscle and liver RNA-seq from 18 individuals was carried out to obtain gene sequences and markers functionally related to growth, resulting in a total of 20,447 genes and 85,344 single nucleotide polymorphisms (SNPs). Many growth-related genes and SNPs were identified and placed in the turbot genome and genetic map to explore their co-localization with growth-QTL markers. Forty-five SNPs on growth-related genes were selected based on QTL co-localization and relevant function for growth traits. Forty-three SNPs were technically feasible and validated in a wild Atlantic population, where 91% were polymorphic. The integration of functional and structural genomic resources in turbot provides a practical approach for QTL mining in this species. Validated SNPs represent a useful set of growth-related gene markers for future association, functional and population studies in this flatfish species.

Molecular characterization and expression of three preprosomatostatin genes and their association with growth in common carp (Cyprinus carpio)

Somatostatins (SSs) are a structurally diverse family of peptides that play important roles in the regulation of growth, development and metabolism in vertebrates. In this study, three preprosomatostatin genes (PSSs) in the common carp, Cyprinus carpio (Cc) were identified and characterized. Based on cloned sequences and genome BLAST, six isoforms of the PSS gene in C. carpio (CcPSS) were identified and included CcPSS1a and CcPSS1b, CcPSS2a and CcPSS2b, and finally, CcPSS3a and CcPSS3b. The open reading frames (ORF) of CcPSS1a, CcPSS2a and CcPSS3a consist of 345, 336 and 363 nucleotides. During embryonic development, the expressions of CcPSS2 and CcPSS3 were first observed at the stage of optic vesicle, and CcPSS1 mRNA was initially detected at the stage of muscular effect. The highest mRNA levels of CcPSS1, CcPSS2 and CcPSS3 were observed at 1-day post-hatch (dph), 2-dph and the stage of heart beating, respectively. In the adult brain, the distributions of three CcPSS mRNAs were differential but overlapping in the hypothalamus, telencephalon and medulla oblongata. For peripheral tissues, all three CcPSS mRNAs were detected in the mid-intestine, and CcPSS1 and CcPSS3 mRNAs were also expressed in the liver. Owing to the importance of somatostatins on regulating growth, functional mutations of CcPSSs were identified in a C. carpio population. A total of 23 polymorphic sites were detected in CcPSS1a and CcPSS3a. Of them, two SNPs (CcPSS1a-g.922C>T, and CcPSS3a-g.1125C>A) were significantly associated with growth traits, indicating their potential applications in gene (marker)-assisted selective breeding in C. carpio.

Molecular characterization of heat shock protein 70 (HSP 70) promoter in Japanese flounder (Paralichthys olivaceus), and the association of Pohsp70 SNPs with heat-resistant trait

Ambient temperature is one of the major abiotic environmental factors determining the main parameters of fish vital activity. HSP70 plays an essential role in heat response. In this investigation, the promoter and structure of Paralichthys olivaceus hsp70 (Pohsp70) gene was cloned and predicted. 2558 bp upstream regulatory region of Pohsp70 was annotated with four potential promoter elements and four putative binding sites of transcription factors heat shock elements (HSE, nGAAn) in the upstream of the transcription start site. In addition, one intron with 454 bp in the 5'-noncoding region was found. Quantitative Real Time PCR analysis indicated that the transcript level of Pohsp70 was raised markedly after 1 h by heat shocked. Furthermore, 25 SNPs were identified in Pohsp70 by resequencing, seven of which was associated with heat resistance. In addition, two of the seven SNPs, namely SNP14 and SNP16, were observed in strong linkage disequilibrium. The haplotype with association analysis showed TAGGAG haplotype was more represented in heat susceptible group while (DEL/T) GAATA haplotype was more frequent in heat resistant group. The heat resistant SNPs and haplotype could be candidate markers potentially serving for selective breeding programs of Japanese flounder aimed at improving anti-stress and production.

So long to genetic diversity, and thanks for all the fish

The world faces a global fishing crisis. Wild marine fisheries comprise nearly 15% of all animal protein in the human diet, but, according to the U.N. Food and Agriculture Organization, nearly 60% of all commercially important marine fish stocks are overexploited, recovering, or depleted (FAO 2012; Fig. 1). Some authors have suggested that the large population sizes of harvested marine fish make even collapsed populations resistant to the loss of genetic variation by genetic drift (e.g. Beverton 1990). In contrast, others have argued that the loss of alleles because of overfishing may actually be more dramatic in large populations than in small ones (Ryman et al. 1995). In this issue, Pinsky & Palumbi (2014) report that overfished populations have approximately 2% lower heterozygosity and 12% lower allelic richness than populations that are not overfished. They also performed simulations which suggest that their estimates likely underestimate the actual loss of rare alleles by a factor of three or four. This important paper shows that the harvesting of marine fish can have genetic effects that threaten the long-term sustainability of this valuable resource.

Targeted approach to identify genetic loci associated with evolved dioxin tolerance in Atlantic killifish (Fundulus heteroclitus)

Background

The most toxic aromatic hydrocarbon pollutants are categorized as dioxin-like compounds (DLCs) to which extreme tolerance has evolved independently and contemporaneously in (at least) four populations of Atlantic killifish (Fundulus heteroclitus). Surprisingly, the magnitude and phenotype of DLC tolerance is similar among these killifish populations that have adapted to varied, but highly aromatic hydrocarbon-contaminated urban/industrialized estuaries of the US Atlantic coast. Multiple tolerant and neighboring sensitive killifish populations were compared with the expectation that genetic loci associated with DLC tolerance would be revealed.

Results

Since the aryl hydrocarbon receptor (AHR) pathway partly or fully mediates DLC toxicity in vertebrates, single nucleotide polymorphisms (SNPs) from 42 genes associated with the AHR pathway were identified to serve as targeted markers. Wild fish (N = 36/37) from four highly tolerant killifish populations and four nearby sensitive populations were genotyped using 59 SNP markers. Similar to other killifish population genetic analyses, strong genetic differentiation among populations was detected, consistent with isolation by distance models. When DLC-sensitive populations were pooled and compared to pooled DLC-tolerant populations, multi-locus analyses did not distinguish the two groups. However, pairwise comparisons of nearby tolerant and sensitive populations revealed high differentiation among sensitive and tolerant populations at these specific loci: AHR 1 and 2, cathepsin Z, the cytochrome P450s (CYP1A and 3A30), and the NADH dehydrogenase subunits. In addition, significant shifts in minor allele frequency were observed at AHR2 and CYP1A loci across most sensitive/tolerant pairs, but only AHR2 exhibited shifts in the same direction across all pairs.

Conclusions

The observed differences in allelic composition at the AHR2 and CYP1A SNP loci were identified as significant among paired sensitive/tolerant populations of Atlantic killifish with multiple statistical tests. The genetic patterns reported here lend support to the argument that AHR2 and CYP1A play a role in the adaptive response to extreme DLC contamination. Additional functional assays are required to isolate the exact mechanism of DLC tolerance.

Development of molecular resources for an intertidal clam, Sinonovacula constricta, using 454 transcriptome sequencing

Background

The razor clam Sinonovacula constricta is a benthic intertidal bivalve species with important commercial value. Despite its economic importance, knowledge of its transcriptome is scarce. Next generation sequencing technologies offer rapid and efficient tools for generating large numbers of sequences, which can be used to characterize the transcriptome, to develop effective molecular markers and to identify genes associated with growth, a key breeding trait.

Results

Total RNA was isolated from the mantle, gill, liver, siphon, gonad and muscular foot tissues. High-throughput deep sequencing of S. constricta using 454 pyrosequencing technology yielded 859,313 high-quality reads with an average read length of 489 bp. Clustering and assembly of these reads produced 16,323 contigs and 131,346 singletons with average lengths of 1,376 bp and 458 bp, respectively. Based on transcriptome sequencing, 14,615 sequences had significant matches with known genes encoding 147,669 predicted proteins. Subsequently, previously unknown growth-related genes were identified. A total of 13,563 microsatellites (SSRs) and 13,634 high-confidence single nucleotide polymorphism loci (SNPs) were discovered, of which almost half were validated.

Conclusion

De novo sequencing of the razor clam S. constricta transcriptome on the 454 GS FLX platform generated a large number of ESTs. Candidate growth factors and a large number of SSRs and SNPs were identified. These results will impact genetic studies of S. constricta.

Microevolution in time and space: SNP analysis of historical DNA reveals dynamic signatures of selection in Atlantic cod

Little is known about how quickly natural populations adapt to changes in their environment and how temporal and spatial variation in selection pressures interact to shape patterns of genetic diversity. We here address these issues with a series of genome scans in four overfished populations of Atlantic cod (Gadus morhua) studied over an 80-year period. Screening of >1000 gene-associated single-nucleotide polymorphisms (SNPs) identified 77 loci that showed highly elevated levels of differentiation, likely as an effect of directional selection, in either time, space or both. Exploratory analysis suggested that temporal allele frequency shifts at certain loci may correlate with local temperature variation and with life history changes suggested to be fisheries induced. Interestingly, however, largely nonoverlapping sets of loci were temporal outliers in the different populations and outliers from the 1928 to 1960 period showed almost complete stability during later decades. The contrasting microevolutionary trajectories among populations resulted in sequential shifts in spatial outliers, with no locus maintaining elevated spatial differentiation throughout the study period. Simulations of migration coupled with observations of temporally stable spatial structure at neutral loci suggest that population replacement or gene flow alone could not explain all the observed allele frequency variation. Thus, the genetic changes are likely to at least partly be driven by highly dynamic temporally and spatially varying selection. These findings have important implications for our understanding of local adaptation and evolutionary potential in high gene flow organisms and underscore the need to carefully consider all dimensions of biocomplexity for evolutionarily sustainable management.

Spatiotemporal SNP analysis reveals pronounced biocomplexity at the northern range margin of Atlantic cod Gadus morhua

Accurate prediction of species distribution shifts in the face of climate change requires a sound understanding of population diversity and local adaptations. Previous modeling has suggested that global warming will lead to increased abundance of Atlantic cod (Gadus morhua) in the ocean around Greenland, but the dynamics of earlier abundance fluctuations are not well understood. We applied a retrospective spatiotemporal population genomics approach to examine the temporal stability of cod population structure in this region and to search for signatures of divergent selection over a 78-year period spanning major demographic changes. Analyzing >900 gene-associated single nucleotide polymorphisms in 847 individuals, we identified four genetically distinct groups that exhibited varying spatial distributions with considerable overlap and mixture. The genetic composition had remained stable over decades at some spawning grounds, whereas complete population replacement was evident at others. Observations of elevated differentiation in certain genomic regions are consistent with adaptive divergence between the groups, indicating that they may respond differently to environmental variation. Significantly increased temporal changes at a subset of loci also suggest that adaptation may be ongoing. These findings illustrate the power of spatiotemporal population genomics for revealing biocomplexity in both space and time and for informing future fisheries management and conservation efforts.

Rapid, economical single-nucleotide polymorphism and microsatellite discovery based on de novo assembly of a reduced representation genome in a non-model organism: a case study of Atlantic cod Gadus morhua

By combining next-generation sequencing technology (454) and reduced representation library (RRL) construction, the rapid and economical isolation of over 25 000 potential single-nucleotide polymorphisms (SNP) and >6000 putative microsatellite loci from c. 2% of the genome of the non-model teleost, Atlantic cod Gadus morhua from the Celtic Sea, south of Ireland, was demonstrated. A small-scale validation of markers indicated that 80% (11 of 14) of SNP loci and 40% (6 of 15) of the microsatellite loci could be amplified and showed variability. The results clearly show that small-scale next-generation sequencing of RRL genomes is an economical and rapid approach for simultaneous SNP and microsatellite discovery that is applicable to any species. The low cost and relatively small investment in time allows for positive exploitation of ascertainment bias to design markers applicable to specific populations and study questions.

Rates of assay success and genotyping error when single nucleotide polymorphism genotyping in non-model organisms: a case study in the Antarctic fur seal

Although single nucleotide polymorphisms (SNPs) are increasingly being recognized as powerful molecular markers, their application to non-model organisms can bring significant challenges. Among these are imperfect conversion rates of assays designed from in silico resources and the enhanced potential for genotyping error relative to pre-validated, highly optimized human SNPs. To explore these issues, we used Illumina's GoldenGate assay to genotype 480 Antarctic fur seal (Arctocephalus gazella) individuals at 144 putative SNPs derived from a 454 transcriptome assembly. One hundred and thirty-five polymorphic SNPs (93.8%) were automatically validated by the program GenomeStudio, and the initial genotyping error rate, estimated from nine replicate samples, was 0.004 per reaction. However, an almost tenfold further reduction in the error rate was achieved by excluding 31 loci (21.5%) that exhibited unclear clustering patterns, manually editing clusters to allow rescoring of ambiguous or incorrect genotypes, and excluding 18 samples (3.8%) with unreliable genotypes. After stringent quality filtering, we also found a counter-intuitive negative relationship between in silico minor allele frequency and the conversion rate, suggesting that some of our assays may have been designed from paralogous loci. Nevertheless, we obtained over 45 000 individual SNP genotypes with a final error rate of 0.0005, indicating that the GoldenGate assay is eminently capable of generating large, high-quality data sets for non-model organisms. This has positive implications for future studies of the evolutionary, behavioural and conservation genetics of natural populations.

Monitoring adaptive genetic responses to environmental change

Widespread environmental changes including climate change, selective harvesting and landscape alterations now greatly affect selection regimes for most organisms. How animals and plants can adapt to these altered environments via contemporary evolution is thus of strong interest. We discuss how to use genetic monitoring to study adaptive responses via repeated analysis of the same populations over time, distinguishing between phenotypic and molecular genetics approaches. After describing monitoring designs, we develop explicit criteria for demonstrating adaptive responses, which include testing for selection and establishing clear links between genetic and environmental change. We then review a few exemplary studies that explore adaptive responses to climate change in Drosophila, selective responses to hunting and fishing, and contemporary evolution in Daphnia using resurrected resting eggs. We further review a broader set of 44 studies to assess how well they meet the proposed criteria, and conclude that only 23% fulfill all criteria. Approximately half (43%) of these studies failed to rule out the alternative hypothesis of replacement by a different, better-adapted population. Likewise, 34% of the studies based on phenotypic variation did not test for selection as opposed to drift. These shortcomings can be addressed via improved experimental designs and statistical testing. We foresee monitoring of adaptive responses as a future valuable tool in conservation biology, for identifying populations unable to evolve at sufficiently high rates and for identifying possible donor populations for genetic rescue. Technological advances will further augment the realization of this potential, especially next-generation sequencing technologies that allow for monitoring at the level of whole genomes.

Novel tools for conservation genomics: comparing two high-throughput approaches for SNP discovery in the transcriptome of the European hake

The growing accessibility to genomic resources using next-generation sequencing (NGS) technologies has revolutionized the application of molecular genetic tools to ecology and evolutionary studies in non-model organisms. Here we present the case study of the European hake (Merluccius merluccius), one of the most important demersal resources of European fisheries. Two sequencing platforms, the Roche 454 FLX (454) and the Illumina Genome Analyzer (GAII), were used for Single Nucleotide Polymorphisms (SNPs) discovery in the hake muscle transcriptome. De novo transcriptome assembly into unique contigs, annotation, and in silico SNP detection were carried out in parallel for 454 and GAII sequence data. High-throughput genotyping using the Illumina GoldenGate assay was performed for validating 1,536 putative SNPs. Validation results were analysed to compare the performances of 454 and GAII methods and to evaluate the role of several variables (e.g. sequencing depth, intron-exon structure, sequence quality and annotation). Despite well-known differences in sequence length and throughput, the two approaches showed similar assay conversion rates (approximately 43%) and percentages of polymorphic loci (67.5% and 63.3% for GAII and 454, respectively). Both NGS platforms therefore demonstrated to be suitable for large scale identification of SNPs in transcribed regions of non-model species, although the lack of a reference genome profoundly affects the genotyping success rate. The overall efficiency, however, can be improved using strict quality and filtering criteria for SNP selection (sequence quality, intron-exon structure, target region score).

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.