Retinoic acid homeostasis through aldh1a2 and cyp26a1 mediates meiotic entry in Nile tilapia (Oreochromis niloticus)

Micronutrient Deficiencies in Pediatric IBD: How Often, Why, and What to Do?

Inflammatory bowel disease (IBDs), including Crohn's disease (CD), and ulcerative colitis (UC) are complex diseases with a multifactorial etiology, associated with genetic, dietetic, and other environmental risk factors. Children with IBD are at increased risk for nutritional inadequacies, resulting from decreased oral intake, restrictive dietary patterns, malabsorption, enhanced nutrient loss, surgery, and medications. Follow-up of IBD children should routinely include evaluation of specific nutritional deficits and dietetic and/or supplementation strategies should be implemented in case deficiencies are detected. This narrative review focuses on the prevalence, risk factors, detection strategy, and management of micronutrient deficiencies in pediatric IBD.

Sex-specific reproductive toxicity of subacute exposure to glufosinate-ammonium in zebrafish (Danio rerio)

Glufosinate-ammonium (GLA) is a common agricultural herbicide used worldwide. It can be transported into water bodies and can persist for long periods, posing a risk to non-target aquatic organisms. In this study, adult zebrafish were exposed to GLA (0, 0.6, 6, and 60 mg/L) for 21 days to evaluate its effect on reproduction. Fecundity, offspring development, sex hormone levels, histological changes, and apoptosis in the gonads were measured, and the expression levels of genes related to the hypothalamic-pituitary-gonadal-liver (HPG-L) axis, the cell cycle, and apoptosis were examined to assess the toxic mechanisms of GLA. Higher GLA concentrations were measured in the ovaries than in the testes. Decreases in spawning count, sperm density, and motility were observed. Meanwhile, the offspring survival rate decreased, and larval offspring swimming behavior was inhibited. GLA exposure significantly increased estradiol levels in females and reduced testosterone levels in males by affecting the expression of HPG-L axis genes. Furthermore, GLA exposure induced apoptosis in gonadal cells by controlling the expression of genes involved in cell cycle regulation and apoptotic pathways. Notably, the smaller effects of GLA concentration on body weight, gonad somatic index value, gonadal cell composition, and gonadal cell apoptosis were observed in male fish than in female fish. Taken together, GLA can accumulate in the gonads and cause sex-specific alterations in the expression of genes involved in the HPG-L axis and subsequent steroidogenesis and gametogenesis, which may be responsible for GLA-induced reproductive and developmental toxicities.

Meiosis and retinoic acid in the mouse fetal gonads: An unforeseen twist

In mammals, differentiation of germ cells is crucial for sexual reproduction, involving complex signaling pathways and environmental cues defined by the somatic cells of the gonads. This review examines the long-standing model positing that all-trans retinoic acid (ATRA) acts as a meiosis-inducing substance (MIS) in the fetal ovary by inducing expression of STRA8 in female germ cells, while CYP26B1 serves as a meiosis-preventing substance (MPS) in the fetal testis by degrading ATRA and preventing STRA8 expression in the male germ cells until postnatal development. Recent genetic studies in the mouse challenge this paradigm, revealing that meiosis initiation in female germ cells can occur independently of ATRA signaling, with key roles played by other intrinsic factors like DAZL and DMRT1, and extrinsic signals such as BMPs and vitamin C. Thus, ATRA can no longer be considered as 'the' long-searched MIS. Furthermore, evidence indicates that CYP26B1 does not prevent meiosis by degrading ATRA in the fetal testis, but acts by degrading an unidentified MIS or synthesizing an equally unknown MPS. By emphasizing the necessity of genetic loss-of-function approaches to accurately delineate the roles of signaling molecules such ATRA in vivo, this chapter calls for a reevaluation of the mechanisms instructing and preventing meiosis initiation in the fetal ovary and testis, respectively. It highlights the need for further research into the molecular identities of the signals involved in these processes.

The pathogenic role of retinoid nuclear receptor signaling in cancer and metabolic syndromes

The retinoid nuclear receptor pathway, activated by the vitamin A metabolite retinoic acid, has been extensively investigated for over a century. This study has resulted in conflicting hypotheses about how the pathway regulates health and how it should be pharmaceutically manipulated. These disagreements arise from a fundamental contradiction: retinoid agonists offer clear benefits to select patients with rare bone growth disorders, acute promyelocytic leukemia, and some dermatologic diseases, yet therapeutic retinoid pathway activation frequently causes more harm than good, both through acute metabolic dysregulation and a delayed cancer-promoting effect. In this review, we discuss controlled clinical, mechanistic, and genetic data to suggest several disease settings where inhibition of the retinoid pathway may be a compelling therapeutic strategy, such as solid cancers or metabolic syndromes, and also caution against continued testing of retinoid agonists in cancer patients. Considerable evidence suggests a central role for retinoid regulation of immunity and metabolism, with therapeutic opportunities to antagonize retinoid signaling proposed in cancer, diabetes, and obesity.

Environmentally relevant estrogens impaired spermatogenesis and sexual behaviors in male and female zebrafish

Environmental estrogens (EEs) are found extensively in natural waters and negatively affect fish reproduction. Research on the reproductive toxicity of EEs mixtures in fish at environmentally relevant concentrations is scarce. In this study, adult male zebrafish were exposed for 60 days to EES (a mixture of EEs), EE2-low (5.55 ng/L, with an estrogenic potency equal to EES), and EE2-high (11.1 ng/L). After exposure, the expression levels of vtg1, vtg3, and esr1 in the livers in EES-treated fish remained unaltered, whereas they were significantly increased in EE2-treated fish. Both EE2-high and EES exposures notably reduced the gonad somatic index and sperm count. A disrupted spermatogenesis was also observed in the testes of EE2-high- and EES-exposed fish, along with an alteration in the expression of genes associated with spermatogonial proliferation (pcna, nanog), cell cycle transition (cyclinb1, cyclind1), and meiosis (aldh1a2, cyp26a1, sycp3). Both EE2 and EES significantly lowered plasma 11-ketotestosterone levels in males, likely by inhibiting the expression level of genes for its synthesis (scc, cyp17a1 and cyp11b2), and increased 17β-estradiol (E2) levels, possibly through upregulating the expression of cyp19a1a. A significant increase in tnfrsf1a expression and the tnfrsf1a/tnfrsf1b ratio in EE2-high and EES-treated males also suggests increased apoptosis via the extrinsic pathway. Further investigation showed that both EE2-high and EES diminished the sexual behavior of male fish, accompanied with reduced E2 levels in the brain and the expression of genes in the kisspeptin/gonadotropin-releasing hormone system. Interestingly, the sexual behavior of unexposed females paired with treated males was also reduced, indicating a synergistic effect. This study suggests that EES have a more severe impact on reproduction than EE2-low, and EEs could interfere not only with spermatogenesis in fish, but also with the sexual behaviors of both exposed males and their female partners, thereby leading to a more significant disruption in fish reproduction.

An Efficient Vector-Based CRISPR/Cas9 System in Zebrafish Cell Line

The clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) system has been widely applied in animals as an efficient genome editing tool. However, the technique is difficult to implement in fish cell lines partially due to the lack of efficient promoters to drive the expression of both sgRNA and the Cas9 protein within a single vector. In this study, it was indicated that the zebrafish U6 RNA polymerase III (ZFU6) promoter could efficiently induce tyrosinase (tyr) gene editing and lead to loss of retinal pigments when co-injection with Cas9 mRNA in zebrafish embryo. Furthermore, an optimized all-in-one vector for expression of the CRISPR/Cas9 system in the zebrafish fibroblast cell line (PAC2) was constructed by replacing the human U6 promoter with ZFU6 promoter, basing on the lentiCRISPRV2 system that widely applied in mammal cells. This new vector could successfully target the cellular communication network factor 2a (ctgfa) gene and demonstrated its function in the PAC2 cell. Notably, the vector could also be used to edit the endogenous EMX1 gene in the mammal 293 T cell line, implying its wide application potential. In conclusion, we established a new gene editing tool for zebrafish cell line, which could be a useful in vitro platform for high-throughput analyzing gene function in fish.

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.

Genome editing in East African cichlids and tilapias: state-of-the-art and future directions

