Role of Alternative Splicing in Sex Determination in Vertebrates

Overexpression of Kdm6b induces testicular differentiation in a temperature-dependent sex determination system

In reptiles, such as the red-eared slider turtle ( Trachemys scripta elegans), gonadal sex determination is highly dependent on the environmental temperature during embryonic stages. This complex process, which leads to differentiation into either testes or ovaries, is governed by the finely tuned expression of upstream genes, notably the testis-promoting gene Dmrt1 and the ovary-promoting gene Foxl2. Recent studies have identified epigenetic regulation as a crucial factor in testis development, with the H3K27me3 demethylase KDM6B being essential for Dmrt1 expression in T. s. elegans. However, whether KDM6B alone can induce testicular differentiation remains unclear. In this study, we found that overexpression of Kdm6b in T. s. elegans embryos induced the male development pathway, accompanied by a rapid increase in the gonadal expression of Dmrt1 at 31°C, a temperature typically resulting in female development. Notably, this sex reversal could be entirely rescued by Dmrt1 knockdown. These findings demonstrate that Kdm6b is sufficient for commitment to the male pathway, underscoring its role as a critical epigenetic regulator in the sex determination of the red-eared slider turtle.

RNA-clique: a method for computing genetic distances from RNA-seq data

BACKGROUND

Although RNA-seq data are traditionally used for quantifying gene expression levels, the same data could be useful in an integrated approach to compute genetic distances as well. Challenges to using mRNA sequences for computing genetic distances include the relatively high conservation of coding sequences and the presence of paralogous and, in some species, homeologous genes.

RESULTS

We developed a new computational method, RNA-clique, for calculating genetic distances using assembled RNA-seq data and assessed the efficacy of the method using biological and simulated data. The method employs reciprocal BLASTn followed by graph-based filtering to ensure that only orthologous genes are compared. Each vertex in the graph constructed for filtering represents a gene in a specific sample under comparison, and an edge connects a pair of vertices if the genes they represent are best matches for each other in their respective samples. The distance computation is a function of the BLAST alignment statistics and the constructed graph and incorporates only those genes that are present in some complete connected component of this graph. As a biological testbed we used RNA-seq data of tall fescue (Lolium arundinaceum), an allohexaploid plant ( 2 n = 14 Gb ), and bluehead wrasse (Thalassoma bifasciatum), a teleost fish. RNA-clique reliably distinguished individual tall fescue plants by genotype and distinguished bluehead wrasse RNA-seq samples by individual. In tests with simulated RNA-seq data, the ground truth phylogeny was accurately recovered from the computed distances. Moreover, tests of the algorithm parameters indicated that, even with stringent filtering for orthologs, sufficient sequence data were retained for the distance computations. Although comparisons with an alternative method revealed that RNA-clique has relatively high time and memory requirements, the comparisons also showed that RNA-clique's results were at least as reliable as the alternative's for tall fescue data and were much more reliable for the bluehead wrasse data.

CONCLUSION

Results of this work indicate that RNA-clique works well as a way of deriving genetic distances from RNA-seq data, thus providing a methodological integration of functional and genetic diversity studies.

Analysis of circRNA-miRNA-mRNA regulatory network of embryonic gonadal development in Mulard duck

The aim of the study was to explore the regulatory mechanism of differences in embryonic gonadal development between intergeneric distance hybrid offspring Mulard ducks and parent ducks. The morphological differences gonadal tissues of Muscovy ducks, Pekin ducks and Mulard ducks at 12.5-day embryonic age were observed by sectioning and hematoxylin-eosin (HE) staining. Then followed by transcriptome sequencing to screen for gonadal development-related differentially expressed circRNAs and mRNAs to construct a competitive endogenous RNA (ceRNA) regulatory network. Finally, qRT-PCR and luciferase reporter system were used to verify the sequencing data and targeting relationship of ceRNA pairs. The results showed that the seminiferous tubule lumen of Mulard ducks was not obvious, while there were obvious seminiferous tubules and tubular structures in testis of Pekin ducks and Muscovy ducks, with number and shape indicating maturity. There were 18 upregulated circRNAs and 16 downregulated circRNAs in Mulard ducks and Pekin ducks, respectively, and 39 upregulated circRNAs and 1 downregulated circRNA in Mulard ducks and Muscovy ducks, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis found that genes involves in dorso-ventral axis formation, for example, neurogenic locus notch homolog protein 1 (NOTCH1), were significantly enriched (P < 0.05). The novel_circ_0002265-gga-miR-122-5p-PAFAH1B2 regulatory network was constructed. The qRT-PCR results showed that the sequencing results were reliable. The dual-luciferase reporter assay showed that gga-miR-122-5p exists binding site of circ_0002265 and PAFAH1B2, indicating circ_0002265-gga-miR-122-5p-PAFAH1B2 targeting relationship. In summary, the embryonic gonadal development of intergeneric hybrid Mulard ducks may be regulated by differentially expressed circRNAs and genes, such as novel_circ_0000519, novel_circ_0003537, NOTCH1, FGFR2, PAFAH1B1, and PAFAH1B2, among which circ_0002265-gga-miR-122-5p-PAFAH1B2 may participate in the targeted regulation of gonadal development in Mulard ducks. The findings of this study are helpful for analyzing the mechanism of embryonic gonadal development differences in avians.

