Transcriptome display during tilapia sex determination and differentiation as revealed by RNA-Seq analysis

A New Investigation to Discriminate Sexes in Alive Nile Tilapia (Oreochromis niloticus) Using Cyp19a1a and Dmrt1 Gene Expression in Tail Fin Tissues

The Nile Tilapia (Oreochromis niloticus), a gonochoristic teleost fish with a XX/XY sex-determination system, is an ideal model for investigating gonadal sex differentiation. During gonadal differentiation, the expression of cyp19a1a in XX gonads and dmrt1 in XY gonads are required for undifferentiated tissues to develop into ovary or testis. In this study, quantitative real-time RT-PCR assessed the expression of cyp19a1a and dmrt1 genes in gonads and tail fin tissues. Differences in gene expression mean among sexually differentiated fish were analyzed using two-way analysis of variance (ANOVA) and validation of mixed model using discriminant analysis (DA) for morphometric traits and the gene expression in gonads and tail fin tissues used to validate and utilize them in discriminating sexes in sex-differentiated Nile Tilapia fish. The results revealed that, cyp19a1a gene expression in female ovaries was more significant than dmrt1 in male testis. In the other hand, the dmrt1 gene expression in the tail fin was higher in males than females. Both, cyp19a1a and dmrt1 genes, can discriminate fish sexes by 100% by using their expression in tail fin tissues. In conclusion, the cyp19a1a and dmrt1 genes could be used as a genetic marker to discriminate between the Nile Tilapia sexes, whereas used as an indicator for ovarian or testis differentiation in sexually differentiated Nile Tilapia using tail fin tissues. It is worth mentioning that this is the first investigation for using cyp19a1a and dmrt1 genes from Nile Tilapia tail fin tissues in sex determination.

Disruption of Zar1 leads to arrested oogenesis by regulating polyadenylation via Cpeb1 in tilapia (Oreochromis niloticus)

Oogenesis is a complex process regulated by precise coordination of multiple factors, including maternal genes. Zygote arrest 1 (zar1) has been identified as an ovary-specific maternal gene that is vital for oocyte-to-embryo transition and oogenesis in mouse and zebrafish. However, its function in other species remains to be elucidated. In the present study, zar1 was identified with conserved C-terminal zinc finger domains in Nile tilapia. zar1 was highly expressed in the ovary and specifically expressed in phase I and II oocytes. Disruption of zar1 led to the failed transition from oogonia to phase I oocytes, with somatic cell apoptosis. Down-regulation and failed polyadenylation of figla, gdf9, bmp15 and wee2 mRNAs were observed in the ovaries of zar1-/- fish. Cpeb1, a gene essential for polyadenylation that interacts with Zar1, was down-regulated in zar1-/- fish. Moreover, decreased levels of serum estrogen and increased levels of androgen were observed in zar1-/- fish. Taken together, zar1 seems to be essential for tilapia oogenesis by regulating polyadenylation and estrogen synthesis. Our study shows that Zar1 has different molecular functions during gonadal development by the similar signaling pathway in different species.

Tilapia, a good model for studying reproductive endocrinology

The Nile tilapia (Oreochromis niloticus), with a system of XX/XY sex determination, is a worldwide farmed fish with a shorter sexual maturation time than that of most cultured fish. Tilapia show a spawning cycle of approximately 14 days and can be artificially propagated in the laboratory all year round to obtain genetically all female (XX) and all male (XY) fry. Its genome sequence has been opened, and a perfect gene editing platform has been established. With a moderate body size, it is convenient for taking enough blood to measure hormone level. In recent years, using tilapia as animal model, we have confirmed that estrogen is crucial for female development because 1) mutation of star2, cyp17a1 or cyp19a1a (encoding aromatase, the key enzyme for estrogen synthesis) results in sex reversal (SR) due to estrogen deficiency in XX tilapia, while mutation of star1, cyp11a1, cyp17a2, cyp19a1b or cyp11c1 affects fertility due to abnormal androgen, cortisol and DHP levels in XY tilapia; 2) when the estrogen receptors (esr2a/esr2b) are mutated, the sex is reversed from female to male, while when the androgen receptors are mutated, the sex cannot be reversed; 3) the differentiated ovary can be transdifferentiated into functional testis by inhibition of estrogen synthesis, and the differentiated testis can be transdifferentiated into ovary by simultaneous addition of exogenous estrogen and androgen synthase inhibitor; 4) loss of male pathway genes amhy, dmrt1, gsdf causes SR with upregulation of cyp19a1a in XY tilapia. Disruption of estrogen synthesis rescues the male to female SR of amhy and gsdf but not dmrt1 mutants; 5) mutation of female pathway genes foxl2 and sf-1 causes SR with downregulation of cyp19a1a in XX tilapia; 6) the germ cell SR of foxl3 mutants fails to be rescued by estrogen treatment, indicating that estrogen determines female germ cell fate through foxl3. This review also summarized the effects of deficiency of other steroid hormones, such as androgen, DHP and cortisol, on fish reproduction. Overall, these studies demonstrate that tilapia is an excellent animal model for studying reproductive endocrinology of fish.

Identification of hub genes in digestive system of mandarin fish (Siniperca chuatsi) fed with artificial diet by weighted gene co-expression network analysis

Mandarin fish (Siniperca chuatsi) is a carnivorous freshwater fish and an economically important species. The digestive system (liver, stomach, intestine, pyloric caecum, esophagus, and gallbladder) is an important site for studying fish domestication. In our previous study, we found that mandarin fish undergoes adaptive changes in histological morphology and gene expression levels of the digestive system when subjected to artificial diet domestication. However, we are not clear which hub genes are highly associated with domestication. In this study, we performed WGCNA on the transcriptomes of 17 tissues and 9 developmental stages and combined differentially expressed genes analysis in the digestive system to identify the hub genes that may play important functions in the adaptation of mandarin fish to bait conversion. A total of 31,657 genes in 26 samples were classified into 23 color modules via WGCNA. The modules midnightblue, darkred, lightyellow, and darkgreen highly associated with the liver, stomach, esophagus, and gallbladder were extracted, respectively. Tan module was highly related to both intestine and pyloric caecum. The hub genes in liver were cp, vtgc, c1in, c9, lect2, and klkb1. The hub genes in stomach were ghrl, atp4a, gjb3, muc5ac, duox2, and chia2. The hub genes in esophagus were mybpc1, myl2, and tpm3. The hub genes in gallbladder were dyst, npy2r, slc13a1, and slc39a4. The hub genes in the intestine and pyloric caecum were slc15a1, cdhr5, btn3a1, anpep, slc34a2, cdhr2, and ace2. Through pathway analysis, modules highly related to the digestive system were mainly enriched in digestion and absorption, metabolism, and immune-related pathways. After domestication, the hub genes vtgc and lect2 were significantly upregulated in the liver. Chia2 was significantly downregulated in the stomach. Slc15a1, anpep, and slc34a2 were significantly upregulated in the intestine. This study identified the hub genes that may play an important role in the adaptation of the digestive system to artificial diet, which provided novel evidence and ideas for further research on the domestication of mandarin fish from molecular level.