African cichlid fishes of the Cichlidae family are a group of teleosts important for aquaculture and research. A thriving research community is particularly interested in the cichlid radiations of the East African Great Lakes. One key goal is to pinpoint genetic variation underlying phenotypic diversification, but the lack of genetic tools has precluded thorough dissection of the genetic basis of relevant traits in cichlids. Genome editing technologies are well established in teleost models like zebrafish and medaka. However, this is not the case for emerging model organisms, such as East African cichlids, where these technologies remain inaccessible to most laboratories, due in part to limited exchange of knowledge and expertise. The Cichlid Science 2022 meeting (Cambridge, UK) hosted for the first time a Genome Editing Workshop, where the community discussed recent advances in genome editing, with an emphasis on CRISPR/Cas9 technologies. Based on the workshop findings and discussions, in this review we define the state-of-the-art of cichlid genome editing, share resources and protocols, and propose new possible avenues to further expand the cichlid genome editing toolkit.

Molecular characterization and potential function of Rxrγ in gonadal differentiation of Chinese soft-shelled turtle (Pelodiscus sinensis)

Retinoid X receptor (RXR) is a member of the ligand-dependent nuclear receptor family. Previous studies revealed that RXRs are involved in reproduction in vertebrates. However, information on the function of RXRs in turtles is scarce. In this study, the Rxrγ cDNA sequence of Pelodiscus sinensis was cloned and analyzed, and a polyclonal antibody was constructed. RXRγ protein showed a positive signal in both mature and differentiated gonads of the turtle. Subsequently, the function of the Rxrγ gene in gonadal differentiation was confirmed using short interfering RNA (RNAi). The full-length cDNA sequence of the Rxrγ gene in P. sinensis was 2152 bp, encoding 407 amino acids and containing typical nuclear receptor family domains, including the DNA-binding domain (DBD), ligand-binding domain (LBD), and activation function 1 (AF1). Moreover, gonadal Ps-Rxrγ showed sexual dimorphism expression patterns in differentiated gonads. Real-time quantitative PCR results revealed that the Rxrγ gene was highly expressed in the turtle ovary. RNAi treatment increased the number of Sertoli cells in ZZ embryonic gonads. Furthermore, RNA interference upregulated Dmrt1 and Sox9 in ZZ and ZW embryonic gonads. However, Foxl2, Cyp19a1, Stra8, and Cyp26b1 were downregulated in embryonic gonads. The results indicated that Rxrγ participated in gonadal differentiation and development in P. sinensis.

Sex Determination and Ovarian Development in Reptiles and Amphibians: From Genetic Pathways to Environmental Influences

BACKGROUND

Reptiles and amphibians provide untapped potential for discovering how a diversity of genetic pathways and environmental conditions are incorporated into developmental processes that can lead to similar functional outcomes. These groups display a multitude of reproductive strategies, and whereas many attributes are conserved within groups and even across vertebrates, several aspects of sexual development show considerable variation.

SUMMARY

In this review, we focus our attention on the development of the reptilian and amphibian ovary. First, we review and describe the events leading to ovarian development, including sex determination and ovarian maturation, through a comparative lens. We then describe how these events are influenced by environmental factors, focusing on temperature and exposure to anthropogenic chemicals. Lastly, we identify critical knowledge gaps and future research directions that will be crucial to moving forward in our understanding of ovarian development and the influences of the environment in reptiles and amphibians.

KEY MESSAGES

Reptiles and amphibians provide excellent models for understanding the diversity of sex determination strategies and reproductive development. However, a greater understanding of the basic biology of these systems is necessary for deciphering the adaptive and potentially disruptive implications of embryo-by-environment interactions in a rapidly changing world.

Genetic Aspects of Micronutrients Important for Inflammatory Bowel Disease

Inflammatory bowel disease (IBD), Crohn’s disease (CD) and ulcerative colitis (UC) are complex diseases whose etiology is associated with genetic and environmental risk factors, among which are diet and gut microbiota. To date, IBD is an incurable disease and the main goal of its treatment is to reduce symptoms, prevent complications, and improve nutritional status and the quality of life. Patients with IBD usually suffer from nutritional deficiency with imbalances of specific micronutrient levels that contribute to the further deterioration of the disease. Therefore, along with medications usually used for IBD treatment, therapeutic strategies also include the supplementation of micronutrients such as vitamin D, folic acid, iron, and zinc. Micronutrient supplementation tailored according to individual needs could help patients to maintain overall health, avoid the triggering of symptoms, and support remission. The identification of individuals’ genotypes associated with the absorption, transport and metabolism of micronutrients can modify future clinical practice in IBD and enable individualized treatment. This review discusses the personalized approach with respect to genetics related to micronutrients commonly used in inflammatory bowel disease treatment.

Developmental changes in gene expression and gonad morphology during sex differentiation in Atlantic salmon (Salmo salar)

The Atlantic salmon (Salmo salar) is a globally important species for its value in fisheries and aquaculture, and as a research model. In order to characterise aspects of sex differentiation at the morphological and mRNA level in this species, the present study examined developmental changes in gonad morphology and gene expression in males and females between 0 and 79 days post hatch (dph). Morphological differentiation of the ovary (indicated by the formation of germ cell cysts) became apparent from 52 dph. By 79 dph, ovarian phenotype was evident in 100% of genotypic females. Testes remained in an undifferentiated-like state throughout the experiment, containing germ cells dispersed singularly within the gonadal region distal to the mesentery. There were no significant sex-related differences in gonad cross-section size, germ cell number or germ cell diameter during the experiment. The expression of genes involved in teleost sex differentiation (anti-müllerian hormone (amh), cytochrome P450, family 19, subfamily A, polypeptide 1a (cyp19a1a), forkhead box L2a (foxl2a), gonadal soma-derived factor (gsdf), r-spondin 1 (rspo1), sexually dimorphic on the Y chromosome (sdY)), retinoic acid-signalling (aldehyde dehydrogenase 1a2 (aldh1a2), cytochrome P450 family 26 a1 (cyp26a1), cytochrome P450 family 26 b1 (cyp26b1), t-box transcription factor 1 (tbx1a)) and neuroestrogen production (cytochrome P450, family 19, subfamily A, polypeptide 1b (cyp19a1b)) was investigated. Significant sex-related differences were observed only for the expression of amh, cyp19a1a, gsdf and sdY. In males, amh, gsdf and sdY were upregulated from 34, 59 and 44 dph respectively. In females, cyp19a1a was upregulated from 66 dph. Independent of sex, foxl2a expression was highest at 0 dph and had reduced ∼ 47-fold by the time of morphological sex differentiation at 52 dph. This study provides new insights into the timing and sequence of some physiological changes associated with sex differentiation in Atlantic salmon. These findings also reveal that some aspects of the mRNA sex differentiation pathways in Atlantic salmon are unique compared to other teleost fishes, including other salmonids.

Single nucleotide replacement in the Atlantic salmon genome using CRISPR/Cas9 and asymmetrical oligonucleotide donors

Background

New breeding technologies (NBT) using CRISPR/Cas9-induced homology directed repair (HDR) has the potential to expedite genetic improvement in aquaculture. The long generation time in Atlantic salmon makes breeding an unattractive solution to obtain homozygous mutants and improving the rates of perfect HDR in founder (F0) fish is thus required. Genome editing can represent small DNA changes down to single nucleotide replacements (SNR). This enables edits such as premature stop codons or single amino acid changes and may be used to obtain fish with traits favorable to aquaculture, e.g. disease resistance. A method for SNR has not yet been demonstrated in salmon.

Results

Using CRISPR/Cas9 and asymmetrical ODNs, we were able to perform precise SNR and introduce a premature stop codon in dnd in F0 salmon. Deep sequencing demonstrated up to 59.2% efficiency in single embryos. In addition, using the same asymmetrical ODN design, we inserted a FLAG element into slc45a2 and dnd, showing high individual perfect HDR efficiencies (up to 36.7 and 32.7%, respectively).

Conclusions

In this work, we demonstrate that precise SNR and knock-in (KI) can be performed in F0 salmon embryos using asymmetrical oligonucleotide (ODN) donors. We suggest that HDR-induced SNR can be applied as a powerful NBT, allowing efficient introgression of favorable alleles and bypassing challenges associated with traditional selective breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07823-8.

Integrated lncRNA and mRNA Transcriptome Analyses in the Ovary of Cynoglossus semilaevis Reveal Genes and Pathways Potentially Involved in Reproduction

