The regulation of prolactin secretion and its targeting function of teleost

Prolactin is involved in regulating various physiological activities of vertebrates and is one of the most momentous pituitary hormones. However, not enough attention is currently paid to prolactin, especially in teleost. This paper aims to gather, organize, and analyze recent studies on the regulation and functions of prolactin. By comparing with other animal groups, it highlights the significant role of prolactin in fish reproduction, immunity, growth, and osmotic pressure regulation, as well as the upstream and downstream factors that may be involved in the regulation of prolactin functions were introduced to provide a theoretical basis for the in-depth study and potential practical application of prolactin.

Genome-Wide Marker Analysis for Traits of Economic Importance in Asian Seabass Lates calcarifer

To date, it is not known whether animal breeding values in Asian seabass (Lates calcarifer) can be estimated using single nucleotide polymorphisms (SNPs) generated from new high-throughput genotyping by sequencing platforms. The principal aim of the present study was to assess the genomic prediction accuracy for growth traits, survival, cannibalism, and disease resistance against Streptococcus iniae in this species L. calcarifer. Additionally, this study attempted to identify markers associated with the five traits studied as well as to understand if the genotype data can be used to estimate genetic parameters for these complex traits. The genomic best linear unbiased prediction (gBLUP) method was used to analyze 11,084 SNPs and showed that the prediction accuracies for growth traits (weight and length) were high (0.67–0.75). By contrast, these estimates for survival were low (0.25). Multi-locus mixed model analyses identified four SNPs significantly associated with body weight (p < 5 × 10−8 or −log10 p ≥ 5). There were, however, no significant associations detected for other traits. Similarly, the SNP heritability was moderate, while the estimates for other traits were approximated to zero and not significant. Genetic correlations between body weight and standard length were close to unity. Collectively, the results obtained from this study suggest that genotyping by sequencing platforms can provide informative DNA markers to conduct genome-wide association analysis, estimation of genetic parameters, and evaluation of genomic prediction accuracy for complex traits in Asian seabass.

Polymorphisms of the growth hormone gene and their association with growth traits and sex in Sarcocheilichthys sinensis

The growth hormone gene (gh) of Sarcocheilichthys sinensis was cloned and characterized in this study. The cDNA length of gh was 973 bp, containing a 5'-UTR of 15 bp, a 3'-UTR of 325 bp and an open reading frame of 633 bp. The genomic DNA of gh was 2135 bp in length containing five exons and four introns. The precursor peptide of gh contained 210 amino acids (aa), including a signal peptide of 22 aa (Met¹-Ala²²) and a mature region of 188 aa (Ser²³-Leu²¹⁰). The similarity and identity ranges of the gh precursor peptide with those of other cyprinids were 88.6%-99.0% and 84.8%-98.6%, respectively. The gh of S. sinensis expressed at the highest level in the pituitary, and its expression was also detected in muscle and brain. Six polymorphic sites were detected in intron 1 (g.51InDel, g.64InDel and g.242InDel), intron 2 (g.864T>C), intron 3 (g.1017InDel) and intron 4 (g.1541A>G). Among these sites, g.242InDel was significantly associated with condition factor, g.1541A>G was associated with all six growth traits, while g.864T>C was associated with sex. The data obtained herein provide useful information for further studies on the regulation mechanisms of growth and sexual growth differences in S. sinensis.

A novel 23 bp indel mutation in PRL gene is associated with growth traits in Luxi Blackhead sheep

Prolactin is a highly versatile pituitary hormone with multiple biological functions. PRL expression is regulated by POU1F1 and the prophet of POU1F1 (PROP1). The aim of this study was to investigate the indel variations in ovine PRL and the directly related (PROP1 and POU1F1) genes, and their associations with growth traits in Luxi Blackhead (LXBH) sheep. A monomorphism in PROP1 and POU1F1 genes, and one novel 23-bp insertion mutation in the PRL gene, were identified in LXBH sheep. The 23 bp insertion mutation within PRL gene was significantly associated with several body measurements (e.g., body weight, body height) in sheep of different ages (p < 0.05). Ram lambs (p = 0.036) of genotype insertion/insertion (II) had significantly higher body weights. Weaners (p = 0.018) of genotypes insertion/insertion (II) and insertion/deletion (ID) also had significantly higher body weights compared with male sheep of deletion/deletion (DD) genotype. Moreover, among ewe lambs, individuals of genotype insertion/insertion (II) had a higher paunch girth compared to those with other genotypes (p = 0.044). These findings indicate that a 23 bp indel variant of the ovine PRL gene is correlated with body measurements in LXBH sheep. The findings have potential utility for sheep breeding programs based on marker-assisted selection.