Genome-wide association study reveals genomic loci of sex differentiation and gonadal development in Plectropomus leopardus

Introduction: Plectropomus leopardus, a commercially significant marine fish, is primarily found in the Western Pacific regions and along the coast of Southeast Asia. A thorough analysis of the molecular mechanisms involved in sex differentiation is crucial for gaining a comprehensive understanding of gonadal development and improving sex control breeding. However, the relevant fundamental studies of P. leopardus are relatively lacking. Methods: In this study, a genome-wide association study (GWAS) was conducted to investigate the genetic basis mechanism of sex differentiation and gonadal developmental traits in P. leopardus utilizing about 6,850,000 high-quality single-nucleotide polymorphisms (SNPs) derived from 168 individuals (including 126 females and 42 males) by the genome-wide efficient mixed-model association (GEMMA) algorithm. Results: The results of these single-trait GWASs showed that 46 SNP loci (-log10 p > 7) significantly associated with sex differentiation, and gonadal development traits were distributed in multiple different chromosomes, which suggested the analyzed traits were all complex traits under multi-locus control. A total of 1,838 potential candidate genes were obtained by considering a less-stringent threshold (-log10 p > 6) and ±100 kb regions surrounding the significant genomic loci. Moreover, 31 candidate genes were identified through a comprehensive analysis of significant GWAS peaks, gene ontology (GO) annotations, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, including taf7, ddx6, apoeb, sgk1, a2m, usf1, hsd3b7, dll4, xbp1, tet3, esr1, and gli3. These trait-associated genes have been shown to be involved in germline development, male sex differentiation, gonad morphogenesis, hormone receptor binding, oocyte development, male gonad development, steroidogenesis, estrogen-synthetic pathway, etc. Discussion: In the present study, multiple genomic loci of P. leopardus associated with sex differentiation and gonadal development traits were identified for the first time by using GWAS, providing a valuable resource for further research on the molecular genetic mechanism and sex control in P. leopardus. Our results also can contribute to understanding the genetic basis of the sex differentiation mechanism and gonadal development process in grouper fish.

Combined effects of rearing temperature regime (thermocycle vs. constant temperature) during early development and thermal treatment on Nile tilapia (Oreochromis niloticus) sex differentiation

In nature, water temperature experiences daily variations known as thermocycles. Temperature is the main environmental factor that influences sex determination in most teleost fish. The purpose of this study was to examine the effects of rearing temperature (thermocycle (TC) vs. constant (CTE)) on development and a posterior thermal shock throughout the period of sex differentiation of Nile tilapia (Oreochromis niloticus). Embryos and larvae were kept under two temperature regimes: TC of 31 °C:25 °C day:night vs. CTE of 28 °C from 0 to 11 dpf. After this period, the larvae from each group were subjected to either heat treatment (HT, 36 °C for 12 days) or kept under the same rearing temperatures until 23 dpf (Control, C). Then all the groups remained at constant temperature until 270 dpf, when blood and gonads were collected. Larval samples were used to examine the expression of genes related to male (amh, ara, sox9a, dmrt1a) and female (cyp19a1a, foxl2, era) sexual differentiation. In juveniles, sex was characterized by histology, the gonadal expression of the genes involved in the sex steroid synthesis was analyzed by qPCR, and plasma testosterone (T) and estradiol (E₂) levels were analyzed by ELISA. In larvae, daily TCs increased the survival rate against HT and up-regulated the expression of ovarian differentiation genes. In juveniles, TC + C induced a higher proportion of females and higher cyp19a1a expression compared to CTE + C. HT induced changes in the CTE group by up-regulating testicular differentiation genes and down-regulating female promoting genes, which did not occur in the TC group. Juveniles from TC + C group presented a higher proportion of females with higher E₂ and cyp19a1a than CTE + HT. Fish from the CTE + HT group showed a higher percentage of males with highest T and amh. These findings indicate that daily TCs during larval development promote ovarian differentiation and diminish the masculinizing effects of HT.

Sex steroid profile during oocyte development and maturation in the intertidal worm Marphysa madrasi (Polychaeta: Eunicidae) from the east coast of India

Marphysa madrasi is a commercially valuable maturation diet in crustacean aquaculture. This study presents the first detailed investigation of oogenesis in the intertidal polychaete worm M. madrasi and reports the steroid profile during oocyte growth and development. Oogenesis is extraovarian type I, originating from coelomic epithelial cells, with four stages of development - primary growth, early vitellogenic, late vitellogenic, and maturation. The primary growth phase contains oogonial cells and previtellogenic oocyte clusters in the early, mid, and late stages of development form a dispersed ovary attached to blood vessels. The late previtellogenic oocytes detach from the ovary at the onset of vitellogenesis. The detached oocytes complete vitellogenesis and final maturation in the coelomic fluid as solitary free-floating cells without any connection with follicle cells. The worms display asynchronous reproduction with a heterogeneous population of developing oocytes. Steroid extracts from the polychaete homogenates in different stages of oogenesis were identified by HPLC and confirmed by LC-MS/MS. In M. madrasi, two vertebrate-type steroids, pregnenolone (P5) and 17α-hydroxyprogesterone (17-OHP) were detected and quantified. The P5 levels were low in immature worms but increased significantly by ∼ 8.3-fold in the previtellogenic stage and peaked during oocyte maturation. 17-OHP levels were low in immature worms but gradually increase as the oogenesis progress to the primary growth and early vitellogenic phase, with a significant increase (p < 0.001) during the late vitellogenic phase. Although an increase in the concentration of P5 and 17-OHP during vitellogenesis and maturation of oocytes points to a possible role in reproduction, the absence of other vertebrate-type steroids in the investigated polychaete signifies a plausible uptake of P5 and 17-OHP from the environment.

Dynamic transcriptome analysis reveals the gene network of gonadal development from the early history life stages in dwarf surfclam Mulinia lateralis

BACKGROUND

Gonadal development is driven by a complex genetic cascade in vertebrates. However, related information remains limited in molluscs owing to the long generation time and the difficulty in maintaining whole life cycle in the lab. The dwarf surfclam Mulinia lateralis is considered an ideal bivalve model due to the short generation time and ease to breed in the lab.

RESULTS

To gain a comprehensive understanding of gonadal development in M. lateralis, we conducted a combined morphological and molecular analysis on the gonads of 30 to 60 dpf. Morphological analysis showed that gonad formation and sex differentiation occur at 35 and 40-45 dpf, respectively; then the gonads go through gametogenic cycle. Gene co-expression network analysis on 40 transcriptomes of 35-60 dpf gonads identifies seven gonadal development-related modules, including two gonad-forming modules (M6, M7), three sex-specific modules (M14, M12, M11), and two sexually shared modules (M15, M13). The modules participate in different biological processes, such as cell communication, glycan biosynthesis, cell cycle, and ribosome biogenesis. Several hub transcription factors including SOX2, FOXZ, HSFY, FOXL2 and HES1 are identified. The expression of top hub genes from sex-specific modules suggests molecular sex differentiation (35 dpf) occurs earlier than morphological sex differentiation (40-45 dpf).

