Development of 101 novel EST-derived single nucleotide polymorphism markers for Zhikong scallop (Chlamys farreri)

Zhikong scallop(Chlamys farreri) is an important maricultured species in China.Many researches on this species,such as population genetics and QTL fine-mapping,need a large number of molecular markers.In this study,based on the expressed sequence tags(EST),a total of 300 putative single nucleotide polymorphisms(SNPs) were selected and validated using high resolution melting(HRM) technology with unlabeled probe.Of them,101(33.7%) were found to be polymorphic in 48 individuals from 4 populations.Further evaluation with 48 individuals from Qingdao population showed that all the polymorphic loci had two alleles with the minor allele frequency ranged from 0.046 to 0.500.The observed and expected heterozygosities ranged from 0.000 to 0.925 and from 0.089 to 0.505,respectively.Fifteen loci deviated significantly from Hardy-Weinberg equilibrium and significant linkage disequilibrate was detected in one pair of markers.BLASTx gave significant hits for 72 of the 101 polymorphic SNP-containing ESTs.Thirty four polymorphic SNP loci were predicted to be non-synonymous substitutions as they caused either the change of codons(33 SNPs) or pretermination of translation(1 SNP).The markers developed can be used for the population studies and genetic improvement on Zhikong scallop.     

Single nucleotide polymorphisms in intron 1 and intron 2 of Larimichthys crocea growth hormone gene are correlated with growth traits

The growth hormone gene (GH) affects animal growth and is a potential target for genetic studies of variation related to growth traits. In this study, we analyzed single nucleotide polymorphisms (SNPs) in GH intron regions and their associations with growth traits in large yellow croaker, Larimichthys crocea, from Zhejiang and Fujian stocks. The results of PCR-single strand conformation polymorphism showed two haplotypes of intron 1, named AA and AB genotypes, in Zhejiang stock. AB exhibited an SNP at position 196 (G A) that was negatively correlated with body height and positively correlated with standard length/body height (P 0.05). Two different genotypes, CC and CD, were identified in intron 2 in Fujian stock, with CD showing an SNP at position 692 (T C). The CD genotype had a significantly positive correlation with both weight and total length (P 0.01). These basic data highlight the potential for using GH as a genetic marker of fish growth in marker assisted selection.     

Influence of white spot syndrome virus infection on hepatopancreas gene expression of ‘Huanghai No. 2’ shrimp (<Emphasis Type="Italic">Fenneropenaeus chinensis</Emphasis>)

To elucidate the molecular response of shrimp hepatopancreas to white spot syndrome virus (WSSV) infection, microarray was applied to investigate the differentially expressed genes in the hepatopancreas of ‘Huanghai No. 2’ (Fenneropenaeus chinensis). A total of 59137 unigenes were designed onto a custom-made 60K Agilent chip. After infection, the gene expression profiles in the hepatopancreas of the shrimp with a lower viral load at early (48–96 h), peak (168–192 h) and late (264–288 h) infection phases were analyzed. Of 18704 differentially expressed genes, 6412 were annotated. In total, 5453 differentially expressed genes (1916 annotated) expressed at all three phases, and most of the annotated were either up- or down-regulated continuously. These genes function diversely in, for example, immune response, cytoskeletal system, signal transduction, stress resistance, protein synthesis and processing, metabolism among others. Some of the immune-related genes, including antilipopolysaccharide factor, Kazal-type proteinase inhibitor, C-type lectin and serine protease encoding genes, were up-regulated after WSSV infection. These genes have been reported to be involved in the anti-WSSV responses. The expression of genes related to the cytoskeletal system, including β-actin and myosin but without tubulin genes, were down-regulated after WSSV infection. Astakine was found for the first time in the WSSV-infected F. chinensis. To further confirm the expression of differentially expressed genes, quantitative real-time PCR was performed to test the expression of eight randomly selected genes and verified the reliability and accuracy of the microarray expression analysis. The data will provide valuable information to understanding the immune mechanism of shrimp’s response to WSSV.     

