Dissecting transposable elements and endogenous retroviruses upregulation by HDAC inhibitors in leiomyosarcoma cells: Implications for the interferon response.

Transposable elements (TEs) are of interest as immunomodulators for cancer therapies. TEs can fold into dsRNAs that trigger the interferon response. Here, we investigated the effect of different HDAC inhibitors (HDACIs) on the expression of TEs in leiomyosarcoma cells. Our data show that endogenous retroviruses (ERVs), especially ERV1 elements, are upregulated after treatment with HDAC1/2/3-specific inhibitors. Surprisingly, the interferon response was not activated. We observed an increase in A-to-I editing of upregulated ERV1. This could have an impact on the stability of dsRNAs and the activation of the interferon response. We also found that H3K27ac levels are increased in the LTR12 subfamilies, which could be regulatory elements controlling the expression of proapoptotic genes such as TNFRSF10B. In summary, we provide a detailed characterization of TEs modulation in response to HDACIs and suggest the use of HDACIs in combination with ADAR inhibitors to induce cell death and support immunotherapy in cancer.

Transcriptomic Module Discovery of Diarrhea-Predominant Irritable Bowel Syndrome: A Causal Network Inference Approach.

Irritable bowel syndrome with diarrhea (IBS-D) is the most prevalent subtype of IBS, characterized by chronic gastrointestinal symptoms in the absence of identifiable pathological findings. This study aims to investigate the molecular mechanisms underlying IBS-D using transcriptomic data. By employing causal network inference methods, we identify key transcriptomic modules associated with IBS-D. Utilizing data from public databases and applying advanced computational techniques, we uncover potential biomarkers and therapeutic targets. Our analysis reveals significant molecular alterations that affect cellular functions, offering new insights into the complex pathophysiology of IBS-D. These findings enhance our understanding of the disease and may foster the development of more effective treatments.

TEspeX: consensus-specific quantification of transposable element expression preventing biases from exonized fragments.

SUMMARY: Transposable elements (TEs) play key roles in crucial biological pathways. Therefore, several tools enabling the quantification of their expression were recently developed. However, many of the existing tools lack the capability to distinguish between the transcription of autonomously expressed TEs and TE fragments embedded in canonical coding/non-coding non-TE transcripts. Consequently, an apparent change in the expression of a given TE may simply reflect the variation in the expression of the transcripts containing TE-derived sequences. To overcome this issue, we have developed TEspeX, a pipeline for the quantification of TE expression at the consensus level. TEspeX uses Illumina RNA-seq short reads to quantify TE expression avoiding counting reads deriving from inactive TE fragments embedded in canonical transcripts. AVAILABILITY AND IMPLEMENTATION: The tool is implemented in python3, distributed under the GNU General Public License (GPL) and available on Github at https://github.com/fansalon/TEspeX (Zenodo URL: https://doi.org/10.5281/zenodo.6800331). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Transcription of endogenous retroviruses in senescent cells contributes to the accumulation of double-stranded RNAs that trigger an anti-viral response that reinforces senescence.

An important epigenetic switch marks the onset and maintenance of senescence. This allows transcription of the genetic programs that arrest the cell cycle and alter the microenvironment. Transcription of endogenous retroviruses (ERVs) is also a consequence of this epigenetic switch. In this manuscript, we have identified a group of ERVs that are epigenetically silenced in proliferating cells but are upregulated during replicative senescence or during various forms of oncogene-induced senescence, by RAS and Akt, or after HDAC4 depletion. In a HDAC4 model of senescence, removal of the repressive histone mark H3K27me3 is the plausible mechanism that allows the transcription of intergenic ERVs during senescence. We have shown that ERVs contribute to the accumulation of dsRNAs in senescence, which can initiate the antiviral response via the IFIH1-MAVS signaling pathway and thus contribute to the maintenance of senescence. This pathway, and MAVS in particular, plays an active role in shaping the microenvironment and maintaining growth arrest, two essential features of the senescence program.

NRF2 interacts with distal enhancer and inhibits nitric oxide synthase 2 expression in KRAS-driven pancreatic cancer cells.

Nitric oxide is a pleiotropic free radical produced by three nitric oxide synthases (NOS1-3), of which inducible NOS2 is involved in tumor initiation and progression. In this study, RNA-seq, ChIP-seq and qRT-PCR experiments combined with bioinformatic analyses showed that NRF2 is a repressor of NOS2 gene by maintaining a distal enhancer located 22 kb downstream of TSS in an inactive state. Deletion of NRF2 leads to activation of the enhancer, which exerts a pioneering function before it is fully activated. Specifically, NRF2 controls the expression of NOS2 in response to intracellular oxidative stress and extracellular oxygen pressure. We found that abrogation of NOS2 expression by siRNAs partially reduced the ability of WT Panc-1 cells to form 3D spheroids, but strongly reduced the formation of 3D spheroids by NRF2-depleted Panc-1 cells. Mechanistically, this effect correlates with the finding that NOS2 and nitric oxide stimulate epithelial-to-mesenchymal transition in NRF2-depleted Panc-1 and MIA PaCa-2 cells. We also found that knockdown of NOS2 leads to blockade of 3D matrigel invasion of NRF2-depleted PDAC cells, demonstrating that a short-circuit in the reciprocal regulation of NOS2 and NRF2 attenuates the malignancy of PDAC cells. In summary, we show for the first time that: (i) NRF2 is a suppressor of NOS2 in pancreatic cancer cells; (ii) NRF2 binds to and inactivates an enhancer located 22 kb downstream of TSS of the NOS2 gene; (iii) activation of NOS2 requires suppression of NRF2; (iv) NOS2 is required for NRF2-depleted Panc-1 cells to maintain their malignancy and invasiveness.