CONCLUSION

This study provides a deep insight into the molecular basis of gonad formation, sex differentiation and gametogenesis in M. lateralis, which will contribute to a comprehensive understanding of the reproductive regulation network in molluscs.

Testosterone nanoemulsion produced masculinized Nile tilapia (Oreochromis niloticus)

The objective of this work was to develop a food additive for the sex reversal of Nile tilapia (Oreochromis niloticus) based on a simple oil in water (O/W) nanoemulsion with testosterone propionate for incorporation into commercial feed. Oil screening and evaluation of the organoleptic and physicochemical characteristics were carried out to determine the best formulation. A palatability test was also performed. Sex reversal test was assayed using 5 experimental groups: negative control - macerated feed without hormone; free testosterone - macerated feed with 60 mg/kg of testosterone propionate diluted in ethanol; and macerated feed with testosterone propionate nanoemulsion at a concentration of 30, 60, and 90 mg/kg. Stable nanoemulsions (size 76-210 nm) with testosterone propionate were produced. All nanoemulsion-added feed was palatable to tilapia. We obtained sex reversal values of ≈65, 75, and 72% in the groups of 30, 60, and 90 mg/kg, respectively. We can conclude that the nanoemulsion showed promising results; it is capable of inducing sex reversal in tilapia, is suitable as a commercial product, and has the potential to promote safety for rural staff and reduce the environmental impact of hormones.

Identification and Characterization of Sex-Biased miRNAs in the Golden Pompano (Trachinotus blochii)

The golden pompano (Trachinotus blochii) is a marine fish of considerable commercial importance in China. It shows notable sexual size dimorphism; the growth rate of females is faster than that of males. Therefore, sex-biased research is of great importance in T. blochii breeding. However, there have been few studies on sex differentiation and mechanisms underlying sex determination in T. blochii. MicroRNAs (miRNAs) play crucial roles in sex differentiation and determination in animals. However, limited miRNA data are available on fish. In this study, two small RNA libraries prepared from the gonads of T. blochii were constructed and sequenced. The RNA-seq analysis yielded 1366 known and 69 novel miRNAs with 289 significantly differentially expressed miRNAs (p < 0.05). Gene ontology (GO) analysis confirmed that the TFIIA transcription factor complex (GO: 0005672) was the most significantly enriched GO term. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the differentially expressed miRNAs and target genes were mainly related to sex determination and gonadal developmental signaling pathways, specifically the Wnt signaling pathway, MAPK signaling pathway, and steroid biosynthetic pathway. MiRNA-mRNA co-expression network analysis strongly suggested a role for sex-biased miRNAs in sex determination/differentiation and gonadal development. For example, gata4, foxo3, wt1, and sf1 genes were found to be regulated by bta-miR-2898; esr2 and foxo3 by novel_176, and ar by oar-let-7b. Quantitative real-time polymerase chain reaction analysis of selected mRNAs and miRNAs validated the integrated analysis. This study established a set of sex-biased miRNAs that are potential regulatory factors in gonadal development in T. blochii. These results provide new insight into the function of miRNAs in sex differentiation and determination in T. blochii and highlight some key miRNAs for future studies.

Developmental roles and molecular mechanisms of Asterix/GTSF1

Maintenance of germline genomic integrity is critical for the survival of animal species. Consequently, many cellular and molecular processes have evolved to ensure genetic stability during the production of gametes. Here, we describe the discovery, characterization, and emerging molecular mechanisms of the protein Asterix/Gametocyte‐specific factor 1 (GTSF1), an essential gametogenesis factor that is conserved from insects to humans. Beyond its broad importance for healthy germline development, Asterix/GTSF1 has more specific functions in the Piwi‐interacting RNA (piRNA)–RNA interference pathway. There, it contributes to the repression of otherwise deleterious transposons, helping to ensure faithful transmission of genetic information to the next generation.

This article is categorized under:

Regulatory RNAs/RNAi/Riboswitches > RNAi: Mechanisms of Action

Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs

RNA Interactions with Proteins and Other Molecules > RNA‐Protein Complexes

RNA Interactions with Proteins and Other Molecules > Protein‐RNA Interactions: Functional Implications

Genome-Wide Identification and Characterization of the BRD Family in Nile Tilapia (Oreochromis niloticus)

Simple Summary

Nile tilapia is a good model for genome-wide identification and examination of the expression and role of gene families. In this study, we identified 54 bromodomain genes (BRDs) divided into eight subfamilies in Nile tilapia. Phylogenetic analysis revealed a high conservation of the BRDs family in vertebrates, with BRDs expansion due to fish-specific duplications. Most of the BRDs displayed sexually dimorphic expression in the gonads at 90 and 180 dah (days after hatching), including 21 testis-dominated genes (brdt, brd4a and brd2b, etc.), and 9 ovary-dominated genes (brd3b, brd2a and kat2a, etc.). Male fish treated with JQ1 (BET subfamily inhibitor) displayed abnormal spermatogenesis. The numbers of germ cells were reduced and the expression of steroidogenic enzyme genes was downregulated, while the expression of apoptosis-promoting genes was elevated in the testes of treated fish.

Abstract

The bromodomain (BRD) proteins specifically recognize the N-acetyllysine motifs, which is a key event in the reading process of epigenetic marks. BRDs are evolutionarily highly conserved. Over recent years, BRDs attracted great interest because of their important roles in biological processes. However, the genome-wide identification of this family was not carried out in many animal groups, in particular, in teleosts. Moreover, the expression patterns were not reported for any of the members in this family, and the role of the BRD family was not extensively studied in fish reproduction. In this study, we identified 16 to 120 BRD genes in 24 representative species. BRDs expanded significantly in vertebrates. Phylogenetic analysis showed that the BRD family was divided into eight subfamilies (I–VIII). Transcriptome analysis showed that BRDs in Nile tilapia (Oreochromis niloticus) exhibited different expression patterns in different tissues, suggesting that these genes may play different roles in growth and development. Gonadal transcriptome analysis showed that most of the BRDs display sexually dimorphic expression in the gonads at 90 and 180 dah (days after hatching), including 21 testis-dominated genes (brdt, brd4a and brd2b, etc.), and nine ovary-dominated genes (brd3b, brd2a and kat2a, etc.). Consistent with transcriptomic data, the results of qRT-PCR and fluorescence in situ hybridization showed that brdt expression was higher in the testis than in the ovary, suggesting its critical role in the spermatogenesis of the tilapia. Male fish treated with JQ1 (BET subfamily inhibitor) displayed abnormal spermatogenesis. The numbers of germ cells were reduced, and the expression of steroidogenic enzyme genes was downregulated, while the expression of apoptosis-promoting genes was elevated in the testis tissue of treated fish. Our data provide insights into the evolution and expression of BRD genes, which is helpful for understanding their critical roles in sex differentiation and gonadal development in teleosts.

