Two sox9 genes on duplicated zebrafish chromosomes: expression of similar transcription activators in distinct sites

The exposure to nandrolone changed the expression of genes associated with sexual differentiation and disrupted the levels of hormones in the hypothalamic-pituitary-gonadal (HPG) axis, ultimately inducing male differentiation in zebrafish (Danio rerio)

Various anabolic androgenic steroids (AAS) are frequently detected in aquatic environments, yet the potential ecotoxicological impacts of their increased load have not received adequate attention. This study investigates the effects of exposure to the androgen nandrolone on sex differentiation, gonadal development, and the expression of hormones, enzymes, and proteins related to the hypothalamic-pituitary-gonadal (HPG) axis in zebrafish. The results indicate that exposure to nandrolone at 60 days post-fertilization promotes the maturation of spermatogonia, inhibits the development of oocytes, and induces their apoptosis. Additionally, there were dose-dependent changes in male differentiation and hormone levels within the HPG axis: in the high-dose group (5 μg/L), males constituted 73 % of the population, and levels of E2 hormones in gonadal tissues decreased by 33 %, while T hormone levels increased by 55 %. Furthermore, the expression of genes essential for ovarian differentiation were downregulated, whereas the expression of testicular differentiation-related genes were upregulated. The molecular docking results show that the drug can form hydrogen bonds with key proteins, which have high affinity, and are more likely to cause conformational changes in the proteins, thereby altering their activity. KEGG analysis suggests that the MAPK signaling pathway is likely the target pathway through which nandrolone exposure induces reproductive toxicity.

Understanding Genetic Regulation of Sex Differentiation in Hermaphroditic Fish

As a fundamental taxonomic group within vertebrates, fish represent an invaluable resource for investigating the mechanisms underlying sex determination and differentiation owing to their extensive geographical distribution and rich biodiversity. Within this biological cohort, the processes of sex determination and differentiation are intricately governed by both genetic factors and the complex interplay of environmental cues. While variations in external environmental factors, particularly temperature, can exert a modulatory influence on sex differentiation in fish to a limited degree, genetic factors remain the primary determinants of sexual traits. Hermaphroditic fish display three distinct types of sexual transitions: protandry (male to female), protogyny (female-to-male), bidirectional sex change (both directions serially). These fish, characterized by their unique reproductive strategies and sexual plasticity, serve as exemplary natural models for elucidating the mechanisms of sex differentiation and sexual transitions in fish. The present review delves into the histological dynamics during gonadal development across three types of sequential hermaphroditic fish, meticulously delineating the pivotal characteristics at each stage, from the inception of primordial gonads to sexual specialization. Furthermore, it examines the regulatory genes and associated signaling pathways that orchestrate sex determination and differentiation. By systematically synthesizing these research advancements, this paper endeavors to offer a comprehensive and profound insight into the intricate mechanisms governing sex differentiation in sequential hermaphroditic fish.

Impact of hexafluoropropylene oxide trimer acid (HFPO-TA) on sex differentiation after exposures during different development stages of zebrafish (Danio rerio)

Hexafluoropropylene oxide trimer acid (HFPO-TA), a novel alternative to perfluorooctanoic acid (PFOA), has been widely used and ubiquitously detected in aquatic environments. However, its potential effects on sex differentiation of aquatic organisms are not well known. Therefore, in this study, zebrafish were exposed to HFPO-TA at different development stages (0-21, 21-42, and 42-63 dpf) to investigate the effects on sex differentiation and its underlying mechanisms. All three exposures to HFPO-TA resulted in the feminization of zebrafish, and the impact of Stage II was most significant. The transcription levels of key genes related to female differentiation (bpm15, cyp19a1a, esr1, vtg1, and sox9b) were up-regulated, while those of key genes related to male differentiation (dmrt1, gata4, amh, and sox9a) were down-regulated, which could lead to the feminization. In addition, it was found that the dysregulations of these genes were prolonged in adult zebrafish even through a long recovery, which could cause sex imbalance in populations. Therefore, HFPO-TA might not be a safe alternative to PFOA, and more evidences from multi- and transgenerational toxicology are warranted.

Adamts9 is required for the development of primary ovarian follicles and maintenance of female sex in zebrafish†

Previous studies have suggested that adamts9 (a disintegrin and metalloprotease with thrombospondin type-1 motifs, member 9), an extracellular matrix (ECM) metalloprotease, participates in primordial germ cell (PGC) migration and is necessary for female fertility. In this study, we found that adamts9 knockout (KO) led to reduced body size, and female-to-male sex conversion in late juvenile or adult zebrafish; however, primary sex determination was not affected in early juveniles of adamts9 KO. Overfeeding and lowering the rearing density rescued growth defects in female adamts9 KO fish but did not rescue defects in ovarian development in adamts9 KO. Delayed PGC proliferation, significantly reduced number and size of Stage IB follicles (equivalent to primary follicles) in early juveniles of adamts9 KO, and arrested development at Stage IB follicles in mid- or late-juveniles of adamts9 KO are likely causes of female infertility and sex conversion. Via RNAseq, we found significant enrichment of differentially expressed genes involved in ECM organization during sexual maturation in ovaries of wildtype fish; and significant dysregulation of these genes in adamts9 KO ovaries. RNAseq analysis also showed enrichment of inflammatory transcriptomic signatures in adult ovaries of these adamts9 KO. Taken together, our results indicate that adamts9 is critical for development of primary ovarian follicles and maintenance of female sex, and loss of adamts9 leads to defects in ovarian follicle development, female infertility, and sex conversion in late juveniles and mature adults. These results show that the ECM and extracellular metalloproteases play major roles in maintaining ovarian follicle development in zebrafish.

Histological observations and transcriptome analyses reveal the dynamic changes in the gonads of the blotched snakehead (Channa maculata) during sex differentiation and gametogenesis

BACKGROUND

Blotched snakehead (Channa maculata) displays significant sexual dimorphism, with males exhibiting faster growth rates and larger body sizes compared to females. The cultivation of the all-male population of snakeheads holds substantial economic and ecological value. Nonetheless, the intricate processes governing the development of bipotential gonads into either testis or ovary in C. maculata remain inadequately elucidated. Therefore, it is necessary to determine the critical time window of sex differentiation in C. maculata, providing a theoretical basis for sex control in production practices.

METHODS

The body length and weight of male and female C. maculata were measured at different developmental stages to reveal when sexual dimorphism in growth initially appears. Histological observations and spatiotemporal comparative transcriptome analyses were performed on ovaries and testes across various developmental stages to determine the crucial time windows for sex differentiation in each sex and the sex-related genes. Additionally, qPCR and MG2C were utilized to validate and locate sex-related genes, and levels of E2 and T were quantified to understand sex steroid synthesis.

RESULTS

Sexual dimorphism in growth became evident starting from 90 dpf. Histological observations revealed that morphological sex differentiation in females and males occurred between 20 and 25 dpf or earlier and 30-35 dpf or earlier, respectively, corresponding to the appearance of the ovarian cavity or efferent duct anlage. Transcriptome analyses revealed divergent gene expression patterns in testes and ovaries after 30 dpf. The periods of 40-60 dpf and 60-90 dpf marked the initiation of molecular sex differentiation in females and males, respectively. Male-biased genes (Sox11a, Dmrt1, Amh, Amhr2, Gsdf, Ar, Cyp17a2) likely play crucial roles in male sex differentiation and spermatogenesis, while female-biased genes (Foxl2, Cyp19a1a, Bmp15, Figla, Er) could be pivotal in ovarian differentiation and development. Numerous biological pathways linked to sex differentiation and gametogenesis were also identified. Additionally, E2 and T exhibited sexual dimorphism during sex differentiation and gonadal development. Based on these results, it is hypothesized that in C. maculata, the potential male sex differentiation pathway, Sox11a-Dmrt1-Sox9b, activates downstream sex-related genes (Amh, Amhr2, Gsdf, Ar, Cyp17a2) for testicular development, while the antagonistic pathway, Foxl2/Cyp19a1a, activates downstream sex-related genes (Bmp15, Figla, Er) for ovarian development.