Long non-coding RNAs (lncRNAs) have been reported to be involved in multiple biological processes. However, the roles of lncRNAs in the reproduction of half-smooth tongue sole (Cynoglossus semilaevis) are unclear, especially in the molecular regulatory mechanism driving ovarian development and ovulation. Thus, to explore the mRNA and lncRNA mechanisms regulating reproduction, we collected tongue sole ovaries in three stages for RNA sequencing. In stage IV vs. V, we identified 312 differentially expressed (DE) mRNAs and 58 DE lncRNAs. In stage V vs. VI, we identified 1,059 DE mRNAs and 187 DE lncRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that DE mRNAs were enriched in ECM-receptor interaction, oocyte meiosis and steroid hormone biosynthesis pathways. Furthermore, we carried out gene set enrichment analysis (GSEA) to identify potential reproduction related-pathways additionally, such as fatty metabolism and retinol metabolism. Based on enrichment analysis, DE mRNAs with a potential role in reproduction were selected and classified into six categories, including signal transduction, cell growth and death, immune response, metabolism, transport and catabolism, and cell junction. The interactions of DE lncRNAs and mRNAs were predicted according to antisense, cis-, and trans-regulatory mechanisms. We constructed a competing endogenous RNA (ceRNA) network. Several lncRNAs were predicted to regulate genes related to reproduction including cyp17a1, cyp19a1, mmp14, pgr, and hsd17b1. The functional enrichment analysis of these target genes of lncRNAs revealed that they were involved in several signaling pathways, such as the TGF-beta, Wnt signaling, and MAPK signaling pathways and reproduction related-pathways such as the progesterone-mediated oocyte maturation, oocyte meiosis, and GnRH signaling pathway. RT-qPCR analysis showed that two lncRNAs (XR_522278.2 and XR_522171.2) were mainly expressed in the ovary. Dual-fluorescence in situ hybridization experiments showed that both XR_522278.2 and XR_522171.2 colocalized with their target genes cyp17a1 and cyp19a1, respectively, in the follicular cell layer. The results further demonstrated that lncRNAs might be involved in the biological processes by modulating gene expression. Taken together, this study provides lncRNA profiles in the ovary of tongue sole and further insight into the role of lncRNA involvement in regulating reproduction in tongue sole.

Identification, Expression and Evolution of Short-Chain Dehydrogenases/Reductases in Nile Tilapia (Oreochromis niloticus)

The short-chain dehydrogenases/reductases (SDR) superfamily is involved in multiple physiological processes. In this study, genome-wide identification and comprehensive analysis of SDR superfamily were carried out in 29 animal species based on the latest genome databases. Overall, the number of SDR genes in animals increased with whole genome duplication (WGD), suggesting the expansion of SDRs during evolution, especially in 3R-WGD and polyploidization of teleosts. Phylogenetic analysis indicated that vertebrates SDRs were clustered into five categories: classical, extended, undefined, atypical, and complex. Moreover, tandem duplication of hpgd-a, rdh8b and dhrs13 was observed in teleosts analyzed. Additionally, tandem duplications of dhrs11-a, dhrs7a, hsd11b1b, and cbr1-a were observed in all cichlids analyzed, and tandem duplication of rdh10-b was observed in tilapiines. Transcriptome analysis of adult fish revealed that 93 SDRs were expressed in more than one tissue and 5 in one tissue only. Transcriptome analysis of gonads from different developmental stages showed that expression of 17 SDRs were sexually dimorphic with 11 higher in ovary and 6 higher in testis. The sexually dimorphic expressions of these SDRs were confirmed by in situ hybridization (ISH) and qPCR, indicating their possible roles in steroidogenesis and gonadal differentiation. Taken together, the identification and the expression data obtained in this study contribute to a better understanding of SDR superfamily evolution and functions in teleosts.

Loss of stra8 Increases Germ Cell Apoptosis but Is Still Compatible With Sperm Production in Atlantic Salmon (Salmo salar)

Entering meiosis strictly depends on stimulated by retinoic acid 8 (Stra8) gene function in mammals. This gene is missing in a number of fish species, including medaka and zebrafish, but is present in the majority of fishes, including Atlantic salmon. Here, we have examined the effects of removing stra8 on male fertility in Atlantic salmon. As in mammals, stra8 expression was restricted to germ cells in the testis, transcript levels increased during the start of puberty, and decreased when blocking the production of retinoic acid. We targeted the salmon stra8 gene with two gRNAs one of these were highly effective and produced numerous mutations in stra8, which led to a loss of wild-type (WT) stra8 expression in F0 salmon testis. In maturing stra8 crispants, the spermatogenetic tubuli were partially disorganized and displayed a sevenfold increase in germ cell apoptosis, in particular among type B spermatogonia and spermatocytes. The production of spermatogenic cysts, on the other hand, increased in maturing stra8 crispants. Gene expression analysis revealed unchanged (lin28a, ret) or reduced levels (egr1, dusp4) of transcripts associated with undifferentiated spermatogonia. Decreased expression was recorded for some genes expressed in differentiating spermatogonia including dmrt1 and ccnd2 or in spermatocytes, such as ccna1. Different from Stra8-deficient mammals, a large number of germ cells completed spermatogenesis, sperm was produced and fertilization rates were similar in WT and crispant males. While loss of stra8 increased germ cell apoptosis during salmon spermatogenesis, crispants compensated this cell loss by an elevated production of spermatogenic cysts, and were able to produce functional sperm. It appears that also in a fish species with a stra8 gene in the genome, the critical relevance this gene has attained for mammalian spermatogenesis is not yet given, although detrimental effects of the loss of stra8 were clearly visible during maturation.

Cloning and characterization of rec8 gene in orange-spotted grouper (Epinephelus coioides) and Dmrt1 regulation of rec8 promoter activity

Meiosis is a specialized type of cell division critical for gamete production during sexual reproduction in eukaryotes. The meiotic recombination protein Rec8 has been identified as an important factor in germ cell meiotic initiation in vertebrates; however, its equivalent role in teleosts is poorly characterized. In this study, we cloned and sequenced the rec8 gene from orange-spotted grouper (Epinephelus coioides). The cDNA sequence consisted of 2244 base pairs (bp), including a 5' untranslated region (UTR) of 198 bp and a 3'UTR of 284 bp. The open reading frame of grouper rec8 was 1752 bp, encoding 584 amino acids. Expression levels of rec8 were higher in the ovary, intersex gonad, and testis. A neighbor-joining phylogenetic tree based on the deduced amino acid sequence indicated a common origin for grouper and other teleost rec8 molecules. Immunohistochemistry using a polyclonal anti-Rec8 antibody localized the protein in the oogonia and primary oocytes in the ovary and in spermatogonia and spermatocytes in the intersex gonad and testis, suggesting that Rec8 may play an important role in the meiotic division and the development of grouper germ cells. In addition, we found that the transcription factor Dmrt1 increased rec8 promoter activity through the second binding site, based on dual-luciferase assays. Together, these results suggest that Rec8 plays a crucial role in meiosis and may be regulated by Dmrt1 to affect meiosis in groupers.

Approach for in vivo delivery of CRISPR/Cas system: a recent update and future prospect

The clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated (Cas) system provides a groundbreaking genetic technology that allows scientists to modify genes by targeting specific genomic sites. Due to the relative simplicity and versatility of the CRISPR/Cas system, it has been extensively applied in human genetic research as well as in agricultural applications, such as improving crops. Since the gene editing activity of the CRISPR/Cas system largely depends on the efficiency of introducing the system into cells or tissues, an efficient and specific delivery system is critical for applying CRISPR/Cas technology. However, there are still some hurdles remaining for the translatability of CRISPR/Cas system. In this review, we summarized the approaches used for the delivery of the CRISPR/Cas system in mammals, plants, and aquacultures. We further discussed the aspects of delivery that can be improved to elevate the potential for CRISPR/Cas translatability.

Differential expression patterns of the two paralogous Rec8 from Nile tilapia and their responsiveness to retinoic acid signaling