Whole genome sequencing reveals the impact of recent artificial selection on red sea bream reared in fish farms

Red sea bream, a popular fish resource in Korea and Japan, is being bred in fish farms of the two countries. It is hypothesized that the genomes of red sea bream are influenced by decades of artificial selection. This study investigates the impact of artificial selection on genomes of red sea bream. Whole genome sequencing was conducted for 40 samples of red sea bream either from Ehime, Nagasaki and Tongyeong fish farms or from the wild. Population stratification based on whole genome data was investigated and the genomic regions of fish farm populations under selection were identified using XP-EHH and relative nucleotide diversity. Gene ontology analysis revealed that different functions were enriched in different fish farms. In conclusion, this study highlights the difference between independently cultured red sea bream populations by showing that influence of artificial selection acted upon completely different genes related to different functions including metabolic and developmental processes.

Genome-Wide Association Study Identifies Loci Associated with Resistance to Viral Nervous Necrosis Disease in Asian Seabass

Viral nervous necrosis disease (VNN), caused by nervous necrosis virus (NNV), is one major threat to mariculture. Identifying loci and understanding the mechanisms associated with resistance to VNN are important in selective breeding programs. We performed a genome-wide association study (GWAS) using genotyping-by-sequencing (GBS) to study the genomic architecture of resistance to NNV infection in Asian seabass. We genotyped 986 individuals from 43 families produced by 15 founders with 44498 bi-allelic genetic variants using GBS. The GWAS identified three genome-wide significant loci on chromosomes 16, 19, and 20, respectively, and six suggestive loci on chromosomes 1, 8, 14, 15, 21, and 24, respectively, associated with resistance to NNV infection measured as binary and quantitative traits. Using the 500 most significant markers in combination with a training population of 800 samples could reach a genomic prediction accuracy of 0.7. Candidate genes significantly associated with resistance to NNV, including lysine-specific demethylase 2A, beta-defensin 1, and cystatin-B, which play important roles in immune responses against virus infection, were identified. Almost all the candidate genes were differentially expressed in different tissues against NNV infection. The significant genetic variants can be used in genomic selection and help understand the mechanism of resistance to VNN. Future studies should use populations of large effective size and whole genome resequencing to identify more useful genetic variants.

BAC-pool sequencing and analysis confirms growth-associated QTLs in the Asian seabass genome

The Asian seabass is an important marine food fish that has been cultured for several decades in Asia Pacific. However, the lack of a high quality reference genome has hampered efforts to improve its selective breeding. A 3D BAC pool set generated in this study was screened using 22 SSR markers located on linkage group 2 which contains a growth-related QTL region. Seventy-two clones corresponding to 22 FPC contigs were sequenced by Illumina MiSeq technology. We co-assembled the MiSeq-derived scaffolds from each FPC contig with error-corrected PacBio reads, resulting in 187 sequences covering 9.7 Mb. Eleven genes annotated within this region were found to be potentially associated with growth and their tissue-specific expression was investigated. Correlation analysis demonstrated that SNPs in ctsb, skp1 and ppp2ca can be potentially used as markers for selecting fast-growing fingerlings. Conserved syntenies between seabass LG2 and five other teleosts were identified. This study i) provided a 10 Mb targeted genome assembly; ii) demonstrated NGS of BAC pools as a potential approach for mining candidates underlying QTLs of this species; iii) detected eleven genes potentially responsible for growth in the QTL region; and iv) identified useful SNP markers for selective breeding programs of Asian seabass.

Transcriptome Analysis of Blunt Snout Bream (Megalobrama amblycephala) Reveals Putative Differential Expression Genes Related to Growth and Hypoxia

