Comparative phylotranscriptomics reveals putative sex differentiating genes across eight diverse bivalve species

Mollusks, especially bivalves, exhibit a great diversity of sex determining mechanisms, including both genetic and environmental sex determination. Some bivalve species can be gonochoristic (separate sexes), while others are hermaphroditic (sequential or simultaneous). Several models have been proposed for specific bivalve species, utilizing information gained from gene expression data, as well as limited RAD-seq data (e.g., from Crassostrea gigas). However, these mechanisms are not as well studied as those in model organisms (e.g., Mus musculus, Drosophila melanogaster, Caenorhabditis elegans) and many genes involved in sex differentiation are not well characterized. We used phylotranscriptomics to better understand which possible sex differentiating genes are in bivalves and how these genes relate to similar genes in diverse phyla. We collected RNAseq data from eight phylogenetically diverse bivalve species: Argopecten irradians, Ensis directus, Geukensia demissa, Macoma tenta, Mercenaria mercenaria, Mya arenaria, Mytilus edulis, and Solemya velum. Using these data, we assembled representative transcriptomes for each species. We then searched for candidate sex differentiating genes using BLAST and confirmed the identity of nine genes using phylogenetics analyses from nine phyla. To increase the confidence of identification, we included ten bivalve genomes in our analyses. From the analysis of doublesex and mab-3 related transcription factor (DMRT) genes, we confirmed the identify of a Mollusk-specific sex determining DMRT gene: DMRT1L. Based on gene expression data from M. edulis and previous research, DMRT1L and FoxL2 are key genes for male and female development, respectively.

H3K27ac chromatin acetylation and gene expression analysis reveal sex- and situs-related differences in developing chicken gonads

BACKGROUND

Birds exhibit a unique asymmetry in terms of gonad development. The female left gonad generates a functional ovary, whereas the right gonad regresses. In males, both left and right gonads would develop into testes. How is this left/right asymmetry established only in females but not in males remains unknown. The epigenetic regulation of gonadal developmental genes may contribute to this sex disparity. The modification of histone tails such as H3K27ac is tightly coupled to chromatin activation and gene expression. To explore whether H3K27ac marked chromatin activation is involved in the asymmetric development of avian gonads, we probed genome-wide H3K27ac occupancy in left and right gonads from both sexes and related chromatin activity profile to the expression of gonadal genes. Furthermore, we validated the effect of chromatin activity on asymmetric gonadal development by manipulating the chromatin histone acetylation levels.

METHODS

The undifferentiated gonads from both sides of each sex were collected and subjected to RNA-Seq and H3K27ac ChIP-Seq experiments. Integrated analysis of gene expression and active chromatin regions were performed to identify the sex- and situs-specific regulation and expression of gonadal genes. The histone deacetylase inhibitor trichostatin A (TSA) was applied to the undifferentiated female right gonads to assess the effect of chromatin activation on gonadal gene expression and cell proliferation.

RESULTS

Even before sex differentiation, the gonads already show divergent gene expression between different sexes and between left/right sides in females. The sex-specific H3K27ac chromatin distributions coincide with the higher expression of male/female specification genes in each sex. Unexpectedly, the H3K27ac marked chromatin activation show a dramatic difference between left and right gonads in both sexes, although the left/right asymmetric gonadal development was observed only in females but not in males. In females, the side-specific H3K27ac occupancy instructs the differential expression of developmental genes between the pair of gonads and contributes to the development of left but not right gonad. However, in males, the left/right discrepancy of H3K27ac chromatin distribution does not drive the side-biased gene expression or gonad development. The TSA-induced retention of chromatin acetylation causes up-regulation of ovarian developmental genes and increases cell proliferation in the female right gonad.

CONCLUSIONS

We revealed that left/right asymmetry in H3K27ac marked chromatin activation exists in both sexes, but this discrepancy gives rise to asymmetric gonadal development only in females. Other mechanisms overriding the chromatin activation would control the symmetric development of male gonads in chicken.

The effect of methyltestosterone (MT) on sex differentiation and growth in juvenile yellow perch (Perca flavescens)