REC8 (meiotic recombination protein 8) is an essential component of meiotic cohesion complexes. Interestingly, two paralogous rec8 genes happen to exist in the stra8 (stimulated by retinoic acid gene 8)-absent fishes but not in stra8-existing fishes. Stra8 is usually considered as the prerequirement during RA (retinoic acid)-mediated meiosis initiation in mammals. However, how RA triggers meiosis in the stra8-absent fishes just like Nile tilapia (Oreochromis niloticus) remains elusive. Here we characterized the two paralogous rec8 genes in Nile tilapia (Onrec8a and Onrec8b), and investigated their expression patterns and responsiveness to RA signaling by treatment of ex vivo testicular culture and promoter luciferase reporter assay. OnRec8a and OnRec8b share 36% identity to each other and are true orthologs of REC8. Their expression was predominantly restricted to meiotic germline cells with differential spatiotemporal patterns. During spermatogenesis, OnRec8b predominantly exhibited nuclear expression in spermatocytes from 60 dah (days after hatching), while OnRec8a exhibited cytoplasmic expression from 90 dah. During oogenesis, OnRec8a was expressed from 30 dah, while OnRec8b from 90 dah. Further study shows that RA signaling could upregulate the expression of both Onrec8a and Onrec8b. Collectively, our data implies that OnRec8a and OnRec8b might have differential function during meiosis and be involved in RA-mediated meiosis program.

Crosstalk Between Retinoic Acid and Sex-Related Genes Controls Germ Cell Fate and Gametogenesis in Medaka

Sex determination (SD) is a highly diverse and complex mechanism. In vertebrates, one of the first morphological differences between the sexes is the timing of initiation of the first meiosis, where its initiation occurs first in female and later in male. Thus, SD is intimately related to the responsiveness of the germ cells to undergo meiosis in a sex-specific manner. In some vertebrates, it has been reported that the timing for meiosis entry would be under control of retinoic acid (RA), through activation of Stra8. In this study, we used a fish model species for sex determination and lacking the stra8 gene, the Japanese medaka (Oryzias latipes), to investigate the connection between RA and the sex determination pathway. Exogenous RA treatments act as a stress factor inhibiting germ cell differentiation probably by activation of dmrt1a and amh. Disruption of the RA degrading enzyme gene cyp26a1 induced precocious meiosis and oogenesis in embryos/hatchlings of female and even some males. Transcriptome analyzes of cyp26a1–/–adult gonads revealed upregulation of genes related to germ cell differentiation and meiosis, in both ovaries and testes. Our findings show that germ cells respond to RA in a stra8 independent model species. The responsiveness to RA is conferred by sex-related genes, restricting its action to the sex differentiation period in both sexes.

Expression of Vasa, Nanos2 and Sox9 during initial testicular development in Nile tilapia (Oreochromis niloticus) submitted to sex reversal

Sexual differentiation and early gonadal development are critical events in vertebrate reproduction. In this study, the initial testis development and expression of the Vasa, Nanos2 and Sox9 proteins were examined in Nile tilapia Oreochromis niloticus submitted to induced sex reversal. To that end, 150O. niloticus larvae at 5 days post-hatching (dph) were kept in nurseries with no hormonal addition (control group) and 150 larvae were kept with feed containing 17α-methyltestosterone to induce male sex reversal (treated group). Morphological sexual differentiation of Nile tilapia occurred between 21 and 25 dph and sex reversal resulted in 94% males, whereas the control group presented 53% males. During sexual differentiation, gonocytes (Gon) were the predominant germ cells, which decreased and disappeared after that stage in both groups. Undifferentiated spermatogonia (Aund) were identified at 21 dph in the control group and at 23 dph in the treated group. Differentiated spermatogonia (Adiff) were found at 23 dph in both groups. Vasa and Nanos2 occurred in Gon, Aund and Adiff and there were no significant differences between groups. Vasa-labelled Adiff increased at 50 dph in both groups and Nanos2 presented a high proportion of labelled germ cells during sampling. Sertoli cells expressed Sox9 throughout the experiment and its expression was significantly greater during sexual differentiation in the control group. The results indicate that hormonal treatment did not alter initial testis development and expression of Vasa and Nanos2 in Nile tilapia, although lower expression of Sox9 and a delay in sexual differentiation was detected in the treated group.

Physiological impact and comparison of mutant screening methods in piwil2 KO founder Nile tilapia produced by CRISPR/Cas9 system

The application of genome engineering techniques to understand the mechanisms that regulate germ cell development opens promising new avenues to develop methods to control sexual maturation and mitigate associated detrimental effects in fish. In this study, the functional role of piwil2 in primordial germ cells (PGCs) was investigated in Nile tilapia using CRISPR/Cas9 and the resultant genotypes were further explored. piwil2 is a gonad-specific and maternally deposited gene in Nile tilapia eggs which is known to play a role in repression of transposon elements and is therefore thought to be important for maintaining germline cell fate. A functional domain of piwil2, PIWI domain, was targeted by injecting Cas9 mRNA and sgRNAs into Nile tilapia embryos at 1 cell stage. Results showed 54% of injected mutant larvae had no or less putative PGCs compared to control fish, suggesting an essential role of piwil2 in survival of PGCs. The genotypic features of the different phenotypic groups were explored by next generation sequencing (NGS) and other mutant screening methods including T7 endonuclease 1 (T7E1), CRISPR/Cas-derived RNA-guided engineered nuclease (RGEN), high resolution melt curve analysis (HRMA) and fragment analysis. Linking phenotypes to genotypes in F0 was hindered by the complex mosacism and wide indel spectrum revealed by NGS and fragment analysis. This study strongly suggests the functional importance of piwil2 in PGCs survival. Further studies should focus on reducing mosaicism when using CRISPR/Cas9 system to facilitate direct functional analysis in F0.

Transcriptome analysis of Clarias magur brain and gonads suggests neuro-endocrine inhibition of milt release from captive GnRH-induced males

Transcriptome analysis of Clarias magur brain and gonads at preparatory, mature, 6 and 16 h post-GnRH injection (hpi) stages yielded 9.5 GB data with 39,738 contigs. Sequences of 45 reproductive genes were identified for the first time in C. magur along with unique and differentially expressed genes. The expression of 20 genes was validated by qRT-PCR. Upregulation of Cyp11A1, Cyp17A1 and FTZF1 genes in the 16hpi testis accompanied by the 17β-HSD3 expression indicates testosterone (T) synthesis in response to LH surge, while reduced expression of CYP11B1 suggests a high T: 11-KT ratio. It is evident by the gene expression analysis that the inhibitory neurotransmitter GABA, altered T: 11-KT, increased testicular bile acids, and oxytocin-like neuropeptide in the male brain, appear to be involved in arresting the pulsatile motion of testicular smooth muscles. The work generates important leads for an effective induced breeding strategy for silurid catfish.

Tributyltin impaired spermatogenesis and reproductive behavior in male zebrafish

Tributyltin (TBT) was reported to affect sexual behavior and gametogenesis in fish. However, the modes of action involved are largely unclear. In order to elucidate the toxicological mechanisms of TBT in reproduction, zebrafish (Danio rerio) males were exposed to TBT at concentrations of 100 and 500 ng/L for 28 days. After exposure, the sperm count of the treated fish was sharply decreased though the testis weight and gonadosomatic index remained unchanged. Moreover, reduced number of spermatogonia and spermatozoa and increased spermatocytes were observed in TBT-treated fish by histological observation and PCNA-immunostaining. Increased number of apoptotic-positive spermatocytes was also present in TBT-treated fish, indicating an enhanced apoptosis in these cells. Consistent to decreased number of spermatogonia, down-regulated expressions of genes responsible for germ cell proliferation (cyclind1 and pcna) were observed in TBT-treated fish. In contrast, TBT elevated the expressions of genes involved in meiotic entry and maintenance (aldhla2, sycp3 and dmc1) while suppressed the mRNA level of gene responsible for terminus of meiotic entry (cyp26a1), in agreement with arrested meiosis and reduced sperm count. Furthermore, TBT significantly elevated the ratios of bax/bcl-2 and tnfrsf1a/tnfrsf1b in testis, which are markers for intrinsic- and extrinsic-apoptotic pathways, consistent with the enhanced TUNEL positive signals in spermatocytes. Moreover, TBT also significantly affected the parameter of reproductive behaviors in treated fish (reflected by decreased frequency of meeting, visits and time spent in spawning area). Consistently, the expressions of genes responsible for the modulation of reproductive behaviors in brain (such as cyp19a1b, kiss2, gnrh3 and ompb) were significantly down-regulated in treated-fish. Interestingly, disrupted reproductive behaviors of untreated female fish were also observed in the present study. The present study indicated that TBT might affect the reproduction of zebrafish male by disrupting the spermatogenesis and reproductive behavior of the fish.

Indel locations are determined by template polarity in highly efficient in vivo CRISPR/Cas9-mediated HDR in Atlantic salmon

