Sex-determining Gene Research Articles

Advance in sex differentiation in cucumber.

Cucumber belongs to the family Cucurbitaceae (melon genus) and is an annual herbaceous vegetable crop. Cucumber is an important cash crop that is grown all over the world. From morphology to cytology, from canonical genetics to molecular biology, researchers have performed much research on sex differentiation and its regulatory mechanism in cucumber, mainly in terms of cucumber sex determination genes, environmental conditions, and the effects of plant hormones, revealing its genetic basis to improve the number of female flowers in cucumber, thus greatly improving the yield of cucumber. This paper reviews the research progress of sex differentiation in cucumber in recent years, mainly focusing on sex-determining genes, environmental conditions, and the influence of phytohormones in cucumber, and provides a theoretical basis and technical support for the realization of high and stable yield cultivation and molecular breeding of cucumber crop traits.

The sex determination gene of Populus deltoides, PdFERR, interacts with ERF96 to promote the development of female flower organs.

The FEMALE-SPECIFICALLY EXPRESSED RESPONSE REGULATOR (PdFERR) gene in Populus deltoides, a sex determination gene (an orthologous gene of ARR17 in P. tremula), was found to promote femaleness in heterologous expression lines of Arabidopsis. None of the genes in the Arabidopsis genome seem to be orthologous to PdFERR. Although originating from two evolutionarily distant plants, the dioecious poplar FERR might promote femaleness in the hermaphroditic Arabidopsis through an evolutionary consistent regulatory pathway. However, there is no molecular evidence to support this viewpoint. In this study, to identify the shared downstream orthologous gene of PdFERR, we used yeast two-hybrid (Y2H) to screen potential interactors of PdFERR in Arabidopsis. We identified the ETHYLENE RESPONSE FACTOR 96 (AtERF96) and confirmed the interaction via in vivo and in vitro assays. The ERF96 orthologous gene in P. deltoides was also experimentally confirmed to interact with PdFERR. PdFERR could then promote femaleness in poplar or Arabidopsis through interactions with ERF96, which provide a new perspective for understanding the PdFERR gene regulating sex differentiation. This article is protected by copyright. All rights reserved.

Diverse Molecular Mechanisms Underlying Microbe-Inducing Male Killing in the Moth Homona magnanima.

Male killing (MK) is a type of reproductive manipulation induced by microbes, where sons of infected mothers are killed during development. MK is a strategy that enhances the fitness of the microbes, and the underlying mechanisms and the process of their evolution have attracted substantial attention. Homona magnanima, a moth, harbors two embryonic MK bacteria, namely, Wolbachia (Alphaproteobacteria) and Spiroplasma (Mollicutes), and a larval MK virus, Osugoroshi virus (OGV; Partitiviridae). However, whether the three distantly related male killers employ similar or different mechanisms to accomplish MK remains unknown. Here, we clarified the differential effects of the three male killers on the sex-determination cascades and development of H. magnanima males. Reverse transcription-PCR demonstrated that Wolbachia and Spiroplasma, but not OGVs, disrupted the sex-determination cascade of males by inducing female-type splice variants of doublesex (dsx), a downstream regulator of the sex-determining gene cascade. We also found that MK microbes altered host transcriptomes in different manners; Wolbachia impaired the host dosage compensation system, whereas Spiroplasma and OGVs did not. Moreover, Wolbachia and Spiroplasma, but not OGVs, triggered abnormal apoptosis in male embryos. These findings suggest that distantly related microbes employ distinct machineries to kill males of the identical host species, which would be the outcome of the convergent evolution. IMPORTANCE Many microbes induce male killing (MK) in various insect species. However, it is not well understood whether microbes adopt similar or different MK mechanisms. This gap in our knowledge is partly because different insect models have been examined for each MK microbe. Here, we compared three taxonomically distinct male killers (i.e., Wolbachia, Spiroplasma, and a partiti-like virus) that infect the same host. We provided evidence that microbes can cause MK through distinct mechanisms that differ in the expression of genes involved in sex determination, dosage compensation, and apoptosis. These results imply independent evolutionary scenarios for the acquisition of their MK ability.

