Genome-wide identification, evolution and expression analysis of nuclear receptor superfamily in Nile tilapia, Oreochromis niloticus

Molecular cloning of estrogen receptor and its function on vitellogenesis in pompano (Trachinotus ovatus)

Estrogen receptors (ERs) play a critical role in vitellogenesis (Vtgs). However, the contribution of each ER for the regulation of vtgs expression was not analyzed clearly in teleosts. In the present study, three ers isoforms (erα, erβ1, and erβ2) were cloned in pompano (Trachinotus ovatus). Real-time PCR and enzyme-linked immunosorbent assay (ELISA) was used to detect the effects of 17β-estradiol (E2) on ERs and Vtgs in the liver of pompano. In vivo injection experiments showed that E2 significantly increased the expressions of ers and vtgs. ER broad spectrum antagonist Fulvestrant significantly attenuated the E2- induced up-regulation of ers and vtgs in a dose-dependent manner. ERα antagonist Methyl-piperidino pyrazole (MPP) significantly attenuated the up-regulation of erα, erβ2, vtg-B and vtg-C, and promoted the expressions of erβ1 and vtg-A. ERβ antagonist Cyclofenil significantly inhibited the expressions of erβ1, erβ2, vtg-A and vtg-C, and promoted the expressions of erα and vtg-B. In addition, E2 significantly increased the protein level of Vtg, while Fulvestrant, MPP and Cyclofenil significantly inhibited the protein level of Vtg in a dose-dependent manner. Our results indicate that E2 may regulate the expression of each vtg with different subtypes of ERs, and shows a distinct compensatory expression effect on the regulation for ers and vtgs, which provides a theoretical basis for reproductive endocrinology study in pompano.

Genome-wide identification, evolution of histone lysine demethylases (KDM) genes and their expression during gonadal development in Nile tilapia

Histone lysine demethylases (KDM) are responsible for histone demethylation and are involved in gene expression regulation. Previous studies have shown that histone lysine demethylation plays an important role in gonadal development of vertebrates. The KDM family consists of eight subfamilies, i.e., kdm1, kdm2, kdm3, kdm4, kdm5, kdm6, kdm7 and JmjC-only subfamily. In this study, 13 to 63 KDM genes in 23 representative species were identified based on the available version of genome assembly. Phylogenetic relationships, domain architecture, and synteny of these genes were comprehensively analyzed and the results suggested KDM genes probably originated from the early diverging metazoan and significantly expanded in vertebrates with multiple whole genome duplication, especially in the third-round whole genome duplication (3R-WGD) and polyploidization of teleosts. The subfamilies of kdm2, kdm3, kdm4, kdm5, kdm6 and kdm7 were duplicated with 1R-2R events, and duplicates of kdm2a, kdm4a, kdm5b and kdm6b were resulted from 3R-WGD. Based on transcriptome data, the KDM genes were found to be dominantly expressed in the ovary and testis. More than 80% of KDM genes displayed sexual dimorphic expression, with 15 genes dominantly expressed in ovaries, and 12 genes dominantly expressed in testes. Importantly, from transcriptome data, qRT-PCR and fluorescence in situ hybridization during sex reversal, genes with higher expression in ovary than testis, such as kdm1b and two JmjC-only subfamily members hspbap1 and riox1, were downregulated, while other genes, such as kdm3c, kdm5bb, kdm6ba, kdm6bb and kdm7b, with higher expression in testis than ovary, were upregulated in ovotestis, indicating these genes play critical roles in the gonadal development and sex reversal. This study provided new insights into the evolution of the KDM genes and a fundamental clue for understanding their important roles in sex differentiation and gonadal development in teleosts.

Characterization, expression, and regulatory effects of nr0b1a and nr0b1b in spotted scat (Scatophagus argus)