Precise gene editing such as CRISPR/Cas9-mediated homology directed repair (HDR) can increase our understanding of gene function and improve traits of importance for aquaculture. This fine-tuned technology has not been developed for farmed fish including Atlantic salmon. We performed knock-in (KI) of a FLAG element in the slc45a2 gene in salmon using sense (S), anti-sense (AS) and double-stranded (ds) oligodeoxynucleotide (ODN) templates with short (24/48/84 bp) homology arms. We show in vivo ODN integration in almost all the gene edited animals, and demonstrate perfect HDR rates up to 27% in individual F0 embryos, much higher than reported previously in any fish. HDR efficiency was dependent on template concentration, but not homology arm length. Analysis of imperfect HDR variants suggest that repair occurs by synthesis-dependent strand annealing (SDSA), as we show for the first time in any species that indel location is dependent on template polarity. Correct ODN polarity can be used to avoid 5'-indels interrupting the reading frame of an inserted sequence and be of importance for HDR template design in general.

Endocrine and local signaling interact to regulate spermatogenesis in zebrafish: follicle-stimulating hormone, retinoic acid and androgens

Retinoic acid (RA) is crucial for mammalian spermatogonia differentiation, and stimulates Stra8 expression, a gene required for meiosis. Certain fish species, including zebrafish, have lost the stra8 gene. While RA still seems important for spermatogenesis in fish, it is not known which stage(s) respond to RA or whether its effects are integrated into the endocrine regulation of spermatogenesis. In zebrafish, RA promoted spermatogonia differentiation, supported androgen-stimulated meiosis, and reduced spermatocyte and spermatid apoptosis. Follicle-stimulating hormone (Fsh) stimulated RA production. Expressing a dominant-negative RA receptor variant in germ cells clearly disturbed spermatogenesis but meiosis and spermiogenesis still took place, although sperm quality was low in 6-month-old adults. This condition also activated Leydig cells. Three months later, spermatogenesis apparently had recovered, but doubling of testis weight demonstrated hypertrophy, apoptosis/DNA damage among spermatids was high and sperm quality remained low. We conclude that RA signaling is important for zebrafish spermatogenesis but is not of crucial relevance. As Fsh stimulates androgen and RA production, germ cell-mediated, RA-dependent reduction of Leydig cell activity may form a hitherto unknown intratesticular negative-feedback loop.

Meiotic initiation in chicken germ cells is regulated by Cyp26b1 and mesonephros

Our knowledge of mechanisms involved in the meiosis of chicken germ cells is very limited. In mammalian fetal ovaries, the onset of meiosis is dependent on retinoic acid and subsequent upregulation of the Stra8 gene. To clarify the mechanism of meiotic initiation in chicken germ cells, we investigated the role of Cyp26b1, a retinoic acid-degrading enzyme. The Cyp26b1-inhibitor, ketoconazole was used to treat the ex vivo-cultured stage 36 gonads/mesonephroi. Then, the progression of meiosis was studied by histological and immunohistochemical analysis and the level of the transcript for Stra8 was evaluated by a quantitative reverse transcription-polymerase chain reaction in individual ketoconazole-treated gonads after 6 days in culture. The results revealed that meiosis was induced in both testes and right ovary upon inhibition of Cyp26b1 in the ex vivo-cultured gonads, despite downregulation of Stra8 messenger RNA in the treated gonads. Also, meiosis was observed only when mesonephros was cultured alongside the left ovary. These findings demonstrate that in chicken, Stra8 is not the only factor for the entrance into meiosis, and Cyp26b1 and mesonephros play critical regulatory roles for the sex-specific timing of meiotic initiation in birds.

The Genetics and Epigenetics of Sex Change in Fish

Fish show extraordinary sexual plasticity, changing sex naturally as part of their life cycle or reversing sex because of environmental stressors. This plasticity shows that sexual fate is not an irreversible process but the result of an ongoing tug-of-war for supremacy between male and female signaling networks. The behavioral, gonadal, and morphological changes involved in this process are well described, yet the molecular events that underpin those changes remain poorly understood. Epigenetic modifications emerge as a critical link between environmental stimuli, the onset of sex change, and subsequent maintenance of sexual phenotype. Here we synthesize current knowledge of sex change, focusing on the genetic and epigenetic processes that are likely involved in the initiation and regulation of sex change. We anticipate that better understanding of sex change in fish will shed new light on sex determination and development in vertebrates and on how environmental perturbations affect sexual fate.

Involvement of the retinoic acid signaling pathway in sex differentiation and pubertal development in the European sea bass Dicentrarchus labrax

Retinoic Acid (RA) is a vitamin A derivative present in many biological processes including embryogenesis, organ development and cell differentiation. The RA signaling pathway is essential for the onset of meiosis in tetrapods, although its role in fish reproduction needs further evidence. This study reports the expression profiles of several genes involved in this pathway during sex differentiation and the first reproductive season in European sea bass (Dicentrarchus labrax) gonads. The assessed genes are representative of several steps of the pathway including retinol transport, RA synthesis, nuclear receptors, RA transport and degradation. The study includes a synteny analysis of stra8, a tetrapod meiosis gatekeeper, in several taxa. The results show that, these genes were overexpressed during early gonad development and their expression decreased during meiosis progression in males and during vitellogenesis in females. Specifically, a decrease of cyp26a1, involved in RA degradation, together with an increase of aldh1a2 and aldh1a3, in charge of RA-synthesis, might ensure the availability of high RA levels at the time of meiosis in males and females. Moreover, the absence of stra8 in the European sea bass genome, as well as the conserved genomic neighbourhood found in other taxa, suggest a stra8 independent signaling for RA during meiosis. Taken together, our results might help to better understand the role of RA signaling in teleost gonad development.

Development of a Bicistronic Vector for the Expression of a CRISPR/Cas9-mCherry System in Fish Cell Lines

The clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) system has been widely used in animals as an efficient genome editing tool. In fish cells, the technique has been difficult to implement due to the lack of proper vectors that use active promoters to drive the expression of both small guide RNA (sgRNA) and the S. pyogenes Cas9 (spCas9) protein within a single expression platform. Until now, fish cells have been modified using co-transfection of the mRNA of both the sgRNA and the spCas9. In the present study, we describe the optimization of a new vector for the expression of a CRISPR/Cas9 system, designed to edit the genome of fish cell lines, that combines a gene reporter (mCherry), sgRNA, and spCas9 in a single vector, facilitating the study of the efficiency of piscine and non-piscine promoters. A cassette containing the zebrafish U6 RNA III polymerase (U6ZF) promoter was used for the expression of the sgRNA. The new plasmid displayed the expression of spCas9, mCherry, and sgRNA in CHSE/F fish cells. The results demonstrate the functionality of the mammalian promoter and the U6ZF promoter in fish cell lines. This is the first approach aimed at developing a unified genome editing system in fish cells using bicistronic vectors, thus creating a powerful biotechnological platform to study gene function.

The Gene Toolkit Implicated in Functional Sex in Sparidae Hermaphrodites: Inferences From Comparative Transcriptomics

Sex-biased gene expression is the mode through which sex dimorphism arises from a nearly identical genome, especially in organisms without genetic sex determination. Teleost fishes show great variations in the way the sex phenotype forms. Among them, Sparidae, that might be considered as a model family displays a remarkable diversity of reproductive modes. In this study, we sequenced and analyzed the sex-biased transcriptome in gonads and brain (the tissues with the most profound role in sexual development and reproduction) of two sparids with different reproductive modes: the gonochoristic common dentex, Dentex dentex, and the protandrous hermaphrodite gilthead seabream, Sparus aurata. Through comparative analysis with other protogynous and rudimentary protandrous sparid transcriptomes already available, we put forward common male and female-specific genes and pathways that are probably implicated in sex-maintenance in this fish family. Our results contribute to the understanding of the complex processes behind the establishment of the functional sex, especially in hermaphrodite species and set the groundwork for future experiments by providing a gene toolkit that can improve efforts to control phenotypic sex in finfish in the ever-increasingly important field of aquaculture.

Long-term exposure of xenoestrogens with environmental relevant concentrations disrupted spermatogenesis of zebrafish through altering sex hormone balance, stimulating germ cell proliferation, meiosis and enhancing apoptosis