Alternative splicing of histone demethylase Kdm6bb mediates temperature-induced sex reversal in the Nile tilapia

Sex determination in many fish species is remarkably plastic and temperature sensitive. Nile tilapia display a genetic sex-determination system (XX/XY). However, high-temperature treatment during critical thermosensitive periods can induce XX females into XXm pseudo-males, and this phenomenon is termed temperature-induced sex reversal (TISR). To investigate the molecular mechanism of TISR in Nile tilapia, we performed Iso-seq analysis and found a dramatic effect of high temperature on gene alternative splicing (AS). Kdm6bb histone demethylase showed a novel AS at intron 5 that generates Kdm6bb_tv1 transcripts without intron 5 and Kdm6bb_tv2 with intron 5. Kdm6bb_tv1 encodes a full-length protein while Kdm6bb_tv2 encodes a truncated protein. Expression analysis revealed that intron 5 splicing of Kdm6bb is male and gonad biased at larval stage, and only gonad biased at adult stage. High-temperature treatment induced intron 5 splicing in the gonads of XX and XY fish, resulting in increased Kdm6bb_tv1 expression. To directly test the role of Kdm6bb_tv1 in Nile tilapia TISR, we knocked out expression of Kdm6bb_tv1. However, Kdm6bb_tv1-/- homozygous mutants showed embryonic lethality. Overexpression of Kdm6bb_tv1, but not Kdm6bb_tv2, induced sex reversal of XX females into pseudo-males. Overexpression of Kdm6bb_tv1, as with high-temperature treatment, modified the promotor region of Gsdf and Dmrt1 by demethylating the trimethylated lysine 27 of histone 3 (H3K27me3), thereby increasing expression. Collectively, these studies demonstrate that AS of Kdm6bb intron 5 increases the expression of Kdm6bb_tv1, which acts as a direct link between high temperature and activation of Gsdf and Dmrt1 expression, leading to male sex determination.

Aquaporin splice variation differentially modulates channel function during marine teleost egg hydration

Aquaporin-mediated oocyte hydration is a developmentally regulated adaptive mechanism that co-occurs with meiosis resumption in marine teleosts. It provides the early embryos with vital water until osmoregulatory systems develop, and in the majority of marine teleosts causes their eggs to float. Recent studies have shown that the subdomains of two water channels (Aqp1ab1 and Aqp1ab2) encoded in a teleost-specific aquaporin-1 cluster (TSA1C) co-evolved with duplicated Ywhaz-like (14-3-3ζ-like) binding proteins to differentially control their membrane trafficking for maximal egg hydration. Here, we report that in species that encode the full TSA1C, in-frame intronic splice variants of Aqp1ab1 result in truncated proteins that cause dominant-negative inhibition of the canonical channel trafficking to the plasma membrane. The inhibition likely occurs through hetero-oligomerization and retention in the endoplasmic reticulum (ER) and ultimate degradation. Conversely, in species that only encode the Aqp1ab2 channel we found an in-frame intronic splice variant that results in an intact protein with an extended extracellular loop E, and an out-of frame intronic splice variant with exon readthrough that results in a truncated protein. Both isoforms cause dominant-negative enhancement of the degradation pathway. However, the extended and truncated Aqp1ab2-type variants can also partially escape from the ER to reach the oocyte plasma membrane, where they dominantly-negatively inhibit water flux. The ovarian follicular expression ratios of the Aqp1ab2 isoforms in relation to the canonical channel are lowest during oocyte hydration, but subsequently highest when the canonical channel is recycled, thus leaving the eggs endowed with >90% water. These findings suggest that the expression of inhibitory isoforms of Aqp1ab1 and Aqp1ab2 may represent a new regulatory mechanism through which the cell-surface expression and the activity of the canonical channels can be physiologically modulated during oocyte hydration in marine teleosts.

High-throughput RNA sequencing reveals differences between the transcriptomes of the five spore forms of Puccinia striiformis f. sp. tritici, the wheat stripe rust pathogen