Nuclear receptor subfamily 0 group B member 1 (Nr0b1) belongs to the nuclear receptor (NR) superfamily. It plays critical roles in sex determination, sex differentiation, and gonadal development in mammals. In this study, the duplicated genes nr0b1a and nr0b1b were identified in spotted scat (Scatophagus argus). Phylogenetic and synteny analyses revealed that, unlike nr0b1a, nr0b1b was retained in several species of teleosts after an nr0b1 gene duplication event but was secondarily lost in other fish species, amphibians, reptiles, birds, and mammals. In a sequence analysis, only 1.5 LXXLL-related repeat motifs were identified in spotted scat Nr0b1a, Nr0b1b, and non-mammalian Nr0b1a/Nr0b1, different from the 3.5 repeat motifs in mammalian Nr0b1. By qPCR, nr0b1a and nr0b1b were highly expressed in testes from stages IV to V and in ovaries from stages II to IV, respectively. Male-to-female sex reversal was induced in XY spotted scat by the administration of exogenous E2. A qPCR analysis showed that nr0b1b mRNA expression was higher in sex-reversed XY fish than in control XY fish, with no difference in nr0b1a. A luciferase assay showed that spotted scat Nr0b1a and Nr0b1b did not individually activate cyp19a1a gene transcription. As in mammals, spotted scat Nr0b1a suppressed Nr5a1-mediated cyp19a1a expression, despite containing only 1.5 LXXLL-related repeat motifs in its N-terminal region, while Nr0b1b stimulated Nr5a1-mediated cyp19a1a transcription. These results demonstrated that nr0b1a and nr0b1b in spotted scat have distinct expression patterns and regulatory effects and further indicate that nr0b1b might be involved in ovarian development by regulating Nr5a1-mediated cyp19a1a expression.

Beta-Hydroxysteroid Dehydrogenase Genes in Orange-Spotted Grouper (Epinephelus coioides): Genome-Wide Identification and Expression Analysis During Sex Reversal

Beta-hydroxysteroid dehydrogenases (β-HSDs) are a group of steroidogenic enzymes that are involved in steroid biosynthesis and metabolism, and play a crucial role in mammalian physiology and development, including sex determination and differentiation. In the present study, a genome-wide analysis identified the numbers of β-hsd genes in orange-spotted grouper (Epinephelus coioides) (19), human (Homo sapiens) (22), mouse (Mus musculus) (24), chicken (Gallus gallus) (16), xenopus (Xenopus tropicalis) (24), coelacanth (Latimeria chalumnae) (17), spotted gar (Lepisosteus oculatus) (14), zebrafish (Danio rerio) (19), fugu (Takifugu rubripes) (19), tilapia (Oreochromis niloticus) (19), medaka (Oryzias latipes) (19), stickleback (Gasterosteus aculeatus) (17) and common carp (Cyprinus carpio) (27) samples. A comparative analysis revealed that the number of β-hsd genes in teleost fish was no greater than in tetrapods due to gene loss followed by a teleost-specific whole-genome duplication event. Based on transcriptome data from grouper brain and gonad samples during sex reversal, six β-hsd genes had relatively high expression levels in the brain, indicating that these genes may be required for neurogenesis or the maintenance of specific biological processes in the brain. In the gonad, two and eight β-hsd genes were up- and downregulated, respectively, indicating their important roles in sex reversal. Our results demonstrated that β-hsd genes may be involved in the sex reversal of grouper by regulating the synthesis and metabolism of sex steroid hormones.

Characterization of Matrix Metalloprotease-9 Gene from Nile tilapia (Oreochromis niloticus) and Its High-Level Expression Induced by the Streptococcus agalactiae Challenge