Screening white spot syndrome virus (WSSV)-resistant molecular markers from <Emphasis Type="Italic">Fenneropenaeus chinensis</Emphasis>

White spot syndrome virus (WSSV)-resistant molecular markers were screened from the selectively bred new variety ‘Huanghai No. 2’ of Fenneropenaeus chinensis using unlabeled-probe high-resolution melting (HRM) technique. After the artificial infection with WSSV, the first 96 dead shrimps and the last 96 surviving shrimps were collected, representing WSSV-susceptible and -resistant populations, respectively. The genotypes at well-developed 39 single nucleotide polymorphisms (SNPs) loci were obtained. As revealed in the Chi-square test, 3 SNPs, genotype A/A of contig C364-89AT, genotype A/A of C2635-527CA and genotype C/T of contig C12355-592CT, were positively correlated with disease-resistance traits. Other 2 SNPs, genotype G/G of contig C283-145AG and genotype C/C of contig C12355-592CT, were negatively correlated. Moreover, analysis with BlastX program for disease-resistant SNPs indicated that 3 contigs, Contig283, Contig364 and Contig12355, matched to the functional genes of effector caspase of Penaeus monodon, peptide transporter family 1-like protein, and 40S ribosomal protein S2 of Perca flavescens with high sequence similarity. The results will be helpful to provide theoretical and technical supports for molecular marker-assisted selective breeding of F. chinensis.     

Large-Scale Screening of Growth-Related Variants in Chinese Tongue Sole(Cynoglossus semilaevis)

Growth is the most valuable economic trait for improving aquaculture fish species, since fast growth can reduce labor cost and make more economic benefits. However, the knowledge about how many and which genes are related to the growth of Chinese tongue sole is limited. High-throughput sequencing screening of variants is a fast, economical and accurate assay to identify genes related to growth in crops, livestock and other aquaculture products. In the present study, genome-wide resequencing of 30 female Chinese tongue sole individuals from large and small groups to screen growth trait-related variations(SNPs, InDels) was carried out. In total, 6545735 SNPs and 1016745 InDels were detected, while 31 genes related to growth traits were identified. Their functions were mainly involved in muscle architecture, post-embryonic development, neurosensory development and hormone regulation. Furthermore, 18 of the 31 genes related to growth trait are located on W chromosome, indicating that W chromosome plays an important role in regulating the body size of female tongue soles. The markers and genes identified in our study can be applied to selective breeding of tongue sole and thus will promote the aquaculture industry and increase economic efficiency.     

A genetic linkage map of marine shrimp Penaeus (Fenneropenaeus) chinensis based on AFLP, SSR, and RAPD markers

The Chinese shrimp Penaeus (Fenneropaeneus) chinensis is an important species in marine fishery and aquaculture in China. A female Chinese shrimp Penaeus (Fenneropaeneus) chinensis was captured from west coast of the Korean peninsula and mated with a “Yellow Sea No. 1” male to produce the first filial generation (F₁) 100 F₂ full-sib progeny from brother-sister crosses between F₁ families was used for the mapping study. A genetic linkage map of the Chinese shrimp was constructed, based on 354 markers, including 300 amplified fragment length polymorphism (AFLP) markers, 42 microsatellite (SSR) markers, and 12 randomly amplified polymorphism (RAPD) markers. Forty-seven linkage groups (LGs) were identified. The total map length was 4 580.5 cM, with an average spacing of 11.3 cM, covering 75.8% of the estimated genome size. The construction of this genetic linkage map was part of a genetic breeding program. This linkage map will contribute to the discovery of genes and quantitative trait loci (QTLs) in Chinese shrimp.     

Relationship Between Single Nucleotide Polymorphism of Glycogen Synthase Gene of Pacific Oyster Crassostrea gigas and Its Glycogen Content

