Sex-specific Splicing Research Articles

Sex-specific mRNA alternative splicing patterns and Dmrt1 isoforms contribute to sex determination and differentiation of oyster.

Alternative splicing (AS) of pre-mRNA is a crucial mechanism that regulates the expression of genes involved in sex determination and differentiation. Despite its importance, AS has been rarely characterized in molluscs. In this study, PacBio Iso-Seq was employed to obtain full-length transcriptome and unveil AS patterns of gonads in the Pacific oyster Crassostrea gigas. A total of 24,783 AS events were identified across 6259 genes, with many enriched in phosphorylation-related processes. Splicing factors were found to drive a high frequency of AS events in gonads. Significant sex-based differences in isoform abundance and the incidence of AS events were observed. Comparative analysis of mature female and male gonads revealed a subset of overlapping differential alternative splicing genes and differentially expressed genes enriched in processes related to microtubule function and cell motility. In addition, the expression levels of sex-biased genes were found correlated with their isoform number in both female and male gonads. A novel isoform of Dmrt1 was identified with male specific expression in mature gonads. This study provides the first comprehensive understanding of full-length transcriptome and AS patterns in molluscan gonads, shedding light on the post-transcriptional regulatory mechanisms underlying sex determination and differentiation in molluscs and potentially across other animals.

Establishing a dominant early larval sex-selection strain in the Asian malaria vector Anopheles stephensi

BackgroundGenetic biocontrol interventions targeting mosquito-borne diseases require the release of male mosquitoes exclusively, as only females consume blood and transmit pathogens. Releasing only males eliminates the risk of increasing mosquito bites and spreading pathogens while enabling effective population control. The aim of this study is to develop robust sex-sorting methods for early larval stages in mosquitoes, enabling scalable male-only releases for genetic biocontrol interventions.MethodsTo address the challenge of sex-sorting in the Asian malaria vector Anopheles stephensi, we engineer Sexing Element Produced by Alternative RNA-splicing of a Transgenic Observable Reporter (SEPARATOR). This dominant fluorescent-based method, previously proven effective in Aedes aegypti, exploits sex-specific alternative splicing of a reporter to ensure exclusive male-specific expression early in development. The sex-specific alternative RNA splicing of the doublesex gene was selected as a target for engineering SEPARATOR due to its evolutionary conservation in insects. To expand SEPARATOR’s applicability for genetic sexing, we assessed the cross-species sex-specific RNA splicing activity of the An. gambiae doublesex (AngDsx) splicing module in An. stephensi. Male-specific enhanced green fluorescent protein (EGFP) expression was verified throughout the mosquito life cycle using a fluorescent stereomicroscope.ResultsOur results confirm that SEPARATOR regulates male-specific EGFP expression in An. stephensi and enables reliable positive male selection from the first instar larval stages. Molecular analysis demonstrates that male-specific EGFP expression is dependent on doublesex sex-specific splicing events. Additionally, the splicing module from An. gambiae operates effectively in An. stephensi, demonstrating evolutionary conservation in sex-specific splicing events between these species.ConclusionsSEPARATOR’s independence from sex-chromosome linkage provides resistance to breakage that could be mediated by meiotic recombination and chromosomal rearrangements, making it highly suitable for mass male releases. By enabling precise male selection from the first instar larval stages, SEPARATOR represents a significant advancement that will aid in the genetic biocontrol for Anopheles mosquitoes.Graphical

Hybridization between Aedes aegypti and Aedes mascarensis mosquitoes leads to disruption of male sex determination

Understanding the sex determination pathway and its disruptions in mosquitoes is critical for the effective control of disease vectors through genetic manipulations based on sex separation. When male hybrids of Aedes aegypti females and Ae. mascarensis males are backcrossed to Ae. aegypti females, a portion of the backcross progeny manifests as males with abnormal sexual differentiation. We discovered a significant correlation between pupal abnormalities and the feminization of subsequent adults exemplified by the relative abundance of ovarian and testicular tissues. All intersex individuals were genetic males as they expressed a male determining factor, Nix. Further, our analysis of the sex-specific splicing of doublesex and fruitless transcripts demonstrated the presence of both male and female splice variants indicating that sex determination is disrupted. A comparative transcriptomic analysis revealed similar expression levels of most female-associated genes in reproductive organs and carcasses between intersexual males and normal females. Moreover, intersexes had largely normal gene expression in testes but significant gene downregulation in male accessory glands when compared with normal males. We conclude that evolving hybrid incompatibilities between Ae. aegypti and Ae. mascarensis involve disruption of sex determination and are accompanied by changes in gene expression associated with sexual differentiation.