The bacterial diseases of tilapia caused by Streptococcus agalactiae have resulted in the high mortality and huge economic loss in the tilapia industry. Matrix metalloproteinase-9 (MMP-9) may play an important role in fighting infection. However, the role of MMP-9 in Nile tilapia against S. agalactiae is still unclear. In this work, MMP-9 cDNA of Nile tilapia (NtMMP-9) has been cloned and characterized. NtMMP-9 has 2043 bp and encodes a putative protein of 680 amino acids. NtMMP-9 contains the conserved domains interacting with decorin and inhibitors via binding forces compared to those in other teleosts. Quantitative real-time-polymerase chain reaction (qPCR) analysis reveals that NtMMP-9 distinctly upregulated following S. agalactiae infection in a tissue- and time-dependent response pattern, and the tissues, including liver, spleen, and intestines, are the major organs against a S. agalactiae infection. Besides, the proteolytic activity of NtMMP-9 is also confirmed by heterologous expression and zymography, which proves the active function of NtMMP-9 interacting with other factors. The findings indicate that NtMMP-9 was involved in immune responses against the bacterial challenge at the transcriptional level. Further work will focus on the molecular mechanisms of NtMMP-9 to respond and modulate the signaling pathways in Nile tilapia against S. agalactiae invasion and the development of NtMMP-9-related predictive biomarkers or vaccines for preventing bacterial infection in the tilapia industry.

Establishment of three estrogen receptors (esr1, esr2a, esr2b) knockout lines for functional study in Nile tilapia

Estrogens play fundamental roles in regulating reproductive activities and they act through estrogen receptors (ESRs) in all vertebrates. To date, distinct roles of estrogen receptors have been characterized only in human and model organisms, including mouse, rat, zebrafish and medaka. Physiological role of estrogen/receptor signaling in reproduction remains poorly defined in non-model organisms. In the present study, we successfully generated esr1, esr2a and esr2b mutant lines in tilapia by CRISPR/Cas9 and examined their phenotypes. Surprisingly, the esr1 mutants showed no phenotypes of reproductive development and function in both females and males. The esr2a mutant females showed significantly delayed ovarian development and follicle growth at 90 and 180 dah, and the development caught up later at 360 dah. The esr2a mutant males showed no phenotypes at 90 dah, and displayed smaller gonads and efferent ducts, less spermatogonia and more abnormal sperms at 180 dah. In contrast, the esr2b mutants displayed abnormal development of ovarian ducts and efferent ducts which failed to connect to the genital orifice, and which in turn, resulted in infertility in female and male, respectively, although they produced gametes in their gonads. Taken together, our study provides evidence for differential functions of esr1, esr2a and esr2b in fish reproduction.

Genome-wide identification, evolution of DNA methyltransferases and their expression during gonadal development in Nile tilapia

DNA methyltransferases (dnmts) are responsible for DNA methylation and play important roles in organism development. In this study, seven dnmts genes (dnmt1, dnmt2, dnmt3aa, dnmt3ab, dnmt3ba, dnmt3bb.1, dnmt3bb.2) were identified in Nile tilapia. Comprehensive analyses of dnmts were performed using available genome databases from representative animal species. Phylogenetic analysis revealed that the dnmts family were highly conserved in teleosts. Based on transcriptome data from eight adult tilapia tissues, the dnmts were found to be dominantly expressed in the head kidney, testis and ovary. Analyses of the gonadal transcriptome data in different developmental stages revealed that all dnmts were expressed in both ovary and testis, and four de novo dnmts (dnmt3aa, dnmt3ab, dnmt3bb.1, dnmt3bb.2) showed higher expression in the testis than in the ovary. Furthermore, during sex reversal induced by Fadrozole, the expression of these four de novo dnmts increased significantly in treated group compared to female control group. By in situ hybridization, the seven dnmts were found to be expressed mainly in phase I and II oocytes of the ovary and spermatocytes of the testis. When gonads were incubated with a methyltransferase inhibitor (5-AzaCdR) in vitro, the expression of dnmts genes were down-regulated significantly, while the expression of cyp19a1a (a key gene in female pathway) and dmrt1 (a key gene in male pathway) increased significantly. Our results revealed the conservation of dnmts during evolution and indicated a potential role of dnmts in epigenetic regulation of gonadal development.

Identification and Evolution of TGF-β Signaling Pathway Members in Twenty-Four Animal Species and Expression in Tilapia

