Gonadal transcriptomes associated with sex phenotypes provide potential male and female candidate genes of sex determination or early differentiation in Crassostrea gigas, a sequential hermaphrodite mollusc

Functional verification of a landmark gene EVM713 involved in spermatogenesis in the marine bivalve Chlamys nobilis

The formation of broodstock gametes is closely linked to the yield and quality in aquaculture production, yet molecular mechanisms underlying this process remain insufficiently understood. The noble scallop Chlamys nobilis, an economically significant dioecious bivalve species, serves as an excellent model for studying gametogenesis. In this study, the adult scallops with testis at different developmental stage were chosen for histological examination and transcriptome analysis to dig genes related gonad development. Totally, 2663 DEGs and their set modules significantly related to spermatogenesis were obtained using WGCNA, including 40 candidate genes represented by EVM713. The gene was specifically expressed in the testis. RNA interference (RNAi) of EVM713 led to impaired testis development, marked by sparse sperm cell arrangement, spermatocytes detaching from the follicle wall, and reduced spermatocyte numbers. Meanwhile, 24 h after RNAi, the expression levels of Bax, and Caspase3 significantly increased (P < 0.05), while those of Bcl2, Dmrt2 and Tssk4 were significantly decreased (P < 0.05). These results indicate that EVM713 is essential for spermatogenesis in bivalves, regulating testis development through the modulation of Dmrt2 and Tssk4 expression. This study provides the first evidence of EVM713 function in mollusks, which is conducive to better understanding molecular mechanisms underlying gametogenesis in marine invertebrates.

Exposure to plastic debris alters expression of biomineralization, immune, and stress-related genes in the eastern oyster (Crassostrea virginica)

The degradation of marine plastic debris poses a threat to organisms by fragmenting into micro- and nano-scale pieces and releasing a complex chemical leachate into the water. Numerous studies have investigated harms from plastic pollution such as microplastic ingestion and exposure to single chemicals. However, few studies have examined the holistic threat of plastic exposure and the synergistic impacts of chemical mixtures. The objective of this study was to measure changes in gene expression of gill and gonadal tissue of the eastern oyster (Crassostrea virginica) in response to plastic debris exposure during their first year, using RNA-seq to explore multiple types of physiological responses. Shell and polyethylene terephthalate plastic were used as substrate for the metamorphosis of larval oysters in a settlement tank. Substrate pieces were then transferred to metal cages and outplanted in pairs - shell cage and plastic cage - onto restoration reefs in the St. Mary's River, Maryland, USA. After 10 months of growth, the oysters were collected, gill and gonadal tissue removed, and sex identified. The tissues of six oysters from each sex and substrate type were then analyzed in RNA-seq. Both gill and gonadal tissue samples had altered expression of immune and stress-response genes in response to plastic exposure. Genes upregulated in response to plastic were enriched for gene ontology functions of proteolysis and fibrinolysis. Downregulated genes were involved in shell biomineralization and growth. One male oyster exposed to plastic had "feminized" gene expression patterns despite developing mature sperm, suggesting plastic leachate can alter gene expression and shift protandric individuals to develop as females. Plastic pollution may therefore reduce shell growth, initiate immune and stress responses, alter sex differentiation, and impact reproductive output of eastern oysters through changes in transcription.

Integrated transcriptomic analysis reveals potential ceRNA network and hub genes involved in sex determination and differentiation of the Pacific oyster (Crassostrea gigas)

Non-coding RNA (ncRNA) and competing endogenous RNA (ceRNA) network play vital roles in gene expression regulation, but their involvement in sex determination and differentiation remains unclear in molluscs. In this study, a comprehensive transcriptomic analysis was performed to investigate ncRNAs and ceRNA network in female and male gonads of Crassostrea gigas. Differential expression analysis identified 3496 mRNAs, 582 lncRNAs, 184 miRNAs, and 644 circRNAs with sex-biased expression. Functional enrichment analyses highlighted key pathways such as the cell cycle, oocyte meiosis, energy metabolism, and lipid metabolism, underscoring their involvement in sex determination and differentiation. A ceRNA network was constructed involving 398 lncRNAs, 119 circRNAs, 140 miRNAs, and 720 mRNAs. Hub genes, such as Fem1c, Spef1, Dgkq, Ppp1ca, Nkd1, Morn3, Dpf2, and Gabarap were identified, with pronounced sex-biased expression and localization in specific gonadal cell types, as revealed by bulk and single-nucleus RNA-seq analysis. These genes are associated with critical processes, including follicle development, spermatogenesis, and hormonal regulation. Collectively, these findings provide novel insights into the ceRNA-mediated regulatory mechanisms underlying sex determination and differentiation in C. gigas, contributing to a deeper understanding of molluscan reproductive biology.