Establishing a Male-Positive Genetic Sexing Strain in the Asian Malaria Vector Anopheles stephensi.

Genetic biocontrol interventions targeting mosquito-borne diseases require the release of male mosquitoes exclusively, as only females consume blood and transmit human pathogens. This reduces the risk of spreading pathogens while enabling effective population control. Robust sex sorting methods to enable early larval sorting in mosquitoes need to be developed to allow for scalable sex sorting for genetic biocontrol interventions. This study applies the SEPARATOR (Sexing Element Produced by Alternative RNA-splicing of A Transgenic Observable Reporter) system, previously developed for Aedes aegypti, to the Asian malaria vector Anopheles stephensi. We hypothesized that the intron from the doublesex gene in Anopheles gambiae would function in An. stephensi due to evolutionary conservation. Our results confirm that the splicing module from An. gambiae operates effectively in An. stephensi, demonstrating evolutionary conservation in sex-specific splicing events between these species. This system enables reliable positive male selection from first instar larval to pupal stages. RT-PCR analysis demonstrates that male-specific EGFP expression is dependent on doublesex sex-specific splicing events. The SEPARATOR system's independence from sex-chromosome linkage confers resistance to meiotic recombination and chromosomal rearrangements. This approach may facilitate the mass release of males, and the cross-species portability of SEPARATOR establishes it as a valuable tool for genetic biocontrol interventions across various pest species.

Sexual dimorphism in the tardigrade Paramacrobiotus metropolitanus transcriptome

BackgroundIn gonochoristic animals, the sex determination pathway induces different morphological and behavioral features that can be observed between sexes, a condition known as sexual dimorphism. While many components of this sex differentiation cascade show high levels of diversity, factors such as the Doublesex-Mab-3-Related Transcription factor (DMRT) are widely conserved across animal taxa. Species of the phylum Tardigrada exhibit remarkable diversity in morphology and behavior between sexes, suggesting a pathway regulating this dimorphism. Despite the wealth of genomic and zoological knowledge accumulated in recent studies, the sexual differences in tardigrades genomes have not been identified. In the present study, we focused on the gonochoristic species Paramacrobiotus metropolitanus and employed omics analyses to unravel the molecular basis of sexual dimorphism.ResultsTranscriptome analysis between sex-identified specimens revealed numerous differentially expressed genes, of which approximately 2,000 male-biased genes were focused on 29 non-male-specific genomic loci. From these regions, we identified two Macrobiotidae family specific DMRT paralogs, which were significantly upregulated in males and lacked sex specific splicing variants. Furthermore, phylogenetic analysis indicated all tardigrade genomes lack the doublesex ortholog, suggesting doublesex emerged after the divergence of Tardigrada. In contrast to sex-specific expression, no evidence of genomic differences between the sexes was found. We also identified several anhydrobiosis genes that exhibit sex-biased expression, suggesting a possible mechanism for protection of sex-specific tissues against extreme stress.ConclusionsThis study provides a comprehensive analysis for analyzing the genetic differences between sexes in tardigrades. The existence of male-biased, but not male-specific, genomic loci and identification of the family specific male-biased DMRT subfamily provides the foundation for understanding the sex determination cascade. In addition, sex-biased expression of several tardigrade-specific genes which are involved their stress tolerance suggests a potential role in protecting sex-specific tissue and gametes.

The evolution of sexual dimorphism in gene expression in response to a manipulation of mate competition.

Many genes are differentially expressed between males and females and patterns of sex-biased gene expression (SBGE) vary among species. Some of this variation is thought to have evolved in response to differences in mate competition among species that cause varying patterns of sex-specific selection. We used experimental evolution to test this by quantifying SBGE and sex-specific splicing in 15 Drosophila melanogaster populations that evolved for 104 generations in mating treatments that removed mate competition via enforced monogamy, or allowed mate competition in either small, simple, or larger, structurally more complex mating environments. Consistent with sex-specific selection affecting SBGE, initially sex-biased genes diverged in expression more among treatments than unbiased genes, and there was greater expression divergence for male- than female-biased genes. It has been suggested the transcriptome should be "feminized" under monogamy because of the removal of sexual selection on males; we did not observe this, likely because selection differs in additional ways between monogamy vs. polygamy. Significant divergence in average expression dimorphism between treatments was observed and, in some treatment comparisons, the direction of the divergence differed across different sex-bias categories. There was not a generalized reduction in expression dimorphism under enforced monogamy.