Transforming growth factor β (TGF-β) signaling controls diverse cellular processes during embryogenesis as well as in mature tissues of multicellular animals. Here we carried out a comprehensive analysis of TGF-β pathway members in 24 representative animal species. The appearance of the TGF-β pathway was intrinsically linked to the emergence of metazoan. The total number of TGF-β ligands, receptors, and smads changed slightly in all invertebrates and jawless vertebrates analyzed. In contrast, expansion of the pathway members, especially ligands, was observed in jawed vertebrates most likely due to the second round of whole genome duplication (2R) and additional rounds in teleosts. Duplications of TGFB2, TGFBR2, ACVR1, SMAD4 and SMAD6, which were resulted from 2R, were first isolated. Type II receptors may be originated from the ACVR2-like ancestor. Interestingly, AMHR2 was not identified in Chimaeriformes and Cypriniformes even though they had the ligand AMH. Based on transcriptome data, TGF-β ligands exhibited a tissue-specific expression especially in the heart and gonads. However, most receptors and smads were expressed in multiple tissues indicating they were shared by different ligands. Spatial and temporal expression profiles of 8 genes in gonads of different developmental stages provided a fundamental clue for understanding their important roles in sex determination and reproduction. Taken together, our findings provided a global insight into the phylogeny and expression patterns of the TGF-β pathway genes, and hence contribute to the greater understanding of their biological roles in the organism especially in teleosts.

Bioinformatic analyses of zona pellucida genes in vertebrates and their expression in Nile tilapia

Zona pellucida (ZP) genes encode ZP glycoproteins which constitute the coat surrounding oocytes and early embryos. Genome-wide identification of ZP genes is still lacking in vertebrates, especially in fish species. Herein, we conducted bioinformatic analyses of the ZP genes of the Nile tilapia and other vertebrates. Totally 16, 9, 17, 27, 21, 20, 26, 19, 14,11, 24, 17, 9, 18, 8, 11, 9, 8, 5, and 4 ZP genes belonging to 5 subfamilies (ZPA, ZPB, ZPC, ZPD, and ZPAX) were found in the sea lamprey, elephant shark, coelacanth, spotted gar, zebrafish, medaka, stickleback, Nile tilapia, Amazon molly, platyfish, seahorse, Northern snakehead, cavefish, tetraodon, clawed frog, turtle, chicken, platypus, kangaroo rat, and human genomes, respectively. The expansion of ZP genes in basal vertebrates was mainly achieved by gene duplication of ZPB, ZPC, and ZPAX subfamilies, while the shrink of ZP gene number in viviparous mammals was achieved by keeping only one copy of the ZP genes in each subfamily or even secondary loss of some subfamilies. The number of ZP gene is related to the environment where the eggs are fertilized and the embryos develop in vertebrates. Transcriptomic analysis showed that 14 ZP genes were expressed in the ovary of Nile tilapia, while two (ZPB2b and ZPC2) were highly expressed in the liver. On the other hand, ZPB1a and ZPB2c were not found to be expressed in any tissue or at any developmental stage of the gonads examined. In the ovary, the expression of ZP genes started from 30 dah (days after hatching), significantly upregulated at 90 dah and maintained this level at 180 dah. The expression of ZPC2 in the liver and ZPC5-2 and ZPAX1 in the ovary was confirmed by in situ hybridization. The ovary- and liver-expressed ZP genes are expressed coordinately with oocyte growth in tilapia.

Duplication and gene expression patterns of β-catenin in Nile tilapia

β-Catenin, a key transcriptional coactivator of the Wnt pathway, plays an important role in animal embryonic development and organogenesis. In our earlier study, we have reported that two types of β-catenin (β-catenin-1 and β-catenin-2) were ubiquitously expressed in almost all the tissues examined in tilapia. However, the immunolocalization of β-catenin in those tissues, especially in extra-gonadal tissues, remains unclear. In the present study, we further confirm that these two types of β-catenin gene exist only in teleosts and are derived from 3R (third round of genome duplication) by phylogenetic and syntenic analyses. Moreover, the transcriptome analysis conducted in this investigation reveals that two β-catenins exhibited similar expression patterns in seven adult tissues and four key developmental stages of XX and XY gonads. Finally, immunohistochemistry analysis was performed to detect the cell localization of β-catenin. A positive signal of β-catenin was observed in various tissues of tilapia, such as the intestine, liver, kidney, spleen, eye, brain, and gonads. The results of our study indicate that tilapia β-catenin might be involved in the organ development and play some specific functions in biological processes; all these data will provide basic reference for understanding the molecular mechanism of β-catenin in regulating of teleost organogenesis.