A study was conducted to evaluate the gonad differentiation of juvenile yellow perch (YP, Perca flavencens) and determine the latest labile period related to hormone treatment. Juvenile fish were subjected to two dietary concentrations of methyltestosterone (MT; 20 and 50 mg/kg feed) for 60 days in three (3) age groups of 38-, 46-, and 67-days post-hatching (dph), where control group were fed with standard commercial feed. Following a 10-month on-growing period, sex phenotypes were determined by gross and histological gonad morphology. Results showed the juvenile YP responded to the exogenous hormone when it was applied at 38 dph for both 20 and 50 mg/kg feed resulting in 100% males. At 46 dph, only 50 mg/kg feed resulted in 100% males. Both MT-treated at 38 and 46 dph significantly differed (P < 0.01) from the expected normal population of male:female (1:1). MT-treated at 67 dph resulted in 37% and 25% intersex fish for both 20 and 50 mg/kg feed dosage groups, respectively. MT-treated at 38 and 46 dph promoted growth and showed significantly heavier mean body weight (P < 0.05) compared to control. The gonadosomatic index (GSI) of MT-treated at 38 and 46 dph was significantly lower than that in control. This study provides the first evidence that juvenile YP can be successfully masculinized when the treatment is initiated at the age of up to 46 dph. The result is important for sex control in aquaculture.

Systematic genome-wide analysis of the ethylene-responsive ACS gene family: Contributions to sex form differentiation and development in melon and watermelon

Ethylene is an important regulatory phytohormone for sex differentiation and flower development. As the rate-limiting enzyme encoding genes in ethylene biosynthesis, ACS gene family has been well studied in cucumber; however, little is known in other cucurbit crops, such as melon and watermelon, which show diverse sex types in the field. Here, we identified and characterized eight ACS genes each in the genomes of melon and watermelon. According to the conserved serine residues at C-terminal, all the ACS genes could be characterized into three groups, which were supported by the exon-intron organizations and conserved motif distributions. ACS genes displayed diverse tissue-specific expression patterns among four melon and three watermelon sex types. Furthermore, a comparative expression analysis in the shoot apex identified orthologous pairs with potential functions in sex determination, e.g., ACS1s and ACS6s. All ACS orthologs in melon and watermelon exhibited similar expression patterns in monoecious and gynoecious genotypes, except for ACS11s and ACS12s. As expected, the majority of ACS genes were responsive to exogenous ethephon; however, some orthologs exhibited opposite expression patterns, such as ACS1s, ACS9s, and ACS10s. Collectively, our findings provide valuable ACS candidates related to flower development in various sex types of melon and watermelon.

Histological evaluation of sex differentiation and early sex identification in hatchery-produced greater amberjack (Seriola dumerili) reared in sea cages

The histological process of gonadal differentiation, together with the endocrine changes of sex steroid hormones and some of their precursors, was studied in hatchery-produced greater amberjack Seriola dumerili from 101 until 408 days post-hatching (dph), with samplings conducted every 50 days. Histological processing showed that sex differentiation began at 101 dph with the formation of the ovarian cavity in females, while the presumptive males did not yet contain any germ cells in their gonad. At 150 dph, we observed the first germ cells in the developing testes. Sex differentiation in almost all sampled individuals was complete at 408 dph. No size dimorphism was observed between the sexes, and the sex ratio was 1:1, suggesting that there was no influence of early rearing in captivity on sex differentiation. Plasma concentrations of adrenosterone (Ad), androstenedione (Δ4), 11-ketotestosterone (11ΚΤ), testosterone (Τ), estradiol (Ε₂), progesterone (P4) and 17,20β-dihydroxy-4-pregnen-3-one (17,20βP) were measured in males and females with the use of liquid chromatography tandem mass spectrometry (LC-MS/MS) to examine their role in the sex differentiation process. From the seven hormones, the only one that exhibited differences between the sexes was 11-KT and the plasma 11-KT concentration was found to be a useful indication of greater amberjack sex. Variations were observed in the mean values of Ad, Δ4, 11-KT, T, P4 and 17,20βP over time in one or both sexes, indicating their involvement in the sex differentiation process.

Profile of gene expression changes during estrodiol-17β-induced feminization in the Takifugu rubripes brain

Background

As the critical tissue of the central nervous system, the brain has been found to be involved in gonad development. Previous studies have suggested that gonadal fate may be affected by the brain. Identifying brain-specific molecular changes that occur during estrodiol-17β (E₂) -induced feminization is crucial to our understanding of the molecular control of sex differentiation by the brains of fish.

Results

In this study, the differential transcriptomic responses of the Takifugu rubripes larvae brain were compared after E₂ treatment for 55 days. Our results showed that 514 genes were differentially expressed between E₂-treated-XX (E-XX) and Control-XX (C-XX) T. rubripes, while 362 genes were differentially expressed between E₂-treated-XY (E-XY) and Control-XY (C-XY). For example, the expression of cyp19a1b, gnrh1 and pgr was significantly up-regulated, while st, sl, tshβ, prl and pit-1, which belong to the growth hormone/prolactin family, were significantly down-regulated after E₂ treatment, in both sexes. The arntl1, bhlbe, nr1d2, per1b, per3, cry1, cipc and ciart genes, which are involved in the circadian rhythm, were also found to be altered. Differentially expressed genes (DEGs), which were identified between E-XX and C-XX, were significantly enriched in neuroactive ligand-receptor interaction, arachidonic acid metabolism, cytokine-cytokine receptor interaction and the calcium signaling pathway. The DEGs that were identified between E-XY and C-XY were significantly enriched in tyrosine metabolism, phenylalanine metabolism, arachidonic acid metabolism and linoleic acid metabolism.