Comprehensive transcriptome analysis reveals lncRNA-mRNA interactions and immune response mechanisms in Mytilus coruscus upon Vibrio alginolyticus infection

Long non-coding RNAs (lncRNAs) play diverse biological roles within cells. Despite not encoding proteins, they are crucial in regulating gene expression, chromatin structure and function, cell differentiation and development, and the occurrence of diseases. Vibrio alginolyticus (V. alginolyticus) is a common bacterium found in marine environments that poses a threat to shellfish by infecting them through filtration feeding. Research has demonstrated the substantial involvement of lncRNAs in the immune response of shellfish. However, the specific mechanism by which lncRNAs participate in the immune regulatory process following infection of Mytilus coruscus (M. coruscus) with V. alginolyticus has not been investigated. Therefore, the transcription profiles of lncRNAs in M. coruscus hemocytes were investigated. A grand total of 48,246 lncRNAs were detected, with 2421 genes that exhibited and 717 lncRNAs that had differential expression. To gain a better understanding of the potential roles of the differentially expressed lncRNAs (DE-lncRNAs), GO and KEGG pathway analyses were performed on their target mRNAs, suggesting that lncRNAs have the ability to control gene expression levels and consequently influence immune-related pathways, hence regulating the immune response in M. coruscus. Additionally, a total of 138 lncRNA-mRNA pairs were identified through the calculation of co-expression relationships between DE-lncRNAs and immune-related DE-mRNAs. These findings provide new insights into the role of lncRNAs in the immune response of M. coruscus and offer important resources for further investigation into the role of lncRNAs in M. coruscus pathogen infection.

Insights from the single-cell level: lineage trajectory and somatic-germline interactions during spermatogenesis in dwarf surfclam Mulinia lateralis

BACKGROUND

Spermatogenesis is a complex process of cellular differentiation that commences with the division of spermatogonia stem cells, ultimately resulting in the production of functional spermatozoa. However, a substantial gap remains in our understanding of the molecular mechanisms and key driver genes that underpin this process, particularly in invertebrates. The dwarf surfclam (Mulinia lateralis) is considered an optimal bivalve model due to its relatively short generation time and ease of breeding in laboratory settings.

RESULTS

In this study, over 4,600 testicular cells from various samples were employed to identify single-cell heterogeneity on a more comprehensive scale. The four germ cell populations (spermatogonia, primary spermatocytes, secondary spermatocytes, and round spermatids/spermatozoa) and three somatic populations (follicle cell, hemocyte, and nerve cell) were characterized. The four types of germ cells exhibited disparate cell cycle statuses and an uninterrupted developmental trajectory, progressing from spermatogonia to spermatids/spermatozoa. Pseudotime analysis indicates that gene expression, translation, ATP metabolic process, and microtubule-based process are involved in the transition of germ cell types. Weighted gene coexpression network analysis (WGCNA) identified four modules corresponding to the four types of germ cells, as well as key transcription factors (e.g., MYC, SREBF1, SOXH) that may play a critical role in these cell types. Furthermore, our findings revealed that there is extensive bidirectional communication between the somatic cells and the germline cells, including the FGF and TGF-β signaling pathways, as well as other ligand-receptor pairs, such as NTN1-NEO1 and PLG-PLGRKT.

CONCLUSIONS

This study provides a comprehensive single-cell transcriptome landscape of the gonad, which will contribute to the understanding of germ cell fate transition during spermatogenesis, and the development of germ cell manipulation technologies in mollusks.

Low pH Means More Female Offspring: A Multigenerational Plasticity in the Sex Ratio of Marine Bivalves