CONCLUSIONS

This study provides a comprehensive overview of gonadal dynamic changes during sex differentiation and gametogenesis in C. maculata, establishing a scientific foundation for sex control in this species.

CRISPR/dCas9-Mediated DNA Methylation Editing on emx2 in Chinese Tongue Sole (Cynoglossus semilaevis) Testis Cells

DNA methylation is a key epigenetic mechanism orchestrating gene expression networks in many biological processes. Nonetheless, studying the role of specific gene methylation events in fish faces challenges. In this study, we validate the regulation of DNA methylation on empty spiracles homeobox 2 (emx2) expression with decitabine treatment in Chinese tongue sole testis cells. We used the emx2 gene as the target gene and developed a new DNA methylation editing system by fusing dnmt3a with catalytic dead Cas9 (dCas9) and demonstrated its ability for sequence-specific DNA methylation editing. Results revealed that utilizing dCas9-dnmt3a to target emx2 promoter region led to increased DNA methylation levels and decreased emx2 expression in Chinese tongue sole testis cells. More importantly, the DNA methylation editing significantly suppressed the expression of MYC proto-oncogene, bHLH transcription factor (myc), one target gene of emx2. Furthermore, we assessed the off-target effects of dCas9-dnmt3a and confirmed no significant impact on the predicted off-target gene expression. Taken together, we developed the first DNA methylation editing system in marine species and demonstrated its effective editing ability in Chinese tongue sole cells. This provides a new strategy for both epigenetic research and molecular breeding of marine species.

Gonadal transcriptome analysis of sex-biased gene and genome-wide investigation of dmrt gene family in Acanthogobius ommaturus

Acanthogobius ommaturus is one of the largest goby fish, and widely distributed in the Northwestern Pacific as an annual benthic fish. This study aims to report the gonadal transcriptome of A. ommaturus and identify differentially expressed genes (DEGs) between sexes. A total of 5460 (27.94 %) DEGs were detected from genome, with 3301 (16.89 %) biased towards males and 2159 (11.05 %) towards females. Analysis of 76 known vertebrate sex-related genes revealed multiple key genes, including the male-biased genes dmrt1 (Doublesex and Mab-3 related transcription factor 1) and amh (Anti-Mullerian Hormone), and the female-biased genes foxl2 (Forkhead Box L2) and cyp19a1a (Cytochrome P450 Aromatase 19 Subfamily A1). Furthermore, a genome-wide gene family analysis focused on the most significantly differentially expressed male-biased gene, dmrt1, was conducted using the chromosomal-level genome. Six Aodmrt genes were identified and subjected to phylogenetic and protein interaction network analyses. To validate the expression pattern, quantitative real-time PCR (qRT-PCR) was performed and compared with gonadal transcriptome data. The results showed that only dmrt1 exhibited significant male-bias, while the expression levels and sex differences of other dmrt genes in the gonads were inconclusive. Interestingly, the other dmrt genes displayed higher expression levels in other tissues, suggesting currently unknown functions. In conclusion, this study provides valuable genetic information contributing to the understanding of the sex determination mechanism of A. ommaturus and bony fish.

Direct male development in chromosomally ZZ zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish (Danio rerio), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome, or fewer than two Z chromosomes, is essential to initiate oocyte development; and without the W factor, or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

Direct Male Development in Chromosomally ZZ Zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish ( Danio rerio ), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB strain fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome or fewer than two Z chromosomes is essential to initiate oocyte development; and without the W factor or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

gjSOX9 Cloning, Expression, and Comparison with gjSOXs Family Members in Gekko japonicus

SOX9 plays a crucial role in the male reproductive system, brain, and kidneys. In this study, we firstly analyzed the complete cDNA sequence and expression patterns for SOX9 from Gekko japonicus SOX9 (gjSOX9), carried out bioinformatic analyses of physiochemical properties, structure, and phylogenetic evolution, and compared these with other members of the gjSOX family. The results indicate that gjSOX9 cDNA comprises 1895 bp with a 1482 bp ORF encoding 494aa. gjSOX9 was not only expressed in various adult tissues but also exhibited a special spatiotemporal expression pattern in gonad tissues. gjSOX9 was predicted to be a hydrophilic nucleoprotein with a characteristic HMG-Box harboring a newly identified unique sequence, "YKYQPRRR", only present in SOXE members. Among the 20 SOX9 orthologs, gjSOX9 shares the closest genetic relationships with Eublepharis macularius SOX9, Sphacrodactylus townsendi SOX9, and Hemicordylus capensis SOX9. gjSOX9 and gjSOX10 possessed identical physicochemical properties and subcellular locations and were tightly clustered with gjSOX8 in the SOXE group. Sixteen gjSOX family members were divided into six groups: SOXB, C, D, E, F, and H with gjSOX8, 9, and 10 in SOXE among 150 SOX homologs. Collectively, the available data in this study not only facilitate a deep exploration of the functions and molecular regulation mechanisms of the gjSOX9 and gjSOX families in G. japonicus but also contribute to basic research regarding the origin and evolution of SOX9 homologs or even sex-determination mode in reptiles.

Exposure to the aryl hydrocarbon receptor agonist dioxin disrupts formation of the muscle, nerves, and vasculature in the developing jaw

Human exposure to environmental pollutants can disrupt embryonic development and impact juvenile and adult health outcomes by adversely affecting cell and organ function. Notwithstanding, environmental contamination continues to increase due to industrial development, insufficient regulations, and the mobilization of pollutants as a result of extreme weather events. Dioxins are a class of structurally related persistent organic pollutants that are highly toxic, carcinogenic, and teratogenic. 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is the most potent dioxin compound and has been shown to induce toxic effects in developing organisms by activating the aryl hydrocarbon receptor (AHR), a ligand activated transcription factor targeted by multiple persistent organic pollutants. Contaminant-induced AHR activation results in malformations of the craniofacial cartilages and neurocranium; however, the mechanisms mediating these phenotypes are not well understood. In this study, we utilized the optically transparent zebrafish model to elucidate novel cellular targets and potential transcriptional targets underlying TCDD-induced craniofacial malformations. To this end, we exposed zebrafish embryos at 4 h post fertilization to TCDD and employed a mixed-methods approach utilizing immunohistochemistry staining, transgenic reporter lines, fixed and in vivo confocal imaging, and timelapse microscopy to determine the targets mediating TCDD-induced craniofacial phenotypes. Our data indicate that embryonic TCDD exposure reduced jaw and pharyngeal arch Sox10+ chondrocytes and Tcf21+ pharyngeal mesoderm progenitors. Exposure to TCDD correspondingly led to a reduction in collagen type II deposition in Sox10+ domains. Embryonic TCDD exposure impaired development of tissues derived from or guided by Tcf21+ progenitors, namely: nerves, muscle, and vasculature. Specifically, TCDD exposure disrupted development of the hyoid and mandibular arch muscles, decreased neural innervation of the jaw, resulted in compression of cranial nerves V and VII, and led to jaw vasculature malformations. Collectively, these findings reveal novel structural targets and potential transcriptional targets of TCDD-induced toxicity, showcasing how contaminant exposures lead to congenital craniofacial malformations.

