Associations between single nucleotide polymorphisms in candidate genes and growth rate in Arctic charr (Salvelinus alpinus L.)

Characterisation of Single Nucleotide Polymorphisms and Haplotypes of MSTN Associated with Growth Traits in European Sea Bass (Dicentrarchus labrax)

The myostatin (MSTN) gene, known as growth differentiation factor-8 (GDF-8), is a member of the transforming growth factor-β (TGF-β) superfamily and plays a specific inhibitory role during the critical phases of skeletal muscle mass development in vertebrates. This study was conducted to investigate MSTN polymorphisms in harvest size European sea bass reared in Turkey. Nine single nucleotide polymorphisms (SNPs) and two indels were identified in exons 1-3 of MSTN in the European sea bass population The associations between the g.16612A indel located in intron 1 and standard length were significant. The MSTN g.15252 T > A locus in intron 2 was significantly related to the total weight, fillet weight and standard length (P < 0.05). The relationship between the g.14873C > T locus in exon 3 of MSTN and standard height, head length, body length, pre-anal length, abdominal length, post-anal length and head width was significant (P < 0.05). According to the results of the haplotype analysis, two haplogroup and eight haplotype combinations were detected in the population. The haplogroup 2 had significant associations with all measured growth traits (P < 0.05). Thus, SNPs and haplotypes identified in this study could be useful for European sea bass breeding and marker-assisted selection.

Transcriptome analyses reveal differentially expressed genes associated with development of the palatal organ in bighead carp (Hypophthalmichthys nobilis)

The palatal organ is a filter-feeding related organ and occupies a considerable proportion of the head of bighead carp (Hypophthalmichthys nobilis), a large cyprinid fish intensive aquaculture in Asia. In this study, we performed RNA-seq of the palatal organ during growth periods of two (M2), six (M6) and 15 (M15) months of age after hatching. The numbers of differentially expressed genes (DEGs) were 1384, 481 and 1837 for M2 VS M6, M6 VS M15 and M2 VS M15 respectively. The following signaling pathways of energy metabolism and cytoskeleton function were enriched, including ECM-receptor interaction, Cardiac muscle contraction, Steroid biosynthesis and PPAR signaling pathway. Several members of collagen family (col1a1, col2a1, col6a2, col6a3, col9a2), Laminin gamma 1 (lamc1), integrin alpha 1 (itga1), Fatty acid binding protein 2 (fads2) and lipoprotein lipase (lpl), and Protein tyrosine kinase 7 (Ptk7) are candidate genes for growth and development of basic tissues of the palatal organ. Furthermore, taste-related genes such as fgfrl1, fgf8a, fsta and notch1a were also identified, which may be involved in the development of taste buds of the palatal organ. The transcriptome data obtained in this study provide insights into the understanding functions and development mechanisms of palatal organ, and potential candidate genes that may be related to the genetic modulation of head size of bighead carp.

SNP within the bovine ASB-3 gene and their association analysis with stature traits in three Chinese cattle breeds

ASB-3 is one of the 18 members of ASB gene family. As a special negative regulation factor of TNF-R2, ASB-3 inhibits the signal transduction of JNK-TNF-R2 and JNK-STAT signaling pathway by TNF-R2 protein. In this study, the genetic polymorphisms of ASB-3 were detected in total of 637 from Qinchuan, Jinnan and Xianan cattle using the sequence of mixed DNA pool, Tetra-primer ARMS-PCR and PCR-RFLP methods. Four mutation sites were detected including the g.C41255T, g.G74754A, and g.T75438C were synonymous mutation, whereas the g.C115213T was missense mutation (Pro > Ser). The associated analysis of four polymorphic loci of ASB-3 gene respectively with growth traits in the three cattle breeds. The result showed that SNP1 site was significantly related with Qinchuan cattle height and TT was the dominant genotype; SNP2 had a significant relationship with body length of Xianan cattle and cross department height of Qinchuan cattle, AA was the dominant genotype; SNP3 was significantly related to cross height of Xianan cattle, TT was the dominant genotype; SNP4 site was significantly correlated with body height of Xianan cattle and cross height of Jinnan cattle. Genotype combinations were only significantly correlated with the hucklebone width in the adult Qinchuan cattle. The combination genotype CTAGCTCC was outperformed other combination genotypes of Qinchuan cattle. The results showed that ASB-3 could be an important candidate gene and the four SNPs in ASB-3 can be used for molecular marker-assisted selection of four beef cattle breeds.

A SNP of the hemocyanin gene (LvHc) is a marker for high growth and ammonia-tolerance in Pacific white shrimp Litopenaeus vannamei

Expression levels of hemocyanin (LvHc), activating transcription factor 4 (LvAtf4), glutathione S-transferase (LvGst), caspase 2 (LvCasp2) and anti-lipopolysaccharide factor (LvAlf) were examined in the hepatopancreas of Pacific white shrimp Litopenaeus vannamei juveniles exposed to a lethal concentration of ammonia-N (32.15 mg/l). The expression levels of all transcripts except LvAlf were significantly greater (P < 0.05) in tolerant shrimp (Lv-AT; N = 30) that survived up to 72 h post treatment (hpt) than in susceptible shrimp (Lv-AS24 and Lv-AS72; N = 45 and 15), that died within 24 h or between 24 and 72 hpt, respectively. Subsequently, effects of non-lethal concentrations of ammonia-N (control, 10 and 20 mg/l) on the expression of LvHc in juvenile shrimp were examined. Compared to the control, expression levels of LvHc transcripts in hemocytes and the hepatopancreas of tested shrimp changed after exposure to ammonia-N. One SNP (C > T₅₄₅) was found in the LvHc₃₂₂ gene segment. Real-time PCR amplification of specific alleles (real-time PASA) was developed for detection of C > T₅₄₅ genotypes. Juveniles in the lethal exposure test that carried a C/T₅₄₅ genotype showed a greater average body weight and total length (8.46 ± 0.36 g and 10.05 ± 0.16 cm) than those with a C/C₅₄₅ genotype (7.48 ± 0.31 g and 9.60 ± 0.13 cm) (P < 0.05). Similar results were found in the second generation (G2) of a growth-improved stock (3 and 4 families of BIOTEC-G2-L1 and BIOTEC-G2-L2) and in commercially farmed shrimp (2 groups). Accordingly, expression levels and SNP of LvHc can serve as markers for selection high growth performance in ammonia-tolerant L. vannamei.