Conclusion

A number of genes and pathways were identified in the brain of E₂-treated T. rubripes larvae by RNA-seq. It provided the opportunity for further study on the possible involvement of networks in the brain-pituitary-gonadal axis in sex differentiation in T. rubripes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08158-0.

Transcriptomic analysis revealed gene expression profiles during the sex differentiation of Chinese tongue sole (Cynoglossus semilaevis)

Sex differentiation in aquatic fish is important both for theoretical study and practical production, as growth dimorphism frequently appears in different sexes, especially in marine fish. The deciphered genome, identification of the male-determining gene dmrt1 and established genotypic sex screening method make Chinese tongue sole (Cynoglossus semilaevis) an ideal model to study sex differentiation in fish. In this study, comparative gonadal transcriptomic analyses were conducted for genetic females and males at 48, 68, and 108 days post hatching (dph), representing pre-, during- and post-gonadal differentiation stages, although the gonad is not completely differentiated and isolable in 48 and 68 dph individuals, while it is in 108 dph individuals. Altogether, 28 libraries were constructed, and a mean of 46.64 M clean reads was obtained. Differentially expressed gene (DEG) analysis revealed that 179 genes had similar expression patterns in males and females in all three stages. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis suggested that the enriched pathways included ubiquitin-mediated proteolysis, lysosomes, and RNA degradation. Moreover, weighted gene coexpression network analyses (WGCNA) identified 14 modules, one of which was closely correlated with female differentiation, exhibiting female-biased expression in all three stages (48, 68, 108 dph). An illustrated core gene interaction network of this module identified 50 genes, most of which are on W chromosomes. Six genes, including two ubiquitin conjugating enzymes, were selected for further investigation, and their female-biased expression was confirmed in even earlier stages, at 10 and 30 dph. These data facilitate our understanding of sex differentiation in fish and provide a genomic rationale for screening candidate genes (preferentially W-linked genes) that could be involved in the female differentiation process.

Non-Coding RNAs: lncRNAs, miRNAs, and piRNAs in Sexual Development

Non-coding RNAs (ncRNAs) are a group of RNAs that do not encode functional proteins, including long non-coding RNAs (lncRNAs), microRNAs (miRNAs), PIWI-interacting RNAs (piRNAs), and short interfering RNAs (siRNAs). In the last 2 decades an effort has been made to uncover the role of ncRNAs during development and disease, and nowadays it is clear that these molecules have a regulatory function in many of the developmental and physiological processes where they have been studied. In this review, we provide an overview of the role of ncRNAs during gonad determination and development, focusing mainly on mammals, although we also provide information from other species, in particular when there is not much information on the function of particular types of ncRNAs during mammalian sexual development.

Characterization of the sex differentiation and gonadal development in small yellow croaker (Larimichthys polyactis) and its hybrid (L. polyactis ♀ × L. crocea ♂)

Interspecific hybridization has been considered as a possible approach to improve biological traits and has been applied in aquaculture practices. In the present study, artificial hybridization was carried out in the small yellow croaker (SYC; Larimichthys polyactis) ♀ × large yellow croaker (LYC; L. crocea) ♂ by artificial insemination, and the processes of sex differentiation and gonadal development in SYC and its hybrid were investigated under controlled conditions. Histological analysis of SYC larvae showed that migrating primordial germ cells (PGCs) were observed at 4 days post-hatching (dph), a genital ridge was formed on the dorsal side of the peritoneum at 6 dph, and a pair of primary gonads was first observed at 10 dph. Signs of the differentiated ovary and ovarian cavity were observed at 45 dph. However, some presumptive testes showed alterations in morphology, including an increase in the number of oocytes and an enhanced basophilia at 50 dph. These presumptive testes seemed to alter again, and numerous gonial cells were arranged in cyst-like groups with several degenerating oocytes that developed into residual body-like structures during 60-90 dph. Compared with SYC, the hybrid had a lower number of PGCs and showed retarded gonadal development at the early stage. Ovarian differentiation in the hybrid was observed at 50 dph, while testicular differentiation occurred at 60 dph. The presence of vitellogenic oocytes and spermatozoa at 360 dph in the hybrids suggested that hybrid individuals can undergo successful gametogenesis in females and males, respectively. Overall, the present results suggest that morphological sex differentiation occurred at 40 and 50 dph in SYC and its hybrid, respectively, both of which have normal gametogenesis. Moreover, some level of heterosis (hybrid vigor) occurred in the growth of the hybrid (total length and body weight) compared with that in the growth of SCY over time. Gonadal development of the hybrid was also found to be advanced at 360 dph. The present information will contribute to the potential use and management of these fish for aquaculture.