Comparing Transcriptomes Reveals Key Metabolic Mechanisms in Superior Growth Performance Nile Tilapia (Oreochromis niloticus)

Metabolic capacity is intrinsic to growth performance. To investigate superior growth performance in Nile tilapia, three full-sib families were bred and compared at the biochemical and transcriptome levels to determine metabolic mechanisms involved in significant growth differences between individuals under the same culture environment and feeding regime. Biochemical analysis showed that individuals in the higher growth group had significantly higher total protein, total triglyceride, total cholesterol, and high- and low-density lipoproteins, but significantly lower glucose, as compared with individuals in the lower growth group. Comparative transcriptome analysis showed 536 differentially expressed genes (DEGs) were upregulated, and 622 DEGs were downregulated. These genes were significantly enriched in three key pathways: the tricarboxylic acid cycle (TCA cycle), fatty acid biosynthesis and metabolism, and cholesterol biosynthesis and metabolism. Conjoint analysis of these key pathways and the biochemical parameters suggests that Nile tilapia with superior growth performance have higher ability to consume energy substrates (e.g., glucose), as well as higher ability to biosynthesize fatty acids and cholesterol. Additionally, the fatty acids biosynthesized by the superior growth performance individuals were less active in the catabolic pathway overall, but were more active in the anabolic pathway, and might be used for triglyceride biosynthesis to store excess energy in the form of fat. Furthermore, the tilapia with superior growth performance had lower ability to convert cholesterol into bile acids, but higher ability to convert it into sterols. We discuss the molecular mechanisms of the three key metabolic pathways, map the pathways, and note key factors that may impact the growth of Nile tilapia. The results provide an important guide for the artificial selection and quality enhancement of superior growth performance in tilapia.

Tissue-Specific Expression Pattern in Ancherythroculter nigrocauda, a Sexually Size Dimorphic Fish

Certain members of the Actinopterygii class are known to exhibit sexual dimorphism (SD) that results in major phenotypic differences between male and female fishes of a species. One of the most common differences between the two sexes is in body weight, a factor with a high economic value in aquaculture. In this study, we used RNA sequencing (RNA-seq) to study the liver and brain transcriptomes of Ancherythroculter nigrocauda, a fish exhibiting SD. Females attain about fourfold body weight of males at sexual maturity. Sample clustering showed that both sexes were grouped well with their sex phenotypes. In addition, 2,395 and 457 differentially expressed genes (DEGs) were identified in the liver and brain tissues, respectively. The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses predicted the association of PPAR signaling, cytochrome P450, and steroid hormone biosynthesis to the differences in sexual size. In addition, weighted gene co-expression network analyses (WGCNA) were conducted, and the green module was identified to be significantly correlated with sexual size dimorphism (SSD). Altogether, these results improve our understanding of the molecular mechanism underlying SSD in A. nigrocauda.

Identification and Characterization of Dimorphic Expression of Sex-Related Genes in Rock Bream, a Fish With Multiple Sex Chromosomes

The rock bream (Oplegnathus fasciatus) is a typical fish with a unique multiple sex chromosome system. In this study, we investigated the gene expression profiling in the gonads and brains of both males and females using RNA-Seq to identify sex-related genes and pathways. In accordance with the dimorphic expression profiles, combined with Gene ontology and KEGG enrichment analyses, a number of potential genes and pathways associated with sex determination were obtained from transcriptional analysis, especially some sex-biased genes and pathways. Next, we selected 18 candidate genes and analyzed their expression in different tissues and developmental stages. We found that the expression levels of Amh, Dmrt1, Sox9, Dmrtb1, and Nanos2 were significantly higher in the testis than those in the ovary or other tissues, whereas the expression levels of ZP4, Bouncer, RNF208, FoxH1, and TOB were significantly higher in the ovary than those in the testis. Furthermore, the expression levels of these genes in different developmental stages of gonads also showed sexually dimorphic patterns, suggesting that they might play important roles during gonadal development. These genes are useful markers for investigating sex determination and differentiation in rock bream. The findings of this study can provide insights into the molecular mechanisms of sex determination and differentiation in fish with multiple sex chromosome systems.

A Chromosome-Level Genome Assembly of Mozambique Tilapia (Oreochromis mossambicus) Reveals the Structure of Sex Determining Regions

The Mozambique tilapia (Oreochromis mossambicus) is a fascinating taxon for evolutionary and ecological research. It is an important food fish and one of the most widely distributed tilapias. Because males grow faster than females, genetically male tilapia are preferred in aquaculture. However, studies of sex determination and sex control in O. mossambicus have been hindered by the limited characterization of the genome. To address this gap, we assembled a high-quality genome of O. mossambicus, using a combination of high coverage of Illumina and Nanopore reads, coupled with Hi-C and RNA-Seq data. Our genome assembly spans 1,007 Mb with a scaffold N50 of 11.38 Mb. We successfully anchored and oriented 98.6% of the genome on 22 linkage groups (LGs). Based on re-sequencing data for male and female fishes from three families, O. mossambicus segregates both an XY system on LG14 and a ZW system on LG3. The sex-patterned SNPs shared by two XY families narrowed the sex determining regions to ∼3 Mb on LG14. The shared sex-patterned SNPs included two deleterious missense mutations in ahnak and rhbdd1, indicating the possible roles of these two genes in sex determination. This annotated chromosome-level genome assembly and identification of sex determining regions represents a valuable resource to help understand the evolution of genetic sex determination in tilapias.

Comparative transcriptomic analysis of gonadal development and renewal in the ovoviviparous black rockfish (Sebastes schlegelii)

BACKGROUND

The black rockfish (Sebastes schlegelii) has an ovoviviparous reproductive pattern and long-term sperm storage, resulting in asynchronous gonadal development between the sexes. However, the comprehensive understanding of gonadal development in black rockfish has not yet been achieved. Here, we studied gonadal development and germ cell renewal using histology and RNA-seq.

RESULTS

In this study, RNA-seq was performed on testes and ovaries to characterize key pathways and genes that are active during development and gamete maturation in black rockfish. Differentially expressed genes (DEGs) were identified and annotated in 4 comparisons (F_III vs. F_IV, F_IV vs. F_V, M_III vs. M_IV and M_IV vs. M_V). Based on analysis of DEGs enriched in the testis, 11 and 14 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were mapped to the M_III vs. M_IV group and the M_IV vs. M_V group, respectively. DEGs in ovarian development were also classified into 10 groups according to their biological functions. The expression patterns of the selected genes determined by qPCR were significantly correlated with the RNA-Seq results, supporting the reliability and accuracy of the RNA-Seq analysis. E2 levels showed down regulation from previtellogenesis to mature stage in female and T level showed down regulation from spermatogenesis to regressed stage in the male.

CONCLUSIONS