Efficient sex separation by exploiting differential alternative splicing of a dominant marker in Aedes aegypti.

Only female mosquitoes consume blood giving them the opportunity to transmit deadly human pathogens. Therefore, it is critical to remove females before conducting releases for genetic biocontrol interventions. Here we describe a robust sex-sorting approach termed SEPARATOR (Sexing Element Produced by Alternative RNA-splicing of A Transgenic Observable Reporter) that exploits sex-specific alternative splicing of an innocuous reporter to ensure exclusive dominant male-specific expression. Using SEPARATOR, we demonstrate reliable sex selection from early larval and pupal stages in Aedes aegypti, and use a Complex Object Parametric Analyzer and Sorter (COPAS) to demonstrate scalable high-throughput sex-selection of first instar larvae. Additionally, we use this approach to sequence the transcriptomes of early larval males and females and find several genes that are sex-specifically expressed. SEPARATOR can simplify mass production of males for release programs and is designed to be cross-species portable and should be instrumental for genetic biocontrol interventions.

The sex-specific factor SOA controls dosage compensation in Anopheles mosquitos.

The Anopheles mosquito is one of thousands of species in which sex differences play a central role in their biology, as only females need a blood meal in order to produce eggs. Sex differentiation is regulated by sex chromosomes, but their presence creates a dosage imbalance between males (XY) and females (XX). Dosage compensation (DC) can re-equilibrate the expression of sex-chromosomal genes, but because DC mechanisms have only been fully characterized in a few model organisms, key questions about its evolutionary diversity and functional necessity remain unresolved 1. Here we report the discovery of a previously uncharacterized gene (SOA, for sex chromosome activation) as a master regulator of DC in the malaria mosquito Anopheles gambiae. Sex-specific alternative splicing prevents functional SOA protein expression in females. The male isoform encodes a DNA-binding protein that binds the promoters of active X chromosomal genes. Expressing male SOA is sufficient to induce DC in female cells. Male mosquitoes lacking SOA or female mosquitos ectopically expressing the male isoform exhibit X chromosome misregulation, which is compatible with viability but causes developmental delay. Thus, our molecular analysis of the first DC master regulator in a non-model organism elucidates the evolutionary steps leading to the establishment of a chromosome-specific fine-tuning mechanism.

The identification and expression pattern of the sex determination genes and their sex-specific variants in the egg parasitoid Trichogramma dendrolimi Matsumura (Hymenoptera: Trichogrammatidae).

Introduction: Trichogramma wasps are egg parasitoids of agricultural lepidopteran pests. The sex of Trichogramma is determined by its ploidy as well as certain sex ratio distorters, such as the endosymbiotic bacteria Wolbachia spp. and the paternal sex ratio (PSR) chromosome. The sex determination systems of hymenopterans, such as Trichogramma spp., involve cascades of the genes transformer (tra), transformer-2 (tra2), and doublesex (dsx) and are associated with sex-specific tra and dsx splicing. First, these genes and their sex-specific variants must be identified to elucidate the interactions between the sex ratio disorders and the sex determination mechanism of Trichogramma. Methods: Here, we characterized the sex determination genes tra, tra2, and dsx in Trichogramma dendrolimi. Sex-specific tra and dsx variants were detected in cDNA samples obtained from both male and female Trichogramma wasps. They were observed in the early embryos (1-10h), late embryos (12-20h), larvae (32h and 48h), pre-pupae (96h), and pupae (144h, 168h, 192h, and 216h) of both male and female T. dendrolimi offspring. Results: We detected female-specific tra variants throughout the entire early female offspring stage. The male-specific variant began to express at 9-10h as the egg was not fertilized. However, we did not find any maternally derived, female-specific tra variant in the early male embryo. This observation suggests that the female-specific tra variant expressed in the female embryo at 1-9h may not have originated from the maternal female wasp. Discussion: The present study might be the first to identify the sex determination genes and sex-specific gene splicing in Trichogramma wasps. The findings of this study lay the foundation for investigating the sex determination mechanisms of Trichogramma and other wasps. They also facilitate sex identification in immature T. dendrolimi and the application of this important egg parasitoid in biological insect pest control programs.