Growth and gonadal development retardations after long-term exposure to estradiol in little yellow croaker, Larimichthys polyactis

Endocrine disrupting chemicals (EDCs) including 17β-estradiol (E2) are widely distributed in the aquatic environment and are known to negatively affect the reproductive system of many animals, including fish. EDCs leading to feminization, altered sex ratio and reduced fecundity, it is possibly posing potential risks to the ecosystems. To investigate the potentially toxic effects of E2 exposure on little yellow croaker (Larimichthys polyactis, L. poliactis) who have a unique gonadal development pattern that males undergo a hermaphroditic stage. An experiment was set up where L. poliactis were maintained in tanks and exposed to E2 concentrations of 10 μg/L or no E2 exposure (the ethanol and control groups) from 30 to 90 days post-hatching (dph). After exposure, the E2 withdrawal and continual cultured to 150 and 365 dph. The morphological and histological analyses were used to compare the changes in the fish body and gonad under E2 exposure. The results showed that E2 exposure caused three major phenotypes at 30 and 60 days after treatment (dat), including ovary, ovotestis and gonadal development retardation compared with the control groups. The average ratio of these three phenotypes is 60.6%, 11.97% and 27.43%, respectively. The body length and weight of E2 exposure groups were repressed during the E2 exposure period, while it can recover after E2 withdrawal. However, the gonadal development (Gonadosomatic Index) of E2 exposure groups testis were retarded at 60 dat and doesn't recover until 365 dph. The sex determination/differentiation-related genes erα, erβI, erβII, fshβ and cyp11b2 were significantly decreased in E2-exposure male fish. This research highlights the E2 leads to feminization, disrupts testis maturation and spermatogenesis, this effect persisted into the stage of sexual maturity. Collectively, our findings provide insights into the molecular mechanisms underlying E2 disturbance of a marine economic fish reproduction.

Developmental aspects of the hypothalamic-pituitary network related to reproduction in teleost fish

The hypothalamic-pituitary-gonadal axis is the main system that regulates reproduction in vertebrates through a complex network that involves different neuropeptides, neurotransmitters, and pituitary hormones. Considering that this axis is established early on life, the main goal of the present work is to gather information on its development and the actions of its components during early life stages. This review focuses on fish because their neuroanatomical characteristics make them excellent models to study neuroendocrine systems. The following points are discussed: i) developmental functions of the neuroendocrine components of this network, and ii) developmental disruptions that may impact adult reproduction. The importance of the components of this network and their susceptibility to external/internal signals that can alter their specific early functions and/or even the establishment of the reproductive axis, indicate that more studies are necessary to understand this complex and dynamic network.

Turning the spotlight on the C11-oxy androgens in human fetal development

Research into the biosynthesis of C11-oxy C₁₉ steroids during human fetal development, specifically fetal adrenal development and during the critical period of sex differentiation, is currently lacking. Cortisol, which possesses a C11-hydroxyl moiety has, however, been firmly established in this context. Compelling questions are whether the C11-oxy C₁₉ steroids (11β-hydroxyandrostenedione, 11β-hydroxytestosterone, 11-ketoandrostenedione and 11-ketotestosterone [11KT]) and the C11-oxy C₂₁ steroids (11β-hydroxyprogesterone and 11-ketoprogesterone) are biosynthesised during gestation, and whether these hormones circulate between the placenta and the developing fetus, and between the placenta and the mother. This review will consider the role of cortisol, 11KT and 11β-hydroxysteroid dehydrogenase type 2 (11βHSD2) in determining the sex of teleost fish, while these hormones and 11βHSD2 will also be discussed with regards to murine mammals. The focus of the review will shift to highlight the potential role of C11-oxy steroids in human fetal development based on the timely expression of steroidogenic enzymes in the adrenal, testes and ovary, as well as in the placenta; summarising reported evidence of C11-oxy steroids in neonatal life.

A sex-linked SNP mutation in amhr2 is responsible for male differentiation in obscure puffer (Takifugu obscurus)

BACKGROUND

Anti-Mullerian hormone receptor type II (Amhr2) is a key receptor of Amh signaling in regulating gonad development. The amhr2 gene has been identified in numerous species, including a few teleost fishes. However, the roles of Amhr2 in Amh signaling in fish are poorly studied.

METHODS AND RESULTS