Global changes can profoundly affect the sex determination and reproductive output of marine organisms, disrupting the population structure and ecosystems. High CO2driven low pH in the context of ocean acidification (OA) has been shown to severely affect various calcifiers, but less is known about the extent to which low pH influences sex determination and reproduction of marine organisms, particularly mollusks. This study is the first to report a biased sex ratio over multiple generations toward females, driven by exposure to high CO2-induced low pH environments, using the ecologically and economically important Portuguese oyster (Crassostrea angulata) as a model. This phenomenon, which we term pH-mediated sex determination (PSD), has no consequences for fecundity, gonadal development, or reproductive function in the offspring. Moreover, PSD persisted into a second year of reproduction and was inherited across multiple generations. Transcriptomic analysis indicates PSD is associated with the activation of the Wnt signaling pathway in females and inhibition of spermiogenesis-related functions in males. This work expands our understanding of environmental sex determination and highlights the possible impact of global changes on reproduction and population dynamics of mollusks and other marine organisms.

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.

DNA methylation correlates with transcriptional noise in response to elevated pCO2 in the eastern oyster (Crassostrea virginica)

Ocean acidification significantly affects marine calcifiers like oysters, warranting the study of molecular mechanisms like DNA methylation that contribute to adaptive plasticity in response to environmental change. However, a consensus has not been reached on the extent to which methylation modules gene expression, and in turn plasticity, in marine invertebrates. In this study, we investigated the impact of pCO2 on gene expression and DNA methylation in the eastern oyster, Crassostrea virginica. After a 30-day exposure to control (572 ppm) or elevated pCO2 (2827 ppm), whole-genome bisulfite sequencing (WGBS) and RNA-seq data were generated from adult female gonad tissue and male sperm samples. Although differentially methylated loci (DMLs) were identified in females (89) and males (2916), there were no differentially expressed genes and only one differentially expressed transcript in females. However, gene body methylation impacted other forms of gene activity in sperm, such as the maximum number of transcripts expressed per gene and changes in the predominant transcript expressed. Elevated pCO2 exposure increased gene expression variability (transcriptional noise) in males but decreased noise in females, suggesting a sex-specific role of methylation in gene expression regulation. Functional annotation of genes with changes in transcript-level expression or containing DMLs revealed several enriched biological processes potentially involved in elevated pCO2 response, including apoptotic pathways and signal transduction, as well as reproductive functions. Taken together, these results suggest that DNA methylation may regulate gene expression variability to maintain homeostasis in elevated pCO2 conditions and could play a key role in environmental resilience in marine invertebrates.

Different metazoan parasites, different transcriptomic responses, with new insights on parasitic castration by digenetic trematodes in the schistosome vector snail Biomphalaria glabrata

BACKGROUND

Gastropods of the genus Biomphalaria (Family Planorbidae) are exploited as vectors by Schistosoma mansoni, the most common causative agent of human intestinal schistosomiasis. Using improved genomic resources, overviews of how Biomphalaria responds to S. mansoni and other metazoan parasites can provide unique insights into the reproductive, immune, and other systems of invertebrate hosts, and their responses to parasite challenges.

RESULTS

Using Illumina-based RNA-Seq, we compared the responses of iM line B. glabrata at 2, 8, and 40 days post-infection (dpi) to single infections with S. mansoni, Echinostoma paraensei (both digenetic trematodes) or Daubaylia potomaca (a nematode parasite of planorbid snails). Responses were compared to unexposed time-matched control snails. We observed: (1) each parasite provoked a distinctive response with a predominance of down-regulated snail genes at all time points following exposure to either trematode, and of up-regulated genes at 8 and especially 40dpi following nematode exposure; (2) At 2 and 8dpi with either trematode, several snail genes associated with gametogenesis (particularly spermatogenesis) were down-regulated. Regarding the phenomenon of trematode-mediated parasitic castration in molluscs, we define for the first time a complement of host genes that are targeted, as early as 2dpi when trematode larvae are still small; (3) Differential gene expression of snails with trematode infection at 40dpi, when snails were shedding cercariae, was unexpectedly modest and revealed down-regulation of genes involved in the production of egg mass proteins and peptide processing; and (4) surprisingly, D. potomaca provoked up-regulation at 40dpi of many of the reproduction-related snail genes noted to be down-regulated at 2 and 8dpi following trematode infection. Happening at a time when B. glabrata began to succumb to D. potomaca, we hypothesize this response represents an unexpected form of fecundity compensation. We also document expression patterns for other Biomphalaria gene families, including fibrinogen domain-containing proteins (FReDs), C-type lectins, G-protein coupled receptors, biomphalysins, and protease and protease inhibitors.