Sex-lethal regulates back-splicing and generation of the sex-differentially expressed circular RNAs.

Conversely to canonical splicing, back-splicing connects the upstream 3' splice site (SS) with a downstream 5'SS and generates exonic circular RNAs (circRNAs) that are widely identified and have regulatory functions in eukaryotic gene expression. However, sex-specific back-splicing in Drosophila has not been investigated and its regulation remains unclear. Here, we performed multiple RNA analyses of a variety sex-specific Drosophila samples and identified over ten thousand circular RNAs, in which hundreds are sex-differentially and -specifically back-spliced. Intriguingly, we found that expression of SXL, an RNA-binding protein encoded by Sex-lethal (Sxl), the master Drosophila sex-determination gene that is only spliced into functional proteins in females, promoted back-splicing of many female-differential circRNAs in the male S2 cells, whereas expression of a SXL mutant (SXLRRM) did not promote those events. Using a monoclonal antibody, we further obtained the transcriptome-wide RNA-binding sites of SXL through PAR-CLIP. After splicing assay of mini-genes with mutations in the SXL-binding sites, we revealed that SXL-binding on flanking exons and introns of pre-mRNAs facilitates back-splicing, whereas SXL-binding on the circRNA exons inhibits back-splicing. This study provides strong evidence that SXL has a regulatory role in back-splicing to generate sex-specific and -differential circRNAs, as well as in the initiation of sex-determination cascade through canonical forward-splicing.

The use of in silico extreme pathway (ExPa) analysis to identify conserved reproductive transcriptional-regulatory networks in humans, mice, and zebrafish

Vertebrate sex determination and differentiation are coordinated by the activations and maintenance of reproductive transcriptional-regulatory networks (TRNs). There is considerable interest in studying the conserved design principles and functions of reproductive TRNs given that their intricate regulation is susceptible to disruption by gene mutations or exposures to exogenous endocrine disrupting chemicals (or EDCs). In this manuscript, the Boolean rules describing reproductive TRNs in humans, mice, and zebrafish, were represented as a pseudo-stoichiometric matrix model. This model mathematically described the interactions of 35 transcription factors with 21 sex determination and differentiation genes across the three species. The in silico approach of Extreme Pathway (ExPa) analysis was used to predict the extent of TRN gene activations subject to the species-specific transcriptomics data, from across various developmental life-stages. A goal of this work was to identify conserved and functional reproductive TRNs across the three species. ExPa analyses predicted the sex differentiation genes, DHH, DMRT1, and AR, to be highly active in male humans, mice, and zebrafish. Whereas FOXL2 was the most active gene in female humans and mice; and CYP19A1A in female zebrafish. These results agree with the expectation that regardless of a lack of sex determination genes in zebrafish, the TRNs responsible for canalizing male vs. female sexual differentiation are conserved with mammalian taxa. ExPa analysis therefore provides a framework with which to study the TRNs that influence the development of sexual phenotypes. And the in silico predicted conservation of sex differentiation TRNs between mammals and zebrafish identifies the piscine species as an effective in vivo model to study mammalian reproductive systems under normal or perturbed pathologies.

Promoter generation for the chimeric sex-determining gene dm-W in Xenopus frogs.

Many sex-determining genes (SDGs) were generated as neofunctionalized genes through duplication and/or mutation of gonadal formation-related genes. We previously identified dm-W as an SDG in the African clawed frog Xenopus laevis and found that a partial duplication of the masculinization gene dmrt1 created the neofunctionalized dm-W after allotetraploidization by interspecific hybridization. The allotetraploid Xenopus species have two dmrt1 genes, dmrt1.L and dmrt1.S. Xenopus laevis dm-W has four exons: two dmrt1.S-derived exons (exons 2 and 3) and two other exons (noncoding exon 1 and exon 4). Our recent work revealed that exon 4 originated from a DNA transposon, hAT-10. Here, to clarify when and how the noncoding exon 1 and its coexisting promoter evolved during the establishment of dm-W after allotetraploidization, we newly determined nucleotide sequences of the dm-W promoter region from two other allotetraploid species, X. largeni and X. petersii, and performed an evolutionary analysis. We found that dm-W acquired a new exon 1 and TATA-type promoter in the common ancestor of the three allotetraploid Xenopus species, resulting in the deletion of the dmrt1.S-derived TATA-less promoter. In addition, we demonstrated that the TATA box contributes to dm-W promoter activity in cultured cells. Collectively, these findings suggest that this novel TATA-type promoter was important for the establishment of dm-W as a sex-determining gene, followed by the degeneration of the preexisting promoter.