The categories "intercellular interaction and cytoskeleton", "molecule amplification" and "repair in the cell cycle" were revealed to be crucial in testis development and spermatogenesis, as was the biosynthesis of a series of metabolites. Our results provide comprehensive insight into black rockfish gonadal development and provide a basis for further study of reproductive physiology and molecular biology in ovoviviparity teleosts.

Clustering of Sex-Biased Genes and Transposable Elements in the Genome of the Medaka Fish Oryzias latipes

Although genes with similar expression patterns are sometimes found in the same genomic regions, almost nothing is known about the relative organization in genomes of genes and transposable elements (TEs), which might influence each other at the regulatory level. In this study, we used transcriptomic data from male and female gonads of the Japanese medaka Oryzias latipes to define sexually biased genes and TEs and analyze their relative genomic localization. We identified 20,588 genes expressed in the adult gonads of O. latipes. Around 39% of these genes are differentially expressed between male and female gonads. We further analyzed the expression of TEs using the program SQuIRE and showed that more TE copies are overexpressed in testis than in ovaries (36% vs. 10%, respectively). We then developed a method to detect genomic regions enriched in testis- or ovary-biased genes. This revealed that sex-biased genes and TEs are not randomly distributed in the genome and a part of them form clusters with the same expression bias. We also found a correlation of expression between TE copies and their closest genes, which increases with decreasing intervening distance. Such a genomic organization suggests either that TEs hijack the regulatory sequences of neighboring sexual genes, allowing their expression in germ line cells and consequently new insertions to be transmitted to the next generation, or that TEs are involved in the regulation of sexual genes, and might therefore through their mobility participate in the rewiring of sex regulatory networks.

Sex Determination and Differentiation in Teleost: Roles of Genetics, Environment, and Brain

The fish reproductive system is a complex biological system. Nonetheless, reproductive organ development is conserved, which starts with sex determination and then sex differentiation. The sex of a teleost is determined and differentiated from bipotential primordium by genetics, environmental factors, or both. These two processes are species-specific. There are several prominent genes and environmental factors involved during sex determination and differentiation. At the cellular level, most of the sex-determining genes suppress the female pathway. For environmental factors, there are temperature, density, hypoxia, pH, and social interaction. Once the sexual fate is determined, sex differentiation takes over the gonadal developmental process. Environmental factors involve activation and suppression of various male and female pathways depending on the sexual fate. Alongside these factors, the role of the brain during sex determination and differentiation remains elusive. Nonetheless, GnRH III knockout has promoted a male sex-biased population, which shows brain involvement during sex determination. During sex differentiation, LH and FSH might not affect the gonadal differentiation, but are required for regulating sex differentiation. This review discusses the role of prominent genes, environmental factors, and the brain in sex determination and differentiation across a few teleost species.

Formation of autotriploid Carassius auratus and its fertility-related genes analysis

Background

Formation of triploid organism is useful in genetics and breeding. In this study, autotriploid Carassius auratus (3nRR, 3n = 150) was generated from Carassius auratus red var. (RCC, 2n = 100) (♀) and autotetraploid Carassius auratus (4nRR, 4n = 200) (♂). The female 3nRR produced haploid, diploid and triploid eggs, whereas the male 3nRR was infertile. The aim of the present study was to explore fertility of potential candidate genes of 3nRR.

Results

Gonadal transcriptome profiling of four groups (3 females RCC (FRCC), 3 males 4nRR (M4nRR), 3 males 3nRR (M3nRR) and 3 females 3nRR (F3nRR)) was performed using RNA-SEq. A total of 78.90 Gb of clean short reads and 24,262 differentially expressed transcripts (DETs), including 20,155 in F3nRR vs. FRCC and 4,107 in M3nRR vs. M4nRR were identified. A total of 106 enriched pathways were identified through KEGG enrichment analysis. Out of the enriched pathways, 44 and 62 signalling pathways were identified in F3nRR vs. FRCC and M3nRR vs. M4nRR, respectively. A total of 80 and 25 potential candidate genes for fertility-related in F3nRR and M3nRR were identified, respectively, through GO, KEGG analyses and the published literature. Moreover, protein-protein interaction (PPI) network construction of these fertility-associated genes were performed. Analysis of the PPI networks showed that 6 hub genes (MYC, SOX2, BMP4, GATA4, PTEN and BMP2) were involved in female fertility of F3nRR, and 2 hub genes (TP53 and FGF2) were involved in male sterility of M3nRR.

Conclusions

Establishment of autotriploid fish offers an ideal model to study reproductive traits of triploid fish. RNA-Seq data revealed 6 genes, namely, MYC, SOX2, BMP4, GATA4, PTEN and BMP2, involved in the female fertility of the F3nRR. Moreover, 2 genes, namely, TP53 and FGF2, were related to the male sterility of the M3nRR. These findings provide information on reproduction and breeding in triploid fish.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07753-5.

Oryzias curvinotus in Sanya Does Not Contain the Male Sex-Determining Gene dmy

Hainan medaka (Oryzias curvinotus) is distributed in the coastal waters of the South China Sea and is able to adapt to a wide range of salinities. In this study, we characterized O. curvinotus in Sanya River (SY-medaka), which lacks dmy (a male sex-determining gene in O. latipes and O. curvinotus). In a comparison of SY-medaka and Gaoqiao medaka (GQ-medaka), the morphological difference between the two populations does not reach the subspecies level and they can be considered two geographic populations of O. curvinotus. A mitochondrial cytochrome oxidase subunit I (CoI) sequence alignment showed that the sequence identities between SY-medaka and other geographic populations of O. curvinotus are as high as 95%. A phylogenetic analysis of the mitochondrial genome also indicated that SY-medaka belongs to O. curvinotus. Molecular marker-based genetic sex assays and whole genome re-sequencing showed that SY-medaka does not contain dmy. Further, in RNA-Seq analyses of the testis and ovaries of sexually mature SY-medaka, dmy expression was not detected. We speculate that high temperatures resulted in the loss of dmy in SY-medaka during evolution, or the lineage has another sex-determining gene. This study provides a valuable dataset for elucidating the mechanism underlying sex determination in Oryzias genus and advances research on functional genomics or reproduction biology in O. curvinotus.

Testis Development and Differentiation in Amphibians

Sex is determined genetically in amphibians; however, little is known about the sex chromosomes, testis-determining genes, and the genes involved in testis differentiation in this class. Certain inherent characteristics of the species of this group, like the homomorphic sex chromosomes, the high diversity of the sex-determining mechanisms, or the existence of polyploids, may hinder the design of experiments when studying how the gonads can differentiate. Even so, other features, like their external development or the possibility of inducing sex reversal by external treatments, can be helpful. This review summarizes the current knowledge on amphibian sex determination, gonadal development, and testis differentiation. The analysis of this information, compared with the information available for other vertebrate groups, allows us to identify the evolutionarily conserved and divergent pathways involved in testis differentiation. Overall, the data confirm the previous observations in other vertebrates-the morphology of the adult testis is similar across different groups; however, the male-determining signal and the genetic networks involved in testis differentiation are not evolutionarily conserved.