Gonad Transcriptome Analysis of High-Temperature-Treated Females and High-Temperature-Induced Sex-Reversed Neomales in Nile Tilapia

Background: Nowadays, the molecular mechanisms governing TSD (temperature-dependent sex determination) or GSD + TE (genotypic sex determination + temperature effects) remain a mystery in fish. Methods: We developed three all-female families of Nile tilapia (Oreochromis niloticus), and the family with the highest male ratio after high-temperature treatment was used for transcriptome analysis. Results: First, gonadal histology analysis indicated that the histological morphology of control females (CF) was not significantly different from that of high-temperature-treated females (TF) at various development stages. However, the high-temperature treatment caused a lag of spermatogenesis in high-temperature-induced neomales (IM). Next, we sequenced the transcriptome of CF, TF, and IM Nile tilapia. 79, 11,117, and 11,000 differentially expressed genes (DEGs) were detected in the CF–TF, CF–IM, and TF–IM comparisons, respectively, and 44 DEGs showed identical expression changes in the CF–TF and CF–IM comparisons. Principal component analysis (PCA) indicated that three individuals in CF and three individuals in TF formed a cluster, and three individuals in IM formed a distinct cluster, which confirmed that the gonad transcriptome profile of TF was similar to that of CF and different from that of IM. Finally, six sex-related genes were validated by qRT-PCR. Conclusions: This study identifies a number of genes that may be involved in GSD + TE, which will be useful for investigating the molecular mechanisms of TSD or GSD + TE in fish.

Annotation of the Nuclear Receptors in an Estuarine Fish species, Fundulus heteroclitus

The nuclear receptors (NRs) are ligand-dependent transcription factors that respond to various internal as well as external cues such as nutrients, pheromones, and steroid hormones that play crucial roles in regulation and maintenance of homeostasis and orchestrating the physiological and stress responses of an organism. We annotated the Fundulus heteroclitus (mummichog; Atlantic killifish) nuclear receptors. Mummichog are a non-migratory, estuarine fish with a limited home range often used in environmental research as a field model for studying ecological and evolutionary responses to variable environmental conditions such as salinity, oxygen, temperature, pH, and toxic compounds because of their hardiness. F. heteroclitus have at least 74 NRs spanning all seven gene subfamilies. F. heteroclitus is unique in that no RXRα member was found within the genome. Interestingly, some of the NRs are highly conserved between species, while others show a higher degree of divergence such as PXR, SF1, and ARα. Fundulus like other fish species show expansion of the RAR (NR1B), Rev-erb (NR1D), ROR (NR1F), COUPTF (NR2F), ERR (NR3B), RXR (NR2B), and to a lesser extent the NGF (NR4A), and NR3C steroid receptors (GR/AR). Of particular interest is the co-expansion of opposing NRs, Reverb-ROR, and RAR/RXR-COUPTF.

Cloning, expression and functional characterization on vitellogenesis of estrogen receptors in Scatophagus argus

Estrogen receptors (Er) play a critical role in vitellogenesis. Three ers (erα, erβ1 and erβ2) and vitellogenins (vtg-A, vtg-B and vtg-C) subtypes were isolated in various fish species, while the contribution of each Er to the regulation of vtgs expression was not analyzed in detail. Here, erα, erβ1 and erβ2 were cloned and all were found to be expressed in female liver in Scatophagus argus. During proteic vitellogenesis stage, erα was simultaneously up-regulated, while erβ1 and erβ2 were not, with three vtgs in female liver. The effects of 17β-estradiol (E₂) alone or combined with Er antagonists on ers, vtgs mRNA expressions and Vtg protein content in incubated male liver were examined by real-time PCR and enzyme-linked immunosorbent assay (ELISA), respectively. The expressions of erα, erβ1, vtgs mRNA and Vtg protein increased significantly after 24h incubation with E₂ (0.1, 1 and 10μM), while Er nonselective antagonist ICI 182 780 (0.01, 0.1 and 1μM) significantly attenuated the up-regulation effects of E₂ on ers, vtgs mRNA and Vtg protein in a dose-dependent manner. Erα selective antagonist Methyl-piperidinopyrazole (MPP) (0.01, 0.1 and 1μM) significantly attenuated the up-regulation effects of E₂ on erα, vtg-B, vtg-C mRNA and Vtg protein, while promoted the expression of erβ1 and vtg-A. Erβ selective antagonist Cyclofenil (0.01, 0.1 and 1μM) attenuated the up-regulation effects of E₂ on erβ1, erβ2, vtg-A, vtg-C mRNA and Vtg protein while promoted the expression of erα and vtg-B. Our results suggest that the regulation of Ers on different vtgs was divergent in S. argus.