Variation in sex determination mechanisms may constrain parthenogenesis-induction by endosymbionts in haplodiploid systems.

Endosymbionts are maternally transmitted, and therefore benefit from maximizing female offspring numbers. Parthenogenesis-induction (PI) is the most effective type of manipulation for transmission, but has solely been detected in haplodiploid species, whereas cytoplasmic incompatibility (CI) is detected frequently across the arthropod phylum, including haplodiploids. This puzzling observation led us to hypothesize that the molecular sex-determination mechanism of the haplodiploid host may be a constraining factor in the ability of endosymbionts to induce parthenogenesis. Recent insights indicate that PI-endosymbionts may be able to directly manipulate sex-determination genes to induce the necessary steps required for PI in haplodiploids. However, sex-determination cascades vary extensively, so PI-induction would require a specialized and host-dependent tool set. Contrastingly, CI-related genes target conserved cell-cycle mechanisms, are located on mobile elements, and spread easily. Finally, endosymbiont-manipulations may have a strong impact on the effectiveness of haplodiploid biocontrol agents, but can also be used to enhance their efficacy.

Expression and functional analysis of transformer-2 in Phytoseiulus persimilis and other genes potentially participating in reproductive regulation

Transformer-2 (tra-2) is an important sex-determining gene in insects. It also plays a role in the reproduction of phytoseiid mites. We performed bioinformatic analyses for the tra-2 ortholog in Phytoseiulus persimilis (termed Pptra-2), measured its expression at different stages and quantitatively identified its function in reproduction. This gene encodes 288 amino acids with a conserved RRM domain. The peak of its expression was observed in adult females, especially ca. 5 days after mating. In addition, expression is also higher in eggs than in other stages and adult males. When Pptra-2 was silenced through RNA interference with oral delivery of dsRNA, 56% of the females had their egg hatching rates decreased in the first 5 days, from ca. 100% to ca. 20%, and maintained at low levels during the rest of the oviposition period. To detect other genes functionally related to Pptra-2, transcriptome analyses were performed on day 5 after mating. We compared mRNA expressions among interfered females with significantly reduced egg hatching rate, interfered females without significant hatching rate and CK. In total 403 differential genes were identified, of which 42 functional genes involved in the regulation of female reproduction and embryonic development were screened and discussed.

Insights into Early Ontogenesis of Salmo salar: RNA Extraction, Housekeeping Gene Validation and Transcriptional Expression of Important Primordial Germ Cell and Sex-Determination Genes

The challenge in extracting high-quality RNA impedes the investigation of the transcriptome of developing salmonid embryos. Furthermore, the mRNA expression pattern of important PGC and SD genes during the initial embryonic development of Salmo salar is yet to be studied. So, in the present study, we aimed to isolate high-quality RNA from eggs and developing embryos to check vasa, dnd1, nanos3a, sdf1, gsdf, amh, cyp19a, dmrt1 and foxl2 expression by qPCR. Additionally, four HKGs (GAPDH, UB2L3, eEf1a and β-actin) were validated to select the best internal control for qPCR. High-quality RNA was extracted, which was confirmed by spectrophotometer, agarose gel electrophoresis and Agilent TapeStation analysis. UB2L3 was chosen as a reference gene because it exhibited lower intra- and inter-sample variation. vasa transcripts were expressed in all the developmental stages, while dnd1 was expressed only up to 40 d°C. Nanos3a was expressed in later stages and remained at its peak for a shorter period, while sdf1 showed an irregular pattern of mRNA expression. The mRNA expression levels of SD genes were observed to be upregulated during the later stages of development, prior to hatching. This study presents a straightforward methodology for isolating high-quality RNA from salmon eggs, and the resulting transcript profiles of significant PGC and SD genes in S. salar could aid in improving our comprehension of reproductive development in this commercially important species.