Pomc Plays an Important Role in Sexual Size Dimorphism in Tilapia

Sexual dimorphism is common across the animal kingdom. Knowledge of the mechanisms of sexual size dimorphism is limited although it is important in biology and aquaculture. Tilapia is the common name for ~ 100 species of cichlid fish. Some are important aquaculture species and males outgrow females. To gain novel insights into the mechanisms underlying sexual size dimorphism, we analyzed the differences of brain transcriptomes between males and females in Mozambique tilapia and studied the function of the pro-opiomelanocortin (Pomc) gene in tilapia and zebrafish. The transcriptome analysis identified 123, 55, and 2706 sex-biased genes at 5, 30, and 90 dph (days post-hatch), respectively, indicating sexual dimorphism of gene expressions in the brain. The expression of Pomc in the tilapia brain was a female-biased at 30, 90, and 120 dph. An analysis of the DNA sequence located upstream of the tilapia Pomc transcriptional start site identified two estrogenic response elements. In vitro luciferase assay of the two elements revealed that β-estradiol significantly enhanced the expression of luciferase activity, suggesting that the expression of Pomc is mediated by estrogen. We knocked out Pomc in zebrafish using Crispr/Cas-9. The Pomc-knockout zebrafish showed faster growth and higher sensitivity to feeding as compared to the wild-type fish. Taken together, our results indicate that Pomc contributes to sexual size dimorphism and suggest that the high estrogen level in females promotes the expression of Pomc and suppresses feeding in female tilapias, which leads to the slower growth of female tilapias.

Sexual Development of the Hermaphroditic Scallop Argopecten irradians Revealed by Morphological, Endocrine and Molecular Analysis

Simultaneous or functional hermaphrodites possessing both ovary and testis at the same time are good materials for studying sexual development. However, previous research on sex determination and differentiation was mainly conducted in gonochoristic species and studies on simultaneous hermaphrodites are still limited. In this study, we conducted a combined morphological, endocrine and molecular study on the gonadal development of a hermaphroditic scallop Argopecten irradians aged 2–10 month old. Morphological analysis showed that sex differentiation occurred at 6 months of age. By examining the dynamic changes of progesterone, testosterone and estradiol, we found testosterone and estradiol were significantly different between the ovaries and testes almost throughout the whole process, suggesting the two hormones may be involved in scallop sex differentiation. In addition, we identified two critical sex-related genes FoxL2 and Dmrt1L, and investigated their spatiotemporal expression patterns. Results showed that FoxL2 and Dmrt1L were female- and male-biased, respectively, and mainly localized in the germ cells and follicular cells, indicating their feasibility as molecular markers for early identification of sex. Further analysis on the changes of FoxL2 and Dmrt1L expression in juveniles showed that significant sexual dimorphic expression of FoxL2 occurred at 2 months of age, earlier than that of Dmrt1L. Moreover, FoxL2 expression was significantly correlated with estradiol/testosterone ratio (E₂/T). All these results indicated that molecular sex differentiation occurs earlier than morphological sex differentiation, and FoxL2 may be a key driver that functions through regulating sex steroid hormones in the scallop. This study will deepen our understanding of the molecular mechanism underlying sex differentiation and development in spiralians.

Identification of a Goat Intersexuality-Associated Novel Variant Through Genome-Wide Resequencing and Hi-C

Background: Polled intersex syndrome (PIS) leads to reproductive disorders in goats and exerts a heavy influence on goat breeding. Since 2001, the core variant of an 11.7 kb deletion at ~129 Mb on chromosome 1 (CHI1) has been widely used as a genetic diagnostic criterion. In 2020, a ~0.48 Mb insertion within the PIS deletion was identified by sequencing in XX intersex goats. However, the suitability of this variation for the diagnosis of intersex goats worldwide and its further molecular genetic mechanism need to be clarified.

Results: The whole-genome selective sweep of intersex goats from China was performed with whole-genome next-generation sequencing technology for large sample populations and a case–control study on interbreeds. A series of candidate genes related to the goat intersexuality phenotype were found. We further confirmed that a ~0.48 Mb duplicated fragment (including ERG and KCNJ15) downstream of the ~20 Mb PIS region was reversely inserted into the PIS locus in intersex Chinese goats and was consistent with that in European Saanen and Valais black-necked goats. High-throughput chromosome conformation capture (Hi-C) technology was then used to compare the 3D structures of the PIS variant neighborhood in CHI1 between intersex and non-intersex goats. A newly found structure was validated as an intrachromosomal rearrangement. This inserted duplication changed the original spatial structure of goat CHI1 and caused the appearance of several specific loop structures in the adjacent ~20 kb downstream region of FOXL2.

Conclusions: Results suggested that the novel complex PIS variant genome was sufficient as a broad-spectrum clinical diagnostic marker of XX intersexuality in goats from Europe and China. A series of private dense loop structures caused by segment insertion into the PIS deletion might affect the expression of FOXL2 or other neighboring novel candidate genes. However, these structures require further in-depth molecular biological experimental verification. In general, this study provided new insights for future research on the molecular genetic mechanism underlying female-to-male sex reversal in goats.

High-quality chromosome-level genomes of two tilapia species reveal their evolution of repeat sequences and sex chromosomes

Tilapias are one of the most farmed fishes that are coined as "aquatic chicken" by the food industry. Nile tilapia and blue tilapia exhibit very recent transition of sex chromosome systems since their divergence approximately five million years ago, making them a great model for elucidating the molecular and evolutionary mechanisms of sex chromosome turnovers. Studies of their sex-determining pathways are also critical for developing genetic sex control in aquaculture. We report here the newly produced genomes of Nile tilapia and blue tilapia that integrate long-read sequencing and chromatin conformation data. The two nearly complete genomes have anchored over 97% of the sequences into linkage groups (LGs), and assembled majorities of complex repetitive regions including telomeres, centromeres and rDNA clusters. In particular, we inferred two episodes of repeat expansion at LG3 respectively in the ancestor of cichlids and that of tilapias. The consequential large heterochromatic region concentrated at one end of LG3 comprises tandem arrays of mRNA and small RNA genes, among which we have identified a candidate female determining gene Paics in blue tilapia. Paics shows female-specific patterns of single-nucleotide variants, copy numbers and expression patterns in gonads during early gonadogenesis. Our work provides a very important genomic resource for functional studies of cichlids, and suggested that unequal distribution of repeat content that impacts the local recombination rate might make some chromosomes more likely to become sex chromosomes.

Genomics to accelerate genetic improvement in tilapia