Sexual dimorphic eyestalk transcriptome of kuruma prawn Marsupenaeus japonicus

Kuruma prawn (Marsupenaeus japonicus) is a benthic decapod crustacean that is widely distributed in the Indo-West Pacific region. It is one of the most important fishery resources in Japan, but its annual catches have declined sharply since the 1990s. To increase stocks, various approaches such as seed production and aquaculture were attempted. Since the demand for important fishery species, including kuruma prawn, is expected to increase worldwide in the future, there is a need to develop new technologies that will make aquaculture more efficient. Historically, the eyestalk endocrine organ is known to consist of the X-organ and sinus gland (XO/SG) complex that synthesizes and secrets various neuropeptide hormones that regulate growth, molt, sexual maturation, reproduction, and changes in body color. In the current study, eyestalk-derived neuropeptides were identified in the transcriptome. In addition, most orthologs of sex-determination genes were expressed in eyestalks. We identified two doublesex genes (MjapDsx1 and MjapDsx2) and found that MjapDsx1 showed male-biased expression in the eyestalk ganglion with no sex-specific splicing, unlike insect species. Therefore, this study will provide an opportunity to advance the research of neuropeptides and sex determination in the kuruma prawn.

Sex-specific splicing occurs genome-wide during early Drosophila embryogenesis.

Sex-specific splicing is an essential process that regulates sex determination and drives sexual dimorphism. Yet, how early in development widespread sex-specific transcript diversity occurs was unknown because it had yet to be studied at the genome-wide level. We use the powerful Drosophila model to show that widespread sex-specific transcript diversity occurs early in development, concurrent with zygotic genome activation. We also present a new pipeline called time2Splice to quantify changes in alternative splicing over time. Furthermore, we determine that one of the consequences of losing an essential maternally deposited pioneer factor called CLAMP (chromatin-linked adapter for MSL proteins) is altered sex-specific splicing of genes involved in diverse biological processes that drive development. Overall, we show that sex-specific differences in transcript diversity exist even at the earliest stages of development..

Female Sex Determination Factors in Ceratitis capitata: Molecular and Structural Basis of TRA and TRA2 Recognition.

In the model system for genetics, Drosophila melanogaster, sexual differentiation and male courtship behavior are controlled by sex-specific splicing of doublesex (dsx) and fruitless (fru). In vitro and in vivo studies showed that female-specific Transformer (TRA) and the non-sex-specific Transformer 2 (TRA2) splicing factors interact, forming a complex promoting dsx and fru female-specific splicing. TRA/TRA2 complex binds to 13 nt long sequence repeats in their pre-mRNAs. In the Mediterranean fruitfly Ceratitis capitata (Medfly), a major agricultural pest, which shares with Drosophila a ~120 million years old ancestor, Cctra and Cctra2 genes seem to promote female-specific splicing of Ccdsx and Ccfru, which contain conserved TRA/TRA2 binding repeats. Unlike Drosophila tra, Cctra autoregulates its female-specific splicing through these putative regulatory repeats. Here, a yeast two-hybrid assay shows that CcTRA interacts with CcTRA2, despite its high amino acid divergence compared to Drosophila TRA. Interestingly, CcTRA2 interacts with itself, as also observed for Drosophila TRA2. We also generated a three-dimensional model of the complex formed by CcTRA and CcTRA2 using predictive approaches based on Artificial Intelligence. This structure also identified an evolutionary and highly conserved putative TRA2 recognition motif in the TRA sequence. The Y2H approach, combined with powerful predictive tools of three-dimensional protein structures, could use helpful also in this and other insect species to understand the potential links between different upstream proteins acting as primary sex-determining signals and the conserved TRA and TRA2 transducers.

The draft genome sequence of the Japanese rhinoceros beetle Trypoxylus dichotomus septentrionalis towards an understanding of horn formation