Glycogen is important not only for the energy supplementary of oysters, but also for human consumption. High glycogen content can improve the stress survival of oyster. A key enzyme in glycogenesis is glycogen synthase that is encoded by glycogen synthase gene GYS. In this study, the relationship between single nucleotide polymorphisms(SNPs) in coding regions of Crassostrea gigas GYS(Cg-GYS) and individual glycogen content was investigated with 321 individuals from five full-sib families. Single-strand conformation polymorphism(SSCP) procedure was combined with sequencing to confirm individual SNP genotypes of Cg-GYS. Least-square analysis of variance was performed to assess the relationship of variation in glycogen content of C. gigas with single SNP genotype and SNP haplotype. As a consequence, six SNPs were found in coding regions to be significantly associated with glycogen content(P 0.01), from which we constructed four main haplotypes due to linkage disequilibrium. Furthermore, the most effective haplotype H2(GAGGAT) had extremely significant relationship with high glycogen content(P 0.0001). These findings revealed the potential influence of Cg-GYS polymorphism on the glycogen content and provided molecular biological information for the selective breeding of good quality traits of C. gigas.     

Mutations in exons of the CYP17-II gene affect sex steroid concentration in male Japanese flounder (Paralichthys olivaceus)

As a specific gene of fish,cytochrome P450c17-Ⅱ(CYP17-Ⅱ) gene plays a key role in the growth,development and reproduction level of fish.In this study,the single-stranded conformational polymorphism(SSCP) technique was used to characterize polymorphisms within the coding region of CYP17-Ⅱ gene in a population of 75 male Japanese flounder(Paralichthys olivaceus).Three single nucleotide polymorphisms(SNPs) were identified in CYP17-Ⅱ gene of Japanese flounder.They were c.G594A(p.G188R),c.G939A and c.G1502A(p.G490D).SNP1(c.G594A),located in exon 4 of CYP17-Ⅱ gene,was significantly associated with gonadosomatic index(GSI).Individuals with genotype GG of SNP1 had significantly lower GSI(P < 0.05) than those with geno-type AA or AG.SNP2(c.G939A) located at the CpG island of CYP17-Ⅱ gene.The mutation changed the methylation of exon 6.Indi-viduals with genotype AA of SNP2 had significantly lower serum testosterone(T) level and hepatosomatic index(HSI) compared to those with genotype GG.The results suggested that SNP2 could influence the reproductive endocrine of male Japanese flounder.How-ever,the SNP3(c.G1502A) located in exon 9 did not affect the four measured reproductive traits.This study showed that CYP17-Ⅱ gene could be a potentially useful candidate gene for the research of genetic breeding and physiological aspects of Japanese flounder.     

Predictive ability of genomic selection models for breeding value estimation on growth traits of Pacific white shrimp Litopenaeus vannamei

Genomic selection(GS)can be used to accelerate genetic improvement by shortening the selection interval.The successful application of GS depends largely on the accuracy of the prediction of genomic estimated breeding value(GEBV).This study is a fi rst attempt to understand the practicality of GS in Litopenaeus vannamei and aims to evaluate models for GS on growth traits.The performance of GS models in L.vannamei was evaluated in a population consisting of 205 individuals,which were genotyped for 6 359 single nucleotide polymorphism(SNP)markers by specifi c length amplifi ed fragment sequencing(SLAF-seq)and phenotyped for body length and body weight.Three GS models(RR-BLUP,Bayes A,and Bayesian LASSO)were used to obtain the GEBV,and their predictive ability was assessed by the reliability of the GEBV and the bias of the predicted phenotypes.The mean reliability of the GEBVs for body length and body weight predicted by the dif ferent models was 0.296 and 0.411,respectively.For each trait,the performances of the three models were very similar to each other with respect to predictability.The regression coeffi cients estimated by the three models were close to one,suggesting near to zero bias for the predictions.Therefore,when GS was applied in a L.vannamei population for the studied scenarios,all three models appeared practicable.Further analyses suggested that improved estimation of the genomic prediction could be realized by increasing the size of the training population as well as the density of SNPs.     

Genome-wide mining,characterization, and development of microsatellite markers in <Emphasis Type="Italic">Marsupenaeus japonicus</Emphasis> by genome survey sequencing