Salmonid Antibacterial Immunity: An Aquaculture Perspective

The aquaculture industry is continuously threatened by infectious diseases, including those of bacterial origin. Regardless of the disease burden, aquaculture is already the main method for producing fish protein, having displaced capture fisheries. One attractive sector within this industry is the culture of salmonids, which are (a) uniquely under pressure due to overfishing and (b) the most valuable finfish per unit of weight. There are still knowledge gaps in the understanding of fish immunity, leading to vaccines that are not as effective as in terrestrial species, thus a common method to combat bacterial disease outbreaks is the use of antibiotics. Though effective, this method increases both the prevalence and risk of generating antibiotic-resistant bacteria. To facilitate vaccine design and/or alternative treatment efforts, a deeper understanding of the teleost immune system is essential. This review highlights the current state of teleost antibacterial immunity in the context of salmonid aquaculture. Additionally, the success of current techniques/methods used to combat bacterial diseases in salmonid aquaculture will be addressed. Filling the immunology knowledge gaps highlighted here will assist in reducing aquaculture losses in the future.

Identification of a Growth-Associated Single Nucleotide Polymorphism (SNP) in Cyclin C of the Giant Tiger Shrimp Penaeus monodon

The full-length cDNA of cyclin C of the giant tiger shrimp Penaeus monodon (PmCyC) was isolated by RACE-PCR. It was 1443 bp in length containing an open reading frame (ORF) of 804 bp and 267 deduced amino acids. Tissue distribution analysis indicated that PmCyC was more abundantly expressed in ovaries and testes than other tissues of female and male juveniles (P < 0.05). A pair of primers was designed, and an amplification product of 403 bp containing an intron of 123 bp was obtained. Polymorphism of amplified PmCyC gene segments of the 5th (3-month-old G5, N = 30) and 7th (5-month-old G7, N = 18) generations of domesticated juveniles was analyzed. Four conserved SNPs (T>C₁₃₄, T>C₁₈₈, G>A₃₇₉, and T>C₃₈₂) were found within the examined sequences. A TaqMan genotyping assay was developed for detection of a T>C₁₃₄ SNP. Association analysis indicated that this SNP displayed significant association with body weight (P < 4.2e-10) and total length (P < 2e-09) of the examined G7 P. monodon (N = 419) with an allele substitution effect of 5.02 ± 0.78 g and 1.41 ± 0.19 cm, respectively. Juveniles with C/C₁₃₄ (22.80 ± 2.51 g and 12.97 ± 0.53 cm, N = 19) and T/C₁₃₄ (20.41 ± 0.93 g and 12.77 ± 0.21 cm, N = 129) genotypes exhibited a significantly greater average body weight and total length than those with a T/T₁₃₄ genotype (14.72 ± 0.53 g and 11.37 ± 0.13 cm, N = 271) (P < 0.05).

A first genetic linage map construction and QTL mapping for growth traits in Larimichthys polyactis

Larimichthys polyactis is a commercially important marine fish species in Eastern Asia, yet very few genetic resources exist. In particular, genetic linkage maps are critical tools for genetic breeding. In this study, we generated a high resolution linkage map from a family of 110 individuals and their parents by resequencing the individuals. 3,802 effective SNPs were mapped to 24 linkage groups (LGs). The map spanned 2,567.39 cm, with an average marker interval of 0.66 cm. We used the map to conduct QTL analysis for growth traits, and found 31 markers were significantly associated with growth-related traits. Specifically, three SNPs were identified for total length, nineteen SNPs for body length, and nine SNPs for body weight. The identified SNPs could explain 15.2-22.6% of the phenotypic variation. SNPs associated with growth traits were distributed on LG6 and LG11, and candidate genes included, kif26b, bat1, gna1, gbgt1, and amfr, which may regulate growth. The linkage map and mapped QTLs would be useful for improving the quality of L. polyactis via marker-assisted selection.

Identification of novel genes significantly affecting growth in catfish through GWAS analysis

Growth is the most important economic trait in aquaculture. Improvements in growth-related traits can enhance production, reduce costs and time to produce market-size fish. Catfish is the major aquaculture species in the United States, accounting for 65% of the US finfish production. However, the genes underlying growth traits in catfish were not well studied. Currently, the majority of the US catfish industry uses hybrid catfish derived from channel catfish female mated with blue catfish male. Interestingly, channel catfish and blue catfish exhibit differences in growth-related traits, and therefore the backcross progenies provide an efficient system for QTL analysis. In this study, we conducted a genome-wide association study for catfish body weight using the 250 K SNP array with 556 backcross progenies generated from backcross of male F1 hybrid (female channel catfish × male blue catfish) with female channel catfish. A genomic region of approximately 1 Mb on linkage group 5 was found to be significantly associated with body weight. In addition, four suggestively associated QTL regions were identified on linkage groups 1, 2, 23 and 24. Most candidate genes in the associated regions are known to be involved in muscle growth and bone development, some of which were reported to be associated with obesity in humans and pigs, suggesting that the functions of these genes may be evolutionarily conserved in controlling growth. Additional fine mapping or functional studies should allow identification of the causal genes for fast growth in catfish, and elucidation of molecular mechanisms of regulation of growth in fish.

Genetic polymorphisms and DNA methylation in exon 1 CpG-rich regions of PACAP gene and its effect on mRNA expression and growth traits in half smooth tongue sole (Cynoglossus semilaevis)

