EST analysis and identification of gonad-related genes from the normalized cDNA library of large yellow croaker, Larimichthys crocea

Genome-Wide Association Study Identifies Genomic Loci of Sex Determination and Gonadosomatic Index Traits in Large Yellow Croaker (Larimichthys crocea)

Larimichthys crocea is one of the traditional marine culture fishes in China, widely distributed in South China Sea, East Sea, and southern Yellow Sea. Sex dimorphism is evident in this species that females present a substantial growth strength than males, suggesting breeding females could obtain more economic benefits in L. crocea aquaculture industry. With the continuous expansion of aquaculture industry, both identifying sex-associated genome region and understanding the genetic basis underlying gonad differentiation and development matter to not only sex control aquaculture but also breeding industry. Thus, genome-wide association analysis (GWAS) of sex determination was conducted with a random breeding population of 905 individuals (including 463 females and 442 males) by ddRAD sequencing. For sex determination, 21 significant single nucleotide polymorphisms (SNPs) in chromosome (Chr) 22 were identified. Surrounding these SNPs, we founded 14 candidate genes, including dmrt1, dmrt3, and piwil2, fam102a, and odf2. The sex-associated region was narrowed down further to 2.4 Mb on Chr22 through F_st scanning and insertion-deletion (InDel) analysis. Besides, 3 SNPs in the supposed sex-determining region on Chr22 were identified as highly associated with gonad differentiation through GWAS on gonadosomatic index (GSI) in 350 males and 231 females. Because of the significant difference of GSI between females and males of L. crocea, GWAS on GSI of different genders was also conducted independently. Finally, we identified a SNP in Chr18 showing genome-wide significant association with male GSI (MGSI) and three genes axl, cyp2a10, and cyp2g1 involved in the gonadal development regulation process of aromatase. Overall, this study explored the genetic basis of sex determination mechanism and provided novel insights into gonad differentiation and development, offering solid genetic support for sex control breeding, marker-assisted selection, and marine resources conservation.

Molecular cloning, characterization and expression of Lc-Sox11a in large yellow croaker Larimichthys crocea

Sox genes play important roles in various developmental processes such as sex determination, embryogenesis, oogenesis, neurogenesis, and larval development. In order to clarify the roles of Sox genes in the developmental process of large yellow croaker, the full-length cDNA of the Sox11a gene (Lc-Sox11a) was cloned for the first time. Bioinformatics analysis indicated that Lc-Sox11a contains a protein of 366 amino acids with a Ser-rich region, a C-terminal conserved region, and a high mobility group box. The expression of Lc-Sox11a in different tissues of both sexes and in different developmental embryonic stages revealed that Lc-Sox11a were expressed with tissue and gender specificity, of which the expression level in female was ovary>brain>eye>gill; in male was brain>testis>gill. The gender differences occurred in the brain and eye with the male brain>female brain, female eye>male eye. Moreover, the expression of Lc-Sox11a in the gonad and brain at different growth stages was detected. Significant up-regulated expression of Lc-Sox11a was found in the ovary and the male brain at 1000dph (days post hatching) compared with 270dph and 635dph. However, significant down-regulated expression of Lc-Sox11a occurred in the testis with growth. Besides, the expression of Lc-Sox11a in the female brain showed a trend of first rising then falling, with the highest peak in 635dph. The results of in situ hybridization displayed that Lc-Sox11a was widely distributed only in cytoplasm of oocytes at each stage in oogenesis. In early stage of oocytes, Lc-Sox11a was expressed weakly and evenly. As the appearance of vacuoles and synthesis of yolks, positive signals of Lc-Sox11a distributed intensively in the residual cytoplasm. In spermatogenesis, Lc-Sox11a was distributed in cytoplasm of all male germ cells except spermatozoon with spermatogonium>spermatocyte>spermatid. During embryogenesis, Lc-Sox11a was expressed in most embryonic stages, the highest expression occurred in the formation-of-eye-lens stage, closely followed by the closure-of-blastopore stage, then the beginning-of-heart-pulsation stage. The results of whole mount in situ hybridization showed that the expression of Lc-Sox11a began to increase beginning with the multiple-cell stage, with the major distribution of Lc-Sox11a in the brain and eye areas in the pre-hatching stage.

Functional marker detection and analysis on a comprehensive transcriptome of large yellow croaker by next generation sequencing