In this study, an amhr2 homolog from obscure puffer (Takifugu obscurus) was identified, and its molecular characteristics were systematically analyzed. Expression analysis revealed that amhr2 was highly expressed in the gonads of adult pufferfish and significantly upregulated during sex differentiation. Significantly, a sex-linked SNP site was verified in obscure puffer amhr2. Females exhibited a homozygous genotype (C/C), while males possessed a heterozygous genotype (C/G), resulting in an amino acid variation (His/Asp384) in the kinase domain of Amhr2. Then, the functions of the different Amhr2 genotypes were further investigated. The male genotype protein (Amhr2^D384) showed an enhanced ability to interact with the type I receptor (Bmpr1a) compared to the female genotype (Amhr2^H384). The phosphorylation levels of Smads and activity of the target gene (id3) induced by the male genotype were also much higher than those induced by the female genotype. These results confirmed that the male genotype had an enhanced effect on the Amh signaling pathway compared with the female genotype.

CONCLUSIONS

This study provides direct experimental evidence for the roles of different Amhr2 genotypes in pufferfish and suggests that amhr2 is responsible for male sex differentiation in obscure puffer.

Gonadal sex differentiation and early ovarian/testicular development in cultured Pacific bluefin tuna, Thunnus orientalis (Temminck et Schlegel)

Pacific bluefin tuna (PBT), Thunnus orientalis, is one of the most important species for aquaculture in Japan. Recently, the reduction in muscle fat content associated with sexual maturation in farmed PBT has become a serious problem. To develop technologies for inducing sterility, detailed and reliable data on gonadal development in PBT are needed. Here, we demonstrated the process of gonadal sex differentiation, and of early ovarian and testicular development during the immature stages in PBT. Gonadal sex differentiation was first characterized by the formation of the ovarian cavity in female and of the efferent ducts in male 57 days post hatching (dph). The gonads then differentiated into ovaries or testes according to the genotypic sex until 83 dph. During this period, primordial germ cells, oogonia, and type-A spermatogonia were solitarily distributed in the gonads, and the number of germ cells did not differ between sexes. After gonadal sex differentiation, gonads of PBTs developed in a sexually dimorphic manner: proliferation and differentiation of germ cells occurred earlier in the ovaries than in the testes. The oogonia in ovaries formed cysts at 185 dph, but the type-A spermatogonia were solitarily distributed in testes at this stage, and cysts of type-A spermatogonia were first observed at 247 dph. Moreover, the oogonia entered meiosis and differentiated into chromatin-nucleolus stage oocytes until 247 dph, and subsequently into peri-nucleolus stage oocytes until 285 dph, whereas the type-A spermatogonia differentiated into type-B spermatogonia, spermatocytes, spermatids, and spermatozoa from 446 dph onwards. We believe the results of this study provide the necessary basis for future studies on sterile PBT production.

Metabolism and Sex Differentiation in Animals from a Starvation Perspective

Animals determine their sex genetically (GSD: genetic sex determination) and/or environmentally (ESD: environmental sex determination). Medaka (Oryzias latipes) employ a XX/XY GSD system, however, they display female-to-male sex reversal in response to various environmental changes such as temperature, hypoxia, and green light. Interestingly, we found that 5 days of starvation during sex differentiation caused female-to-male sex reversal. In this situation, the metabolism of pantothenate and fatty acid synthesis plays an important role in sex reversal. Metabolism is associated with other biological factors such as germ cells, HPG axis, lipids, and epigenetics, and supplys substances and acts as signal transducers. In this review, we discuss the importance of metabolism during sex differentiation and how metabolism contributes to sex differentiation.

Differential gene expression patterns during gametophyte development provide insights into sex differentiation in the dioicous kelp Saccharina japonica

BACKGROUND

In brown algae, dioicy is the prevalent sexual system, and phenotypic differences between male and female gametophytes have been found in many dioicous species. Saccharina japonica show remarkable sexual dimorphism in gametophytes before gametogenesis. A higher level of phenotypic differentiation was also found in female and male gametes after gametogenesis. However, the patterns of differential gene expression throughout gametophyte development and how these changes might relate to sex-specific fitness at the gamete stage in S. japonica are not well known.

RESULTS

In this study, differences in gene expression between male and female gametophytes in different developmental stages were investigated using comparative transcriptome analysis. Among the 20,151 genes expressed in the haploid gametophyte generation, 37.53% were sex-biased. The abundance of sex-biased genes in mature gametophytes was much higher than that in immature gametophytes, and more male-biased than female-biased genes were observed in the mature stage. The predicted functions of most sex-biased genes were closely related to the sex-specific characteristics of gametes, including cell wall biosynthesis, sperm motility, and sperm and egg recognition. In addition, 51 genes were specifically expressed in males in both stages, showing great potential as candidate male sex-determining region (SDR) genes.