The pituitary adenylate cyclase activating polypeptide (PACAP) is a new type of hypophysiotropic hormone and plays an important role in regulating the synthesis and secretion of growth hormone and gonadotropin. The research on the relationship between PACAP and different growth traits would contribute to explain its function during the process of growth. Moreover, epigenetic modifications, especially DNA methylation at the CpG sites of the SNPs, play important roles in regulating gene expression. The results suggest that a SNP mutation (c.C151G) in the PACAP gene of male half smooth tongue sole (Cynoglossus semilaevis) is significantly associated with growth traits and serum physiological and biochemical parameters such as inorganic phosphorus (P < 0.05). The SNP is located in a CpG-rich region of exon 1. Intriguingly, the transition (C→G) added a new methylation site of PACAP gene. This SNP was also significantly related to the expression and methylation level of PACAP (P < 0.05). Individuals with GG genotype had faster growth rates than those of CG and CC genotypes. Moreover, GG genotype had significantly higher PACAP expression level and lower methylation level than CG and CC genotypes. In the serum indexes, only inorganic phosphorus content within GG genotypes was significantly higher than CC genotypes. This implied that the mutation and methylation status of PACAP gene could influence growth traits and this locus could be considered as a candidate genetic or epigenetic marker for Cynoglossus semilaevis molecular breeding.

Transcriptome Analysis for Identification of Genes Related to Gonad Differentiation, Growth, Immune Response and Marker Discovery in The Turbot (Scophthalmus maximus)

Background

Turbot Scophthalmus maximus is an economically important species extensively aquacultured in China. The genetic selection program is necessary and urgent for the sustainable development of this industry, requiring more and more genome background knowledge. Transcriptome sequencing is an excellent alternative way to identify transcripts involved in specific biological processes and exploit a considerable quantity of molecular makers when no genome sequences are available. In this study, a comprehensive transcript dataset for major tissues of S. maximus was produced on basis of an Illumina platform.

Results

Total RNA was isolated from liver, spleen, kidney, cerebrum, gonad (testis and ovary) and muscle. Equal quantities of RNA from each type of tissues were pooled to construct two cDNA libraries (male and female). Using the Illumina paired-end sequencing technology, nearly 44.22 million clean reads in length of 100 bp were generated and then assembled into 106,643 contigs, of which 71,107 were named unigenes with an average length of 892 bp after the elimination of redundancies. Of these, 24,052 unigenes (33.83% of the total) were successfully annotated. GO, KEGG pathway mapping and COG analysis were performed to predict potential genes and their functions. Based on our sequence analysis and published documents, many candidate genes with fundamental roles in sex determination and gonad differentiation (dmrt1), growth (ghrh, myf5, prl/prlr) and immune response (TLR1/TLR21/TLR22, IL-15/IL-34), were identified for the first time in this species. In addition, a large number of credible genetic markers, including 21,192 SSRs and 8,642 SNPs, were identified in the present dataset.

Conclusion

This informative transcriptome provides valuable new data to increase genomic resources of Scophthalmus maximus. The future studies of corresponding gene functions will be very useful for the management of reproduction, growth and disease control in turbot aquaculture breeding programs. The molecular markers identified in this database will aid in genetic linkage analyses, mapping of quantitative trait loci, and acceleration of marker assisted selection programs.

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.

GHRH, PRP-PACAP and GHRHR Target Sequencing via an Ion Torrent Personal Genome Machine Reveals an Association with Growth in Orange-Spotted Grouper (Epinephelus coioides)

Growth hormone-releasing hormone (GHRH) and the receptor, GHRHR, constitute important components of the hypothalamus-pituitary growth axis and act on the downstream growth hormone (GH). PACAP-related peptide/pituitary adenylate cyclase activating polypeptide (PRP-PACAP) is a paralog of GHRH. These genes all play key roles in development and growth patterns. To improve the quality of cultured fish strains, natural genetic variation must be examined and understood. A mixed linear model has been widely used in association mapping, taking the population structures and pairwise kinship patterns into consideration. In this study, a mass cross population of orange-spotted grouper (Epinephelus coioides) was examined. These candidate genes were found to harbor low nucleotide diversity (θ_w from 0.00154 to 0.00388) and linkage disequilibrium levels (delay of 50% within 2 kbp). Association mapping was employed, and two single-nucleotide polymorphisms (KR269823.1:g.475A>C and KR269823.1:g.2143T>C) were found to be associated with growth (false discovery rate Q < 0.05), explaining 9.0%–17.0% of the phenotypic variance. The association of KR269823.1:g.2143T>C was also found via haplotype-based association (p < 0.05). The identified associations offer new insights into gene functions, and the associated single-nucleotide polymorphisms (SNPs) may be used for breeding purposes.

Genetic and genomic analyses for economically important traits and their applications in molecular breeding of cultured fish

The traits of cultured fish must continually be genetically improved to supply high-quality animal protein for human consumption. Economically important fish traits are controlled by multiple gene quantitative trait loci (QTL), most of which have minor effects, but a few genes may have major effects useful for molecular breeding. In this review, we chose relevant studies on some of the most intensively cultured fish and concisely summarize progress on identifying and verifying QTLs for such traits as growth, disease and stress resistance and sex in recent decades. The potential applications of these major-effect genes and their associated markers in marker-assisted selection and molecular breeding, as well as future research directions are also discussed. These genetic and genomic analyses will be valuable for elucidating the mechanisms modulating economically important traits and to establish more effective molecular breeding techniques in fish.

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.

Novel single nucleotide polymorphisms of the insulin-like growth factor-I gene and their associations with growth traits in common carp (Cyprinus carpio L.)