Pioneer statoacoustic neurons guide neuroblast behaviour during otic ganglion assembly

Cranial ganglia are aggregates of sensory neurons that mediate distinct types of sensation. The statoacoustic ganglion (SAG) develops into several lobes that are spatially arranged to connect appropriately with hair cells of the inner ear. To investigate the cellular behaviours involved in the 3D organization of the SAG, we use high-resolution confocal imaging of single-cell, labelled zebrafish neuroblasts (NBs), photoconversion, photoablation, and genetic perturbations. We show that otic NBs delaminate out of the otic epithelium in an epithelial-mesenchymal transition-like manner, rearranging apical polarity and primary cilia proteins. We also show that, once delaminated, NBs require RhoGTPases in order to perform active migration. Furthermore, tracking of recently delaminated NBs revealed their directed migration and coalescence around a small population of pioneer SAG neurons. These pioneer SAG neurons, not from otic placode origin, populate the coalescence region before otic neurogenesis begins and their ablation disrupts delaminated NB migratory pathways, consequentially affecting SAG shape. Altogether, this work shows for the first time the role of pioneer SAG neurons in orchestrating SAG development.

Ovary of Zebrafish during Spawning Season: Ultrastructure and Immunohistochemical Profiles of Sox9 and Myostatin

This study sought to examine the ovarian cellular and stromal components of the zebrafish (Danio rerio) throughout the spawning season using light and electron microscopic tools. The ovaries of zebrafish showed oocytes in all stages of follicular development and degeneration (atresia). Six stages of oogenesis were demonstrated: oogonia, early oocytes, late oocytes, vacuolated follicles, the yolk globule stage (vitellogenesis), and mature follicles. The SOX9 protein was expressed in the ooplasm of the primary and previtellogenic oocytes and the theca cell layer of the mature follicles. Myostatin was expressed in the granulosa and theca cells. Many stem cells in the ovarian stroma expressed myostatin and SOX9. During the spawning season, the EM results indicated that the zona radiata increased in thickness and was crossed perpendicularly by pore canals that contained processes from both oocytes and zona granulosa. The granulosa cells contained many mitochondria, rER, sER, and vesicles. Meanwhile, the thecal layer consisted of fibroblast-like cells. Atretic follicles could be demonstrated that involved both oocytes and their follicular walls. Several types of cells were distinguished in the ovarian stroma, including mast cells, telocytes, lymphocytes, fibroblasts, endocrine cells, macrophages, adipocytes, dendritic cells, and steroidogenic (stromal) cells. The ovary of the zebrafish serves as a model to investigate follicular development.

Loss of dlx5a/dlx6a Locus Alters Non-Canonical Wnt Signaling and Meckel's Cartilage Morphology

The dlx genes encode transcription factors that establish a proximal-distal polarity within neural crest cells to bestow a regional identity during craniofacial development. The expression regions of dlx paralogs are overlapping yet distinct within the zebrafish pharyngeal arches and may also be involved in progressive morphologic changes and organization of chondrocytes of the face. However, how each dlx paralog of dlx1a, dlx2a, dlx5a and dlx6a affects craniofacial development is still largely unknown. We report here that the average lengths of the Meckel's, palatoquadrate and ceratohyal cartilages in different dlx mutants were altered. Mutants for dlx5a-/- and dlx5i6-/-, where the entire dlx5a/dlx6a locus was deleted, have the shortest lengths for all three structures at 5 days post fertilization (dpf). This phenotype was also observed in 14 dpf larvae. Loss of dlx5i6 also resulted in increased proliferation of neural crest cells and expression of chondrogenic markers. Additionally, altered expression and function of non-canonical Wnt signaling were observed in these mutants suggesting a novel interaction between dlx5i6 locus and non-canonical Wnt pathway regulating ventral cartilage morphogenesis.

De novo assembly, characterization and comparative transcriptome analysis of gonads reveals sex-biased genes in Coreoperca whiteheadi

The wild Coreoperca whiteheadi is considered as the primordial species in sinipercine fish, which has valuable genetic information. Unfortunately, C. whiteheadi was listed as a near-threatened species because of the environmental pollution, over-exploitation and species invasion. Therefore, more genetic information is needed to have a better understanding of gonadal development in C. whiteheadi. Here, the first gonadal transcriptomes analysis of C. whiteheadi was conducted and 277.14 million clean reads were generated. A total of 96,753 unigenes were successfully annotated. By comparing ovary and testis transcriptomes, a total of 21,741 differentially expressed genes (DEGs) were identified, of which 12,057 were upregulated and 9684 were downregulated in testes. Among them, we also identified about 53 differentially expressed sex-biased genes. Subsequently, the expression of twenty-four DEGs were confirmed by real-time fluorescence quantitative PCR. Furthermore, the histological analysis was conducted on ovaries and testes of one-year-old C. whiteheadi. Our results provided basic support for further studies on the function of sex-biased genes and the molecular mechanism of sex determination and reproduction in C. whiteheadi.

Exposure to the persistent organic pollutant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD, dioxin) disrupts development of the zebrafish inner ear

Dioxins are a class of highly toxic and persistent environmental pollutants that have been shown through epidemiological and laboratory-based studies to act as developmental teratogens. 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), the most potent dioxin congener, has a high affinity for the aryl hydrocarbon receptor (AHR), a ligand activated transcription factor. TCDD-induced AHR activation during development impairs nervous system, cardiac, and craniofacial development. Despite the robust phenotypes previously reported, the characterization of developmental malformations and our understanding of the molecular targets mediating TCDD-induced developmental toxicity remains limited. In zebrafish, TCDD-induced craniofacial malformations are produced, in part, by the downregulation of SRY-box transcription factor 9b (sox9b), a member of the SoxE gene family. sox9b, along with fellow SoxE gene family members sox9a and sox10, have important functions in the development of the otic placode, the otic vesicle, and, ultimately, the inner ear. Given that sox9b is a known target of TCDD and that transcriptional interactions exist among SoxE genes, we asked whether TCDD exposure impaired the development of the zebrafish auditory system, specifically the otic vesicle, which gives rise to the sensory components of the inner ear. Using immunohistochemistry, in vivo confocal imaging, and time-lapse microscopy, we assessed the impact of TCDD exposure on zebrafish otic vesicle development. We found exposure resulted in structural deficits, including incomplete pillar fusion and altered pillar topography, leading to defective semicircular canal development. The observed structural deficits were accompanied by reduced collagen type II expression in the ear. Together, our findings reveal the otic vesicle as a novel target of TCDD-induced toxicity, suggest that the function of multiple SoxE genes may be affected by TCDD exposure, and provide insight into how environmental contaminants contribute to congenital malformations.

Exposure to Tetrabutylammonium Bromide Impairs Cranial Neural Crest Specification, Neurogenic Program, and Brain Morphogenesis