CONCLUSIONS

Our study is relevant in identifying several genes involved in reproduction that are targeted by parasites in the vector snail B. glabrata and that might be amenable to manipulation to minimize their ability to serve as vectors of schistosomes.

The role of DNA methylation reprogramming during sex determination and sex reversal in the Pacific oyster Crassostrea gigas

DNA methylation is instrumental in vertebrate sex reversal. However, the mechanism of DNA methylation regulation regarding sex reversal in invertebrates is unclear. In this study, we used whole genome bisulfite sequencing (WGBS) to map single-base resolution methylation profiles of the Pacific oyster, including female-to-male (FMa-to-FMb) and male-to-female (MFa-to-MFb) sex reversal, as well as sex non-reversed males and females (MMa-to-MMb and FFa-to-FFb). The results showed that global DNA methylation levels increase during female-to-male sex reversals, with a particular increase in the proportion of high methylation levels (mCGs >0.75) and a decrease in the proportion of intermediate methylation levels (0.25 < mCGs <0.75). This increase in DNA methylation was mainly associated with the elevated expression of DNA methylase genes. Genome-wide methylation patterns of females were accurately remodeled to those of males after sex reversal, while the opposite was true for the male-to-female reversal. Those findings directly indicate that alterations in DNA methylation play a significant role in sex reversal in Pacific oysters. Comparative analysis of the DNA methylomes of pre- and post- sex reversal gonadal tissues (FMb-vs-FMa or MFb-vs-MFa) revealed that differentially methylated genes were mainly involved in the biological processes of sex determination or gonadal development. However critical genes such as Dmrt1, Foxl2 and Sox-like, which are involved in the putative sex determination pathway in Pacific oysters, showed almost an absence of methylation modifications, varying greatly from vertebrates. Additionally, comparative analysis of the DNA methylomes of sexual reversal and sex non-reversal (FMa-vs-FFa or MFa-vs-MMa) revealed that heat shock protein genes, such as Hsp68-like and Hsp70B, were important for the occurrence of sex reversal. These findings shed light on the epigenetic mechanisms underlying the maintenance of gonadal plasticity and the reversal of organ architecture in oysters.

Development of Nanopore sequencing-based full-length transcriptome database toward functional genome annotation of the Pacific oyster, Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is a widely cultivated shellfish in the world, while its transcriptome diversity remains less unexplored due to the limitation of short reads. In this study, we used Oxford Nanopore sequencing to develop the full-length transcriptome database of C. gigas. We identified 77,920 full-length transcripts from 21,523 genes, and uncovered 9668 alternative splicing events and 87,468 alternative polyadenylation sites. Notably, a total of 16,721 novel transcripts were annotated in this work. Furthermore, integrative analysis of 25 publicly available RNA-seq datasets revealed the transcriptome diversity involved in post-transcriptional regulation in C. gigas. We further developed a Drupal based webserver, Cgtdb, which can be used for transcriptome visualization, sequence alignment, and functional genome annotation analyses. This work provides valuable resources and a useful tool for integrative analysis of various transcriptome datasets in C. gigas, which will serve as an essential reference for functional annotation of the oyster genome.

Starvation-induced changes in sex ratio involve alterations in sex-related gene expression and methylation in Pacific oyster Crassostrea gigas