Insulin-like growth factor-I (IGF-I) plays an important role in the growth and development of vertebrates. To study polymorphisms of IGF-I, we screened a total of 4555 bp of genomic sequences in four exons and partial introns for the discovery of single nucleotide polymorphism (SNP) in common carp (Cyprinus carpio). Three SNPs (g.3759T>G, g.7627T>A and g.7722T>C) in intron 2 and a nonsynonymous SNP (g.7892C>T) in exon 3 were identified in a pilot population including random parents and their progenies. 289 progenies were further genotyped for studying possible associations between genotypes or combined genotypes and growth traits. The results showed that the locus g.7627T>A was significantly associated with body weight and body length, and fish with genotype AA had a mean body weight 5.9% higher than those with genotype TT. No significant associations were observed between genotypes of other loci and growth traits. However, when both g.7627T>A and g.7722T>C were considered, the combined genotype TT/TT was extremely associated with the lowest values of body length and body weight and the highest K value in comparison with other diplotypes (p < 0.01). These results suggest that genotype AA at g.7627T>A and its combined genotypes with alleles from another locus have positive effects on growth traits, which would be a candidate molecular marker for further studies in marker-assisted selection in common carp.

Single nucleotide polymorphisms in the insulin-like growth factor 1 (IGF1) gene are associated with growth-related traits in farmed Atlantic salmon

Understanding the genetic basis of variation in traits related to growth and fillet quality in Atlantic salmon is of importance to the aquaculture industry. Several growth-related QTL have been identified via the application of genetic markers. The IGF1 gene is considered a highly conserved and crucial growth-regulating gene in salmonid species. However, the association between polymorphisms in the IGF1 gene and growth-related traits in Atlantic salmon is unknown. Therefore, in this study, regions of the Atlantic salmon IGF1 gene were sequenced, aligned and compared across individuals. Three SNPs were identified in the putative promoter (SNP1, g.5763G>T; GenBank no. AGKD01012745), intron 1 (SNP2, g.7292C>T; GenBank no. AGKD01012745) and intron 3 (SNP3, g.4671A>C; GenBank no. AGKD01133398) regions respectively. These SNPs were genotyped in a population of 4800 commercial Atlantic salmon with data on several weight and fillet traits measured at harvest (at approximately 3 years of age). In a mixed model, association analysis of individual SNPs, SNP1 and SNP3 were both significantly associated with several weight traits (P < 0.05). The estimated additive effect on overall harvest weight was approximately 35 and 110 g for SNPs 1 and 3 respectively. A haplotype analysis confirmed the association between genetic variation in the IGF1 gene with overall body weight (P < 0.05) and fillet component traits (P < 0.05). Our findings suggest the identified nucleotide polymorphisms of the IGF1 gene may either affect farmed Atlantic salmon growth directly or be in population-wide linkage disequilibrium with causal variation, highlighting their possible utility as candidates for marker-assisted selection in the aquaculture industry.

A new anti-lipopolysaccharide factor (ALF) gene with its SNP polymorphisms related to WSSV-resistance of Litopenaeus vannamei

Anti-lipopolysaccharide factors (ALFs) of crustacean play an important role against bacteria or virus infection. In this study, the cDNA sequence and genomic sequence of one new isoform of ALF designated as nLvALF1 were reported. The open reading frame (ORF) of nLvALF1 consisted of 369 bp encoding 123 amino acids and the genomic structure of nLvALF1 comprised four introns and three exons. The predicted pI of the deduced protein was 8.82 and the molecular weight (MW) was 13.72 KDa. The deduced amino acid sequence of nLvALF1 contained a typical functional domain of ALF: LPS-binding domain. Phylogenetic analysis showed that nLvALF1 had the closest relationship with FcALF1 from Fenneropenaeus chinensis. The nLvALF1 was specifically expressed in lymphoid organ (Oka) of shrimp. Its transcriptional level was significantly up-regulated after white spot syndrome virus (WSSV) challenge, suggesting that nLvALF1 might participate in defense against WSSV in Litopenaeus vannamei. In order to search potential genetic markers associated with WSSV-resistance, we scanned the polymorphisms of the genomic fragment with 397 bp where the LPS-binding domain encoding sequence located and 18 SNPs were found. The distribution frequency of these SNPs was analyzed in WSSV susceptible shrimp and resistant shrimp separately. Significant differences existed in allelic frequencies at loci g.1361-T > C, g.1370-T > C, g.1419-T > A between the WSSV-resistant group and the WSSV-susceptible/susceptible group (P < 0.05). The specific haplotype CT consisted of g.1415-C > A and g.1419-T > A was associated with susceptibility to WSSV (P < 0.05). These findings provide theoretical support for selection of WSSV-resistant varieties of L. vannamei.

Genetic diversity of upper Lake Constance whitefish Coregonus spp. under the influence of fisheries: a DNA study based on archived scale samples from 1932, 1975 and 2006

This investigation examined changes in the genetic diversity of pelagic upper Lake Constance (ULC) whitefish Coregonus wartmanni population before and after the alteration of fishery methods and management from 1932 to 2006. The study spans a period of pronounced changes in trophic status of the lake and transitions from traditional relatively unselective pelagic seine (Klusgarn) fishing to highly size-selective nylon gillnet techniques. In addition, supportive breeding and stocking became most popular during the phase of eutrophication in the 1970s. The main hypothesis is that size-selective fisheries and breeding lead to an overall decrease in genetic variability over time. A total of 215 archived C. wartmanni scale samples from 1932, 1975 and 2006 were analysed by genotyping 11 microsatellite loci. A comparison of population genetic parameters, including allelic richness, observed and expected heterozygosities, and estimates of effective population sizes, suggests that the genetic diversity of C. wartmanni population has not decreased. The appearance of new alleles in the gene pool in 1975 and 2006 may be indicative of admixture with other forms in the lake or with stocked allochthonous forms. Overall, the fisheries management practice in ULC, including the effects of size-selective fisheries, supportive breeding and stocking, have not significantly altered the genetic diversity of Coregonus spp. over an 80 year period.

A candidate gene association study for growth performance in an improved giant freshwater prawn (Macrobrachium rosenbergii ) culture line