Tetrabutylammonium bromide (TBAB) is a widely used industrial reagent and is commonly found in our aquatic ecosystem as an industrial byproduct. In humans, the ingestion of TBAB causes severe neurological impairments and disorders such as vertigo, hallucinations, and delirium. Yet, the extent of environmental risk and TBAB toxicity to human health is poorly understood. In this study, we aim to determine the developmental toxicity of TBAB using zebrafish embryos as a model and provide novel insights into the mechanism of action of such chemicals on neurodevelopment and the overall embryonic program. Our results show that exposure to TBAB results in impaired development of the brain, inner ear, and pharyngeal skeletal elements in the zebrafish embryo. TBAB treatment resulted in aberrations in the specification of the neural crest precursors, hindbrain segmentation, and otic neurogenesis. TBAB treatment also induced a surge in apoptosis in the head, tail, and trunk regions of the developing embryo. Long-term TBAB exposure resulted in cardiac edema and craniofacial defects. Further, in silico molecular docking analysis indicated that TBAB binds to AMPA receptors and modulates neural developmental genes such as olfactomedin and acetylcholinesterase in the embryonic brain. To summarize, our study highlights the novel effects of TBAB on embryonic brain formation and segmentation, ear morphogenesis, and craniofacial skeletal development.

Identification and Functional Analysis of foxo Genes in Chinese Tongue Sole (Cynoglossus semilaevis)

The Chinese tongue sole (Cynoglossus semilaevis) is a traditional, precious fish in China. Due to the large growth difference between males and females, the investigation of their sex determination and differentiation mechanisms receives a great deal of attention. Forkhead Box O (FoxO) plays versatile roles in the regulation of sex differentiation and reproduction. Our recent transcriptomic analysis has shown that foxo genes may participate in the male differentiation and spermatogenesis of Chinese tongue sole. In this study, six Csfoxo members (Csfoxo1a, Csfoxo3a, Csfoxo3b, Csfoxo4, Csfoxo6-like, and Csfoxo1a-like) were identified. Phylogenetic analysis indicated that these six members were clustered into four groups corresponding to their denomination. The expression patterns of the gonads at different developmental stages were further analyzed. All members showed high levels of expression in the early stages (before 6 months post-hatching), and this expression was male-biased. In addition, promoter analysis found that the addition of C/EBPα and c-Jun transcription factors enhanced the transcriptional activities of Csfoxo1a, Csfoxo3a, Csfoxo3b, and Csfoxo4. The siRNA-mediated knockdown of the Csfoxo1a, Csfoxo3a, and Csfoxo3b genes in the testicular cell line of Chinese tongue sole affected the expression of genes related to sex differentiation and spermatogenesis. These results have broadened the understanding of foxo's function and provide valuable data for studying the male differentiation of tongue sole.

The non-canonical Wnt receptor Ror2 is required for cartilage cell polarity and morphogenesis of the craniofacial skeleton in zebrafish

Non-canonical/β-catenin-independent Wnt signaling plays crucial roles in tissue/cell polarity in epithelia, but its functions have been less well studied in mesenchymal tissues, such as the skeleton. Mutations in non-canonical Wnt signaling pathway genes cause human skeletal diseases such as Robinow syndrome and Brachydactyly Type B1, which disrupt bone growth throughout the endochondral skeleton. Ror2 is one of several non-canonical Wnt receptor/co-receptors. Here, we show that ror2-/- mutant zebrafish have craniofacial skeletal defects, including disruptions of chondrocyte polarity. ror1-/- mutants appear to be phenotypically wild type, but loss of both ror1 and ror2 leads to more severe cartilage defects, indicating partial redundancy. Skeletal defects in ror1/2 double mutants resemble those of wnt5b-/- mutants, suggesting that Wnt5b is the primary Ror ligand in zebrafish. Surprisingly, the proline-rich domain of Ror2, but not its kinase domain, is required to rescue its function in mosaic transgenic experiments in ror2-/- mutants. These results suggest that endochondral bone defects in ROR-related human syndromes reflect defects in cartilage polarity and morphogenesis.

Enhancement of zebrafish sperm production via a large body-sized surrogate with germ cell transplantation

Zebrafish (Danio rerio) is a commonly-used vertebrate model species for many research areas. However, its low milt volume limits effective cryopreservation of sperm from a single individual and often precludes dividing a single semen sample to conduct multiple downstream procedures such as genomic DNA/RNA extraction and in-vitro fertilization. Here, we apply germ stem cell transplantation to increase zebrafish sperm production in a closely related larger species from the same subfamily, giant danio Devario aequipinnatus. The endogenous germ cell of the host is depleted by dead-end morpholino antisense oligonucleotide. Histology of the sterile gonad and quantitative PCR of gonadal tissue reveals all sterile giant danio develop the male phenotype. Spermatogonial cells of Tg(ddx4:egfp) transgenic zebrafish are transplanted into sterile giant danio larvae, and 22% of recipients (germline chimera) produce donor-derived sperm at sexual maturation. The germline chimera produce approximately three-fold the volume of sperm and 10-fold the spermatozoon concentration of the donor. The donor-derived sperm is functional and gives rise to viable progeny upon fertilization of donor oocytes. We show that the issue of low milt volume can be effectively addressed by employing a larger surrogate parent.

Exposure to the persistent organic pollutant 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD, dioxin) disrupts development of the zebrafish inner ear

Dioxins are a class of highly toxic and persistent environmental pollutants that have been shown through epidemiological and laboratory-based studies to act as developmental teratogens. 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), the most potent dioxin congener, has a high affinity for the aryl hydrocarbon receptor (AHR), a ligand activated transcription factor. TCDD-induced AHR activation during development impairs nervous system, cardiac, and craniofacial development. Despite the robust phenotypes previously reported, the characterization of developmental malformations and our understanding of the molecular targets mediating TCDD-induced developmental toxicity remains limited. In zebrafish, TCDD-induced craniofacial malformations are produced, in part, by the downregulation of SRY-box transcription factor 9b ( sox9b ), a member of the SoxE gene family. sox9b , along with fellow SoxE gene family members sox9a and sox10 , have important functions in the development of the otic placode, the otic vesicle, and, ultimately, the inner ear. Given that sox9b in a known target of TCDD and that transcriptional interactions exist among SoxE genes, we asked whether TCDD exposure impaired the development of the zebrafish auditory system, specifically the otic vesicle, which gives rise to the sensory components of the inner ear. Using immunohistochemistry, in vivo confocal imaging, and time-lapse microscopy, we assessed the impact of TCDD exposure on zebrafish otic vesicle development. We found exposure resulted in structural deficits, including incomplete pillar fusion and altered pillar topography, leading to defective semicircular canal development. The observed structural deficits were accompanied by reduced collagen type II expression in the ear. Together, our findings reveal the otic vesicle as a novel target of TCDD-induced toxicity, suggest that the function of multiple SoxE genes may be affected by TCDD exposure, and provide insight into how environmental contaminants contribute to congenital malformations.

Highlights

The zebrafish ear is necessary to detect changes in motion, sound, and gravity.Embryos exposed to TCDD lack structural components of the developing ear.TCDD exposure impairs formation of the fusion plate and alters pillar topography.The semicircular canals of the ear are required to detect changes in movement.Following TCDD exposure embryos fail to establish semicircular canals.

Denovo RNA-Seq analysis of ovary and testis reveals potential differentially expressed transcripts associated with gonadal unsynchronization development in Onychostoma macrolepis