Aquatic animals are subject to varying degrees of starvation stress in their natural habitats due to food limitations. Consequently, starvation is a crucial environmental factor for sex determination in many species; however, limited research has been conducted on the effects of starvation on sex determination in shellfish. Here, four full sibling families of Pacific oysters were established and subjected to starvation stress. The results demonstrated that starvation caused the sex ratio (female to male) to change from 1:0.78 to 1:1.44 and resulted in a delay in gonadal development. Further studies revealed that the expression levels of DNA methylation-related genes Dnmt1 (DNA methyltransferase 1), Dnmt3a/b (DNA methyltransferase 3a/b) and Tet3 (tet methylcytosine dioxygenase 3) decreased under starvation stress. Conversely, the upregulation of Dmrt1 (doublesex and mab-3 related transcription factor 1), a gene typically associated with males, in females suggests that the increased proportion of males may be linked to starvation-induced high expression of this particular gene. In addition, the gene Dgkd (diacylglycerol kinase delta), which is involved in the regulation of second messenger protein kinase C, was differentially methylated between males and females, with the methylation level of this gene gradually increasing with male development, while the methylation level of this gene decreased under starvation stress. Correlation analysis of Dgkd methylation levels with expression levels showed a negative correlation between DNA methylation and gene expression. Finally dual fluorescence reporter experiments confirmed that DNA methylation suppressed Dgkd expression in vitro. In summary, the results suggest that starvation may regulate Dgkd gene expression through DNA methylation and thus affect Dmrt1 expression, thereby causing sex reversal in the Pacific oyster. The outcomes resolved how environmental factors are involved in sex reversal from an epigenetic perspective and provided a theoretical basis for sex control breeding in the Pacific oyster.

Early gonadal differentiation is associated with the antagonistic action of Foxl2 and Dmrt1l in the Pacific oyster

As the second largest phylum in the zoological kingdom next to arthropods, the mechanism of gonadal differentiation in mollusca is quite complex. Currently, although much has been carried out on gonadal differentiation in the Pacific oyster, there is still unknown information that needs to be further explored. Here, analysis of the Foxl2 and Dmrt1l expression in samples at different development periods of male and female gonads as well as in annual gonad samples revealed that Log₁₀ (Foxl2/Dmrt1l) values were an effective method for sex identification in oysters. In differentiated gonadal tissue, Log₁₀ (Foxl2/Dmrt1l) values greater than 2 were females and less than 1 for males. Subsequent sequential sampling of the same individuals verified that Log₁₀ (Foxl2/Dmrt1l) values greater than 2 for resting gonads would develop as females and less than 1 would develop as males in the future. Relative expression analysis of Foxl2 and Dmrt1l in the annual samples revealed a negative correlation between Log₁₀ (Foxl2) and Log₁₀ (Dmrt1l). Double fluorescence reporter validation results showed that DMRT1L protein was able to bind the Foxl2 promoter and repress its activity with a weak dosage effect. Antagonism between Dmrt1l and Foxl2 is therefore not restricted to vertebrates, and the competing regulatory networks are of great significance in the maintenance of gonadal sex in oysters after sexual differentiation. This study provides novel ideas and insights into the study of early gonadal differentiation in the adult oyster.

Multi-tissue RNA-Seq Analysis and Long-read-based Genome Assembly Reveal Complex Sex-specific Gene Regulation and Molecular Evolution in the Manila Clam

The molecular factors and gene regulation involved in sex determination and gonad differentiation in bivalve molluscs are unknown. It has been suggested that doubly uniparental inheritance (DUI) of mitochondria may be involved in these processes in species such as the ubiquitous and commercially relevant Manila clam, Ruditapes philippinarum. We present the first long-read-based de novo genome assembly of a Manila clam, and a RNA-Seq multi-tissue analysis of 15 females and 15 males. The highly contiguous genome assembly was used as reference to investigate gene expression, alternative splicing, sequence evolution, tissue-specific co-expression networks, and sexual contrasting SNPs. Differential expression (DE) and differential splicing (DS) analyses revealed sex-specific transcriptional regulation in gonads, but not in somatic tissues. Co-expression networks revealed complex gene regulation in gonads, and genes in gonad-associated modules showed high tissue specificity. However, male gonad-associated modules showed contrasting patterns of sequence evolution and tissue specificity. One gene set was related to the structural organization of male gametes and presented slow sequence evolution but high pleiotropy, whereas another gene set was enriched in reproduction-related processes and characterized by fast sequence evolution and tissue specificity. Sexual contrasting SNPs were found in genes overrepresented in mitochondrial-related functions, providing new candidates for investigating the relationship between mitochondria and sex in DUI species. Together, these results increase our understanding of the role of DE, DS, and sequence evolution of sex-specific genes in an understudied taxon. We also provide resourceful genomic data for studies regarding sex diagnosis and breeding in bivalves.