A candidate gene approach using type I single nucleotide polymorphism (SNP) markers can provide an effective method for detecting genes and gene regions that underlie phenotypic variation in adaptively significant traits. In the absence of available genomic data resources, transcriptomes were recently generated in Macrobrachium rosenbergii to identify candidate genes and markers potentially associated with growth. The characterisation of 47 candidate loci by ABI re-sequencing of four cultured and eight wild samples revealed 342 putative SNPs. Among these, 28 SNPs were selected in 23 growth-related candidate genes to genotype in 200 animals selected for improved growth performance in an experimental GFP culture line in Vietnam. The associations between SNP markers and individual growth performance were then examined. For additive and dominant effects, a total of three exonic SNPs in glycogen phosphorylase (additive), heat shock protein 90 (additive and dominant) and peroxidasin (additive), and a total of six intronic SNPs in ankyrin repeats-like protein (additive and dominant), rolling pebbles (dominant), transforming growth factor-β induced precursor (dominant), and UTP-glucose-1-phosphate uridylyltransferase 2 (dominant) genes showed significant associations with the estimated breeding values in the experimental animals (P =0.001-0.031). Individually, they explained 2.6-4.8 % of the genetic variance (R²=0.026-0.048). This is the first large set of SNP markers reported for M. rosenbergii and will be useful for confirmation of associations in other samples or culture lines as well as having applications in marker-assisted selection in future breeding programs.

Polymorphisms in Myostatin Gene and associations with growth traits in the common carp (Cyprinus carpio L.)

Myostatin (MSTN) is a member of the transforming growth factor-β superfamily that negatively regulates skeletal muscle development and growth. In the present study, partial genomic fragments of MSTN were screened for single nucleotide polymorphisms (SNPs) in selected common carp individuals from wild populations, and two SNPs in intron 2 (c.371 + 749A > G, c.371 + 781T > C) and two synonymous SNPs in exon 3 (c.42A > G, c.72C > T) were identified. Genotyping by direct sequencing of polymerase chain reaction (PCR) products for these four SNPs were performed in 162 individuals from a commercial hatchery population. Association analysis showed that two SNPs in exon 3 were significantly associated with body weight (BW) and condition factor (K), and haplotype analyses revealed that haplotype H7H8 showed better growth performance. Our results demonstrated that some of the SNPs in MSTN may have positive effects on growth traits and suggested that MSTN could be a candidate gene for growth and marker-assisted selection in common carp.

The genetic architecture of growth rate in juvenile Takifugu species

Closely related species have often evolved dramatic differences in body size. Takifugu rubripes (fugu) is a large marine pufferfish whose genome has been sequenced, whereas T. niphobles is the smallest species among Takifugu. We show that, unsurprisingly, the juvenile growth rate of T. rubripes is higher than that of T. niphobles in a laboratory setting. We produced F(2) progenies of their F(1) hybrids and found one quantitative trait locus (QTL) significantly associated with variation in juvenile body size. This QTL region (3.5 Mb) contains no known genes directly related to growth phenotype (such as IGFs) except Fgf21, which inhibits growth hormone signaling in mouse. The QTL in Takifugu spp. is distinct from the region previously known to control body size variations in stickleback or tilapia. Our results suggest that in the fish tested herein, genomic regions underlying body size evolution might have different genetic origins. They also suggest that many diverse traits in Takifugu spp. are amenable to genetic mapping.

Genetic diversity of the sulfotransferase-like gene and one nonsynonymous SNP associated with growth traits of clam, Meretrix meretrix

The sulfotransferase (SULT) is an important factor in phase II metabolism. Three single nucleotide polymorphisms (SNPs; MmeSUl-696, MmeSUl-714 and MmeSUl-806) were discovered in a SULT-like gene of Meretrix meretrix. A modified allele-specific PCR detection assay was developed to perform polymorphism analysis of the three SNPs in two families and two populations of M. meretrix. Results showed that the genotype of the three SNPs is accordance to the expected Mendelian segregation ratios and conforms to Hardy-Weinberg Equilibrium (HWE) in families and populations, respectively. Moreover, sampled clams (n = 151) from the two populations were tested for the associations between SNPs and growth traits. The single SNP association analysis indicated MmeSUl-806 significantly (P < 0.05) correlated with growth traits of M. meretrix. A further haplotype-trait association analysis confirmed that the three SNPs were in linkage disequilibrium and showed that the haplotype TTG and TCG associated with low growth in the two populations, respectively, which supported that MmeSUl-806 was a growth related marker. This study suggested SULT gene might be a candidate gene related to growth in bivalves and MmeSUl-806 could be a gene marker in linkage map construction and breeding programs of M. meretrix.

Robust linear regression methods in association studies

MOTIVATION

It is well known that data deficiencies, such as coding/rounding errors, outliers or missing values, may lead to misleading results for many statistical methods. Robust statistical methods are designed to accommodate certain types of those deficiencies, allowing for reliable results under various conditions. We analyze the case of statistical tests to detect associations between genomic individual variations (SNP) and quantitative traits when deviations from the normality assumption are observed. We consider the classical analysis of variance tests for the parameters of the appropriate linear model and a robust version of those tests based on M-regression. We then compare their empirical power and level using simulated data with several degrees of contamination.

RESULTS

Data normality is nothing but a mathematical convenience. In practice, experiments usually yield data with non-conforming observations. In the presence of this type of data, classical least squares statistical methods perform poorly, giving biased estimates, raising the number of spurious associations and often failing to detect true ones. We show through a simulation study and a real data example, that the robust methodology can be more powerful and thus more adequate for association studies than the classical approach.

AVAILABILITY

The code of the robustified version of function lmekin() from the R package kinship is provided as Supplementary Material.

Analysis of single nucleotide polymorphisms in three chromosomes of European sea bass Dicentrarchus labrax

