Whole Genome Sequencing of the Giant Grouper (Epinephelus lanceolatus) and High-Throughput Screening of Putative Antimicrobial Peptide Genes

Giant groupers, the largest grouper type in the world, are of economic importance in marine aquaculture for their rapid growth. At the same time, bacterial and viral diseases have become the main threats to the grouper industry. Here, we report a high-quality genome of a giant grouper sequenced by an Illumina HiSeq X-Ten and PacBio Bioscience Sequel platform. A total of 254 putative antimicrobial peptide (AMP) genes were identified, which can be divided into 34 classes according to the annotation of the Antimicrobial Peptides Database (APD3). Their locations in pseudochromosomes were also determined. Thrombin-, lectin-, and scolopendin-derived putative AMPs were the three largest parts. In addition, expressions of putative AMPs were measured by our transcriptome data. Two putative AMP genes (gapdh1 and gapdh2) were involved in glycolysis, which had extremely high expression levels in giant grouper muscle. As it has been reported that AMPs inhibit the growth of a broad spectrum of microbes and participate in regulating innate and adaptive immune responses, genome sequencing of this study provides a comprehensive cataloging of putative AMPs of groupers, supporting antimicrobial research and aquaculture therapy. These genomic resources will be beneficial to further molecular breeding of this economically important fish

Chromosome-Level Genome Assembly and Transcriptome Comparison Analysis of Cephalopholis sonnerati and Its Related Grouper Species

The tomato hind, Cephalopholis sonnerati, is a bottom-dwelling coral reef fish, which is widely distributed in the Indo-Pacific and Red Sea. C. sonnerati also features complex social structures and behaviour mechanisms. Here, we present a high-quality, chromosome-level genome assembly for C. sonnerati that was derived using PacBio sequencing and Hi-C technologies. A 1043.66 Mb genome with an N50 length of 2.49 Mb was assembled, produced containing 795 contigs assembled into 24 chromosomes. Overall, 97.2% of the complete BUSCOs were identified in the genome. A total of 26,130 protein-coding genes were predicted, of which 94.26% were functionally annotated. Evolutionary analysis revealed that C. sonnerati diverged from its common ancestor with E. lanceolatus and E. akaara approximately 41.7 million years ago. In addition, comparative genome analyses indicated that the expanded gene families were highly enriched in the sensory system. Finally, we found the tissue-specific expression of 8108 genes. We found that these tissue-specific genes were highly enriched in the brain. In brief, the high-quality, chromosome-level reference genome will provide a valuable genome resource for studies of the genetic conservation, resistance breeding, and evolution of C. sonnerati

An SNP-Based Genetic Map and QTL Mapping for Growth Traits in the Red-Spotted Grouper (Epinephelus akaara)

The red-spotted grouper (Epinephelus akaara) is one of the most commercially important aquatic species in China. However, its seedstock has low larval survival rates, and its stability is confronted with the danger of overexploitation. In this study, a high-density genetic map was constructed using 3435 single nucleotide polymorphisms (SNPs) from 142 first generation (F1) full-sib offspring and two parents of a red-spotted grouper population. The total genetic length of the map was 2300.12 cM with an average intermarker distance of 0.67 cM. Seventeen genome-wide significant quantitative trait loci (QTLs) for growth-related traits were detected on 24 linkage groups, including 5 QTLs for full length, 7 QTLs for body length, and 5 QTLs for body weight. The contribution values of explained phenotypic variance ranged from 10.7% to 12.9%. Moreover, 13 potential candidate genes for growth-related traits were identified. Collectively, these findings will be useful for conducting marker-assisted selection of the red-spotted grouper in future studies

Transcriptome analysis of the Trachinotus ovatus: identification of reproduction, growth and immune-related genes and microsatellite markers.

The Trachinotus ovatus (Teleostei, Carangidae) is an economically important marine fish species in the world. However, the lack of genomic information regarding this species limits our understanding of the genetics and biological mechanisms in Trachinotus ovatus. In this study, high throughput transcriptome sequencing was used to obtain comprehensive genomic information in Trachinotus ovatus.Transcriptome sequencing was performed by using Illumina paired-end sequencing technology. The 98,534,862 high quality reads were yielded, and were de novo assembled into 156,094 unigenes with an average sequence length of 1179 bp. Transcriptome annotation revealed that 75,586 and 67,923 unigenes were functionally annotated in the NCBI non-redundant database and Swiss-Prot protein database, respectively. Functional analysis demonstrated that 67,923 unigenes were grouped into 25 Cluster of Orthologous Groups (COG) functional categories, 37,976 unigenes were clustered into 61 Gene Ontology (GO) terms, and 38,172 unigenes were assigned to 275 different Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Based on the transcriptome dataset, a large number of unigenes associated with reproduction, growth and immunity were identified. Furthermore, a total number of 38,794 simple sequence repeats (SSRs) were discovered and 16 polymorphic loci were characterized in Trachinotus ovatus.The present study is the first transcriptome analysis of a fish species belonging to the genus Trachinotus and provides a valuable genomic resource for novel gene discovery, gene expression and regulation studies, and the identification of genetic markers in Trachinotus ovatus and the other fish of the genus Trachinotus

MT-Feeding-Induced Impermanent Sex Reversal in the Orange-Spotted Grouper during Sex Differentiation

In this study, we systematically investigated the process of sex reversal induced by 17-methyltestosterone (MT) feeding and MT-feeding withdrawal at the ovary differentiation stage in orange-spotted groupers, Epinephelus coioides. Gonadal histology showed that MT feeding induced a precocious sex reversal from immature ovaries to testes, bypassing the formation of an ovarian cavity, and MT-feeding withdrawal led to an ovarian fate. In both the MT feeding and MT-feeding withdrawal phases, cytochrome P450 family 11 subfamily B (cyp11b) gene expression and serum 11-KT levels were not significantly changed, suggesting that the MT-treated fish did not generate endogenous steroids, even though active spermatogenesis occurred. Finally, by tracing doublesex-expressing and Mab-3-related transcription factor 1 (dmrt1)-expressing cells and TUNEL (terminal deoxynucleotidyl transferase 2-deoxyuridine, 5-triphosphate nick end labeling) assays, we found that the efferent duct formed first, and then, the germ cells and somatic cells of the testicular tissue were generated around the efferent duct during MT-feeding-induced precocious sex reversal. Collectively, our findings provide insights into the molecular and cellular mechanisms underlying sex reversal induced by exogenous hormones during sex differentiation in the protogynous orange-spotted grouper

Statistics of the assembled transcripts and unigenes.

<p>Statistics of the assembled transcripts and unigenes.</p

Venn diagram of annotation results against Nr, SwissProt, COG and KEGG databases.

<p>The number in each color block indicates the number of unigenes that is annotated by single or multiple databases.</p

Frequency of classified repeat types of SSRs.

<p>Frequency of classified repeat types of SSRs.</p

Length distribution of all unigenes of <i>Trachinotus ovatus</i>.

<p>Length distribution of all unigenes of <i>Trachinotus ovatus</i>.</p

Comparison of <i>Trachinotus ovatus</i> unigenes to <i>Danio rerio</i> orthologous coding sequences.

<p>(A) The ratio of <i>Trachinotus ovatus</i> unigene length to <i>Danio rerio</i> ortholog length was plotted against <i>Trachinotus ovatus</i> unigene coverage depth. (B) The total percentage of <i>Danio rerio</i> ortholog coding sequence covered by all <i>Trachinotus ovatus</i> unigenes.</p