The Onychostoma macrolepis (O. macrolepis) is a rare and endangered wild species. Their endangered extinction might be due to their low fertility. To further illustrate the molecular mechanism of gonad development of the male and female O. macrolepis, the present study carried out de novo testicular and ovarian transcriptome sequencing. By comparing ovary and testis, 30,869 differentially expressed unigenes (9870 in female, 20999 in male) were identified. In addition, KEGG and GO analysis suggested that the Hedgehog signaling pathway have important roles in testis maintenance and spermatogenesis, whereas the Hippo signaling pathway and Wnt signaling pathway are likely to participate in ovary maintenance. RT-qPCR analysis results were consistent with transcriptome sequencing that all of gender differentiation-related genes (FOXL2, GDF9, WNT4, CYP19A1, SOX9 and GATA4), temperature-enriched genes (NOVA1, CTGF and NR4A1), clock-related genes (PER2, PER3, CRY1, CRY2, BMAL1 and CIPC) were significantly differential expression in testis compared with ovaries. Taken together, these results revealed a potential molecular mechanism that low fertility of the O. macrolepis might strong correlate with the gonadal dyssynchrony development of the male and female, which might provide theoretical basis and technical support for artificial reproduction and germplasm resource protection of the O. macrolepis.

Metzincin metalloproteases in PGC migration and gonadal sex conversion

Development of a functional gonad includes migration of primordial germ cells (PGCs), differentiations of somatic and germ cells, formation of primary follicles or spermatogenic cysts with somatic gonadal cells, development and maturation of gametes, and subsequent releasing of mature germ cells. These processes require extensive cellular and tissue remodeling, as well as broad alterations of the surrounding extracellular matrix (ECM). Metalloproteases, including MMPs (matrix metalloproteases), ADAMs (a disintegrin and metalloproteinases), and ADAMTS (a disintegrin and metalloproteinase with thrombospondin motifs), are suggested to have critical roles in the remodeling of the ECM during gonad development. However, few research articles and reviews are available on the functions and mechanisms of metalloproteases in remodeling gonadal ECM, gonadal development, or gonadal differentiation. Moreover, most studies focused on the roles of transcription and growth factors in early gonad development and primary sex determination, leaving a significant knowledge gap on how differentially expressed metalloproteases exert effects on the ECM, cell migration, development, and survival of germ cells during the development and differentiation of ovaries or testes. We will review gonad development with focus on the evidence of metalloprotease involvements, and with an emphasis on zebrafish as a model for studying gonadal sex differentiation and metalloprotease functions.

Knock-down of amh transcription by antisense RNA reduces FSH and increases follicular atresia in female Oreochromis niloticus

Anti-Müllerian hormone (Amh) plays an important role in regulating gonad development in teleosts. However, little is known about the effects of Amh on follicle development. In this study, we transfected the vector containing antisense RNA fragments of the amh gene to produce Nile tilapia, Oreochromis niloticus, with knocked-down Amh function in vivo. The results confirmed that the antisense RNA effectively inhibited amh transcription and Amh protein expression in female tilapia ovarian tissue. At 180 days of age, compared with control fish, female tilapia with knocked-down Amh function showed significantly increased growth and significantly decreased ovary weight and gonadosomatic index (P < 0.05). Female fish in the control group had ruddy-colored external genitalia, eggs extruded from the abdomen when gently squeezed, and most oocytes were developmental stage V. In contrast, the external genitalia of female fish with knocked-down Amh function did not have the ruddy color, no eggs extruded from the abdomen when squeezed, most oocytes were at developmental stages II and III, and considerable follicular atresia was apparent. At 180 days of age, the transcript levels of amhrII, cyp19a1a, foxl2 and sox9b in ovarian tissue, and the titers of luteinizing hormone, follicle stimulating hormone, and estradiol in the serum, were significantly lower in fish with knocked-down Amh function than in control fish (P < 0.05). We concluded that decreased serum hormone levels and an abnormal AMH signal delayed development and caused follicular degeneration in Nile tilapia with knocked-down Amh function. These findings show that antisense RNA is a feasible approach for gene silencing in fish, and represents an accurate and effective strategy to study gene function.

Potential Involvement of ewsr1-w Gene in Ovarian Development of Chinese Tongue Sole, Cynoglossus semilaevis

Simple Summary

Sexual dimorphism is a phenomenon commonly existing in animals. Chinese tongue sole Cynoglossus semilaevis is an economical marine fish with obvious female-biased size dimorphism. So, it is important to explore the molecular mechanism beyond gonadal development for sex control in aquaculture industry. RNA-binding protein Ewing Sarcoma protein-like (ewsr1) gene is important for mouse gonadal development and reproduction, however there are limited studies on this gene in teleost. In this study, two ewsr1 genes were cloned and characterized from C. semilaevis. The ewsr1-w gene, located in W chromosomes, showed female-biased expression during C. semilaevis gonadal development. In addition, knock-down effect and transcriptional regulation of Cs-ewsr1-w further suggested its essential role in ovarian development. This study broadened our understanding on ewsr1 function in teleost, and provided genetic resources for the further development of sex control breeding techniques in C. semilaevis aquaculture.

Abstract

Ewsr1 encodes a protein that acts as a multifunctional molecule in a variety of cellular processes. The full-length of Cs-ewsr1-w and Cs-ewsr1-z were cloned in Chinese tongue sole (Cynoglossus semilaevis). The open reading frame (ORF) of Cs-ewsr1-w was 1,767 bp that encoded 589 amino acids, while Cs-ewsr1-z was 1,794 bp that encoded 598 amino acids. Real-time PCR assays showed that Cs-ewsr1-w exhibited significant female-biased expression and could be hardly detected in male. It has the most abundant expression in ovaries among eight healthy tissues. Its expression in ovary increased gradually from 90 d to 3 y with C. semilaevis ovarian development and reached the peak at 3 y. After Cs-ewsr1-w knockdown with siRNA interference, several genes related to gonadal development including foxl2, sox9b and pou5f1 were down-regulated in ovarian cell line, suggesting the possible participation of Cs-ewsr1-w in C. semilaevis ovarian development. The dual-luciferase reporter assay revealed that the -733/-154 bp Cs-ewsr1-w promoter fragment exhibited strong transcription activity human embryonic kidney (HEK) 293T cell line. The mutation of a MAF BZIP Transcription Factor K (Mafk) binding site located in this fragment suggested that transcription factor Mafk might play an important role in Cs-ewsr1-w basal transcription. Our results will provide clues on the gene expression level, transcriptional regulation and knock-down effect of ewsr1 gene during ovarian development in teleost.

De novo gonad transcriptome analysis of elongate loach (Leptobotia elongata) provides novel insights into sex-related genes

Elongate loach (Leptobotia elongata) is endemic in middle and upper reaches of the Yangtze River in China. Because of many anthropogenic factors such as overfishing and dam construction, the loach has become a highly endangered species. So far, the genomic resources which benefit for species conservation and utilization are still lacking in elongate loach. Therefore, the first gonad transcriptome analysis of the loach was conducted in this study, providing novel insights into sex-related genes. A total of 286,800,660 clean reads with a total length of 42.02 Gb were obtained. 18,975 differentially expressed genes (DEGs) were identified, where 12,976 DEGs, especially Sox9a, Sox9b, Igf2 and Fgfr2, were upregulated in the testis, and 5999 DEGs, especially Zp3, Eg2, Plk1, Ccnb1, Cdc20 and Mos, were upregulated in the ovary. Meanwhile, some testis-specific genes (i.e. Cald1, Atp1a, Muc2 and Ca2) and ovary-specific genes (i.e. Ca4, Tuba, Acp5, Ccna, Larp6 and Nop4) were identified and verified. According to Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of DEGs, we found a series of enrichment pathways related to reproduction in elongate loach, such as the MAPK signaling pathway, oxytocin signaling pathway and oocyte meiosis pathway. Twelve DEGs were randomly selected to verify RNA-seq results by qPCR. In conclusion, this study provides a data source to study the molecular characteristics and regulatory mechanisms of sex-related genes in elongate loach, which has a potential to improve the resource protection and aquaculture production of the loach.