Environmental estrogens are capable of interfering with the spermatogenesis and fertility of fish. However in natural waters, these chemicals are more likely to occur as a combination rather than a single stressor. Whether and how the mixture of xenoestrogens with environmental relevant concentrations may affect fish spermatogenesis remains largely unknown. In this study, male zebrafish adults were administered to 17alpha-ethinylestradiol (EE2) and a mixture of xenoestrogens (Mix (E2, EE2, DES, 4-t-OP, 4-NP and BPA)), with the estrogenic potency equivalent to EE2. After a 60-day exposures, elevated mRNA levels of vitellogenin 1 (vtg1) and estrogen receptor 1 (esr1) in the liver of fish in both treated groups were observed. Moreover, the plasma level of E2 declined significantly in the Mix group and the ratio of 11-KT/E2 was significantly elevated in both treated groups. Consistently, the mRNA level of P450 side-chain cleavage (scc) in the EE2 group and ovarian type aromatase (cyp19a1a) in the Mix group was significantly suppressed. In addition, decreased gonadosomatic index and sperm count in the fish of Mix group were present. Furthermore, increased number of the proliferating germ cells (such as spermatogonia and spermatocytes) was observed in the fish of both groups, suggesting a stimulated germ cell proliferation and meiosis. Accordingly, both exposures significantly up-regulated the mRNA levels of genes in mitosis (cyclinb1) and meiosis (cyp26a1 in EE2 group, aldh1a2, cyp26a1, sycp3 and spo11 in Mix). In addition, decreased number of spermatozoa and increased number of TUNEL-positive signals were present in the testis of fish in the Mix group, indicating an enhanced apoptosis. Further analyses demonstrated the significant elevated expressions of tnfrsf1a and the ratio of tnfrsf1a/tnfrsf1b in the Mix group, suggesting an elevated apoptosis in the testis of fish in the Mix group via extrinsic pathway. The present study greatly extends our understanding of the underlying mechanisms of the reproductive toxicity of xenoestrogens on fish.

Fsh stimulates Leydig cell Wnt5a production, enriching zebrafish type A spermatogonia

Follicle-stimulating hormone (Fsh) modulates vertebrate spermatogenesis by regulating somatic cell functions in the testis. We have found previously that zebrafish Fsh stimulated the differentiating proliferation of type A undifferentiated spermatogonia (Aund) in an androgen-independent manner by regulating the production of growth factors and other signaling molecules in both Sertoli (SCs) and Leydig cells (LCs). For example, Fsh triggered the release of Igf3 that subsequently activated β-catenin signaling to promote the differentiating proliferation of Aund. In the present study, we report that Fsh moreover uses the non-canonical Wnt pathway to promote the proliferation and accumulation of Aund. Initially, we found that the stimulatory effect of Fsh on the proliferation activity of Aund was further strengthened when β-catenin signaling was inhibited, resulting in an accumulation of Aund. We then showed that this Fsh-induced accumulation of Aund was associated with increased transcript levels of the non-canonical Wnt ligand, wnt5a. In situ hybridization of insl3 mRNA, a gene expressed in LCs, combined with Wnt5a immunocytochemistry identified LCs as the cellular source of Wnt5a in the adult zebrafish testis. Addition of an antagonist of Wnt5a to incubations with Fsh decreased both the proliferation activity and the relative section area occupied by Aund, while an agonist of Wnt5a increased these same parameters for Aund. Taken together, our data suggest that Fsh triggered LCs to release Wnt5a, which then promoted the proliferation and accumulation of Aund. Hence, Fsh uses non-canonical Wnt signaling to ensure the production of Aund, while also triggering β-catenin signaling via Igf3 to ensure spermatogonial differentiation.

Full genome survey and dynamics of gene expression in the greater amberjack Seriola dumerili

Background

Teleosts of the genus Seriola, commonly known as amberjacks, are of high commercial value in international markets due to their flesh quality and worldwide distribution. The Seriola species of interest to Mediterranean aquaculture is the greater amberjack (Seriola dumerili). This species holds great potential for the aquaculture industry, but in captivity, reproduction has proved to be challenging, and observed growth dysfunction hinders their domestication. Insights into molecular mechanisms may contribute to a better understanding of traits like growth and sex, but investigations to unravel the molecular background of amberjacks have begun only recently.

Findings

Illumina HiSeq sequencing generated a high-coverage greater amberjack genome sequence comprising 45 909 scaffolds. Comparative mapping to the Japanese yellowtail (Seriola quinqueriadiata) and to the model species medaka (Oryzias latipes) allowed the generation of in silico groups. Additional gonad transcriptome sequencing identified sex-biased transcripts, including known sex-determining and differentiation genes. Investigation of the muscle transcriptome of slow-growing individuals showed that transcripts involved in oxygen and gas transport were differentially expressed compared with fast/normal-growing individuals. On the other hand, transcripts involved in muscle functions were found to be enriched in fast/normal-growing individuals.

Conclusion

The present study provides the first insights into the molecular background of male and female amberjacks and of fast- and slow-growing fish. Therefore, valuable molecular resources have been generated in the form of a first draft genome and a reference transcriptome. Sex-biased genes, which may also have roles in sex determination or differentiation, and genes that may be responsible for slow growth are suggested.

Identification of critical sex-biased genes in Andrias davidianus by de novo transcriptome

The Chinese giant salamander Andrias davidianus is a protected amphibian with high nutritional and economic value. Understanding its sex determination mechanism is important for improving culture techniques and sex control in breeding. However, little information on the characterization of critical genes involved in sex is available. Herein, sequencing of ovary and test produced 40,783,222 and 46,128,902 raw reads, respectively, which were jointly assembled into 80,497 unigenes. Of these, 36,609 unigenes were annotated, of which 8907 were female-biased and 10,385 were male-biased. Several sex-related pathways were observed, including the Wnt signaling pathway. After elevated temperature and estrogen exposure, neomale and neofemale specimens were identified by a female-specific marker for the first time. RT-qPCR analysis showed the expression profile of ten selected sex-biased genes to be exhibited consistently in male and neomale and in female and neofemale, with the exception of the Amh and TfIIIa genes. Results suggested that these genes may play important roles in A. davidianus sex determination and gonad development. This provides a basis for further investigation of the molecular mechanisms of sex determination in amphibians.

Transcriptome Dynamics During Turbot Spermatogenesis Predicting the Potential Key Genes Regulating Male Germ Cell Proliferation and Maturation

Spermatogenesis is a dynamic developmental process in which spermatogonial stem cells proliferate, differentiate and mature into functional spermatozoa. These processes require an accurate gene regulation network. Here, we investigated the dynamic changes that occur during spermatogenesis through a combination of histological and transcriptome analyses of different developmental stages of the testis. We constructed 18 testis transcriptome libraries, and the average length, N50, and GC content of the unigenes were 1,795 bp; 3,240 bp and 49.25%, respectively. Differentially expressed genes (DEGs) that were related to germ cell proliferation and maturation, such as NANOS3, RARs, KIFs, steroid hormone synthesis-related genes and receptor genes, were identified between pairs of testis at different developmental stages. Gene ontology annotation and pathway analyses were conducted on DEGs with specific expression patterns involved in the regulation of spermatogenesis. Nine important pathways such as steroid hormone biosynthesis related to spermatogenesis were identified. A total of 21 modules that ranged from 49 to 7,448 genes were designed by a weighted gene co-expression network analysis. Furthermore, a total of 83 candidate miRNA were identified by computational methods. Our study provides the first transcriptomic evidence for differences in gene expression between different developmental stages of spermatogenesis in turbot (Scophthalmus maximus).

Identifying SNPs and candidate genes for three litter traits using single-step GWAS across six parities in Landrace and Large White pigs