Single nucleotide polymorphisms (SNPs) are believed to contain relevant information and have been therefore extensively used as genetic markers in population and conservation genetics, and molecular ecology studies. This study reports on the identification of potential SNPs in a diploid European sea bass Dicentrarchus labrax genome by using reference sequences from three assembled chromosomes and mapping all WGS datasets onto them (3× Sanger, 3× 454 and 20× SOLEXA). A total of 20,779 SNPs were identified over the 1469 gene loci and intergenic space analysed. Within chromosomes the occurrence of SNPs was the lowest in exons and higher in introns and intergenic regions, which may be explained by the fact, that coding regions are under strong selective pressure to maintain their biological function. The ratio of nonsynonymous to synonymous mutations was smaller than one for all the chromosomes, suggesting that most of deleterious nonsynonymous mutations were eliminated by negative selection. SNPs were not uniformly distributed over the chromosomes. Two chromosomes exhibited large regions with extremely low SNP density, which might represent homozygous regions in the diploid genome. The results of this study show how SNP detection can take profit from sequencing a single diploid individual, but also uncover the limits of such an approach. SNPs that have been identified will support marker development for genetic linkage mapping, population genetics and aquaculture related questions in general.

Novel methodologies in marine fish larval nutrition

Major gaps in knowledge on fish larval nutritional requirements still remain. Small larval size, and difficulties in acceptance of inert microdiets, makes progress slow and cumbersome. This lack of knowledge in fish larval nutritional requirements is one of the causes of high mortalities and quality problems commonly observed in marine larviculture. In recent years, several novel methodologies have contributed to significant progress in fish larval nutrition. Others are emerging and are likely to bring further insight into larval nutritional physiology and requirements. This paper reviews a range of new tools and some examples of their present use, as well as potential future applications in the study of fish larvae nutrition. Tube-feeding and incorporation into Artemia of (14)C-amino acids and lipids allowed studying Artemia intake, digestion and absorption and utilisation of these nutrients. Diet selection by fish larvae has been studied with diets containing different natural stable isotope signatures or diets where different rare metal oxides were added. Mechanistic modelling has been used as a tool to integrate existing knowledge and reveal gaps, and also to better understand results obtained in tracer studies. Population genomics may assist in assessing genotype effects on nutritional requirements, by using progeny testing in fish reared in the same tanks, and also in identifying QTLs for larval stages. Functional genomics and proteomics enable the study of gene and protein expression under various dietary conditions, and thereby identify the metabolic pathways which are affected by a given nutrient. Promising results were obtained using the metabolic programming concept in early life to facilitate utilisation of certain nutrients at later stages. All together, these methodologies have made decisive contributions, and are expected to do even more in the near future, to build a knowledge basis for development of optimised diets and feeding regimes for different species of larval fish.

Single nucleotide polymorphisms within the estrogen receptor beta gene are linked with reproductive indices in Japanese flounder, Paralichthys olivaceus

The objectives of this study were to characterize polymorphisms within the coding region of estrogen receptor beta (ERbeta) gene in a population of 57 female Japanese flounder (Paralichthys olivaceus) and to analyze the association of ERbeta polymorphisms with reproductive indices by polymerase chain reaction single-strand conformation polymorphism (PCR-SSCP). Two single nucleotide polymorphisms (SNPs), SNP1 (c.577delC) and SNP2 [c.A891T (p.Gln114Leu)], were identified in the ERbeta gene. A one-way ANOVA revealed that SNP1 was significantly associated with the gonadosomatic index (GSI) in female Japanese flounder (P < 0.05). And SNP2 was significantly associated with the serum 17beta-estradiol (E2) level and GSI (P < 0.05). Individuals with genotype AB of SNP2 had significantly higher serum E2 level and GSI than those of genotype AA (P < 0.05). Moreover, the hepatosomatic index (HSI), a marker for genetic effects, was significantly higher for diplotype D2 compared with the other three diplotypes (P < 0.05). These results obtained in this study suggested that SNP2 could influence reproductive endocrinology of female Japanese flounder and be useful as a potential candidate genetic marker for the selection of reproductive indices in female Japanese flounder.

A new method for SNP discovery

Single nucleotide polymorphisms (SNPs) are high-density natural sequence variations in genomes. They are considered to be the major genetic source of phenotypic variability within a given species and serve as excellent genetic markers. SNPs are useful in identifying candidate genes that contribute to disease and phenotypic traits. In non-model organisms, the application of SNPs has been limited, because of the expense and technical difficulties entailed in currently available SNP isolation techniques. In the present study, we have developed a rapid and effective method to isolate SNPs throughout the genome randomly. The DNA fragments containing SNPs could be isolated efficiently from background DNA. We analyzed ten isolated DNA fragments with this method in half-smooth tongue sole (Cynoglossus semilaevis)—a newly exploited and commercially important cultured marine flatfish in China—and found that nine of the fragments contained SNPs. The findings were confirmed successfully in different individuals. The method presented here is cost-effective and applicable to essentially any organism.

Persistence of duplicated PAC1 receptors in the teleost, Sparus auratus

Background:

Duplicated genes are common in vertebrate genomes. Their persistence is assumed to be either a consequence of gain of novel function (neofunctionalisation) or partitioning of the function of the ancestral molecule (sub-functionalisation). Surprisingly few studies have evaluated the extent of such modifications despite the numerous duplicated receptor and ligand genes identified in vertebrate genomes to date. In order to study the importance of function in the maintenance of duplicated genes, sea bream (Sparus auratus) PAC₁ receptors, sequence homologues of the mammalian receptor specific for PACAP (Pituitary Adenylate Cyclase-Activating Polypeptide), were studied. These receptors belong to family 2 GPCRs and most of their members are duplicated in teleosts although the reason why both persist in the genome is unknown.

Results:

Duplicate sea bream PACAP receptor genes (sbPAC₁A and sbPAC₁B), members of family 2 GPCRs, were isolated and share 77% amino acid sequence identity. RT-PCR with specific primers for each gene revealed that they have a differential tissue distribution which overlaps with the distribution of the single mammalian receptor. Furthermore, in common with mammals, the teleost genes undergo alternative splicing and a PAC₁Ahop1 isoform has been characterised. Duplicated orthologous receptors have also been identified in other teleost genomes and their distribution profile suggests that function may be species specific. Functional analysis of the paralogue sbPAC₁s in Cos7 cells revealed that they are strongly stimulated in the presence of mammalian PACAP₂₇ and PACAP₃₈ and far less with VIP (Vasoactive Intestinal Peptide). The sbPAC₁ receptors are equally stimulated (LOG_EC50 values for maximal cAMP production) in the presence of PACAP₂₇ (-8.74 ± 0.29 M and -9.15 ± 0.21 M, respectively for sbPAC₁A and sbPAC₁B, P > 0.05) and PACAP₃₈ (-8.54 ± 0.18 M and -8.92 ± 0.24 M, respectively for sbPAC₁A and sbPAC₁B, P > 0.05). Human VIP was found to stimulate sbPAC₁A (-7.23 ± 0.20 M) more strongly than sbPAC₁B (-6.57 ± 0.14 M, P < 0.05) and human secretin (SCT), which has not so far been identified in fish genomes, caused negligible stimulation of both receptors.

Conclusion:

The existence of functionally divergent duplicate sbPAC₁ receptors is in line with previously proposed theories about the origin and maintenance of duplicated genes. Sea bream PAC₁ duplicate receptors resemble the typical mammalian PAC₁, and PACAP peptides were found to be more effective than VIP in stimulating cAMP production, although sbPAC₁A was more responsive for VIP than sbPAC₁B. These results together with the highly divergent pattern of tissue distribution suggest that a process involving neofunctionalisation occurred after receptor duplication within the fish lineage and probably accounts for their persistence in the genome. The characterisation of further duplicated receptors and their ligands should provide insights into the evolution and function of novel protein-protein interactions associated with the vertebrate radiation.

Genomic organization of the IGF1, IGF2, MYF5, MYF6 and GRF/PACAP genes across Salmoninae genera

Whole-genome duplication in the ancient ray-finned fish and subsequent tetraploidization in the ancestor to the salmonids have complicated genomic and candidate gene studies in these organisms as many genes with multiple copies are present throughout their genomes. In an attempt to identify genes with a potential influence on growth and development, we investigated the genomic positions of insulin-like growth factors 1 and 2 (IGF1, IGF2), myogenic factors 5 and 6 (MYF5, MYF6) and growth hormone-releasing factor/pituitary adenylate cyclase-activating polypeptide (GRF/PACAP) in three salmonid species: rainbow trout (Oncorhynchus mykiss), Atlantic salmon (Salmo salar) and Arctic charr (Salvelinus alpinus). Our results suggest a tight association between the IGF1, MYF5 and MYF6 genes in all three species. We further localized the duplicated copies of IGF1 to the homeologous linkage groups RT-7/15 in rainbow trout and AC-3/24 in Arctic charr, and the two copies of MYF6 to homeologous linkage groups AS-22/24 in Atlantic salmon. Localization of GRF/PACAP to RT-7, AS-31 and AC-27 and IGF2 to RT-27, AS-2 and AC-4 in rainbow trout, Atlantic salmon and Arctic charr respectively is consistent with previously reported homologies among these chromosomal segments identified using other genetic markers. However, localization of the second copy of GRF/PACAP to RT-19 and AC-14 and the duplicated copy of IGF2 to AC-19 suggest a possible new homology/homeology between these chromosomes. These results might also be an indication of a more ancient polyploidization event that occurred deep in the ray-finned fish lineage.

The role of the IGF-I system for vitellogenesis in maturing female sterlet, Acipenser ruthenus Linnaeus, 1758

Transition from previtellogeneic to vitellogenic oocyte growth is a critical phase for folliculogenesis in sturgeon and may often be postponed for several years. Recent findings on the involvement of insulin-like growth factor I (IGF-I) in cell differentiation processes of oocyte follicle and ovarian steroidogenesis of teleosts in vitro led to the hypothesis that paracrine IGF-I could function as a potential trigger in vivo. For the first time, IGF-I and its corresponding receptor (IGF-IR) were identified in a non-teleostean fish. Real-time PCR assays for IGF-I and IGF-IR mRNA were established, normalising mRNA expression of the target genes to beta-microglobulin (beta2m). We clearly show that expression of IGF-I in the gonad is a substantial source for IGF-I-mediated effects in follicles compared to liver, brain, muscle and adipose tissue. Among these tissues, IGF-IR mRNA was highest in the gonad. With regard to different cohorts of coexisting follicles, highest expression of IGF-I and IGF-IR were met in developing follicles, indicating that IGF-I functions as an intraovarian modulator of follicle faith. Comparing previtellogenic follicles in females that matured within two years with non-maturing females f the same age, revealed an increases of 2.3-fold for IGF-I and 2.8-fold for IGF-IR mRNA expression in maturing females. These findings implicate an important role of paracrine IGF-I in early vitellogenesis and identify it as candidate vitellogenesis inducing factor (VIF), determining the faith of the follicle.

Single nucleotide polymorphism (SNP) discovery in duplicated genomes: intron-primed exon-crossing (IPEC) as a strategy for avoiding amplification of duplicated loci in Atlantic salmon (Salmo salar) and other salmonid fishes

Background

Single nucleotide polymorphisms (SNPs) represent the most abundant type of DNA variation in the vertebrate genome, and their applications as genetic markers in numerous studies of molecular ecology and conservation of natural populations are emerging. Recent large-scale sequencing projects in several fish species have provided a vast amount of data in public databases, which can be utilized in novel SNP discovery in salmonids. However, the suggested duplicated nature of the salmonid genome may hamper SNP characterization if the primers designed in conserved gene regions amplify multiple loci.

Results