Large yellow croaker (Larimichthys crocea) is an important economic fish in China and Eastern Asia. Because of the exhaustive fishing and overdense aquaculture, the wild population and the mariculture of the species are facing serious challenges on germplasm degeneration and susceptibility to infectious disease agents. However, a comprehensive transcriptome from multi-tissues of the species has not been reported and functional molecular markers have not yet been detected and analyzed. In this work, we applied RNA-seq with the Illumina Hiseq2000 platform for a multi-tissue sample of large yellow croaker and assembled the transcriptome into 88,103 transcripts. Of them, 52,782 transcripts have been successfully annotated by nt/nr, InterPro, GO and KEGG database. Comparing with public fish proteins, we have found that 34,576 protein coding transcripts are shared in large yellow croaker with zebrafish, medaka, pufferfish, and stickleback. For functional markers, we have discovered 1,276 polymorphic SSRs and 261, 000 SNPs. The functional impact analysis of SNPs showed that the majority (~75%) of small variants cause synonymous mutations in proteins, followed by variations in 3' UTR region. The functional enrichment analysis illuminated that transcripts involved in DNA bindings, enzyme activities, and signal pathways prominently exhibit less single-nucleotide variants but genes for the constituent of the muscular tissue, the cytoskeleton, and the immunity system contain more frequent SNP mutations, which may reflect the structural and functional selections of the translated proteins. This is the first work for the high-throughput detection and analysis of functional polymorphic SSR and SNP markers in a comprehensive transcriptome of large yellow croaker. Our study provides valuable transcript sequence and functional marker resources for the quantitative trait locus (QTL) identification and molecular selection of the species in the research community.

Liver proteomic analysis of the large yellow croaker (Pseudosciaena crocea) following polyriboinosinic:polyribocytidylic acid induction

In the present study, we examined the liver protein profiles of the large yellow croaker (Pseudosciaena crocea) exposed to polyriboinosinic:polyribocytidylic acid [poly(I:C)], a viral mimic, using the differential proteomic approach. Sixteen altered protein spots were identified by matrix-assisted laser desorption ionization time of flight mass spectrometry or matrix-assisted laser desorption ionization time of flight/time of flight mass spectrometry, including eight upregulated proteins and eight downregulated proteins. These altered host proteins were classified into six categories based on their biological function: cellular process, metabolic process, biological regulation, binding, and catabolic process, highlighting the fact that response to poly(I:C) induction in fish seems to be complex and diverse. Moreover, four corresponding genes of the differentially expressed proteins were validated by relative quantitative real-time PCR. Western blot analysis further demonstrated the changes in protein abundance of natural killer enhancing factor and peroxiredoxin 6. These results will be helpful in furthering our understanding of the changes of physiological processes in liver of fish during virus infection.

Molecular cloning, characterization, and gene expression of the androgen receptor in the large yellow croaker, Larimichthys crocea

Androgens mediate a wide range of physiological responses and developmental processes in vertebrates, involving both reproductive and nonreproductive systems. The activity of androgens is mediated by the androgen receptor (AR), a member of the nuclear receptor superfamily. In this study, an AR gene was cloned from the large yellow croaker (Larimichthys crocea) for the first time. qRT-PCR revealed ubiquitous expression of AR in all adult tissues examined, with higher expression in the gonad and liver of both sexes and highest expression in the blastula stage of embryonic development. Using in situ hybridization, we detected positive signals of AR in the spermatogonium, spermatocyte, spermatid, and spermatozoon during spermatogenesis, in the cytoplasm of all oocytes during oogenesis and in the follicle cells of stage IV oocytes. Our findings support the important role that AR plays in gametogenesis, gonadal development, and the early stages of embryonic development.

Multiple isoforms of immune-related genes from hemocytes and eyestalk cDNA libraries of swimming crab Portunus trituberculatus

Expressed sequence tags (ESTs) analysis has been shown to be an efficient approach not only for gene discovery, but also for gene expression profiles performance. Two full-length enriched cDNA libraries were constructed from hemocytes and eyestalk of Portunus trituberculatus, respectively, and randomly sequenced to collect genomic information and identify genes involved in immune defense response. A total of 99 unigenes including 64 unigenes (6.00% of 1066 unigenes) in hemocytes library and 35 unigens (6.86% of 510 unigenes) in eyestalk library are identified to be immune genes. These genes are categorized into six classes, viz. antimicrobial peptides, redox proteins, melanization related proteins, chaperone proteins, clottable proteins and other immune factors. The content and category of immune genes in eyestalk library indicate eyestalk might have unrecognized role in crab immunity. Five immune genes containing multiple protein isoforms are identified and characterized, including anti-lipopolysaccharide factor (PtALF1-7), crustin (PtCrustin1-3), thioredoxin (PtTrx1-2), clip domain serine proteinase (PtcSP1-5) and kazal-type proteinase inhibitor (PtKPI1-4). Sequence alignment and phylogenetic analysis reveal PtALF1-7 contain two conserved cysteine residues and might be encoded by multiple genomic loci. PtCrustin1-3 share the consensus cysteine motif and are considered as Type I crustins. PtTrx1 possesses the critical structural cysteine residue C⁷³ of Trx-1, while PtTrx2 has the N-terminal mitochondrial translocation signal of Trx-2. Sequence analysis shows PtcSP1-5 contain one clip domain and one partial SP catalytic triad domain. PtKPI1-4 present one typical Kazal domain consisting of six conserved cysteine residues. Some protein isoforms are tissue-specific, which might suggest they have different origins and perform diverse functions. Except PtALF1-3 and PtCrustin1, the other isoformes in this study are firstly identified from P. trituberculatus. Especially, PtTrx2 are firstly identified from crustaceans. Our research will provide useful genomic information of P. trituberculatus and be helpful in understanding the molecular mechanisms of crab immunity.