Total number born (TNB), number born alive (NBA), and litter weight born alive (LWB) are critically important traits in pig production. The sow's parity is one of the major factors influencing litter traits. Because of monogenic or polygenic contributions and the presence of temporal gene effects in different sows' parities, it is difficult to clarify the biological and genetic background. To systematically explore the genetic mechanism of litter traits, we conducted 18 GWASs using single-step GWAS (ssGWAS) based on two breeds (908 Landrace and 1,130 Large White sow litter records) for each litter trait in different parities. A total of 300 Landrace and 300 Large White sows were genotyped by sequencing (GBS). ssGWAS was performed separately for each breed and each parity due to population stratification and temporal gene effect. In summary, we identified 80 (15 for Landrace and 65 for Large White), 227 (52 for Landrace, 175 for Large White), and 187 (34 for Landrace, 153 for Large White) single nucleotide polymorphisms (SNPs) affecting TNB, NBA, and LWB, respectively. Of them, we suggest that a total of 22 loci (SSC1: 125098202, SSC1: 117560058, SSC14: 147794697, SSC8: 84823302, SSC9: 143554876, and SSC9: 138766097 for Landrace; SSC1: 4023577, SSC1: 3859573, SSC1: 4891063, SSC16: 5197665, SSC10: 32050819, SSC13: 13552924, SSC13: 92819, SSC17: 3579607, SSC13: 196698221, SSC7: 30918403, SSC16: 46221484, SSC16: 46169204, SSC2: 41988642, SSC2: 44475457, SSC2: 42521875, and SSC7: 58411951 for Large White) are shared by TNB, NBA, and LWB. These results indicate the existence of gene temporal effect in each parity. Furthermore, our findings suggest four interesting candidate genes (FBXL7, ALDH1A2, LEPR, and DDX1) associated with litter traits in different parities that have a major effect on embryonic development progression. In conclusion, 22 crucial SNPs and four interesting candidate genes were identified for three litter traits across six parities. These findings advance our understanding of the genetic architecture of litter traits and confirm the presence of temporal gene effects in different parities. Importantly, functional validation studies for findings of particular interest are recommended in litter traits.

Efficient genome editing using CRISPR/Cas9 ribonucleoprotein approach in cultured Medaka fish cells

Gene editing with CRISPR/Cas9 is a powerful tool to study the function of target genes. Although this technology has demonstrated wide efficiency in many species, including fertilized zebrafish and medaka fish embryos when microinjected, its application to achieve efficient gene editing in cultured fish cells have met some difficulty. Here, we report an efficient and reliable approach to edit genes in cultured medaka (Oryzias latipes) fish cells using pre-formed gRNA-Cas9 ribonucleoprotein (RNP) complex. Both medaka fish haploid and diploid cells were transfected with the RNP complex by electroporation. Efficient gene editing was demonstrated by polymerase chain reaction (PCR) amplification of the target gene from genomic DNA and heteroduplex mobility assay carried out with polyacrylamide gel electrophoresis (PAGE). The heteroduplex bands caused by RNP cleavage and non-homologous end joining could be readily detected by PAGE. DNA sequencing confirmed that these heteroduplex bands contains the mutated target gene sequence. The average gene editing efficiency in haploid cells reached 50%, enabling us to generate a clonal cell line with ntrk3b gene mutation for further study. This RNP transfection method also works efficiently in diploid medaka cells, with the highest mutation efficiency of 61.5%. The specificity of this synthetic RNP CRISPR/Cas9 approach was verified by candidate off-target gene sequencing. Our result indicated that transfection of pre-formed gRNA-Cas9 RNP into fish cells is efficient and reliable to edit target genes in cultured medaka fish cells. This method will be very useful for gene function studies using cultured fish cells.

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.

The transcriptomic signature of different sexes in two protogynous hermaphrodites: Insights into the molecular network underlying sex phenotype in fish

Sex differentiation is a puzzling problem in fish due to the variety of reproductive systems and the flexibility of their sex determination mechanisms. The Sparidae, a teleost family, reflects this remarkable diversity of sexual mechanisms found in fish. Our aim was to capture the transcriptomic signature of different sexes in two protogynous hermaphrodite sparids, the common pandora Pagellus erythrinus and the red porgy Pagrus pagrus in order to shed light on the molecular network contributing to either the female or the male phenotype in these organisms. Through RNA sequencing, we investigated sex-specific differences in gene expression in both species' brains and gonads. The analysis revealed common male and female specific genes/pathways between these protogynous fish. Whereas limited sex differences found in the brain indicate a sexually plastic tissue, in contrast, the great amount of sex-biased genes observed in gonads reflects the functional divergence of the transformed tissue to either its male or female character. Α common "crew" of well-known molecular players is acting to preserve either sex identity of the gonad in these fish. Lastly, this study lays the ground for a deeper understanding of the complex process of sex differentiation in two species with an evolutionary significant reproductive system.

Diethylstilbestrol impaired oogenesis of yellow catfish juveniles through disrupting hypothalamic-pituitary-gonadal axis and germ cell development

Diethylstilbestrol (DES), a non-steroidal estrogen, has been found to cause altered germ cell development and disordered ovarian development in fish females. However, the mechanisms that might be involved are poorly understood. In this study, female juveniles of yellow catfish (Pelteobagrus fulvidraco) (120 days post-hatching) were exposed to two doses (10 and 100 ng l^-1 ) of DES for 28 days. After the endpoint of exposure, decreased ovary weight and gonadosomatic index, as well as various ovarian impairments were observed in response to DES. Besides, DES elevated the mRNA levels of vitellogenin 1 (vtg 1) and estrogen receptor 1 (esr 1) in liver and decreased 17β-estradiol level in plasma. Correspondingly, suppressed mRNA levels of the key genes in the hypothalamic-pituitary-gonadal axis (such as cyp19a1b, gnrh-II, fshβ and lhβ in brain and fshr, lhr and cyp19a1a in ovary) after DES exposure were also observed. The declined level of plasma 17β-estradiol and altered gene expressions of genes in the hypothalamic-pituitary-gonadal axis were thus supposed to be closely related to the disrupted oogenesis in DES-treated fish. Analyses further demonstrated that, higher concentration of DES elevated the expression ratio of bax/bcl-2, indicating the enhanced apoptosis occurred in ovary. Moreover, DES upregulated the expressions of genes involved in proliferation (cyclin d1 and pcna), meiotic entry (cyp26a1 and scp3) and meiotic maintenance (dmc1), resulting in arrested oogenesis in catfish. The present study greatly extended our understanding on the mechanisms underlying of reproductive toxicity of DES on fish oogenesis.

Transcriptomic analysis of the differentiating ovary of the protogynous ricefield eel Monopterus albus

BACKGROUND

The ricefield eel is a protogynous hermaphroditic Synbranchiform species that changes sex naturally from female to male, which offers an interesting model for studying gonadal (particularly ovarian) differentiation in vertebrates. In the present study, transcriptome sequencing of the gonad of ricefield eel larvae was performed to explore the molecular mechanisms underlying the ovarian differentiation and development.

RESULTS

A total of 301,267,988 clean reads were generated from cDNA libraries of gonadal tissues of ricefield eel larvae at 6, 9, 12, and 20 days post hatching (dph), which contained undifferentiated gonads, differentiating ovaries, ovaries with oogonia, and ovaries with meiotic oocytes, respectively. De-novo assembly of all the clean reads generated a total of 265,896 unigenes with a mean size of 720 bp and a N50 of 1107 bp. RT-qPCR analysis of the developmental expression of 13 gonadal development-related functional genes indicated that RNA-seq data are reliable. Transcriptome data suggest that high expression of female development-related genes and low expression of male development-related genes in the early gonads of ricefield eel larvae participate in the cascade of sex differentiation leading to the final female phenotype. The contrasting expression patterns of genes involved in retinoid acid (RA) synthesis and degradation might result in peak production of RA at 12 dph in the gonad of ricefield eel larvae, and induce molecular events responsible for the initiation of meiosis before the meiotic signs could be observed at 20 dph. In addition, only stra6 but not stra8 could be identified in gonadal transcriptome data of ricefield eel larvae, and the expression pattern of stra6 paralleled those of genes involved in RA synthesis, suggesting that stra6 may be a downstream target of RA and play a role in RA metabolism and/or meiotic initiation in the gonad of ricefield eel larvae.

CONCLUSIONS

The present study depicted the first large-scale RNA sequencing of the gonad of ricefield eel larvae, and identified many important functional genes, GO terms and KEGG pathways involved in gonadal development and germ cell meiosis. Results of the present study will facilitate future study on the ovarian differentiation of ricefield eels and other teleosts as well.

Generation of biallelic F0 mutants in medaka using the CRISPR/Cas9 system

Several animal models generated by genome editing methods develop somatic mosaic mutations including wild-type genome sequence in F0 generation because it is difficult to use editing tools at the one-cell stage. Producing complete knockout animals quickly is a great advantage in determining the function of target genes. This study investigated the generation of F0 knockout medaka using the CRISPR/Cas9 system. To determine whether this editing system induced mutations in the medaka genome at the one-cell stage, recombinant Cas9 protein, tracrRNA and crRNA for dead end (dnd), which is essential for germ cell development, were injected into one-cell stage embryos of olvas-DsRedExpress transgenic medaka. This allowed germ cells to be visualized by DsRed fluorescence. Genomic DNA extracted from embryos at the one-cell stage was analyzed by sequencing. Predictably, biallelic mutated sequence patterns in the target sites of dnd were found in the injected embryos. To investigate the phenotypes of the mutated fish, fluorescent and histological observations of germ cells were carried out using fry and adults. The mutations resulted in a complete loss of germ cells, suggesting loss of function of dnd in the injected embryos. Therefore, this system appears to be extremely effective for the production of F0 knockout medaka.