Blockage of androgen and administration of estrogen induce transdifferentiation of testis into ovary

Induction of sex reversal of XY fish has been restricted to the sex undifferentiated period. In the present study, differentiated XY tilapia were treated with trilostane (TR), metopirone (MN) and glycyrrhetinic acid (GA) (inhibitor of 3β-HSD, Cyp11b2 and 11β-HSD, respectively) alone or in combination with 17β-estradiol (E2) from 30 to 90 dah (days after hatching). At 180 dah, E2 alone resulted in 8.3%, and TR, MN and GA alone resulted in no secondary sex reversal (SSR), whereas TR + E2, MN + E2 and GA + E2 resulted in 88.3, 60.0 and 46.7% of SSR, respectively. This sex reversal could be rescued by simultaneous administration of 11-ketotestosterone (11-KT). Compared with the control XY fish, decreased serum 11-KT and increased E2 level were detected in SSR fish. Immunohistochemistry analyses revealed that Cyp19a1a, Cyp11b2 and Dmrt1 were expressed in the gonads of GA + E2, MN + E2 and TR + E2 SSR XY fish at 90 dah, but only Cyp19a1a was expressed at 180 dah. When the treatment was applied from 60 to 120 dah, TR + E2 resulted in 3.3% of SSR, MN + E2 and GA + E2 resulted in no SSR. These results demonstrated that once 11-KT was synthesized, it could antagonize E2-induced male-to-female SSR, which could be abolished by simultaneous treatment with the inhibitor of steroidogenic enzymes. The upper the enzyme was located in the steroidogenic pathway, the higher SSR rate was achieved when it was inhibited as some of the precursors, such as androstenedione, testosterone and 5α-dihydrotestosterone, could act as androgens. These results highlight the key role of androgen in male sex maintenance.

Genome-wide identification, evolution of chromobox family genes and their expression in Nile tilapia

Chromobox (Cbx) family proteins are transcriptional repressors that involved in epigenetic and developmental processes. In this study, comprehensive analyses of Cbxs were performed using available genome databases from representative animal species. The Cbx family were originated from one Polycomb (Pc) gene like the yeast Pc, which duplicated into two and gave rise to the Pc and the Heterochromatin protein 1 (Hp1) identified in invertebrates from protozoon to lancelet. Rapid expansion of Cbx family members was observed in vertebrates as ~8 (5 Pc and 3 Hp1) were identified in spotted gar, coelacanth and tetrapods. Further expansion of the members to ~14 (9 Pc and 5 Hp1) was observed in teleosts due to the third round genome duplication (3R). Based on transcriptome data from eight adult tilapia tissues, most of the Cbxs were found to be dominantly expressed in the brain, testis, ovary and heart. Analyses of the gonadal transcriptome data from four developmental stages revealed that all Cbxs were expressed in both ovary and testis except Cbx7b, with significant increase of the total and average RPKM from 5 to 90dah (days after hatching). By in situ hybridization, the three most highly and sexual dimorphically expressed Cbx genes in gonads, Cbx1b, Cbx3a and Cbx5, were found to be expressed in phase I and II oocytes of the ovary, and in secondary spermatocytes (Cbx1b and Cbx3a) and spermatids (Cbx5) of the testis. Our results revealed the evolution of Cbx genes and indicated a potential role of Cbxs in epigenetic regulation of gametogenesis.