Dynamic transcriptome analysis reveals the gene network of gonadal development from the early history life stages in dwarf surfclam Mulinia lateralis

BACKGROUND

Gonadal development is driven by a complex genetic cascade in vertebrates. However, related information remains limited in molluscs owing to the long generation time and the difficulty in maintaining whole life cycle in the lab. The dwarf surfclam Mulinia lateralis is considered an ideal bivalve model due to the short generation time and ease to breed in the lab.

RESULTS

To gain a comprehensive understanding of gonadal development in M. lateralis, we conducted a combined morphological and molecular analysis on the gonads of 30 to 60 dpf. Morphological analysis showed that gonad formation and sex differentiation occur at 35 and 40-45 dpf, respectively; then the gonads go through gametogenic cycle. Gene co-expression network analysis on 40 transcriptomes of 35-60 dpf gonads identifies seven gonadal development-related modules, including two gonad-forming modules (M6, M7), three sex-specific modules (M14, M12, M11), and two sexually shared modules (M15, M13). The modules participate in different biological processes, such as cell communication, glycan biosynthesis, cell cycle, and ribosome biogenesis. Several hub transcription factors including SOX2, FOXZ, HSFY, FOXL2 and HES1 are identified. The expression of top hub genes from sex-specific modules suggests molecular sex differentiation (35 dpf) occurs earlier than morphological sex differentiation (40-45 dpf).

CONCLUSION

This study provides a deep insight into the molecular basis of gonad formation, sex differentiation and gametogenesis in M. lateralis, which will contribute to a comprehensive understanding of the reproductive regulation network in molluscs.

Differential DNA methylation in Pacific oyster reproductive tissue in response to ocean acidification

BACKGROUND

There is a need to investigate mechanisms of phenotypic plasticity in marine invertebrates as negative effects of climate change, like ocean acidification, are experienced by coastal ecosystems. Environmentally-induced changes to the methylome may regulate gene expression, but methylome responses can be species- and tissue-specific. Tissue-specificity has implications for gonad tissue, as gonad-specific methylation patterns may be inherited by offspring. We used the Pacific oyster (Crassostrea gigas) - a model for understanding pH impacts on bivalve molecular physiology due to its genomic resources and importance in global aquaculture- to assess how low pH could impact the gonad methylome. Oysters were exposed to either low pH (7.31 ± 0.02) or ambient pH (7.82 ± 0.02) conditions for 7 weeks. Whole genome bisulfite sequencing was used to identify methylated regions in female oyster gonad samples. C- > T single nucleotide polymorphisms were identified and removed to ensure accurate methylation characterization.

RESULTS

Analysis of gonad methylomes revealed a total of 1284 differentially methylated loci (DML) found primarily in genes, with several genes containing multiple DML. Gene ontologies for genes containing DML were involved in development and stress response, suggesting methylation may promote gonad growth homeostasis in low pH conditions. Additionally, several of these genes were associated with cytoskeletal structure regulation, metabolism, and protein ubiquitination - commonly-observed responses to ocean acidification. Comparison of these DML with other Crassostrea spp. exposed to ocean acidification demonstrates that similar pathways, but not identical genes, are impacted by methylation.

CONCLUSIONS

Our work suggests DNA methylation may have a regulatory role in gonad and larval development, which would shape adult and offspring responses to low pH stress. Combined with existing molluscan methylome research, our work further supports the need for tissue- and species-specific studies to understand the potential regulatory role of DNA methylation.

Screening of Genes Related to Sex Determination and Differentiation in Mandarin Fish (Siniperca chuatsi)

Mandarin fish has an XX/XY sex-determination system. The female mandarin fish is typically larger than the male. Sex identification and the discovery of genes related to sex determination in mandarin fish have important theoretical significance in the elucidation of the regulation and evolutionary mechanism of animal reproductive development. In this study, the chromosome-level genome of a female mandarin fish was assembled, and we found that LG24 of the genome was an X chromosome. A total of 61 genes on the X chromosome showed sex-biased expression. Only six gonadal genes (LG24G00426, LG24G003280, LG24G003300, LG24G003730, LG24G004200, and LG24G004770) were expressed in the testes, and the expression of the other gene LG24G003870 isoform 1 in the ovaries was significantly higher than that in the testes (p < 0.01). Five (except LG24G003280 and LG24G003300) of the seven aforementioned genes were expressed at the embryonic development stage, suggesting their involvement in early sex determination. The expression of LG24G004770 (encoding HS6ST 3-B-like) was also significantly higher in female muscles than in male muscles (p < 0.01), indicating other functions related to female growth. ZP3 encoded by LG24G003870 isoform 1 increased the C-terminal transmembrane domain, compared with that encoded by other fish zp3 isoforms, indicating their different functions in sex determination or differentiation. This study provides a foundation for the identification of sex-determining genes in mandarin fish.

