Comprehensive Schistosoma mansoni Hierarchical Transcriptome Assembly Points to Novel lncRNAs Associated with Sexual...

Background/Objectives: Schistosomiasis is a neglected tropical disease affecting >200 million people worldwide. Praziquantel is the sole recommended drug against Schistosoma mansoni; however, it lacks activity against juvenile forms and cannot prevent reinfection. Thus, there is an urgent need to identify novel therapeutic targets. Long noncoding RNAs (lncRNAs) are known to regulate various biological processes in S. mansoni, including parasite pairing and fertility; therefore, screening for novel lncRNAs could reveal new potential targets. Methods: We compiled all publicly available RNA-seq data from the Sequence Read Archive (SRA) and performed a hierarchical transcriptome assembly using the multi-sample assembler Ryūtō, combined with version 10 of the S. mansoni genome. We applied HOMER for peak-calling and identification of histone marks and used weighted gene co-expression network analysis (WGCNA) to infer putative functions of lncRNAs in sexual dimorphism. Results: Using a robust pipeline, we identified 10,170 novel lncRNA genes comprising 16,990 novel lncRNA transcripts, including 8783 intergenic, 7918 antisense, and 289 intronic lncRNA transcripts. Most (78.7%) have histone regulatory marks (H3K4me3, H3K27me3, H3K27ac, or H4K20me1) near their transcription start sites, indicating potential expression regulation. Comparing male and female samples, we identified 1991 differentially expressed genes (FDR < 5%, |log₂FC| ≥ 1.5), including 296 known lncRNAs and 339 novel lncRNAs. WGCNA identified hub lncRNAs within co-expression modules, and Gene Ontology enrichment analyses (FDR ≤ 5%) suggest that these lncRNAs are involved in cell differentiation and morphogenesis pathways. Conclusions: We provide a comprehensive catalog of S. mansoni lncRNAs. These findings offer opportunities to discover potential new therapeutic targets, advancing the future development of anti-schistosome therapies.