The kuruma prawn, Marsupenaeus japonicus, is one of the most cultivated and consumed species of shrimp. However, very few molecular genetic/genomic resources are publically available for it. Thus, the characterization and distribution of simple sequence repeats (SSRs) remains ambiguous and the use of SSR markers in genomic studies and marker-assisted selection is limited. The goal of this study is to characterize and develop genome-wide SSR markers in M. japonicus by genome survey sequencing for application in comparative genomics and breeding. A total of 326 945 perfect SSRs were identified, among which dinucleotide repeats were the most frequent class (44.08%), followed by mononucleotides (29.67%), trinucleotides (18.96%), tetranucleotides (5.66%), hexanucleotides (1.07%), and pentanucleotides (0.56%). In total, 151 541 SSR loci primers were successfully designed. A subset of 30 SSR primer pairs were synthesized and tested in 42 individuals from a wild population, of which 27 loci (90.0%) were successfully amplified with specific products and 24 (80.0%) were polymorphic. For the amplified polymorphic loci, the alleles ranged from 5 to 17 (with an average of 9.63), and the average PIC value was 0.796. A total of 58 256 SSR-containing sequences had significant Gene Ontology annotation; these are good functional molecular marker candidates for association studies and comparative genomic analysis. The newly identified SSRs significantly contribute to the M. japonicus genomic resources and will facilitate a number of genetic and genomic studies, including high density linkage mapping, genome-wide association analysis, marker-aided selection, comparative genomics analysis, population genetics, and evolution.     

Performance comparison of two efficient genomic selection methods (gsbay &amp; MixP) applied in aquacultural organisms

Genomic selection is more and more popular in animal and plant breeding industries all around the world, as it can be applied early in life without impacting selection candidates. The objective of this study was to bring the advantages of genomic selection to scallop breeding. Two different genomic selection tools MixP and gsbay were applied on genomic evaluation of simulated data and Zhikong scallop (Chlamys farreri) field data. The data were compared with genomic best linear unbiased prediction (GBLUP) method which has been applied widely. Our results showed that both MixP and gsbay could accurately estimate single-nucleotide polymorphism (SNP) marker effects, and thereby could be applied for the analysis of genomic estimated breeding values (GEBV). In simulated data from different scenarios, the accuracy of GEBV acquired was ranged from 0.20 to 0.78 by MixP; it was ranged from 0.21 to 0.67 by gsbay; and it was ranged from 0.21 to 0.61 by GBLUP. Estimations made by MixP and gsbay were expected to be more reliable than those estimated by GBLUP. Predictions made by gsbay were more robust, while with MixP the computation is much faster, especially in dealing with large-scale data. These results suggested that both algorithms implemented by MixP and gsbay are feasible to carry out genomic selection in scallop breeding, and more genotype data will be necessary to produce genomic estimated breeding values with a higher accuracy for the industry.     

Genetic variation assessed with microsatellites in mass selection lines of the Pacific oyster (<Emphasis Type="Italic">Crassostrea gigas</Emphasis>) in China

Four successive mass selection lines of the Pacific oyster, Crassostrea gigas, selected for faster growth in breeding programs in China were examined at ten polymorphic microsatellite loci to assess the level of allelic diversity and estimate the effective population size. These data were compared with those of their base population. The results showed that the genetic variation of the four generations were maintained at high levels with an average allelic richness of 18.8–20.6, and a mean expected heterozygosity of 0.902–0.921. They were not reduced compared with those of their base population. Estimated effective population sizes based on temporal variances in microsatellite frequencies were smaller to that of sex ratio-corrected broodstock count estimates. Using a relatively large number of broodstock and keeping an equal sex ratio in the broodstock each generation may have contributed to retaining the original genetic diversity and maintaining relatively large effective population size. The results obtained in this study showed that the genetic variation was not affected greatly by mass selection progress and high genetic variation still existed in the mass selection lines, suggesting that there is still potential for increasing the gains in future generations of C. gigas. The present study provided important information for future genetic improvement by selective breeding, and for the design of suitable management guidelines for genetic breeding of C. gigas.