Chromosome-level assembly reveals a putative Y-autosomal fusion in the sex determination system of the Greenland Halibut (Reinhardtius hippoglossoides)

Despite the commercial importance of Greenland Halibut (Reinhardtius hippoglossoides), important gaps still persist in our knowledge of this species, including its reproductive biology and sex determination mechanism. Here, we combined single-molecule sequencing of long reads (Pacific Sciences) with chromatin conformation capture sequencing (Hi-C) data to assemble the first chromosome-level reference genome for this species. The high-quality assembly encompassed more than 598 Megabases (Mb) assigned to 1 594 scaffolds (scaffold N50 = 25 Mb) with 96% of its total length distributed among 24 chromosomes. Investigation of the syntenic relationship with other economically important flatfish species revealed a high conservation of synteny blocks among members of this phylogenetic clade. Sex determination analysis revealed that, similar to other teleost fishes, flatfishes also exhibit a high level of plasticity and turnover in sex-determination mechanisms. A low-coverage whole-genome sequence analysis of 198 individuals revealed that Greenland Halibut possesses a male heterogametic XY system and several putative candidate genes implied in the sex determination of this species. Our study also suggests for the first time in flatfishes that a putative Y-autosomal fusion could be associated with a reduction of recombination typical of the early steps of sex chromosome evolution.

Involvement of Transforming Growth Factor Beta Family Genes in Gonadal Differentiation in Japanese Eel, Anguilla japonica, According to Sex-Related Gene Expressions

The gonochoristic feature with environmental sex determination that occurs during the yellow stage in the eel provides an interesting model to investigate the mechanisms of gonadal development. We previously studied various sex-related genes during gonadal sex differentiation in Japanese eels. In the present study, the members of transforming growth factor beta (TGF-β) superfamily were investigated. Transcript levels of anti-Müllerian hormone, its receptor, gonadal soma-derived factor (amh, amhr2, and gsdf, respectively) measured by real-time polymerase chain reaction (qPCR) showed a strong sexual dimorphism. Transcripts were dominantly expressed in the testis, and their levels significantly increased with testicular differentiation. In contrast, the expressions of amh, amhr2, and gsdf transcripts were low in the ovary of E2-feminized female eels. In situ hybridization detected gsdf (but not amh) transcript signals in undifferentiated gonads. amh and gsdf signals were localized to Sertoli cells and had increased significantly with testicular differentiation. Weak gsdf and no amh signals were detected in early ovaries of E2-feminized female eels. Transcript levels of amh and gsdf (not amhr2) decreased during human chorionic gonadotropin (HCG)-induced spermatogenesis in males. This study suggests that amh, amhr2, and especially gsdf might be involved in the gene pathway regulating testicular differentiation of Japanese eels.

Comparative assessment of Fgf's diverse roles in inner ear development: A zebrafish perspective

Progress in understanding mechanisms of inner ear development has been remarkably rapid in recent years. The research community has benefited from the availability of several diverse model organisms, including zebrafish, chick, and mouse. The complexity of the inner ear has proven to be a challenge, and the complexity of the mammalian cochlea in particular has been the subject of intense scrutiny. Zebrafish lack a cochlea and exhibit a number of other differences from amniote species, hence they are sometimes seen as less relevant for inner ear studies. However, accumulating evidence shows that underlying cellular and molecular mechanisms are often highly conserved. As a case in point, consideration of the diverse functions of Fgf and its downstream effectors reveals many similarities between vertebrate species, allowing meaningful comparisons the can benefit the entire research community. In this review, I will discuss mechanisms by which Fgf controls key events in early otic development in zebrafish and provide direct comparisons with chick and mouse.

Autophagy coordinates chondrocyte development and early joint formation in zebrafish

Autophagy is a catabolic process responsible for the removal of waste and damaged cellular components by lysosomal degradation. It plays a key role in fundamental cell processes, including ER stress mitigation, control of cell metabolism, and cell differentiation and proliferation, all of which are essential for cartilage cell (chondrocyte) development and survival, and for the formation of cartilage. Correspondingly, autophagy dysregulation has been implicated in several skeletal disorders such as osteoarthritis and osteoporosis. To test the requirement for autophagy during skeletal development in zebrafish, we generated an atg13 CRISPR knockout zebrafish line. This line showed a complete loss of atg13 expression, and restricted autophagic activity in vivo. In the absence of autophagy, chondrocyte maturation was accelerated, with chondrocytes exhibiting signs of premature hypertrophy. Focussing on the jaw element, autophagy disruption affected joint articulation causing restricted mouth opening. This gross behavioural phenotype corresponded with a failure to thrive, and death in homozygote atg13 nulls within 17 days. Taken together, our results are consistent with autophagy contributing to the timely regulation of chondrocyte maturation and for extracellular matrix formation.

Y705 and S727 are required for the mitochondrial import and transcriptional activities of STAT3, and for regulation of stem cell proliferation

The STAT3 transcription factor, acting both in the nucleus and mitochondria, maintains embryonic stem cell pluripotency and promotes their proliferation. In this work, using zebrafish, we determined in vivo that mitochondrial STAT3 regulates mtDNA transcription in embryonic and larval stem cell niches and that this activity affects their proliferation rates. As a result, we demonstrated that import of STAT3 inside mitochondria requires Y705 phosphorylation by Jak, whereas its mitochondrial transcriptional activity, as well as its effect on proliferation, depends on the MAPK target S727. These data were confirmed using mouse embryonic stem cells: although the Y705-mutated STAT3 cannot enter mitochondria, the S727 mutation does not affect import into the organelle and is responsible for STAT3-dependent mitochondrial transcription. Surprisingly, STAT3-dependent increase of mitochondrial transcription appears to be independent from STAT3 binding to STAT3-responsive elements. Finally, loss-of-function experiments, with chemical inhibition of the JAK/STAT3 pathway or genetic ablation of stat3 gene, demonstrated that STAT3 is also required for cell proliferation in the intestine of zebrafish.

Summary: Mitochondrial import of STAT3 requires Y705 phosphorylation by Jak, whereas STAT3 mitochondrial transcriptional activity and its effect on proliferation depend on the MAPK target S727.

High temperature induced masculinization of zebrafish by down-regulation of sox9b and esr1 via DNA methylation