CONCLUSIONS

This study describes a thorough investigation into differential gene expression between male and female gametophytes in the dioicous kelp S. japonica. A large number of sex-biased genes in mature gametophytes may be associated with the divergence of phenotypic traits and physiological functions between female gametes (eggs) and male gametes (sperm) during sexual differentiation. These genes may mainly come from new sex-biased genes that have recently evolved in the S. japonica lineage. The duplication of sex-biased genes was detected, which may increase the number of sex-biased genes after gametogenesis in S. japonica to some extent. The excess of male-biased genes over female-biased genes in the mature stage may reflect the different levels of sexual selection across sexes. This study deepens our understanding of the regulation of sex development and differentiation in the dioicous kelp S. japonica.

Effects of synthetic sex steroid hormone exposures on gonadal sex differentiation and dynamics of a male-related gene, Gonadal soma-derived factor (Gsdf) and an estrogen up-regulated gene, Choriogenine-H (ChgH) gene expression in the euryhaline Javafish medaka, Oryzias javanicus, based on genetic sexes

To clarify the basal aspects of sexual development in Javafish medaka, Oryzias javanicus (ZZ/ZW), a model marine species for ecotoxicity testing, we examined the details of gonadal sex differentiation and exogenous sex hormone-dependent sex reversals using genetic sex-linked DNA markers. Sex differences in germ cell numbers were observed at 5 days post hatching (dph), in which there was a significant increase in the germ cells of ZW. In ZW, diplotene oocytes and the ovarian cavity appeared at approximately 10, and 30 dph, respectively. In ZZ, spermatogonial proliferation was observed at approximately 20 dph. A ZZ-dominant expression of Gonadal soma-derived factor (Gsdf) mRNA was detected before hatching. The exposure of embryos to 17α-ethinylestradiol (EE2; 0.1, 1, 10 ng/mL) did not cause sex reversals in most cases. However, EE2 exposures led to significant Choriogenin-H (ChgH) mRNA expression, an estrogen up-regulated gene, in all fry; these exposures did not suppress Gsdf expression in ZZ fry. The exposure of embryos to 17α-methyltestosterone (MT; 0.1, 1, 10 ng/mL) caused sex reversals but only at low frequencies in ZW and ZZ fish. Although the 10 ng/mL MT exposure was accompanied by induction of significant Gsdf expression in ZW fry, induction of ChgH expression was also observed in several fry. Together, the present study indicates for the first time that male-dominant sexual dimorphic expression of Gsdf precedes the first morphological sex difference, i.e., the sex difference in germ cell number, and results strongly suggest that exogenous sex hormone-dependent sex reversal is not induced easily in O. javanicus.

Where the Ends Meet: An Overview of Sex Determination in Atheriniform Fishes

Atheriniform fishes have recently emerged as attractive models for evolutionary, ecological, and molecular/physiological studies on sex determination. Many species in this group have marked temperature-dependent sex determination (TSD) and yet many species also have a sex determinant gene that provides a strong drive for male differentiation. Thus, in these species the 2 forms of sex determination that were once considered to be mutually exclusive, environmental (ESD) and genotypic (GSD) sex determination, can coexist at environmentally relevant conditions. Here, we review the current knowledge on sex determination in atheriniform fishes with emphasis on the molecular and physiological mechanisms of ESD and GSD, the coexistence and cross-talk between these 2 mechanisms, the possibility of extragonadal transduction of environmental information and/or extragonadal onset of sex determination, and the results of field studies applying novel tools such as otolith increment analysis and molecular markers of genetic sex developed for selected New World and Old World atheriniform species. We also discuss the existence of molecular and histological mechanisms to prevent the discrepant differentiation in parts of the gonads because of ambiguous or conflicting environmental and genetic signals and particularly the possibility that the female is the default state in these species.

Comparative proteomics analysis reveals differentially accumulated proteins associated with male and female A. chinensis var. chinensis bud development

Background

Kiwifruit (Actinidia chinensis var. Chinensis) is abundant with vitamin C and is a rapidly developing crop in China, New Zealand, and other countries. It has been widely used as a raw material for food and kiwifruit wine. Among these, A. chinensis var. chinensis and A. chinensis var. deliciosa are the most valuable kiwifruit in production. Kiwifruit is a typical dioecious plant and its female and male plants have different economic values. Therefore, sex identification, especially at the seedling stage, has important implications for the scientific planning of its production and economic benefits. However, the kiwifruit sex regulation mechanism is very complex and molecular studies are in the initial stages. Currently, there is not a universal and effective sex identification method for A. chinensis.

Methods

In this study, we used a label-free quantitative proteomics approach to investigate differentially accumulated proteins, including their presence/absence and significantly different levels of abundances during A. chinensis var. chinensis male and female flower bud development.