The Japanese rhinoceros beetle Trypoxylus dichotomus is a giant beetle with distinctive exaggerated horns present on the head and prothoracic regions of the male. T. dichotomus has been used as a research model in various fields such as evolutionary developmental biology, ecology, ethology, biomimetics, and drug discovery. In this study, de novo assembly of 615 Mb, representing 80% of the genome estimated by flow cytometry, was obtained using the 10 × Chromium platform. The scaffold N50 length of the genome assembly was 8.02 Mb, with repetitive elements predicted to comprise 49.5% of the assembly. In total, 23,987 protein-coding genes were predicted in the genome. In addition, de novo assembly of the mitochondrial genome yielded a contig of 20,217 bp. We also analyzed the transcriptome by generating 16 RNA-seq libraries from a variety of tissues of both sexes and developmental stages, which allowed us to identify 13 co-expressed gene modules. We focused on the genes related to horn formation and obtained new insights into the evolution of the gene repertoire and sexual dimorphism as exemplified by the sex-specific splicing pattern of the doublesex gene. This genomic information will be an excellent resource for further functional and evolutionary analyses, including the evolutionary origin and genetic regulation of beetle horns and the molecular mechanisms underlying sexual dimorphism.

Two forms of sexual dimorphism in gene expression in Drosophila melanogaster: their coincidence and evolutionary genetics.

Phenotypic sexual dimorphism can be mediated by sex differences in gene expression. We examine two forms of sexual dimorphism in gene expression in Drosophila melanogaster: (i) sex-biased gene expression (SBGE) in which the sexes differ in the amount a gene is expressed and (ii) sexual dimorphism in isoform usage, i.e., sex-specific splicing (SSS). In whole body (but not head) expression, we find a negative association between SBGE and SSS, possibly suggesting these are alternate routes to resolving sexual antagonistic selection. Next, we evaluate whether expression dimorphism contributes to the heterogeneity among genes in rmf, the intersexual genetic correlation in body expression that constrains the extent to which a gene's expression can evolve independently between the sexes. We find lower rmf values for genes with than without SSS. We find higher rmf values for male- than female-biased genes (except genes with extreme male-bias), even though male-biased genes are known to have greater evolutionary divergence in expression. Finally, we examine population genetic patterns in relation to SBGE and SSS because genes with expression dimorphism have likely experienced a history of sex differences in selection. SSS is associated with reduced values of Tajima's D and elevated Direction of Selection (DoS) values, suggestive of higher rates of adaptive evolution. Though DoS is highly elevated for genes with extreme male bias, DoS otherwise tends to decline from female-biased to unbiased to male-biased genes. Collectively, the results indicate that SBGE and SSS are differentially distributed across the genome and are associated with different forms of selection.

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.

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.

Cell-based analysis reveals that sex-determining gene signals in Ostrinia are pivotally changed by male-killing Wolbachia.

Wolbachia, a maternally transmitted bacterium, shows male-killing, an adaptive phenotype for cytoplasmic elements, in various arthropod species during the early developmental stages. In lepidopteran insects, lethality of males is accounted for by improper dosage compensation in sex-linked genes owing to Wolbachia-induced feminization. Herein, we established Ostrinia scapulalis cell lines that retained sex specificity per the splicing pattern of the sex-determining gene doublesex (Osdsx). We found that Wolbachia transinfection in male cell lines enhanced the female-specific splice variant of Osdsx (OsdsxF ) while suppressing the male-specific variant (OsdsxM ), indicating that Wolbachia affects sex-determining gene signals even in vitro. Comparative transcriptome analysis isolated only two genes that behave differently upon Wolbachia infection. The two genes were respectively homologous to Masculinizer (BmMasc) and zinc finger-2 (Bmznf-2), male-specifically expressed sex-determining genes of the silkworm Bombyx mori that encode CCCH-type zinc finger motif proteins. By using cultured cells and organismal samples, OsMasc and Osznf-2 were found to be sex-determining genes of O. scapulalis that are subjected to sex-specific alternative splicing depending upon the chromosomal sex, developmental stage, and infection status. Overall, our findings expound the cellular autonomy in insect sex determination and the mechanism through which sex is manipulated by intracellular selfish microbes.

Identification of sex-specific splicing via comparative transcriptome analysis of gonads from sea cucumber Apostichopus japonicus