Identification of conserved genes triggering puberty in European sea bass males (Dicentrarchus labrax) by microarray expression profiling

BACKGROUND

Spermatogenesis is a complex process characterized by the activation and/or repression of a number of genes in a spatio-temporal manner. Pubertal development in males starts with the onset of the first spermatogenesis and implies the division of primary spermatogonia and their subsequent entry into meiosis. This study is aimed at the characterization of genes involved in the onset of puberty in European sea bass, and constitutes the first transcriptomic approach focused on meiosis in this species.

RESULTS

European sea bass testes collected at the onset of puberty (first successful reproduction) were grouped in stage I (resting stage), and stage II (proliferative stage). Transition from stage I to stage II was marked by an increase of 11ketotestosterone (11KT), the main fish androgen, whereas the transcriptomic study resulted in 315 genes differentially expressed between the two stages. The onset of puberty induced 1) an up-regulation of genes involved in cell proliferation, cell cycle and meiosis progression, 2) changes in genes related with reproduction and growth, and 3) a down-regulation of genes included in the retinoic acid (RA) signalling pathway. The analysis of GO-terms and biological pathways showed that cell cycle, cell division, cellular metabolic processes, and reproduction were affected, consistent with the early events that occur during the onset of puberty. Furthermore, changes in the expression of three RA nuclear receptors point at the importance of the RA-signalling pathway during this period, in agreement with its role in meiosis.

CONCLUSION

The results contribute to boost our knowledge of the early molecular and endocrine events that trigger pubertal development and the onset of spermatogenesis in fish. These include an increase in 11KT plasma levels and changes in the expression of several genes involved in cell proliferation, cell cycle progression, meiosis or RA-signalling pathway. Moreover, the results can be applied to study meiosis in this economically important fish species for Mediterranean countries, and may help to develop tools for its sustainable aquaculture.

Identification and expression of cytochrome P450 genes in the Chinese giant salamander Andrias davidianus

Cytochrome P450 (cyp) genes play vital roles in biological processes, including in metabolism of endogenous and exogenous compounds. Exogenous hormone influences on gene expression that leads to disruption of gonad development have been studied in several species, but whether exogenous hormones affect cyp genes that show sexually dimorphic expression remains to be determined. Here, we identified and characterized cyp genes from the widely-cultured Chinese giant salamander Andrias davidianus. We obtained 20 cyp genes including 11 genes with complete sequences. Phylogenetic analyses supported the classification of cyp genes similar to other vertebrates. Expression profile of female and male salamanders showed multiple cyp genes to exhibit higher expression in ovary than in testis, including cyp26a, cyp19a, cyp1a1, cyp4v2, cyp3a24, cyp2c20, cyp2d14, cyp2d15, and cyp4b; while cyp11a, cyp2b11, cyp11b1, cyp4f22, cyp2j6, cyp2k1, cyp2k4, cyp1a5 was higher in testis than in ovary. Seven sex-biased genes were detected after aromatase inhibitor injection and with exposure to high water temperatures. Cyp2k1, cyp11b1, and cyp2j6 expression were increased, while cyp26a, cyp2c20, cyp4b genes were decreased with aromatase inhibition. In ovary, cyp26a and cyp2c20 expression were significantly reduced; whereas cyp2k1, cyp2j6, and cyp2k4 expression were significantly elevated with no changes in cyp11a and cyp11b1 expression after temperature treatment. These findings provide valuable information for further study of sex differentiation mechanisms and cyp gene evolution.

Diethylstilbestrol, flutamide and their combination impaired the spermatogenesis of male adult zebrafish through disrupting HPG axis, meiosis and apoptosis

Both diethylstilbestrol (DES, an environmental estrogen) and flutamide (FLU, an anti-androgen) are found to impair spermatogenesis by disrupting hypothalamic-pituitary-gonadal (HPG) axis and altering androgen levels through different mechanisms/modes of action in fish with poorly understood underlying mechanisms. Furthermore, it is not known whether and how a combined exposure of DES and FLU has a stronger effect than the compounds alone. In this study, male zebrafish adults were exposed to DES, FLU and their combination (DES+FLU) for 30days, and their effects on histological structure and sperm count in testis, androgen level in plasma, as well as the mRNA levels of genes involved in HPG axis, meiotic regulation and apoptosis were analyzed. After exposure, DES and FLU disrupted spermatogenesis in zebrafish, and their combination resulted in even more severe impairment, indicating the inhibitory roles of these chemicals on spermatogenesis and their additive effects on zebrafish. The different regulation of vtg1 expression in the liver in response to DES and FLU further confirmed the different modes of action of these drugs. Gene expression and plasma steroid level analyses demonstrated the suppressed mRNA levels of the key genes (such as gnrh3, fshβ and lhβ in brain and dmrt1, sf1, cyp17a1 and cyp11b2 in testis) in HPG axis and decreased 11-ketotestosterone (11-KT) levels in plasma. The declined level of 11-KT was thus supposed to be closely related to the down-regulation of cyp26a1 (encoding the catabolic enzyme of retinoic acid) and suppression of genes involved in meiotic regulation (nanos1, dmc1 and sycp3). In fish exposed to DES and DES+FLU, enhanced apoptosis (elevated bax/bcl-2 expression ratio) was also observed. The suppression of meiotic regulation in response to all the exposures and enhanced apoptosis in response to DES were thus supposed to result in the spermatogenic impairment in zebrafish. The present study greatly extends our understanding on the mechanisms underlying of reproductive toxicity of environment estrogens and anti-androgens in fish.

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.

Gene editing nuclease and its application in tilapia

Gene editing nucleases including zinc-finger nucleases (ZFNs), transcription activator like effector nucleases (TALENs) and clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated (Cas) system (CRISPR/Cas9) provide powerful tools that improve our ability to understand the physiological processes and their underlying mechanisms. To date, these approaches have already been widely used to generate knockout and knockin models in a large number of species. Fishes comprise nearly half of extant vertebrate species and provide excellent models for studying many aspects of biology. In this review, we present an overview of recent advances in the use of gene editing nucleases for studies of fish species. We focus particularly on the use of TALENs and CRISPR/Cas9 genome editing for studying sex determination in tilapia.

The Reversible Sex of Gonochoristic Fish: Insights and Consequences

Fish sex reversal is a means to understand sex determination and differentiation, but it is also used to control sex in aquaculture. This review discusses sex reversal in gonochoristic fish, with the coexistence of genetic and environmental influences. The different periods of fish sensitivity to sex reversal treatments are presented with the mechanisms implicated. The old players of sex differentiation are revisited with transcriptome data and loss of function studies following hormone- or temperature-induced sex reversal. We also discuss whether cortisol is the universal mediator of sex reversal in fish due to its implication in ovarian meiosis and 11KT increase. The large plasticity in fish for sex reversal is also evident in the brain, with a reversibility existing even in adulthood. Studies on epigenetics are presented, since it links the environment, gene expression, and sex reversal, notably the association of DNA methylation in sex reversal. Manipulations with exogenous factors reverse the primary sex in many fish species under controlled conditions, but several questions arise on whether this can occur under wild conditions and what is the ecological significance. Cases of sex reversal in wild fish populations are shown and their fitness and future perspectives are discussed.

Retinoic acid and meiosis induction in adult versus embryonic gonads of medaka

In vertebrates, one of the first recognizable sex differences in embryos is the onset of meiosis, known to be regulated by retinoic acid (RA) in mammals. We investigated in medaka a possible meiotic function of RA during the embryonic sex determination (SD) period and in mature gonads. We found RA mediated transcriptional activation in germ cells of both sexes much earlier than the SD stage, however, no such activity during the critical stages of SD. In adults, expression of the RA metabolizing enzymes indicates sexually dimorphic RA levels. In testis, RA acts directly in Sertoli, Leydig and pre-meiotic germ cells. In ovaries, RA transcriptional activity is highest in meiotic oocytes. Our results show that RA plays an important role in meiosis induction and gametogenesis in adult medaka but contrary to common expectations, not for initiating the first meiosis in female germ cells at the SD stage.