Spliced leader dependent trans-splicing (SLTS) has been described as an important RNA regulatory process that occurs in different organisms, including the trematode Schistosoma mansoni. We identified more than seven thousand putative SLTS sites in the parasite, comprising genes with a wide spectrum of functional classes, which underlines the SLTS as a ubiquitous mechanism in the parasite. Also, SLTS gene expression levels span several orders of magnitude, showing that SLTS frequency is not determined by the expression level of the target gene, but by the presence of particular gene features facilitating or hindering the trans-splicing mechanism. Our in-depth investigation of SLTS events demonstrates widespread alternative trans-splicing (ATS) acceptor sites occurring in different regions along the entire gene body, highlighting another important role of SLTS generating alternative RNA isoforms in the parasite, besides the polycistron resolution. Particularly for introns where SLTS directly competes for the same acceptor substrate with cis-splicing, we identified for the first time additional and important features that might determine the type of splicing. Our study substantially extends the current knowledge of RNA processing by SLTS in S. mansoni, and provide basis for future studies on the trans-splicing mechanism in other eukaryotes.

Riboswitches are RNA sensors that affect post-transcriptional processes through their ability to bind to small molecules. Thiamine pyrophosphate (TPP) riboswitch class is the most widespread riboswitch occurring in all three domains of life. Even though it controls different genes involved in the synthesis or transport of thiamine and its phosphorylated derivatives in bacteria, archaea, fungi, and plants, the TPP aptamer has a conserved structure. In this study, we aimed at understanding differences in the structural dynamics of TPP riboswitches from Escherichia coli and Arabidopsis thaliana, based on their crystallographic structures (TPPswec and TPPswat, respectively) and dynamics in aqueous solution, both in apo and holo states. A combination of Molecular Dynamics Simulations and Network Analysis empowered to find out slight differences in the dynamical behavior of TPP riboswitches, although relevant for their dynamics in bacteria and plants species. Our results suggest that distinct interactions in the microenvironment surrounding nucleotide U36 of TPPswec (and U35 in TPPswat) are related to different responses to TPP. The network analysis showed that minor structural differences in the aptamer enable enhanced intramolecular communication in the presence of TPP in TPPswec, but not in TPPswat. TPP riboswitches of plants present subtler and slower regulation mechanisms than bacteria do.

RNA molecules are essential players in many fundamental biological processes. Prokaryotes and eukaryotes have distinct RNA classes with specific structural features and functional roles. Computational prediction of protein structures is a research field in which high confidence three-dimensional protein models can be proposed based on the sequence alignment between target and templates. However, to date, only a few approaches have been developed for the computational prediction of RNA structures. Similar to proteins, RNA structures may be altered due to the interaction with various ligands, including proteins, other RNAs, and metabolites. A riboswitch is a molecular mechanism, found in the three kingdoms of life, in which the RNA structure is modified by the binding of a metabolite. It can regulate multiple gene expression mechanisms, such as transcription, translation initiation, and mRNA splicing and processing. Due to their nature, these entities also act on the regulation of gene expression and detection of small metabolites and have the potential to helping in the discovery of new classes of antimicrobial agents. In this review, we describe software and web servers currently available for riboswitch aptamer identification and secondary and tertiary structure prediction, including applications.

Background: Interaction between malignant cells and immune cells that reside within the tumor microenvironment (TME) modulate different aspects of tumor development and progression. Recent works showed the importance of miRNA-containing extracellular vesicles in this crosstalk. Methods: Interested in understanding the interplay between melanoma and immune-related TME cells, we characterized the TCGA's metastatic melanoma samples according to their tumor microenvironment profiles, HLA-I neoepitopes, transcriptome profile and classified them into three groups. Moreover, we combined our results with melanoma single-cell gene expression and public miRNA data to better characterize the regulatory network of circulating miRNAs and their targets related to immune evasion and microenvironment response. Results: The group associated with a worse prognosis showed phenotypic characteristics that favor immune evasion, including a strong signature of suppressor cells and less stable neoantigen:HLA-I complexes. Conversely, the group with better prognosis was marked by enrichment in lymphocyte and MHC signatures. By analyzing publicly available melanoma single-cell RNA and microvesicle microRNAs sequencing data we identified circulating micro-RNAs potentially involved in the crosstalk between tumor and TME cells. Candidate miRNA/target gene pairs with previously reported roles in tumor progression and immune escape mechanisms were further investigated and demonstrated to impact patient's overall survival not only in melanoma but across different tumor types. Conclusion: Our results underscore the impact of tumor-microenvironment interactions on disease outcomes and reveal potential non-invasive biomarkers of prognosis and treatment response.