Results

A total of 6485 proteins were identified, among which, 203 were identified in male buds, which were mainly associated with phenylalanine metabolism, tyrosine metabolism, and plant hormone signal transduction. In female buds, 241 were identified, which were mainly associated with the ErbB signaling pathway, growth hormone synthesis, secretion and action, and mRNA surveillance pathway. A total of 373 proteins were significantly differentially accumulated proteins (fold change > 2; P < 0.05), of which, 168 were upregulated and 205 were downregulated. Significant differences between proteins involved 13 signaling pathways, most of which were involved in flavonoid biosynthesis, phenylpropanoid biosynthesis, and starch and sucrose metabolism. Protein interaction analysis showed that enriched protein nodes included cell division cycle 5-like protein, 40S ribosomal protein S8, ribosomal protein, and 40S ribosomal protein like, which interact with 35, 25, 22, and 22 proteins, respectively.

Conclusions

This study provide valuable information for cloning key genes that control sex traits and functionally analyze their roles, which lay a foundation to the development of molecular markers for male and female kiwifruit identification.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12953-021-00176-w.

Phytohormone and integrated mRNA and miRNA transcriptome analyses and differentiation of male between hermaphroditic floral buds of andromonoecious Diospyros kaki Thunb

BACKGROUND

Persimmon (Diospyros kaki Thunb.) has various labile sex types, and studying its sex differentiation can improve breeding efficiency. However, studies on sexual regulation patterns in persimmon have focused mainly on monoecy and dioecy, whereas little research has been published on andromonoecy. In order to reveal the sex differentiation regulation mechanism of andromonoecious persimmon, we performed histological and cytological observations, evaluated OGI and MeGI expression and conducted phytohormones assays and mRNA and small RNA transcriptome analyses of the male and hermaphroditic floral buds of the andromonoecious persimmon 'Longyanyeshi 1'.

RESULTS

Stages 2 and 4 were identified as the critical morphological periods for sex differentiation of 'Longyanyeshi 1' by histological and cytological observation. At both stages, OGI was differentially expressed in male and hermaphroditic buds, but MeGI was not. This was different from their expressions in dioecious and monoecious persimmons. Meantime, the results of phytohormones assays showed that high IAA, ABA, GA3, and JA levels at stage 2 may have promoted male floral bud differentiation. However, high JA levels at stage 4 and high ZT levels at stages 2 and 4 may have promoted hermaphroditic floral bud differentiation. In these phytohormone biosynthesis and signaling pathways, 52 and 54 differential expression genes (including Aux/IAA, ARFs, DELLA, AHP, A-ARR, B-ARR, CYP735A, CRE1, PP2C, JAZ, MYC2, COI1, CTR1, SIMKK, ACO, and MPK6) were identified, respectively. During the development of male floral buds, five metacaspases genes may have been involved in pistil abortion. In addition, MYB, FAR1, bHLH, WRKY, and MADS transcription factors might play important roles in persimmon floral bud sex differentiation. Noteworthy, miR169v_1, miR169e_3, miR319_1, and miR319 were predicted to contribute to phytohormone biosynthesis and signaling pathways and floral organogenesis and may also regulate floral bud sex differentiation.

CONCLUSION

The present study revealed the differences in morphology and phytohormones content between male and hermaphroditic floral buds of 'Longyanyeshi 1' during the process of sex differentiation, and identified a subset of candidate genes and miRNAs putatively associated with its sex differentiation. These findings can provide a foundation for molecular regulatory mechanism researching on andromonoecious persimmon.

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.

Spin1z induces the male pathway in the chicken by down-regulating Tcf4

SPINDLIN1-Z (SPIN1Z), a member of the Spin/Ssty(Y-linked spermiogenesis specific transcript) protein family, participates in the early embryonic development process. Our previous RNA-seq analysis indicates that the level of Spin1z was abundantly expressed in male embryonic stem cells (ESCs) and primitive germ cells (PGCs), we speculate that Spin1z may play an important role in chicken male differentiation. Therefore, the loss- and gain-of-function experiments provide solid evidence that Spin1z is both necessary and sufficient to initiate male development in chicken. Furthermore, chromatin immunoprecipitation (ChIP) assay and the dual-luciferase assay was performed to further confirm that Spin1z contributed to chicken male differentiation by inhibiting the Tcf4 transcription. Our findings provide a novel insight into the molecular mechanism for chicken male differentiation.

Transcriptome Profiling Reveals the Sexual Dimorphism of Gene Expression Patterns during Gonad Differentiation in the Half-Smooth Tongue Sole (Cynoglossus semilaevis)