The blunt snout bream (Megalobrama amblycephala) is an important freshwater aquaculture species, but it is sensitive to hypoxia. No transcriptome data related to growth and hypoxia response are available for this species. In this study, we performed de novo transcriptome sequencing for the liver and gills of the fast-growth family and slow-growth family derived from ‘Pujiang No.1’ F10 blunt snout bream that were under hypoxic stress and normoxia, respectively. The fish were divided into the following 4 groups: fast-growth family under hypoxic stress, FH; slow-growth family under hypoxic stress, SH; fast-growth family under normoxia, FN; and slow-growth family under normoxia, SN. A total of 185 million high-quality reads were obtained from the normalized cDNA of the pooled samples, which were assembled into 465,582 contigs and 237,172 transcripts. A total of 31,338 transcripts from the same locus (unigenes) were annotated and assigned to 104 functional groups, and 23,103 unigenes were classified into seven main categories, including 45 secondary KEGG pathways. A total of 22,255 (71%) known putative unigenes were found to be shared across the genomes of five model fish species and mammals, and a substantial number (9.4%) of potentially novel genes were identified. When 6,639 unigenes were used in the analysis of differential expression (DE) genes, the number of putative DE genes related to growth pathways in FH, SH, SN and FN was 159, 118, 92 and 65 in both the liver and gills, respectively, and the number of DE genes related to hypoxic response was 57, 33, 23 and 21 in FH, FN, SH and SN, respectively. Our results suggest that growth performance of the fast-growth family should be due to complex mutual gene regulatory mechanisms of these putative DE genes between growth and hypoxia.

Signatures of selection in tilapia revealed by whole genome resequencing

Natural selection and selective breeding for genetic improvement have left detectable signatures within the genome of a species. Identification of selection signatures is important in evolutionary biology and for detecting genes that facilitate to accelerate genetic improvement. However, selection signatures, including artificial selection and natural selection, have only been identified at the whole genome level in several genetically improved fish species. Tilapia is one of the most important genetically improved fish species in the world. Using next-generation sequencing, we sequenced the genomes of 47 tilapia individuals. We identified a total of 1.43 million high-quality SNPs and found that the LD block sizes ranged from 10–100 kb in tilapia. We detected over a hundred putative selective sweep regions in each line of tilapia. Most selection signatures were located in non-coding regions of the tilapia genome. The Wnt signaling, gonadotropin-releasing hormone receptor and integrin signaling pathways were under positive selection in all improved tilapia lines. Our study provides a genome-wide map of genetic variation and selection footprints in tilapia, which could be important for genetic studies and accelerating genetic improvement of tilapia.

Characterization of the LECT2 gene and its associations with resistance to the big belly disease in Asian seabass

Leukocyte cell-derived chemotaxin-2 (LECT2) is an important protein of the innate immune system for the defense against bacterial infection. We cloned and characterized the LECT2 gene from Asian seabass (Lates calcarifer). Its complete cDNA consisted of an open reading frame of 459 bp encoding a protein of 152 amino acids. The genomic DNA sequence of this gene consists of four exons and three introns. Quantitative real-time PCR revealed that the LECT2 gene was expressed predominantly in liver while its expression was moderate in spleen and heart, and weak in other tissues. The LECT2 transcript was up-regulated in the kidney, spleen and liver in response to a challenge with a pathogenic bacterium Vibrio harveyi. In addition, we identified three single nucleotide polymorphisms (SNPs) in the LECT2 gene, and found significant associations between these polymorphisms and resistance to the big belly disease. These results suggest that the LECT2 gene play an important role in resistance to bacterial pathogens in fish. The SNP markers in the gene associated with the resistance to bacterial pathogens may facilitate selecting Asian seabass resistant to bacterial diseases.

Molecular characterization and mapping of Fgf21 gene in a foodfish species asian seabass