Elevated temperature could influence the sex differentiation by altering the expression of sex-related genes in fish. However, the underlying mechanisms by which the gene expression is altered remain poorly understood. Here, we aimed to explore the role of DNA methylation in sex differentiation of zebrafish (Danio rerio) in response to elevated temperature. The results showed that high temperature (33°C) exposure of fish from 20 to 30 days post fertilization (dpf), compared to normal temperature (28°C), resulted in male-biased sex ratio and decreased expression of female-related genes including cyp19a1a, sox9b and esr1. Meanwhile, the expressions of DNA methyltransferases dnmt3a1 and dnmt3a2, and the DNA methylation levels in sox9b and esr1 promoter were significantly increased by high temperature, strongly implying that DNA methylation is involved in high temperature-induced masculinization of zebrafish. Co-treatment with 5-aza-2'-deoxycytidine (a DNA methylation inhibitor) attenuated the high temperature-induced masculinizing effect, recovered the expression of esr1 and sox9b, suppressed the transcription of dnmt3a1 and dnmt3a2, and decreased the methylation of esr1 and sox9b promoter, further confirming that DNA methylation plays an important role in high temperature-induced masculinization of zebrafish. Furthermore, the methylation of sox9b promoter decreased the enrichment of transcription factor CREB (cAMP-responsive element binding proteins). Overall, these findings suggest that high temperature induce masculinization of zebrafish by down-regulation of female-related genes via DNA methylation, providing a new insight in understanding the epigenetic mechanism of thermal-mediated sex differentiation in fish.

Histological and transcriptional effects of androstenedione in adult zebrafish (Danio rerio)

As a natural androgen, androstenedione (AED) may pose potential risks to aquatic organisms due to its ubiquitousness in aquatic environments. Here we assessed the adverse effects of AED on histology of gonads, as well as mRNA expression levels of 34 genes concerned with hypothalamic-pituitary-gonadal (HPG) axis, germ-cell differentiation and sex differentiation in zebrafish (Danio rerio). Adult zebrafish were exposed to solvent control and three measured concentrations of 0.2, 2.3 and 23.7 μg/L AED for 60 days. The results showed that AED did not induce any obvious histological effects in the ovaries and testes. Of the investigated genes, transcriptional expression levels of amh and cyp11c1 genes in the ovaries of females were significantly increased by AED at 2.3 or 23.7 μg/L. However, different exposure concentrations of AED significantly inhibited mRNA expression of gnrh3, atf4b1 and cyp19a1b in the brain of males. In the testes of males, AED at 2.3 μg/L led to a significant induction of sox9b gene, but it at 23.7 μg/L down-regulated nr5a1b gene. These observed transcriptional changes indicated that AED could pose potential androgenic effects in zebrafish.

Systematic identification and expression analysis of the Sox gene family in spotted sea bass (Lateolabrax maculatus)

The Sox gene family encodes a set of transcription factors characterized by a conserved Sry-related high mobility group (HMG)-box domain, which performs a series of essential biological functions in diverse tissues and developmental processes. In this study, the Sox gene family was systematically characterized in spotted sea bass (Lateolabrax maculatus). A total of 26 Sox genes were identified and classified into eight subfamilies, namely, SoxB1, SoxB2, SoxC, SoxD, SoxE, SoxF, SoxH and SoxK. The phylogenetic relationship, exon-intron and domain structure analyses supported their annotation and classification. Comparison of gene copy numbers and chromosome locations among different species indicated that except tandem duplicated paralogs of Sox17/Sox32, duplicated Sox genes in spotted sea bass were generated from teleost-specific whole genome duplication during evolution. In addition, qRT-PCR was performed to detect the expression profiles of Sox genes during development and adulthood. The results showed that the expression of 16 out of 26 Sox genes was induced dramatically at different starting points after the multicellular stage, which is consistent with embryogenesis. At the early stage of sex differentiation, 9 Sox genes exhibited sexually dimorphic expression patterns, among which Sox3, Sox19 and Sox6b showed the most significant ovary-biased expression. Moreover, the distinct expression pattern of Sox genes was observed in different adult tissues. Our results provide a fundamental resource for further investigating the functions of Sox genes in embryonic processes, sex determination and differentiation as well as controlling the homeostasis of adult tissues in spotted sea bass.

Shedding new light on early sex determination in zebrafish

In contrast to established zebrafish gene annotations, the question of sex determination has still not been conclusively clarified for developing zebrafish, Danio rerio, larvae, 28 dpf or earlier. Recent studies indicate polygenic sex determination (PSD), with the genes being distributed throughout the genome. Early genetic markers of sex in zebrafish help unravel co-founding sex-related differences to apply to human health and environmental toxicity studies. A qPCR-based method was developed for six genes: cytochrome P450, family 17, subfamily A, polypeptide 1 (cyp17a1); cytochrome P450, family 19, subfamily A, polypeptide 1a (cyp19a1a); cytochrome P450, family 19, subfamily A, polypeptides 1b (cyp19a1b); vitellogenin 1 (vtg1); nuclear receptor subfamily 0, group B, member 1 (nr0b1), sry (sex-determining region Y)-box 9b (sox9b) and actin, beta 1 (actb1), the reference gene. Sry-box 9a (Sox9a), insulin-like growth factor 3 (igf3) and double sex and mab-3 related transcription factor 1 (dmrt1), which are also known to be associated with sex determination, were used in gene expression tests. Additionally, Next-Generation-Sequencing (NGS) sequenced the genome of two adult female and male and two juveniles. PCR analysis of adult zebrafish revealed sex-specific expression of cyp17a1, cyp19a1a, vtg1, igf3 and dmrt1, the first four strongly expressed in female zebrafish and the last one highly expressed in male conspecifics. From NGS, nine female and four male-fated genes were selected as novel for assessing zebrafish sex, 28 dpf. Differences in transcriptomes allowed allocation of sex-specific genes also expressed in juvenile zebrafish.

Sex determination, gonadal sex differentiation, and plasticity in vertebrate species

A diverse array of sex determination (SD) mechanisms, encompassing environmental to genetic, have been found to exist among vertebrates, covering a spectrum from fixed SD mechanisms (mammals) to functional sex change in fishes (sequential hermaphroditic fishes). A major landmark in vertebrate SD was the discovery of the SRY gene in 1990. Since that time, many attempts to clone an SRY ortholog from nonmammalian vertebrates remained unsuccessful, until 2002, when DMY/dmrt1by was discovered as the SD gene of a small fish, medaka. Surprisingly, however, DMY/dmrt1by was found in only 2 species among more than 20 species of medaka, suggesting a large diversity of SD genes among vertebrates. Considerable progress has been made over the last 3 decades, such that it is now possible to formulate reasonable paradigms of how SD and gonadal sex differentiation may work in some model vertebrate species. This review outlines our current understanding of vertebrate SD and gonadal sex differentiation, with a focus on the molecular and cellular mechanisms involved. An impressive number of genes and factors have been discovered that play important roles in testicular and ovarian differentiation. An antagonism between the male and female pathway genes exists in gonads during both sex differentiation and, surprisingly, even as adults, suggesting that, in addition to sex-changing fishes, gonochoristic vertebrates including mice maintain some degree of gonadal sexual plasticity into adulthood. Importantly, a review of various SD mechanisms among vertebrates suggests that this is the ideal biological event that can make us understand the evolutionary conundrums underlying speciation and species diversity.

A Mutation in VWA1, Encoding von Willebrand Factor A Domain-Containing Protein 1, Is Associated With Hemifacial Microsomia

Background

Hemifacial microsomia (HFM) is a type of rare congenital syndrome caused by developmental disorders of the first and second pharyngeal arches that occurs in one out of 5,600 live births. There are significant gaps in our knowledge of the pathogenic genes underlying this syndrome.

Methods

Whole exome sequencing (WES) was performed on five patients, one asymptomatic carrier, and two marry-in members of a five-generation pedigree. Structure of WARP (product of VWA1) was predicted using the Phyre2 web portal. In situ hybridization and vwa1-knockdown/knockout studies in zebrafish using morpholino and CRISPR/Cas9 techniques were performed. Cartilage staining and immunofluorescence were carried out.