Here we introduce a new intron-primed exon-crossing (IPEC) method in an attempt to overcome this duplication problem, and also evaluate different priming methods for SNP discovery in Atlantic salmon (Salmo salar) and other salmonids. A total of 69 loci with differing priming strategies were screened in S. salar, and 27 of these produced ~13 kb of high-quality sequence data consisting of 19 SNPs or indels (one per 680 bp). The SNP frequency and the overall nucleotide diversity (3.99 × 10^-4) in S. salar was lower than reported in a majority of other organisms, which may suggest a relative young population history for Atlantic salmon. A subset of primers used in cross-species analyses revealed considerable variation in the SNP frequencies and nucleotide diversities in other salmonids.

Conclusion

Sequencing success was significantly higher with the new IPEC primers; thus the total number of loci to screen in order to identify one potential polymorphic site was six times less with this new strategy. Given that duplication may hamper SNP discovery in some species, the IPEC method reported here is an alternative way of identifying novel polymorphisms in such cases.

Insulin-like growth factor signaling in fish

The insulin-like growth factor (IGF) system plays a central role in the neuroendocrine regulation of growth in all vertebrates. Evidence from studies in a variety of vertebrate species suggest that this growth factor complex, composed of ligands, receptors, and high-affinity binding proteins, evolved early during vertebrate evolution. Among nonmammalian vertebrates, IGF signaling has been studied most extensively in fish, particularly teleosts of commercial importance. The unique life history characteristics associated with their primarily aquatic existence has fortuitously led to the identification of novel functions of the IGF system that are not evident from studies in mammals and other tetrapod vertebrates. Furthermore, the emergence of the zebrafish as a preferred model for development genetics has spawned progress in determining the requirements for IGF signaling during vertebrate embryonic development. This review is intended as a summary of our understanding of IGF signaling, as revealed through research into the expression, function, and evolution of IGF ligands, receptors, and binding proteins in fish.

Distribution of genetic variation in the growth hormone 1 gene in Atlantic salmon (Salmo salar) populations from Europe and North America

The level and hierarchical distribution of genetic variation in complete sequences of the Atlantic salmon (Salmo salar) growth hormone (GH1) gene were investigated in populations from Europe and North America with a view to inferring the major evolutionary forces affecting genetic variation at this locus. Seventeen polymorphic sites were identified in complete sequences from nine populations, with levels of noncoding (intron and untranslated region sequences) nucleotide diversity being similar to those observed in other species. No variation, however, was observed in exonic sequences, indicating that nucleotide diversity in the Atlantic salmon GH1 gene is three and 25 times less than that estimated for human and Drosophila coding sequences, respectively. This suggests that purifying selection is the predominant contemporary force controlling the molecular evolution of GH1 coding sequences. Comparison of haplotype relationships within and between populations indicated that differentiation between populations from Europe and North America was greater than within-continent comparisons. However, several haplotypes observed in the northernmost European populations were more similar to those observed in North American than to any other haplotypes observed in Europe. This is most likely to be a result of historical, rather than contemporary, gene flow. Neutrality test statistics, such as Tajima's D, were significantly positive in the European populations in which North American-like haplotypes were observed. Although a positive Tajima's D is commonly interpreted as the signal of balancing selection, a more likely explanation in this case is that either historical migration or ascertainment bias, rather than within population local adaptation, has given rise to an excess of intermediate frequency alleles.

Expressed Sequence Tag-Linked Microsatellites as a Source of Gene-Associated Polymorphisms for Detecting Signatures of Divergent Selection in Atlantic Salmon (Salmo salar L.)

The prediction that selection affects the genome in a locus-specific way also affecting flanking neutral variation, known as genetic hitchhiking, enables the use of polymorphic markers in noncoding regions to detect the footprints of selection. However, as the strength of the selective footprint on a locus depends on the distance from the selected site and will decay with time due to recombination, the utilization of polymorphic markers closely linked to coding regions of the genome should increase the probability of detecting the footprints of selection as more gene-containing regions are covered. The occurrence of highly polymorphic microsatellites in the untranslated regions of expressed sequence tags (ESTs) is a potentially useful source of gene-associated polymorphisms which has thus far not been utilized for genome screens in natural populations. In this study, we searched for the genetic signatures of divergent selection by screening 95 genomic and EST-derived mini- and microsatellites in eight natural Atlantic salmon, Salmo salar L., populations from different spatial scales inhabiting contrasting natural environments (salt-, brackish, and freshwater habitat). Altogether, we identified nine EST-associated microsatellites, which exhibited highly significant deviations from the neutral expectations using different statistical methods at various spatial scales and showed similar trends in separate population samples from different environments (salt-, brackish, and freshwater habitats) and sea areas (Barents vs. White Sea). We consider these ESTs as the best candidate loci affected by divergent selection, and hence, they serve as promising genes associated with adaptive divergence in Atlantic salmon. Our results demonstrate that EST-linked microsatellite genome scans provide an efficient strategy for discovering functional polymorphisms, especially in nonmodel organisms.

Quantitative trait locus analyses and the study of evolutionary process

The past decade has seen a proliferation of studies that employ quantitative trait locus (QTL) approaches to diagnose the genetic basis of trait evolution. Advances in molecular techniques and analytical methods have suggested that an exact genetic description of the number and distribution of genes affecting a trait can be obtained. Although this possibility has met with some success in model systems such as Drosophila and Arabidopsis, the pursuit of an exact description of QTL effects, i.e. individual gene effect, in most cases has proven problematic. We discuss why QTL methods will have difficulty in identifying individual genes contributing to trait variation, and distinguish between the identification of QTL (or marker intervals) and the identification of individual genes or nucleotide differences within genes (QTN). This review focuses on what ecologists and evolutionary biologists working with natural populations can realistically expect to learn from QTL studies. We highlight representative issues in ecology and evolutionary biology and discuss the range of questions that can be addressed satisfactorily using QTL approaches. We specifically address developing approaches to QTL analysis in outbred populations, and discuss practical considerations of experimental (cross) design and application of different marker types. Throughout this review we attempt to provide a balanced description of the benefits of QTL methodology to studies in ecology and evolution as well as the inherent assumptions and limitations that may constrain its application.