In Vitro Comparison of Sex-Specific Splicing Efficiencies of fem Pre-mRNA under Monoallelic and Heteroallelic Conditions of csd, a Master Sex-Determining Gene in the Honeybee.

The sexual fate of honeybees is determined by the complementary sex determination (CSD) model: heterozygosity at a single locus (the CSD locus) determines femaleness, while hemizygosity or homozygosity at the CSD locus determines maleness. The csd gene encodes a splicing factor that regulates sex-specific splicing of the downstream target gene feminizer (fem), which is required for femaleness. The female mode of fem splicing occurs only when csd is present in the heteroallelic condition. To gain insights into how Csd proteins are only activated under the heterozygous allelic composition, we developed an in vitro assay system to evaluate the activity of Csd proteins. Consistent with the CSD model, the co-expression of two csd alleles, both of which lack splicing activity under the single-allele condition, restored the splicing activity that governs the female mode of fem splicing. RNA immunoprecipitation quantitative PCR analyses demonstrated that the CSD protein was specifically enriched in several exonic regions in the fem pre-mRNA, and enrichment in exons 3a and 5 was significantly greater under the heterozygous allelic composition than the single-allelic condition. However, in most cases csd expression under the monoallelic condition was capable of inducing the female mode of fem splicing contrary to the conventional CSD model. In contrast, repression of the male mode of fem splicing was predominant under heteroallelic conditions. These results were reproduced by real-time PCR of endogenous fem expression in female and male pupae. These findings strongly suggest that the heteroallelic composition of csd may be more important for the repression of the male splicing mode than for the induction of the female splicing mode of the fem gene.

Recurrent neo-sex chromosome evolution in kiwifruit.

Sex chromosome evolution is thought to be tightly associated with the acquisition and maintenance of sexual dimorphisms. Plant sex chromosomes have evolved independently in many lineages1,2 and can provide a powerful comparative framework to study this. We assembled and annotated genome sequences of three kiwifruit species (genus Actinidia) and uncovered recurrent sex chromosome turnovers in multiple lineages. Specifically, we observed structural evolution of the neo-Y chromosomes, which was driven via rapid bursts of transposable element insertions. Surprisingly, sexual dimorphisms were conserved in the different species studied, despite the fact that the partially sex-linked genes differ between them. Using gene editing in kiwifruit, we demonstrated that one of the two Y-chromosome-encoded sex-determining genes, Shy Girl, shows pleiotropic effects that can explain the conserved sexual dimorphisms. These plant sex chromosomes therefore maintain sexual dimorphisms through the conservation of a single gene, without a process involving interactions between separate sex-determining genes and genes for sexually dimorphic traits.

The Role of the Heat Shock Cognate Protein 70 Genes in Sex Determination and Differentiation of Chinese Tongue Sole (Cynoglossus semilaevis).

Fish sex determination can be affected by environmental temperature. This process relies on temperature-sensitive proteins such as heat shock proteins (HSPs). Our previous work found that heat shock cognate proteins (HSCs) may participate in high-temperature associated sex reversal of Chinese tongue sole (Cynoglossus semilaevis). However, the role of hsc genes in responding to high temperature and affecting sex determination/differentiation remains unclear. Here, by using C. semilaevis as model, we identified hsc70 and hsc70-like. hsc70 was abundant in the gonads with a testicular-higher expression at all gonadal development stages except for 6 months post fertilization (mpf). Intriguingly, hsc70-like showed higher expression in testes from 6 mpf on. Both long-term heat treatment during the temperature-sensitive sex-determining period and short-term heat stress at the end of this period caused different expression of hsc70/hsc70-like between sexes. The dual-luciferase assay results also suggested that these genes can respond to high temperature rapidly in vitro. Heat treatment of C. semilaevis testis cells overexpressed with hsc70/hsc70-like could affect the expression of sex-related genes sox9a and cyp19a1a. Our results indicated that hsc70 and hsc70-like were key regulators linking external high-temperature signals with sex differentiation in vivo and provide a new idea for understanding the mechanism by which high temperature affects sex determination/differentiation in teleosts.