Results

Through WES and a set of filtration, we identified a c.G905A:p.R302Q point mutation in a novel candidate pathogenic gene, VWA1. The Phyre2 web portal predicted alterations in secondary and tertiary structures of WARP, indicating changes in its function as well. Predictions of protein-to-protein interactions in five pathways related to craniofacial development revealed possible interactions with four proteins in the FGF pathway. Knockdown/knockout studies of the zebrafish revealed deformities of pharyngeal cartilage. A decrease of the proliferation of cranial neural crest cells (CNCCs) and alteration of the structure of pharyngeal chondrocytes were observed in the morphants as well.

Conclusion

Our data suggest that a mutation in VWA1 is functionally linked to HFM through suppression of CNCC proliferation and disruption of the organization of pharyngeal chondrocytes.

Zebrafish as an emerging model to study gonad development

The zebrafish (Danio rerio) has emerged as a popular model organism in developmental biology and pharmacogenetics due to its attribute of pathway conservation. Coupled with the availability of robust genetic and transgenic tools, transparent embryos and rapid larval development, studies of zebrafish allow detailed cellular analysis of many dynamic processes. In recent decades, the cellular and molecular mechanisms involved in the process of gonad development have been the subject of intense research using zebrafish models. In this mini-review, we give a brief overview of these studies, and highlight the essential genes involved in sex determination and gonad development in zebrafish.

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.

Sox Gene Family Revealed Genetic Variations in Autotetraploid Carassius auratus

The Sox gene family encoded transcription factors that played key roles in developmental processes in vertebrates. To further understand the evolutionary fate of the Sox gene family in teleosts, the Sox genes were comprehensively characterized in fish of different ploidy levels, including blunt snout bream (2n = 48, Megalobrama amblycephala, BSB), goldfish (2n = 100, Carassius auratus red var., 2nRCC), and autotetraploid C. auratus (4n = 200, 4nRCC). The 4nRCC, which derived from the whole genome duplication (WGD) of 2nRCC, were obtained through the distant hybridization of 2nRCC (♀) × BSB (♂). Compared with the 26 Sox genes in zebrafish (2n = 50, Danio rerio), 26, 47, and 92 putative Sox genes were identified in the BSB, 2nRCC, and 4nRCC genomes, respectively, and classified into seven subfamilies (B1, B2, C, D, E, F, and K). Comparative analyses showed that 89.36% (42/47) of Sox genes were duplicated in 2nRCC compared with those in BSB, while 97.83% (90/92) of Sox genes were duplicated in 4nRCC compared with those in 2nRCC, meaning the Sox gene family had undergone an expansion in BSB, 2nRCC, and 4nRCC, respectively, following polyploidization events. In addition, potential gene loss, genetic variations, and paternal parent SNP locus insertion occurred during the polyploidization events. Our data provided new insights into the evolution of the Sox gene family in polyploid vertebrates after several rounds of WGD events.

Spatial and temporal expression pattern of sex-related genes in ovo-testis of the self-fertilizing mangrove killifish (Kryptolebias marmoratus)

In vertebrates, sex determination and differentiation comprehend a fine balance between female and male factors, leading the bipotential anlage to develop towards ovary or testis, respectively. Nevertheless, the mangrove killifish, (Kryptolebias marmoratus) a simultaneous hermaphroditic species, could overcome those antagonistic pathways and evolved to develop and maintain reproductively active ovarian and testicular tissues in the same organ. Morphological and mRNA localization analyzes of developing and adult gonads demonstrate that genes related to testis (dmrt1 and amh) and ovary differentiation (foxl2 and sox9a) follow the same expression pattern observed in gonochoristic species, thus functioning as two independent organs. In addition, Amh expression patterns make it a strong candidate for initiation of the formation and maintenance of the testicular tissue in the hermaphroditic gonad. Differently from described so far, foxl3 seems to have an important role in oogenesis as well as spermatogenesis and gonadal structure.

Cbx2, a PcG Family Gene, Plays a Regulatory Role in Medaka Gonadal Development

Chromobox homolog 2 (CBX2), a key member of the polycomb group (PcG) family, is essential for gonadal development in mammals. A functional deficiency or genetic mutation in cbx2 can lead to sex reversal in mice and humans. However, little is known about the function of cbx2 in gonadal development in fish. In this study, the cbx2 gene was identified in medaka, which is a model species for the study of gonadal development in fish. Transcription of cbx2 was abundant in the gonads, with testicular levels relatively higher than ovarian levels. In situ hybridization (ISH) revealed that cbx2 mRNA was predominately localized in spermatogonia and spermatocytes, and was also observed in oocytes at stages I, II, and III. Furthermore, cbx2 and vasa (a marker gene) were co-localized in germ cells by fluorescent in situ hybridization (FISH). After cbx2 knockdown in the gonads by RNA interference (RNAi), the sex-related genes, including sox9 and foxl2, were influenced. These results suggest that cbx2 not only plays a positive role in spermatogenesis and oogenesis but is also involved in gonadal differentiation through regulating the expression levels of sex-related genes in fish.

Nucleoporin 62-Like Protein is Required for the Development of Pharyngeal Arches through Regulation of Wnt/β-Catenin Signaling and Apoptotic Homeostasis in Zebrafish

We have previously observed the predominant expression of nucleoporin 62-like (Nup62l) mRNA in the pharyngeal region of zebrafish, which raises the question whether Nup62l has important implications in governing the morphogenesis of pharyngeal arches (PA) in zebrafish. Herein, we explored the functions of Nup62l in PA development. The disruption of Nup62l with a CRISPR/Cas9-dependent gene knockout approach led to defective PA, which was characterized by a thinned and shortened pharyngeal region and a significant loss of pharyngeal cartilages. During pharyngeal cartilage formation, prechondrogenic condensation and chondrogenic differentiation were disrupted in homozygous nup62l-mutants, while the specification and migration of cranial neural crest cells (CNCCs) were unaffected. Mechanistically, the impaired PA region of nup62l-mutants underwent extensive apoptosis, which was mainly dependent on activation of p53-dependent apoptotic pathway. Moreover, aberrant activation of a series of apoptotic pathways in nup62l-mutants is closely associated with the inactivation of Wnt/β-catenin signaling. Thus, these findings suggest that the regulation of Wnt/β-catenin activity by Nup62l is crucial for PA formation in zebrafish.

A case report of acampomelic campomelic dysplasia and operative difficulties in cleft palate reconstruction

Acampomelic campomelic dysplasia (CD) is a type of CD (CD; OMIM #114290), a rare form of congenital short-limbed dwarfism and is due to mutations in SOX9 gene family. Characteristic phenotypes of CD include bowing of the lower limbs, a narrow thoracic cage, 11 pairs of ribs, hypoplastic scapulae, macrocephaly, flattened supraorbital ridges and nasal bridge, cleft palate and micrognathia. The bending of the long bones is not an obligatory feature and is absent in about 10% of cases, referred to as acampomelic CD. A child previously diagnosed with acampomelic CD was brought to our outpatient clinic for cleft palate reconstruction. Our neurosurgeon cautioned us against performing surgery with extension of the neck in view of the possibility of producing quadriparesis, due to narrowing of the spinal canal as part of the osseous anomaly noted in the magnetic resonance imaging study of the spine, thus making the anaesthesia, surgical and post-operative procedures difficult. The cleft palate reconstruction was performed with all precautions and was uneventful.