The DNA-repair protein APE1 participates with hnRNPA2B1 to motif-enriched and prognostic miRNA secretion.

The base excision repair (BER) Apurinic/apyrimidinic endonuclease 1 (APE1) enzyme is endowed with several non-repair activities including miRNAs processing. APE1 is overexpressed in many cancers but its causal role in the tumorigenic processes is largely unknown. We recently described that APE1 can be actively secreted by mammalian cells through exosomes. However, APE1 role in EVs or exosomes is still unknown, especially regarding a putative regulatory function on vesicular small non-coding RNAs. Through dedicated transcriptomic analysis on cellular and vesicular small RNAs of different APE1-depleted cancer cell lines, we found that miRNAs loading into EVs is a regulated process, dependent on APE1, distinctly conveying RNA subsets into vesicles. We identified APE1-dependent secreted miRNAs characterized by enriched sequence motifs and possible binding sites for APE1. In 33 out of 34 APE1-dependent-miRNA precursors, we surprisingly found EXO-motifs and proved that APE1 cooperates with hnRNPA2B1 for the EV-sorting of a subset of miRNAs, including miR-1246, through direct binding to GGAG stretches. Using TCGA-datasets, we showed that these miRNAs identify a signature with high prognostic significance in cancer. In summary, we provided evidence that the ubiquitous DNA-repair enzyme APE1 is part of the EV protein cargo with a novel post-transcriptional role for this ubiquitous DNA-repair enzyme that could explain its role in cancer progression. These findings could open new translational perspectives in cancer biology.

Unveiling the Impact of Gene Presence/Absence Variation in Driving Inter-Individual Sequence Diversity within the CRP-I Gene Family in Mytilus spp.

Mussels (Mytilus spp.) tolerate infections much better than other species living in the same marine coastal environment thanks to a highly efficient innate immune system, which exploits a remarkable diversification of effector molecules involved in mucosal and humoral responses. Among these, antimicrobial peptides (AMPs) are subjected to massive gene presence/absence variation (PAV), endowing each individual with a potentially unique repertoire of defense molecules. The unavailability of a chromosome-scale assembly has so far prevented a comprehensive evaluation of the genomic arrangement of AMP-encoding loci, preventing an accurate ascertainment of the orthology/paralogy relationships among sequence variants. Here, we characterized the CRP-I gene cluster in the blue mussel Mytilus edulis, which includes about 50 paralogous genes and pseudogenes, mostly packed in a small genomic region within chromosome 5. We further reported the occurrence of widespread PAV within this family in the Mytilus species complex and provided evidence that CRP-I peptides likely adopt a knottin fold. We functionally characterized the synthetic peptide sCRP-I H1, assessing the presence of biological activities consistent with other knottins, revealing that mussel CRP-I peptides are unlikely to act as antimicrobial agents or protease inhibitors, even though they may be used as defense molecules against infections from eukaryotic parasites.

Meta-Analysis Suggests That Intron Retention Can Affect Quantification of Transposable Elements from RNA-Seq Data.

Transposable elements (TEs), also known as "jumping genes", are repetitive sequences with the capability of changing their location within the genome. They are key players in many different biological processes in health and disease. Therefore, a reliable quantification of their expression as transcriptional units is crucial to distinguish between their independent expression and the transcription of their sequences as part of canonical transcripts. TEs quantification faces difficulties of different types, the most important one being low reads mappability due to their repetitive nature preventing an unambiguous mapping of reads originating from their sequences. A large fraction of TEs fragments localizes within introns, which led to the hypothesis that intron retention (IR) can be an additional source of bias, potentially affecting accurate TEs quantification. IR occurs when introns, normally removed from the mature transcript by the splicing machinery, are maintained in mature transcripts. IR is a widespread mechanism affecting many different genes with cell type-specific patterns. We hypothesized that, in an RNA-seq experiment, reads derived from retained introns can introduce a bias in the detection of overlapping, independent TEs RNA expression. In this study we performed meta-analysis using public RNA-seq data from lymphoblastoid cell lines and show that IR can impact TEs quantification using established tools with default parameters. Reads mapped on intronic TEs were indeed associated to the expression of TEs and influence their correct quantification as independent transcriptional units. We confirmed these results using additional independent datasets, demonstrating that this bias does not appear in samples where IR is not present and that differential TEs expression does not impact on IR quantification. We concluded that IR causes the over-quantification of intronic TEs and differential IR might be confused with differential TEs expression. Our results should be taken into account for a correct quantification of TEs expression from RNA-seq data, especially in samples in which IR is abundant.

Transposons Acting as Competitive Endogenous RNAs: In-Silico Evidence from Datasets Characterised by L1 Overexpression.

LINE L1 are transposable elements that can replicate within the genome by passing through RNA intermediates. The vast majority of these element copies in the human genome are inactive and just between 100 and 150 copies are still able to mobilize. During evolution, they could have been positively selected for beneficial cellular functions. Nonetheless, L1 deregulation can be detrimental to the cell, causing diseases such as cancer. The activity of miRNAs represents a fundamental mechanism for controlling transcript levels in somatic cells. These are a class of small non-coding RNAs that cause degradation or translational inhibition of their target transcripts. Beyond this, competitive endogenous RNAs (ceRNAs), mostly made by circular and non-coding RNAs, have been seen to compete for the binding of the same set of miRNAs targeting protein coding genes. In this study, we have investigated whether autonomously transcribed L1s may act as ceRNAs by analyzing public dataset in-silico. We observed that genes sharing miRNA target sites with L1 have a tendency to be upregulated when L1 are overexpressed, suggesting the possibility that L1 might act as ceRNAs. This finding will help in the interpretation of transcriptomic responses in contexts characterized by the specific activation of transposons.