Genomic basis of Y-linked dwarfism in cichlids pursuing alternative reproductive tactics.

Sexually antagonistic selection, which favours different optima in males and females, is predicted to play an important role in the evolution of sex chromosomes. Body size is a sexually antagonistic trait in the shell-brooding cichlid fish Lamprologous callipterus, as "bourgeois" males must be large enough to carry empty snail shells to build nests whereas females must be small enough to fit into shells for breeding. In this species, there is also a second male morph: smaller "dwarf" males employ an alternative reproductive strategy by wriggling past spawning females into shells to fertilize eggs. L.callipterus male morphology is passed strictly from father to son, suggesting Y-linkage. However, sex chromosomes had not been previously identified in this species, and the genomic basis of size dimorphism was unknown. Here we used whole-genome sequencing to identify a 2.4-Mb sex-linked region on scaffold_23 with reduced coverage and single nucleotide polymorphism density in both male morphs compared to females. Within this sex region, distinct Y-haplotypes delineate the two male morphs, and candidate genes for body size (GHRHR, a known dwarfism gene) and sex determination (ADCYAP1R1) are in high linkage disequilibrium. Because differences in body size between females and males are under strong selection in L.callipterus, we hypothesize that sexual antagonism over body size initiated early events in sex chromosome evolution, followed by Y divergence to give rise to bourgeois and dwarf male reproductive strategies. Our results are consistent with the hypothesis that sexually antagonistic traits should be linked to young sex chromosomes.

De Novo Assembly and Annotation of 11 Diverse Shrub Willow (Salix) Genomes Reveals Novel Gene Organization in Sex-Linked Regions.

Poplar and willow species in the Salicaceae are dioecious, yet have been shown to use different sex determination systems located on different chromosomes. Willows in the subgenus Vetrix are interesting for comparative studies of sex determination systems, yet genomic resources for these species are still quite limited. Only a few annotated reference genome assemblies are available, despite many species in use in breeding programs. Here we present de novo assemblies and annotations of 11 shrub willow genomes from six species. Copy number variation of candidate sex determination genes within each genome was characterized and revealed remarkable differences in putative master regulator gene duplication and deletion. We also analyzed copy number and expression of candidate genes involved in floral secondary metabolism, and identified substantial variation across genotypes, which can be used for parental selection in breeding programs. Lastly, we report on a genotype that produces only female descendants and identified gene presence/absence variation in the mitochondrial genome that may be responsible for this unusual inheritance.

Transcriptomic thermal plasticity underlying gonadal development in a turtle with ZZ/ZW sex chromosomes despite canalized genotypic sex determination.

Understanding genome-wide responses to environmental conditions during embryogenesis is essential for discerning the evolution of developmental plasticity and canalization, two processes generating phenotypic variation targeted by natural selection. Here, we present the first comparative trajectory analysis of matched transcriptomic developmental time series from two reptiles incubated under identical conditions, a turtle with a ZZ/ZW system of genotypic sex determination (GSD), Apalone spinifera, and a turtle with temperature-dependent sex determination (TSD), Chrysemys picta. Results from our genome-wide, hypervariate gene expression analysis of sexed embryos across five developmental stages revealed that substantial transcriptional plasticity in the developing gonads can persist for >145 Myr, long after the canalization of sex determination via the evolution of sex chromosomes, while some gene-specific thermal sensitivity drifts or evolves anew. Such standing thermosensitivity represents an underappreciated evolutionary potential harbored by GSD species that may be co-opted during future adaptive shifts in developmental programing, such as a GSD to TSD reversal, if favored by ecological conditions. Additionally, we identified novel candidate regulators of vertebrate sexual development in GSD reptiles, including sex-determining candidate genes in a ZZ/ZW turtle.