Fgf21 is a newly discovered fibroblast growth factor. It is typically induced by fasting and plays important roles in the regulation of glucose and lipid metabolisms and energy balance in mammals, whereas potential functions of this gene in teleosts are still unknown. We identified the Fgf21 gene and studied its functions in Asian seabass (Lates calcarifer). The cDNA of the Fgf21 encoded a protein with 206 amino acids. Analysis of DNA and amino acid sequences of Fgf21 genes revealed that the sequences and structure of the Fgf21 genes were highly conserved in vertebrates. Real-time PCR revealed that Fgf21 was exclusively expressed in the intestine and kidney, which was different from the expression profiles of mammals. Fgf21 was down-regulated under fasting, whereas it was significantly increased during the LPS challenge. Exogenous recombinant FGF21 significantly suppressed the appetite of Asian seabass. Our data suggest that Fgf21 plays a role in energy regulation and acute phase response in Asian seabass, and may have different functions in fish and mammals. In addition, we identified one SNP in Fgf21. By using this SNP, the gene was mapped on the linkage group 23, where a suggestive QTL for growth was mapped previously. Association mapping identified significant associations between Fgf21 genotypes at the SNP and growth traits. These results not only provide important information of the functions of Fgf21, but also suggest that the SNP in this gene can be used as a marker in selecting fast-growing individuals of Asian seabass.

Analysis of two lysozyme genes and antimicrobial functions of their recombinant proteins in Asian seabass

Lysozymes are important proteins of the innate immune system for the defense against bacterial infection. We cloned and analyzed chicken-type (c-type) and goose-type (g-type) lysozymes from Asian seabass (Lates calcarifer). The deduced amino acid sequence of the c-type lysozyme contained 144 residues and possessed typical structure residues, conserved catalytic residues (Glu(50) and Asp(67)) and a "GSTDYGIFQINS" motif. The deduced g-type lysozyme contained 187 residues and possessed a goose egg white lysozyme (GEWL) domain containing three conserved catalytic residues (Glu(71), Asp(84), Asp(95)) essential for catalytic activity. Real time quantitative PCR (qRT-PCR) revealed that the two lysozyme genes were constitutively expressed in all the examined tissues. The c-type lysozyme was most abundant in liver, while the g-type lysozyme was predominantly expressed in intestine and weakly expressed in muscle. The c-type and g-type transcripts were up-regulated in the kidney, spleen and liver in response to a challenge with Vibrio harveyi. The up-regulation of the c-type lysozyme was much stronger than that of the g-type lysozyme in kidney and spleen. The recombinant proteins of the c-type and g-type lysozymes showed lytic activities against the bacterial pathogens Vibrio harveyi and Photobacterium damselae in a dosage-dependent manner. We identified single nucleotide polymorphisms (SNPs) in the two lysozyme genes. There were significant associations of these polymorphisms with resistance to the big belly disease. These results suggest that the c- and g-type genes play an important role in resistance to bacterial pathogens in fish. The SNP markers in the two genes associated with the resistance to bacterial pathogens may facilitate the selection of Asian seabass resistant to bacterial diseases.

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.

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.

Microsatellite length variation in candidate genes correlates with habitat in the gilthead sea bream Sparus aurata

The genetic basis and evolutionary implications of local adaptation in high gene flow marine organisms are still poorly understood. In several Mediterranean fish species, alternative migration patterns exist between individuals entering coastal lagoons that offer favourable conditions for growth and those staying in the sea where environmental conditions are less subject to rapid and stressful change. Whether these coexisting strategies are phenotypically plastic or include a role for local adaptation through differential survival needs to be determined. Here, we explore the genetic basis of alternate habitat use in western Mediterranean populations of the gilthead sea bream (Sparus aurata). Samples from lagoonal and open-sea habitats were typed for three candidate gene microsatellite loci, seven anonymous microsatellites and 44 amplified fragment length polymorphism markers to test for genotype-environment associations. While anonymous markers globally indicated high levels of gene flow across geographic locations and habitats, non-neutral differentiation patterns correlated with habitat type were found at two candidate microsatellite loci located in the promoter region of the growth hormone and prolactin genes. Further analysis of these two genes revealed that a mechanism based on habitat choice alone could not explain the distribution of genotype frequencies at a regional scale, thus implying a role for differential survival between habitats. We also found an association between allele size and habitat type, which, in the light of previous studies, suggests that polymorphisms in the proximal promoter region could influence gene expression by modulating transcription factor binding, thus providing a potential explanatory link between genotype and growth phenotype in nature.