The devastating wheat stripe (yellow) rust pathogen, Puccinia striiformis f. sp. tritici (Pst), is a macrocyclic and heteroecious fungus. Pst produces urediniospores and teliospores on its primary host, wheat, and pycniospores and aeciospores are produced on its alternate hosts, barberry (Berberis spp.) or mahonia (Mahonia spp.). Basidiospores are developed from teliospores and infect alternate hosts. These five spore forms play distinct roles in Pst infection, disease development, and fungal survival, etc. However, the specific genes and mechanisms underlying these functional differences are largely unknown. In this study, we performed, for the first time in rust fungi, the deep RNA sequencing to examine the transcriptomic shift among all five Pst spore forms. Among a total of 29,591 identified transcripts, 951 were specifically expressed in basidiospores, whereas 920, 761, 266, and 110 were specific for teliospores, pycniospores, aeciospores, and urediniospores, respectively. Additionally, transcriptomes of sexual spores, namely pycniospores and basidiospores, showed significant differences from those of asexual spores (urediniospores, teliospores, and aeciospores), and transcriptomes of urediniospores and aeciospores were more similar to each other than to the three other spore forms. Especially, the basidiospores and pycniospores which infected the berberis shows wide differences in the cell wall degrading-enzymes and mating and pheromone response genes. Besides, we also found that there are 6234 differential expressed genes between the urediniospores and pycniospores, while only have 3 genes have alternative splicing enents, suggesting that differential genes expression may make more contribution than AS. This comprehensive transcriptome profiling can substantially improve our understanding of the developmental biology of the wheat stripe rust fungus.

Virus usurps alternative splicing to clear the decks for infection

Since invasion, there will be a tug-of-war between host and virus to scramble cellular resources, for either restraining or facilitating infection. Alternative splicing (AS) is a conserved and critical mechanism of processing pre-mRNA into mRNAs to increase protein diversity in eukaryotes. Notably, this kind of post-transcriptional regulatory mechanism has gained appreciation since it is widely involved in virus infection. Here, we highlight the important roles of AS in regulating viral protein expression and how virus in turn hijacks AS to antagonize host immune response. This review will widen the understandings of host-virus interactions, be meaningful to innovatively elucidate viral pathogenesis, and provide novel targets for developing antiviral drugs in the future.

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.

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.

Temporal Dynamic Analysis of Alternative Splicing During Embryonic Development in Zebrafish

Alternative splicing is pervasive in mammalian genomes and involved in embryo development, whereas research on crosstalk of alternative splicing and embryo development was largely restricted to mouse and human and the alternative splicing regulation during embryogenesis in zebrafish remained unclear. We constructed the alternative splicing atlas at 18 time-course stages covering maternal-to-zygotic transition, gastrulation, somitogenesis, pharyngula stages, and post-fertilization in zebrafish. The differential alternative splicing events between different developmental stages were detected. The results indicated that abundance alternative splicing and differential alternative splicing events are dynamically changed and remarkably abundant during the maternal-to-zygotic transition process. Based on gene expression profiles, we found splicing factors are expressed with specificity of developmental stage and largely expressed during the maternal-to-zygotic transition process. The better performance of cluster analysis was achieved based on the inclusion level of alternative splicing. The biological function analysis uncovered the important roles of alternative splicing during embryogenesis. The identification of isoform switches of alternative splicing provided a new insight into mining the regulated mechanism of transcript isoforms, which always is hidden by gene expression. In conclusion, we inferred that alternative splicing activation is synchronized with zygotic genome activation and discovered that alternative splicing is coupled with transcription during embryo development in zebrafish. We also unveiled that the temporal expression dynamics of splicing factors during embryo development, especially co-orthologous splicing factors. Furthermore, we proposed that the inclusion level of alternative splicing events can be employed for cluster analysis as a novel parameter. This work will provide a deeper insight into the regulation of alternative splicing during embryogenesis in zebrafish.

Detecting signatures of selection on gene expression

A substantial amount of phenotypic diversity results from changes in gene expression levels and patterns. Understanding how the transcriptome evolves is therefore a key priority in identifying mechanisms of adaptive change. However, in contrast to powerful models of sequence evolution, we lack a consensus model of gene expression evolution. Furthermore, recent work has shown that many of the comparative approaches used to study gene expression are subject to biases that can lead to false signatures of selection. Here we first outline the main approaches for describing expression evolution and their inherent biases. Next, we bridge the gap between the fields of phylogenetic comparative methods and transcriptomics to reinforce the main pitfalls of inferring selection on expression patterns and use simulation studies to show that shifts in tissue composition can heavily bias inferences of selection. We close by highlighting the multi-dimensional nature of transcriptional variation and identifying major unanswered questions in disentangling how selection acts on the transcriptome.