New insights from GWAS on BMI-related growth traits in a longitudinal cohort of admixed children with Native American and European ancestry

Body-mass index (BMI) is a hallmark of adiposity. In contrast with adulthood, the genetic architecture of BMI during childhood is poorly understood. The few genome-wide association studies (GWAS) on children have been performed almost exclusively in Europeans and at single ages. We performed cross-sectional and longitudinal GWAS for BMI-related traits on 904 admixed children with mostly Mapuche Native American and European ancestries. We found regulatory variants of the immune gene HLA-DQB3 strongly associated with BMI at years old. A variant in the sex-determining gene DMRT1 was associated with the age at adiposity rebound (Age-AR) in girls (P ). BMI was significantly higher in Mapuche than in Europeans between 5.5 and 16.5 years old. Finally, Age-AR was significantly lower (P ) by 1.94 years and BMI at AR was significantly higher (P ) by 1.2kg/ , in Mapuche children compared with Europeans.

The function and evolution of a genetic switch controlling sexually dimorphic eye differentiation in honeybees

Animals develop sex-specific morphological structures that are diverse between organisms. However, understanding the developmental and evolutionary mechanisms governing these traits is still limited and largely restricted to DM domain genes, which are conserved, sex-specific developmental regulators identified in genetic models. Here, we report a sex-specific developmental regulator gene, glubschauge (glu) that selectively regulates sexually dimorphic eye differentiation in honeybees. We found that the sex determination gene feminizer (fem) controls sex-specific splicing of glu transcripts, establishing a genetic switch in which Glu proteins with a zinc finger (ZnF) domain are only expressed in females. We showed that female coding sequence was essential and sufficient for partial feminization. Comparative sequence and functional studies revealed that the evolutionary origination of the genetic switch was followed by the mutational origin of the essential ZnF domain. Our results demonstrate that glu is a newly evolved sex-specific genetic switch for region-specific regulation of a dimorphic character.

The fate of recessive deleterious or overdominant mutations near mating-type loci under partial selfing

Large regions of suppressed recombination having extended over time occur in many organisms around genes involved in mating compatibility (sex-determining or mating-type genes). The sheltering of deleterious alleles has been proposed to be involved in such expansions. However, the dynamics of deleterious mutations partially linked to genes involved in mating compatibility are not well understood, especially in finite populations. In particular, under what conditions deleterious mutations are likely to be maintained for long enough near mating-compatibility genes remains to be evaluated, especially under selfing, which generally increases the purging rate of deleterious mutations. Using a branching process approximation, we studied the fate of a new deleterious or overdominant mutation in a diploid population, considering a locus carrying two permanently heterozygous mating-type alleles, and a partially linked locus at which the mutation appears. We obtained analytical and numerical results on the probability and purging time of the new mutation. We investigated the impact of recombination between the two loci and of the mating system (outcrossing, intra and inter-tetrad selfing) on the maintenance of the mutation. We found that the presence of a fungal-like mating-type locus (i.e. not preventing diploid selfing) always sheltered the mutation under selfing, i.e. it decreased the purging probability and increased the purging time of the mutations. The sheltering effect was higher in case of automixis (intra-tetrad selfing). This may contribute to explain why evolutionary strata of recombination suppression near the mating-type locus are found mostly in automictic (pseudo-homothallic) fungi. We also showed that rare events of deleterious mutation maintenance during strikingly long evolutionary times could occur, suggesting that deleterious mutations can indeed accumulate near the mating-type locus over evolutionary time scales. In conclusion, our results show that, although selfing purges deleterious mutations, these mutations can be maintained for very long times near a mating-type locus, which may contribute to promote the evolution of recombination suppression in sex-related chromosomes.

Why and how do Y chromosome stop recombining?

Why and how do Y chromosome stop recombining?

Correction to: 'Identification of sex chromosome and sex-determining gene of southern catfish based on XX, XY and YY genome sequencing' (2022) by Zheng et al.

Correction to: 'Identification of sex chromosome and sex-determining gene of southern catfish based on XX, XY and YY genome sequencing' (2022) by Zheng et al.