Alternative splicing (AS) is an essential post-transcriptional regulation mechanism for sex differentiation and gonadal development, which has rarely been reported in marine invertebrates. Sea cucumber (Apostichopus japonicus) is an economically important marine benthic echinoderm with a potential XX/XY sex determination mechanism, whose molecular mechanism in the gonadal differentiation has not been clearly understood. In this study, we analyzed available RNA-seq datasets of male and female gonads to explore if AS mechanism exerts an essential function in sex differentiation and gonadal development of A. japonicus. In our results, a total of 20,338 AS events from 7219 alternatively spliced genes, and 189 sexually differential alternative splicing (DAS) events from 156 genes were identified in gonadal transcriptome of sea cucumber. Gene Ontology analysis indicated that these DAS genes were significantly enriched in spermatogenesis-related GO terms. Maximal Clique Centrality (MCC) was then applied for protein-protein interaction (PPI) analysis to search for protein interactions and hub DAS gene. Among all DAS genes, we identified 10 DAS genes closely related to spermatogenesis and (or) sperm motility and a hub gene dnah1. Thus, this study revealed that alternative isoforms were generated from certain genes in female and male gonads through alternative splicing, which may provide direct evidence that alternative splicing mechanisms participate in female and male gonads. These results suggested a novel perspective for explaining the molecular mechanisms underlying gonadal differentiation between male and female sea cucumbers.

The Sex-Specific Splicing of Doublesex in Brine Shrimp Artemia franciscana.

The understanding of sex determination and differentiation in animals has recently made remarkable strides through the use of advanced research tools. At the gene level, the Mab-3-related transcription factor (Dmrt) gene family, which encodes for the typical DNA-binding doublesex/Mab-3 (DM) domain in their protein, is known for its contribution to sex determination and differentiation in insects. In this study, DNA-binding DM domain screening has identified eight transcripts from Artemia franciscana transcriptomic that encode proteins containing one conserved DNA-binding DM domain. The genome mapping confirmed that these eight transcripts are transcribed from six different loci on the A. franciscana genome assembly. One of those loci, the Af.dsx-4 locus, is closely related to Doublesex, a gene belonging to the Dmrt gene family. This locus could be transcribed into three alternative transcripts, namely Af.dsx4, Af.dsxF and Af.dsxM. While Af.dsx4 and Af.dsxF could putatively be translated to form an identical Af.dsxF protein of 186 aa long, Af.dsxM translates for an Af.dsxM protein of 289 aa long but shares a DNA-binding DM domain. Interestingly, Af.dsxF and Af.dsxM are confirmed as sex-specific transcripts, Af.dsxF is only present in females, and Af.dsxM is only present in male individuals. The results suggest that the sex-specific splicing mechanism of the doublesex described in insects is also present in A. franciscana. Af.dxs-4 locus can be used in further studies to clarify the sex determination pathways in A. fracnciscana.

Comprehensive Transcriptome Analysis Reveals Sex-Specific Alternative Splicing Events in Zebrafish Gonads.

Alternative splicing is an important way of regulating gene functions in eukaryotes. Several key genes involved in sex determination and gonadal differentiation, such as nr5a1 and ddx4, have sex-biased transcripts between males and females, suggesting a potential regulatory role of alternative splicing in gonads. Currently, the sex-specific alternative splicing events and genes have not been comprehensively studied at the genome-wide level in zebrafish. In this study, through global splicing analysis on three independent sets of RNA-seq data from matched zebrafish testes and ovaries, we identified 120 differentially spliced genes shared by the three datasets, most of which haven’t been reported before. Functional enrichment analysis showed that the GO terms of mRNA processing, mRNA metabolism and microtubule-based process were strongly enriched. The testis- and ovary-biased alternative splicing genes were identified, and part of them (tp53bp1, tpx2, mapre1a, kif2c, and ncoa5) were further validated by RT-PCR. Sequence characteristics analysis suggested that the lengths, GC contents, and splice site strengths of the alternative exons or introns may have different influences in different types of alternative splicing events. Interestingly, we identified an unexpected high proportion (over 70%) of non-frameshift exon-skipping events, suggesting that in these cases the two protein isoforms derived from alternative splicing may both have functions. Furthermore, as a representative example, we found that the alternative splicing of ncoa5 causes the loss of a conserved RRM domain in the short transcript predominantly produced in testes. Our study discovers novel sex-specific alternative splicing events and genes with high reliabilities in zebrafish testes and ovaries, which would provide attractive targets for follow-up studies to reveal the biological significances of alternative splicing events and genes in sex determination and gonadal differentiation.