The half-smooth tongue sole (Cynoglossus semilaevis), one of the most economically-important fish species in China, exhibits sexually dimorphic growth. An understanding of sex-related gene expression patterns in the tongue sole may inform sex regulation and breeding processes that increase fish production. However, the gene expression patterns during gonad development in the tongue sole remain unknown. In this study, transcriptome sequencing analyses were performed during gonad differentiation in the tongue sole, namely, at 62 days post-hatching (dph), 100 dph, 120 dph, and 150 dph. Differentially expressed genes associated with sex differentiation and gonad development were identified at each time point. Trend analysis showed that gene expression patterns varied over time. These expression patterns either explained common, non-sexually-dimorphic features or indicated significant sexual dimorphism. Transcript structure analyses identified both sex and time differences among samples. This study investigated the time-dependent expression patterns of several sex-related genes, including Dmrt1, Amh, Foxl2, aromatase encoding gene, Esr, and the Sox gene family, during gonad differentiation in the tongue sole. These results might clarify the significant sexual differences during early development in the tongue sole and might provide insight into the mechanisms controlling sex differentiation and development.

R-spondin1 in loach (Misgurnus anguillicaudatus): Identification, characterization, and analysis of its expression patterns and DNA methylation in response to high-temperature stress

With a well-understood function in mammals, R-spondin1 (Rspo1) is an important regulator of ovarian development via the Wnt/β-catenin pathway. Rspo1 deficiency causes retardation of ovarian development in XX fish, and increases Rspo1 function induces femininity and sex reversal in XY fish. In this study, Rspo1 was successfully cloned from loach (Misgurnus anguillicaudatus), and its expression profile was analyzed. The full-length cDNA of Misgurnus anguillicaudatus Rspo1 (MaRspo1) comprised 1322 bp and included an open reading frame (ORF) of 795 bp, which encoded a predicted polypeptide measuring 264 amino acids in length. Phylogenetic and gene structure analyses showed a highly conserved sequence of MaRspo1 (identical to the Rspo1 genes of other species), consisting of an N-terminal signal peptide (SP), two furin-like cysteine-rich domains (FU1 and FU2), a thrombospondin type 1 repeat (TSP1) and a C-terminal region. Real-time PCR revealed the female-biased expression profile of MaRspo1, with the highest expression level among tested tissues detected in ovary. Investigation of MaRspo1 expression levels throughout the early development stage (10-60 days post hatching) under three temperature treatments (25 °C, 28 °C, and 31 °C) revealed significantly differential expression of MaRspo1 among the three temperature groups, with decreased MaRspo1 expression in the high-temperature (31 °C) group. The results of DNA methylation analysis indicated that exposure to high temperature during early development can increase the average promoter methylation level of MaRspo1 in both females and males. Taken together, the results of this study provide the basis for the further investigation of the molecular mechanism of Rspo1 in response to temperature.

Developmental reproductive toxicity and endocrine activity of propiconazole in the Xenopus tropicalis model

Environmental pollutants and especially endocrine disrupting chemicals (EDCs) are implicated as one of the drivers of the amphibian declines. To advance the understanding of the risks of EDCs to amphibians, methods to determine endocrine-linked adverse effects are needed. The aims were to 1) develop a partial life-cycle assay with the model frog Xenopus tropicalis to determine endocrine perturbation and adverse developmental effects, and 2) determine effects of propiconazole in this assay. Propiconazole is a pesticide with multiple endocrine modes of action in vitro. Its potential endocrine activity and adverse effects in amphibians remain to be elucidated. Tadpoles were exposed to 0, 33 and 384 μg propiconazole/L during critical developmental windows until completed metamorphosis. At metamorphosis, a sub-sample of animals was analysed for endpoints for disruption of estrogen/androgen (sex ratio, brain aromatase activity) and thyroid pathways (time to metamorphosis). The remaining individuals were kept unexposed for 2 months post-metamorphosis to analyze effects on sexual development including gonadal and Müllerian duct maturity and gametogenesis. At metamorphosis, brain aromatase activity was significantly increased in the high-dose group compared to control. In both propiconazole groups, an increased proportion of individuals reached metamorphosis faster than the mean time for controls, suggesting a stimulatory effect on the thyroid system. At 2 months post-metamorphosis, testis size, sperm and Müllerian duct maturity were reduced in the low-dose males, and the liver somatic index in males was increased in both propiconazole groups, compared with controls. In conclusion, our results show that propiconazole exposure caused endocrine perturbations and subsequent hepatic and reproductive effects evident at puberty, indicating persistent disruption of metabolism and male reproductive function. Our findings advance the development of methodology to determine endocrine and adverse effects of EDCs. Moreover, they increase the understanding of endocrine perturbations and consequent risk of adverse effects of azoles in amphibians.