Lung cancer is one of the most frequent types of cancer worldwide. Most patients are diagnosed at advanced stage and thus have poor prognosis. Smoking is a risk factor for lung cancer, however most smokers do not develop lung cancer while 20% of women with lung adenocarcinoma are non-smokers. Therefore, it is possible that these two groups present differences besides the smoking status, including differences in their gene expression signature. The altered expression patterns of non-coding RNAs in complex diseases make them potential biomarkers for diagnosis and treatment. We analyzed data from differentially and constitutively expressed PIWI-interacting RNAs and small nucleolar RNAs from publicly available small RNA high-throughput sequencing data in search of an expression pattern of non-coding RNA that could differentiate these two groups. Here, we report two sets of differentially expressed small non-coding RNAs identified in normal and tumoral tissues of women with lung adenocarcinoma, that discriminate between smokers and non-smokers. Our findings may offer new insights on metabolic alterations caused by tobacco and may be used for early diagnosis of lung cancer.

Oral mucositis is an acute toxicity that occurs in patients submitted to chemoradiotherapy to treat head and neck squamous cell carcinoma. In this study, we evaluated differences in gene expression in the keratinocytes of the oral mucosa of patients treated with photobiomodulation therapy and tried to associate the molecular mechanisms with clinical findings. From June 2009 to December 2010, 27 patients were included in a randomized double-blind pilot study. Buccal smears from 13 patients were obtained at days 1 and 10 of chemoradiotherapy, and overall gene expression of samples from both dates were analyzed by complementary DNA (cDNA) microarray. In addition, samples from other 14 patients were also collected at D1 and D10 of chemoradiotherapy for subsequent validation of cDNA microarray findings by qPCR. The expression array analysis identified 105 upregulated and 60 downregulated genes in our post-treatment samples when compared with controls. Among the upregulated genes with the highest fold change, it was interesting to observe the presence of genes related to keratinocyte differentiation. Among downregulated genes were observed genes related to cytotoxicity and immune response. The results indicate that genes known to be induced during differentiation of human epidermal keratinocytes were upregulated while genes associated with cytotoxicity and immune response were downregulated in the laser group. These results support previous clinical findings indicating that the lower incidence of oral mucositis associated with photobiomodulation therapy might be correlated to the activation of genes involved in keratinocyte differentiation.

Transcriptional profiling of hematopoietic cell subpopulations has helped characterize the developmental stages of the hematopoietic system and the molecular bases of malignant and non-malignant blood diseases for the past three decades. Previously, only the genes targeted by expression microarrays could be profiled genome wide. High-throughput RNA sequencing (RNA-seq), however, encompasses a broader repertoire of RNA molecules, without restriction to previously annotated genes. We analysed the BLUEPRINT consortium RNA-seq data for mature hematopoietic cell types. The data comprised 90 total RNA-seq samples, each composed of one of 27 cell types, and 32 small RNA-seq samples, each composed of one of 11 cell types. We estimated gene and isoform expression levels for each cell type using existing annotations from Ensembl. We then used guided transcriptome assembly to discover unannotated transcripts. We identified hundreds of novel non-coding RNA genes and showed that the majority have cell type dependent expression. We also characterized the expression of circular RNAs and found that these are also cell type specific. These analyses refine the active transcriptional landscape of mature hematopoietic cells, highlight abundant genes and transcriptional isoforms for each blood cell type, and provide a valuable resource for researchers of hematological development and diseases. Finally, we made the data accessible via a web-based interface: https://blueprint.haem.cam.ac.uk/bloodatlas/.

Squamous cell carcinoma (SCC) and adenocarcinoma (ADC) are the most common histological types of cervical cancer (CC). The worse prognosis of ADC cases highlights the need for better molecular characterization regarding differences between these CC types. RNA-Seq analysis of seven SCC and three ADC human papillomavirus 16-positive samples and the comparison with public data from non-tumoral human papillomavirus-negative cervical tissue samples revealed pathways exclusive to each histological type, such as the epithelial maintenance in SCC and the maturity-onset diabetes of the young (MODY) pathway in ADC. The transcriptional regulatory network analysis of cervical SCC samples unveiled a set of six transcription factor (TF) genes with the potential to positively regulate long non-coding RNA genes DSG1-AS1, CALML3-AS1, IGFL2-AS1, and TINCR. Additional analysis revealed a set of MODY TFs regulated in the sequence predicted to be repressed by miR-96-5p or miR-28-3p in ADC. These microRNAs were previously described to target LINC02381, which was predicted to be positively regulated by two MODY TFs upregulated in cervical ADC. Therefore, we hypothesize LINC02381 might act by decreasing the levels of miR-96-5p and miR-28-3p, promoting the MODY activation in cervical ADC. The novel TF networks here described should be explored for the development of more efficient diagnostic tools.