Selective breeding of tilapia populations started in the early 1990s and over the past three decades tilapia has become one of the most important farmed freshwater species, being produced in more than 125 countries around the globe. Although genome assemblies have been available since 2011, most of the tilapia industry still depends on classical selection techniques using mass spawning or pedigree information to select for growth traits with reported genetic gains of up to 20% per generation. The involvement of international breeding companies and research institutions has resulted in the rapid development and application of genomic resources in the last few years. GWAS and genomic selection are expected to contribute to uncovering the genetic variants involved in economically relevant traits and increasing the genetic gain in selective breeding programs, respectively. Developments over the next few years will probably focus on achieving a deep understanding of genetic architecture of complex traits, as well as accelerating genetic progress in the selection for growth-, quality- and robustness-related traits. Novel phenotyping technologies (i.e. phenomics), lower-cost whole-genome sequencing approaches, functional genomics and gene editing tools will be crucial in future developments for the improvement of tilapia aquaculture.

Comparative transcriptome analysis on four types of gonadal tissues of blotched snakehead (Channa maculata)

Blotched snakehead (Channa maculata) is an economically important freshwater fish in China, of which males grow much faster than females. To illuminate the molecular mechanism of sex differentiation and gonad development, RNA-Sequencing was performed to identify sex-related genes and pathway in gonads of 6-month-old normal XX females (XX-F), normal XY males (XY-M), XY sex reversal females (XY-F) and YY super-males (YY-M). The analysis showed that many differentially expressed genes (DEGs) had similar expression patterns in XY-F and XX-F, which were different from XY-M and YY-M. qRT-PCR indicated that Amh, Dmrt1, and Sox9 had relatively high expression in testes of XY-M and YY-M. Taking Amh as an example, there was a relative fold change of 1.0 in XX-F, 2.1 fold change in XY-F, 36.1 fold change in XY-M, and 26.0 fold change in YY-M. Cyp19a1a, Figla, and Foxl2 were highly expressive in ovaries of XX-F and XY-F. Taking Figla as an example, there was a relative fold change of 557 in XX-F, 304.5 fold change in XY-F, 5.6 fold change in XY-M, and 4.4 fold change in YY-M. KEGG analysis revealed many DEGs distributed in pathways related to sex differentiation, steroid hormone synthesis and growth, etc. Significant variation and trends in relative expression levels tested by qRT-PCR were consistent with those recorded by RNA-Sequencing. This is the first time that transcriptome of snakehead has been investigated systematically and in an integrated way. Large quantities of candidate genes involved in sex differentiation, gonad development and growth dimorphism were identified. The study provides useful resources for understanding sex differentiation and growth dimorphism, potentially assisting mono-sex production of snakehead in aquaculture.

Weighted gene correlation network analysis reveals novel biomarkers associated with mesenchymal stromal cell differentiation in early phase

Osteoporosis is a major public health problem that is associated with high morbidity and mortality, and its prevalence is increasing as the world’s population ages. Therefore, understanding the molecular basis of the disease is becoming a high priority. In this regard, studies have shown that an imbalance in adipogenic and osteogenic differentiation of bone marrow mesenchymal stem cells (MSCs) is associated with osteoporosis. In this study, we conducted a Weighted Gene Co-Expression Network Analysis to identify gene modules associated with the differentiation of bone marrow MSCs. Gene Ontology and Kyoto Encyclopedia of Genes and Genome enrichment analysis showed that the most significant module, the brown module, was enriched with genes involved in cell cycle regulation, which is in line with the initial results published using these data. In addition, the Cytoscape platform was used to identify important hub genes and lncRNAs correlated with the gene modules. Furthermore, differential gene expression analysis identified 157 and 40 genes that were upregulated and downregulated, respectively, after 3 h of MSCs differentiation. Interestingly, regulatory network analysis, and comparison of the differentially expressed genes with those in the brown module identified potential novel biomarker genes, including two transcription factors (ZNF740, FOS) and two hub genes (FOXQ1, SGK1), which were further validated for differential expression in another data set of differentiation of MSCs. Finally, Gene Set Enrichment Analysis suggested that the two most important candidate hub genes are involved in regulatory pathways, such as the JAK-STAT and RAS signaling pathways. In summary, we have revealed new molecular mechanisms of MSCs differentiation and identified novel genes that could be used as potential therapeutic targets for the treatment of osteoporosis.

De novo transcriptome assembly and sex-biased gene expression in the gonads of Amur catfish (Silurus asotus)

Amur catfish is extensively distributed and cultured in Asian countries. Despite of economic importance, the genomic information of this species remains limited. A reference transcriptome of Amur catfish was assembled and the sex-biased gene expression in the gonads was characterized using RNA-sequencing. The assembled transcriptome of Amur catfish consisted of 74,840 transcripts. The N50, mean length and max length of transcripts are 1970, 1235 and 16,748 bp. Putative sex-specific transcripts were identified and sex-specific expression of the representative genes was verified by RT-PCR. Differential expression analysis identified 5401 ovary-biased and 5618 testis-biased genes. The ovary-biased genes were mainly enriched in pathways such as RNA transport and ribosome biogenesis in eukaryotes. The testis-biased genes were enriched in calcium signaling and cytokine-cytokine receptor interaction, etc. Our data provide a valuable genomic resource for further investigating the genetic basis of sex determination, sex differentiation and sexual dimorphism of catfish.

Comparative Transcriptome Analysis of Gonads for the Identification of Sex-Related Genes in Giant Freshwater Prawns (Macrobrachium Rosenbergii) Using RNA Sequencing

The giant freshwater prawn (Macrobrachium rosenbergii) exhibits sex dimorphism between the male and female individuals. To date, the molecular mechanism governing gonadal development was unclear, and limited data were available on the gonad transcriptome of M. rosenbergii. Here, we conducted comprehensive gonadal transcriptomic analysis of female (ZW), super female (WW), and male (ZZ) M. rosenbergii for gene discovery. A total of 70.33 gigabases (Gb) of sequences were generated. There were 115,338 unigenes assembled with a mean size of 1196 base pair (bp) and N50 of 2195 bp. Alignment against the National Center for Biotechnology Information (NCBI) non-redundant nucleotide/protein sequence database (NR and NT), the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, SwissProt database, Protein family (Pfam), Gene ontology (GO), and the eukaryotic orthologous group (KOG) database, 36,282 unigenes were annotated at least in one database. Comparative transcriptome analysis observed that 10,641, 16,903, and 3393 genes were significantly differentially expressed in ZW vs. ZZ, WW vs. ZZ, and WW vs. ZW samples, respectively. Enrichment analysis of differentially expressed genes (DEGs) resulted in 268, 153, and 42 significantly enriched GO terms, respectively, and a total of 56 significantly enriched KEGG pathways. Additionally, 23 putative sex-related genes, including Gtsf1, IR, HSP21, MRPINK, Mrr, and other potentially promising candidate genes were identified. Moreover, 56,241 simple sequence repeats (SSRs) were identified. Our findings provide a valuable archive for further functional analyses of sex-related genes and future discoveries of underlying molecular mechanisms of gonadal development and sex determination.

Comparative transcriptome analysis of male and female gonads reveals sex-biased genes in spotted scat (Scatophagus argus)