Genomic identification, rapid evolution, and expression of Argonaute genes in the tilapia, Oreochromis niloticus

Argonaute proteins are key components of the small RNA-induced silencing complex and have multiple roles in RNA-directed regulatory pathways. Argonaute genes can be divided into two subfamilies: the Ago (interacting with microRNA/small interfering RNA) and Piwi subfamilies (interacting with piwi-interacting RNAs (piRNAs)). In the present study, genome-wide analyses firstly yielded the identification of different members of Agos and Piwis in the tilapia, coelacanth, spotted gar, and elephant shark. The additional teleost Ago3b was generated following the fish-specific genome duplication event. Selective pressure analysis on Agos and Piwis between cichlids and other teleosts showed an accelerated evolution of Piwil1 in the cichlid lineages, and the positive selected sites were located in the region of PIWI domain, suggesting that these amino acid substitutions are adapt to targeted cleavage of messenger RNA (mRNA) in cichlids. Ago1 and Ago4 were detected at higher levels at 5 days after hatching (dah) in both ovaries and testes compared with other stages, supporting the previously reported requirement of Ago-mediated pathways to clear the maternal mRNAs during the early embryogenesis. The Piwis were abundantly expressed in tilapia testes, indicating their essential roles in male germline, especially in spermatogenesis. Notable expression of Piwis was also detected in skeletal muscle, indicating that piRNA pathway may not only be confined to development and maintenance of the germline but may also play important roles in somatic tissues. The expression of Piwil1 and Piwil2 was examined by quantitative PCR (qPCR) and in situ hybridization (ISH) to validate the spatial and temporal expression profiles. Taken together, these results present a thorough overview of tilapia Argonaute family and provide a new perspective on the evolution and function of this family in teleosts.

Genome-Wide Identification and Transcriptome-Based Expression Profiling of the Sox Gene Family in the Nile Tilapia (Oreochromis niloticus)

The Sox transcription factor family is characterized with the presence of a Sry-related high-mobility group (HMG) box and plays important roles in various biological processes in animals, including sex determination and differentiation, and the development of multiple organs. In this study, 27 Sox genes were identified in the genome of the Nile tilapia (Oreochromis niloticus), and were classified into seven groups. The members of each group of the tilapia Sox genes exhibited a relatively conserved exon-intron structure. Comparative analysis showed that the Sox gene family has undergone an expansion in tilapia and other teleost fishes following their whole genome duplication, and group K only exists in teleosts. Transcriptome-based analysis demonstrated that most of the tilapia Sox genes presented stage-specific and/or sex-dimorphic expressions during gonadal development, and six of the group B Sox genes were specifically expressed in the adult brain. Our results provide a better understanding of gene structure and spatio-temporal expression of the Sox gene family in tilapia, and will be useful for further deciphering the roles of the Sox genes during sex determination and gonadal development in teleosts.

Draft genome sequence of the silver pomfret fish, Pampus argenteus

Silver pomfret, Pampus argenteus, is a fish species from coastal waters. Despite its high commercial value, this edible fish has not been sequenced. Hence, its genetic and genomic studies have been limited. We report the first draft genome sequence of the silver pomfret obtained using a Next Generation Sequencing (NGS) technology. We assembled 38.7 Gb of nucleotides into scaffolds of 350 Mb with N50 of about 1.5 kb, using high quality paired end reads. These scaffolds represent 63.7% of the estimated silver pomfret genome length. The newly sequenced and assembled genome has 11.06% repetitive DNA regions, and this percentage is comparable to that of the tilapia genome. The genome analysis predicted 16 322 genes. About 91% of these genes showed homology with known proteins. Many gene clusters were annotated to protein and fatty-acid metabolism pathways that may be important in the context of the meat texture and immune system developmental processes. The reference genome can pave the way for the identification of many other genomic features that could improve breeding and population-management strategies, and it can also help characterize the genetic diversity of P. argenteus.