Full-Length Transcriptome Sequences Provide Insight Into Hermaphroditism of Freshwater Pearl Mussel Hyriopsis schlegelii

The freshwater mussel Hyriopsis schlegelii is a cultured bivalve in China, and the quality of the pearls produced is affected by the type of gonads. However, because of the lack of a published genome and the complexity of sex determination, research on sex reversal and development of this species is limited. In this study, Illumina RNA-seq and PacBio Isoform Sequencing (Iso-Seq) were combined to analyze the gonads of H. schlegelii. A total of 201,481 high-quality transcripts were generated. The study identified 7,922 differentially expressed genes in three comparison group (females versus males, hermaphrodites versus females, and hermaphrodites versus males). Twenty-four genes were identified as potential sex-related genes, including sox9 and wnt4 involved in sex determination, and vtg, cyp17a1 and 17β-hsd2 involved in gonadal development. We also speculated a possible pathways for the formation of hermaphroditism in H. schlegelii. The data provide a clear view of the transcriptome for H. schlegelii gonads and will be valuable in elucidating the mechanisms of gonad development.

Transcription Analysis for Core Networks of lncRNAs–mRNAs: Implication for Potential Role in Sterility of Crassostrea gigas

Simple Summary

This study reveals the expression profiles of lncRNA in the gonads of the Pacific oyster Crassostrea gigas. The potential function of lncRNAs was predicted in the case of antisense and cis-regulatory mechanisms based on their physical positions and their coexpression relationships in the case of trans regulation. Sterility-related DEGs and DELs were chosen for subsequent analysis, demonstrating that trans-regulatory lncRNAs might play a vital role in the gametogenesis of C. gigas. We constructed core networks of lncRNAs–mRNAs for triploid sterile females and hermaphrodites based on pathway results, in which 28 lncRNAs and their 54 trans-regulatory genes were detected. Among 28 sterility-specific lncRNAs, MSTRG.79882.3 and MSTRG.79882.4 for triploid sterile females and MSTRG.33704.1, MSTRG.63844.1, and MSTRG.5675.1 for hermaphrodites play the most significant role.

Abstract

Long noncoding RNA (lncRNA), a type of non-protein-coding transcript, is emerging as a crucial regulator of gene expression. However, few roles of lncRNA in the reproductive process of the Pacific oyster (Crassostrea gigas) have been defined, especially in the regulatory mechanism of sterile triploids gametogenesis. To uncover the potential role of lncRNA, the gonads of diploids, sterile triploids, and partially sterile triploids underwent RNA sequencing. A total of 9618 reliable lncRNAs were identified. The target relationship between lncRNA and mRNA was predicted based on cis, trans, and antisense regulation with bioinformatic software. We chose differentially expressed lncRNAs and mRNAs when sterile triploids were compared to partially sterile triploids and diploids for subsequent functional enrichment analysis. Findings revealed that trans-regulatory lncRNAs might play a significant role in the gametogenesis of C. gigas. Combining pathway results, we constructed core networks of lncRNAs–mRNAs for triploid sterile females and hermaphrodites. Fifty-four genes related to cell division, germline-cell maintenance, and glycogen metabolism were found to be associated with sterility. A total of 28 candidate lncRNAs were predicted to trans-regulate these genes. We speculated that MSTRG.79882.3 and MSTRG.79882.4 for triploid sterile females and MSTRG.33704.1, MSTRG.63844.1, and MSTRG.5675.1 for hermaphrodites were highly important as they were predicted to regulate more sterility-specific genes than others. Our work collectively identified sterility-related lncRNAs and implicated the potential mechanism of lncRNA-mediated regulation in the gametogenesis of sterile triploid oysters.