Scatophagus argus is a new emerging aquaculture fish in East and Southeast Asia. To date, research on reproductive development and regulation in S. argus is lacking. Additionally, genetic and genomic information about reproduction, such as gonadal transcriptome data, is also lacking. Herein, we report the first gonadal transcriptomes of S. argus and identify genes potentially involved in reproduction and gonadal development. A total of 136,561 unigenes were obtained by sequencing of testes (n = 3) and ovaries (n = 3) at stage III. Genes upregulated in males and females known to be involved in gonadal development and gametogenesis were identified, including male-biased dmrt1, amh, gsdf, wt1a, sox9b, and nanos2, and female-biased foxl2, gdf9, bmp15, sox3, zar1, and figla. Serum estradiol-17β and 11-ketotestosterone levels were biased in female and male fish, respectively. Sexual dimorphism of serum steroid hormone levels were interpreted after expression analysis of 20 steroidogenesis-related genes, including cyp19a1a and cyp11b2. This gonadal transcript dataset will help investigate functional genes related to reproduction in S. argus.

Transcriptome of Gonads From High Temperature Induced Sex Reversal During Sex Determination and Differentiation in Chinese Tongue Sole, Cynoglossus semilaevis

The sex of Chinese tongue sole (Cynoglossus semilaevis) is determined by both genetic sex determination (GSD) and environmental sex determination (ESD), making it an ideal model to study the relationship between sex-determination and temperature. In the present study, transcriptomes of undifferentiated gonads from genetic females and males, as well as differentiated gonads from males, females, and pseudomales under high and normal temperature treatments were generated for comparative transcriptomic analysis. A mean of 68.24 M high-quality clean reads was obtained for each library. Differentially expressed genes (DEGs) between different sexes and environmental treatments were identified, revealing that the heat shock protein gene family was involved in the high temperature induced sex reversal. The Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways that were enriched in pseudomale and genetic female comparison included neuroactive ligand-receptor interaction, cortisol synthesis and secretion, and steroid hormone biosynthesis. Furthermore, weighted gene co-expression network analyses were conducted on all samples, and two modules were positive correlated with pseudomale under high temperature. An illustrated protein-protein interaction map of the module identified a hub gene, hsc70. These findings provide insights into the genetic network that is involved in sex determination and sexual differentiation, and improve our understanding of genes involved in sex reversal under high temperature.

Gonadal transcriptome analysis of the common carp, Cyprinus carpio: Identification of differentially expressed genes and SSRs

Common carp (Cyprinus carpio) is a world-wide freshwater fish of eutrophic waters. C. carpio, have various reproductive traits, including early sexual maturity, that may make them excellent, large, realistic, aquaculture model species. In the present work, de novo assembly of gonadal (testicular and ovarian) transcriptomes from juvenile common carp was performed to identify genes involved in gonadal development. A total of 81,757 and 43,257 transcripts with average lengths of 769 and 856 bp, were obtained from the immature testicular and ovarian transcriptomes, respectively. About 84,367 unigenes were constructed after removing redundancy involving representation of transcripts in both gonadal transcriptomes. Gene ontology (39,171 unigenes), clusters of orthologous group's analysis (6651 unigenes) and Kyoto encyclopedia of genes, and genomes automatic annotation server analysis (4783 unigenes) were performed to identify potential genes along with their functions. Furthermore, 18,342 (testis) and 8693 (ovary) simple sequence repeats were identified. About 298 differentially expressed genes were identified, of which 171 and 127 genes were up-regulated in testis and ovary, respectively. Quantitative real-time reverse transcription PCR was performed to validate differential expression of selected genes in testis and ovary. Nearly 809 genes related to reproduction were identified, sex-wise expression pattern of genes related to steroid synthesis, endocrine regulation, germ cell maintenance and others factors related to gonadal differentiation was observed, and expression analysis of nanos, ad4bp/sf-1, and gdf9 was performed. The present study identified certain important genes/factors involved in the gonadal development of C. carpio which may provide insights into the understanding of sex-differentiation and gonadal development processes.

Transcriptome Analysis Reveals Common and Differential Response to Low Temperature Exposure Between Tolerant and Sensitive Blue Tilapia (Oreochromis aureus)

Tilapias are very important to the world's aquaculture. As befitting fish of their tropical origin, their distribution, and culture practices are highly affected by low temperatures. In this study, we used genetic and genomic methodologies to reveal pathways involved in the response and tolerance of blue tilapia (Oreochromis aureus) to low temperature stress. Cold tolerance was characterized in 66 families of blue tilapia. Fish from cold-tolerant and cold-sensitive families were sampled at 24 and 12°C, and the transcriptional responses to low-temperature exposure were measured in the gills and liver by high-throughput mRNA sequencing. Four hundred and ninety four genes displayed a similar temperature-dependent expression in both tolerant and sensitive fish and in the two tissues, representing the core molecular response to low temperature exposure. KEGG pathway analysis of these genes revealed down-regulation of focal-adhesion and other cell-extracellular matrix (ECM) interactions, and up-regulation of proteasome and various intra-cellular proteolytic activities. Differential responses between cold-tolerant and cold-sensitive fish were found with genes and pathways that were up-regulated in one group and down-regulated in the other. This reverse response was characterized by genes involved in metabolic pathways such as glycolysis/gluconeogenesis in the gills and biosynthesis of amino-acids in the liver, with low temperature down-regulation in tolerant fish and up-regulation in sensitive fish.

FOXL2 and DMRT1L Are Yin and Yang Genes for Determining Timing of Sex Differentiation in the Bivalve Mollusk Patinopecten yessoensis

Sex determination and differentiation have long been a research hotspot in metazoans. However, little is known about when and how sex differentiation occurs in most mollusks. In this study, we conducted a combined morphological and molecular study on sex differentiation in the Yesso scallop Patinopecten yessoensis. Histological examination on gonads from 5- to 13-month-old juveniles revealed that the morphological sex differentiation occurred at 10 months of age. To determine the onset of molecular sex differentiation, molecular markers were screened for early identification of sex. The gonadal expression profiles of eight candidate genes for sex determination or differentiation showed that only two genes displayed sexually dimorphic expression, with FOXL2 being abundant in ovaries and DMRT1L in testes. In situ hybridization revealed that both of them were detected in germ cells and follicle cells. We therefore developed LOG₁₀(DMRT1L/FOXL2) for scallop sex identification and confirmed its feasibility in differentiated individuals. By tracing its changes in 5- to 13-month-old juveniles, molecular sex differentiation time was determined: some scallops differentiate early in September when they are 7 months old, and some do late in December when they are 10 months old. Two kinds of coexpression patterns were found between FOXL2 and DMRT1L: expected antagonism after differentiation and unexpected coordination before differentiation. Our results revealed that scallop sex differentiation co-occurs with the formation of follicles, and molecular sex differentiation is established prior to morphological sex differentiation. Our study will assist in a better understanding of the molecular mechanism underlying bivalve sex differentiation.