GAIA: G-quadruplexes in alive creature database.

G-quadruplexes (G4) are 3D structures that are found in both DNA and RNA. Interest in this structure has grown over the past few years due to both its implication in diverse biological mechanisms and its potential use as a therapeutic target, to name two examples. G4s in humans have been widely studied; however, the level of their study in other species remains relatively minimal. That said, progress in this field has resulted in the prediction of G4s structures in various species, ranging from bacteria to eukaryotes. These predictions were analysed in a previous study which revealed that G4s are present in all living kingdoms. To date, eleven different databases have grouped the various G4s depending on either their structures, on the proteins that might bind them, or on their location in the various genomes. However, none of these databases contains information on their location in the transcriptome of many of the implicated species. The GAIA database was designed so as to make this data available online in a user-friendly manner. Through its web interface, users can query GAIA to filter G4s, which, we hope, will help the research in this field. GAIA is available at: https://gaia.cobius.usherbrooke.ca.

Reovirus µ2 protein modulates host cell alternative splicing by reducing protein levels of U5 snRNP core components.

Mammalian orthoreovirus (MRV) is a double-stranded RNA virus from the Reoviridae family presenting a promising activity as an oncolytic virus. Recent studies have underlined MRV's ability to alter cellular alternative splicing (AS) during infection, with a limited understanding of the mechanisms at play. In this study, we investigated how MRV modulates AS. Using a combination of cell biology and reverse genetics experiments, we demonstrated that the M1 gene segment, encoding the µ2 protein, is the primary determinant of MRV's ability to alter AS, and that the amino acid at position 208 in µ2 is critical to induce these changes. Moreover, we showed that the expression of µ2 by itself is sufficient to trigger AS changes, and its ability to enter the nucleus is not required for all these changes. Moreover, we identified core components of the U5 snRNP (i.e. EFTUD2, PRPF8, and SNRNP200) as interactors of µ2 that are required for MRV modulation of AS. Finally, these U5 snRNP components are reduced at the protein level by both MRV infection and µ2 expression. Our findings identify the reduction of U5 snRNP components levels as a new mechanism by which viruses alter cellular AS.

Guanine Nucleotide-Binding Protein-Like 1 (GNL1) binds RNA G-quadruplex structures in genes associated with Parkinson's disease.

RNAs are highly regulated at the post-transcriptional level in neurodegenerative diseases and just a few mutations can significantly affect the fate of neuronal cells. To date, the impact of G-quadruplex (G4) regulation in neurodegenerative diseases like Parkinson's disease (PD) has not been analysed. In this study, in silico potential G4s located in deregulated genes related to the nervous system were initially identified and were found to be significantly enriched. Several G4 sequences found in the 5' untranslated regions (5'UTR) of mRNAs associated with Parkinson's disease were demonstrated to in fact fold in vitro by biochemical assays. Subcloning of the full-length 5'UTRs of these candidates upstream of a luciferase reporter system led to the demonstration that the G4s of both Parkin RBR E3 Ubiquitin Protein Ligase (PRKN) and Vacuolar Protein Sorting-Associated Protein 35 (VPS35) significantly repressed the translation of both genes in SH-SY5Y cells. Subsequently, a strategy of using label-free RNA affinity purification assays with either of these two G4 sequences as bait isolated the Guanine Nucleotide-Binding Protein-Like 1 (GNL1). The latter was shown to have a higher affinity for the G4 sequences than for their mutated version. This study sheds light on new RNA G-quadruplexes located in genes dysregulated in Parkinson disease and a new G4-binding protein, GNL1.

G-quadruplex located in the 5'UTR of the BAG-1 mRNA affects both its cap-dependent and cap-independent translation through global secondary structure maintenance.

The anti-apoptotic BAG-1 protein isoforms are known to be overexpressed in colorectal tumors and are considered to be potential therapeutic targets. The isoforms are derived from alternative translation initiations occuring at four in-frame start codons of a single mRNA transcript. Its 5'UTR also contains an internal ribosome entry site (IRES) regulating the cap-independent translation of the transcript. An RNA G-quadruplex (rG4) is located at the 5'end of the BAG-1 5'UTR, upstream of the known cis-regulatory elements. Herein, we observed that the expression of BAG-1 isoforms is post-transcriptionally regulated in colorectal cancer cells and tumors, and that stabilisation of the rG4 by small molecules ligands reduces the expression of endogenous BAG-1 isoforms. We demonstrated a critical role for the rG4 in the control of both cap-dependent and independent translation of the BAG-1 mRNA in colorectal cancer cells. Additionally, we found an upstream ORF that also represses BAG-1 mRNA translation. The structural probing of the complete 5'UTR showed that the rG4 acts as a steric block which controls the initiation of translation at each start codon of the transcript and also maintains the global 5'UTR secondary structure required for IRES-dependent translation.

Citrus exocortis viroid causes ribosomal stress in tomato plants.

Viroids are naked RNAs that do not code for any known protein and yet are able to infect plants causing severe diseases. Because of their RNA nature, many studies have focused on the involvement of viroids in RNA-mediated gene silencing as being their pathogenesis mechanism. Here, the alterations caused by the Citrus exocortis viroid (CEVd) on the tomato translation machinery were studied as a new aspect of viroid pathogenesis. The presence of viroids in the ribosomal fractions of infected tomato plants was detected. More precisely, CEVd and its derived viroid small RNAs were found to co-sediment with tomato ribosomes in vivo, and to provoke changes in the global polysome profiles, particularly in the 40S ribosomal subunit accumulation. Additionally, the viroid caused alterations in ribosome biogenesis in the infected tomato plants, affecting the 18S rRNA maturation process. A higher expression level of the ribosomal stress mediator NAC082 was also detected in the CEVd-infected tomato leaves. Both the alterations in the rRNA processing and the induction of NAC082 correlate with the degree of viroid symptomatology. Taken together, these results suggest that CEVd is responsible for defective ribosome biogenesis in tomato, thereby interfering with the translation machinery and, therefore, causing ribosomal stress.

First report of Grapevine yellow speckle viroid 1 infecting grapevine (Vitis vinifera L.) in Canada.

Quebec is the third largest wine grape producer in Canada in acreage, tonnage, and wine grape sales (Carisse et al. 2017; Ben Moussa et al. 2019). To evaluate the diversity of viruses infecting grapevine in Quebec, a total of 77 leaf tissue samples (cv. Vidal) were collected from July to October in 2020 in three different vineyards located in Frelighsburg, Hemmingford and Saint-Jacques-le-Mineur in Quebec, Canada. Double-stranded RNA was extracted from each sample and used for cDNA library preparation with the Nextera XT DNA Library Preparation Kit (Illumina) as described previously (Kesanakurti et al. 2016). High-throughput sequencing (HTS, 2x300 bp) was conducted on dual-indexed libraries in a v3 flow cell using the Illumina MiSeq platform (Adkar-Purushothama et al. 2020). The obtained raw FASTQ data was de-multiplexed into 154 separate sequence files, and the adapters and barcode sequences were trimmed. The quality of the sequences was verified using Trimmomatic V.0.32 and the "clean" sequences were analyzed using Virtool and VirFind virus detection pipelines described elsewhere (Ho and Tzanetakis 2014; Rott et al. 2017) to screen for all possible viruses in the databases. Over 100,000 reads per sample were obtained with a percentage of mapped viral reads ranging from 1.47 to 19.43% of total number of reads. Out of 77 samples, 16 revealed the sequence of grapevine yellow speckle viroid 1 (GYSVd-1), for which the length coverage ranged from 98.5 to 99.1%; the depth ranged from 2X to 856X. The GYSVd-1 positive sequence files were subjected to whole genome assembly on CLC genomics Workbench v20.0.4 with the isolate SY-BR from Brazil (KU880715) used as reference. Seven complete genomes of GYSVd-1 of 366-368 nucleotides (nt) in size were deposited (GenBank Acc. MW732682 to MW732688). BLASTN analysis of the sequences showed 98-100% nt identities with isolate SY-BR. Other viruses and viroids such as Grapevine fleck virus, Grapevine rupestris stem pitting-associated virus, Grapevine rupestris vein feathering virus and Hop stunt viroid were also detected. To confirm GYSVd-1 presence in Quebec vineyards, seven of the 16 HTS-positive grapevine leaf tissue samples were subjected to total RNA extraction, followed by RT-PCR assay as before (Adkar-Purushothama et al. 2015; Sahana et al. 2013); all were positive by RT-PCR. The PCR products were directly Sanger-sequenced, and they showed 100% nt identity to the HTS derived sequences. Three of the seven GYSVd-1 positive grapevines exhibited yellow leaf spots and flecks and tiny yellow leaves, but their mixed infection status makes definitive symptoms association difficult to determine. Previously, Hop stunt viroid was reported from grapevines in Canada (Xiao et al. 2019; Fall et al. 2020) but to the best of our knowledge, this is the first report of GYSVd-1 infecting grapevines in Canada, specifically in the province of Quebec. Further research is required to assess the GYSVd-1 related yield loss. Monitoring and testing for GYSVd-1 infection is necessary to prevent propagation of infected materials, spread, and potential negative impact for the Canadian grapevine industry.

Impact of Nucleic Acid Sequencing on Viroid Biology.

The early 1970s marked two breakthroughs in the field of biology: (i) The development of nucleotide sequencing technology; and, (ii) the discovery of the viroids. The first DNA sequences were obtained by two-dimensional chromatography which was later replaced by sequencing using electrophoresis technique. The subsequent development of fluorescence-based sequencing method which made DNA sequencing not only easier, but many orders of magnitude faster. The knowledge of DNA sequences has become an indispensable tool for both basic and applied research. It has shed light biology of viroids, the highly structured, circular, single-stranded non-coding RNA molecules that infect numerous economically important plants. Our understanding of viroid molecular biology and biochemistry has been intimately associated with the evolution of nucleic acid sequencing technologies. With the development of the next-generation sequence method, viroid research exponentially progressed, notably in the areas of the molecular mechanisms of viroids and viroid diseases, viroid pathogenesis, viroid quasi-species, viroid adaptability, and viroid-host interactions, to name a few examples. In this review, the progress in the understanding of viroid biology in conjunction with the improvements in nucleotide sequencing technology is summarized. The future of viroid research with respect to the use of third-generation sequencing technology is also briefly envisaged.

3' UTR G-quadruplexes regulate miRNA binding.

MicroRNAs (miRNAs) are small noncoding RNAs that repress the translation of their target genes. It has previously been shown that a target's availability to miRNA can be affected by its structure. G-quadruplexes (G4) are noncanonical structures adopted by G-rich nucleic acids that have been shown to have multiple biological functions. In this study, whether or not G4 structures' presence in the 3' UTRs of mRNAs can hinder miRNA binding was investigated. Putative G4 overlapping with predicted miRNAs' binding sites was searched for, and 44,294 hits were found in humans. The FADS2 mRNA/mir331-3p pair was selected as a model example. In-line probing and G4-specific fluorescent ligand experiments binding were performed and confirmed the presence of a G4 near the predicted miRNA binding site. Subsequent luciferase assays showed that the presence of the G4 prevents the binding of mir331-3p in cellulo. Together, these results served as proof of concept that a G4 structure present in a 3' UTR sequence should be taken into consideration when predicting miRNA binding sites.

The Epstein-Barr virus EBNA1 protein modulates the alternative splicing of cellular genes.

BACKGROUND: Alternative splicing (AS) is an important mRNA maturation step that allows increased variability and diversity of proteins in eukaryotes. AS is dysregulated in numerous diseases, and its implication in the carcinogenic process is well known. However, progress in understanding how oncogenic viruses modulate splicing, and how this modulation is involved in viral oncogenicity has been limited. Epstein-Barr virus (EBV) is involved in various cancers, and its EBNA1 oncoprotein is the only viral protein expressed in all EBV malignancies. METHODS: In the present study, the ability of EBNA1 to modulate the AS of cellular genes was assessed using a high-throughput RT-PCR approach to examine AS in 1238 cancer-associated genes. RNA immunoprecipitation coupled to RNA sequencing (RIP-Seq) assays were also performed to identify cellular mRNAs bound by EBNA1. RESULTS: Upon EBNA1 expression, we detected modifications to the AS profiles of 89 genes involved in cancer. Moreover, we show that EBNA1 modulates the expression levels of various splicing factors such as hnRNPA1, FOX-2, and SF1. Finally, RNA immunoprecipitation coupled to RIP-Seq assays demonstrate that EBNA1 immunoprecipitates specific cellular mRNAs, but not the ones that are spliced differently in EBNA1-expressing cells. CONCLUSION: The EBNA1 protein can modulate the AS profiles of numerous cellular genes. Interestingly, this modulation protein does not require the RNA binding activity of EBNA1. Overall, these findings underline the novel role of EBNA1 as a cellular splicing modulator.

Motif independent identification of potential RNA G-quadruplexes by G4RNA screener.

MOTIVATION: G-quadruplex structures in RNA molecules are known to have regulatory impacts in cells but are difficult to locate in the genome. The minimal requirements for G-quadruplex folding in RNA (G≥3N1-7 G≥3N1-7 G≥3N1-7 G≥3) is being challenged by observations made on specific examples in recent years. The definition of potential G-quadruplex sequences has major repercussions on the observation of the structure since it introduces a bias. The canonical motif only describes a sub-population of the reported G-quadruplexes. To address these issues, we propose an RNA G-quadruplex prediction strategy that does not rely on a motif definition. RESULTS: We trained an artificial neural network with sequences of experimentally validated G-quadruplexes from the G4RNA database encoded using an abstract definition of their sequence. This artificial neural network, G4NN, evaluates the similarity of a given sequence to known G-quadruplexes and reports it as a score. G4NN has a predictive power comparable to the reported G richness and G/C skewness evaluations that are the current state-of-the-art for the identification of potential RNA G-quadruplexes. We combined these approaches in the G4RNA screener, a program designed to manage and evaluate the sequences to identify potential G-quadruplexes. AVAILABILITY AND IMPLEMENTATION: G4RNA screener is available for download at http://gitlabscottgroup.med.usherbrooke.ca/J-Michel/g4rna_screener. CONTACT: jean-michel.garant@usherbrooke.ca or jean-pierre.perreault@usherbrooke.ca or michelle.scott@usherbrooke.ca. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Structurexplor: a platform for the exploration of structural features of RNA secondary structures.

SUMMARY: Discovering function-related structural features, such as the cloverleaf shape of transfer RNA secondary structures, is essential to understand RNA function. With this aim, we have developed a platform, named Structurexplor, to facilitate the exploration of structural features in populations of RNA secondary structures. It has been designed and developed to help biologists interactively search for, evaluate and select interesting structural features that can potentially explain RNA functions. AVAILABILITY AND IMPLEMENTATION: Structurxplor is a web application available at http://structurexplor.dinf.usherbrooke.ca. The source code can be found at http://jpsglouzon.github.io/structurexplor/. CONTACT: shengrui.wang@usherbrooke.ca. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The super-n-motifs model: a novel alignment-free approach for representing and comparing RNA secondary structures.

MOTIVATION: Comparing ribonucleic acid (RNA) secondary structures of arbitrary size uncovers structural patterns that can provide a better understanding of RNA functions. However, performing fast and accurate secondary structure comparisons is challenging when we take into account the RNA configuration (i.e. linear or circular), the presence of pseudoknot and G-quadruplex (G4) motifs and the increasing number of secondary structures generated by high-throughput probing techniques. To address this challenge, we propose the super-n-motifs model based on a latent analysis of enhanced motifs comprising not only basic motifs but also adjacency relations. The super-n-motifs model computes a vector representation of secondary structures as linear combinations of these motifs. RESULTS: We demonstrate the accuracy of our model for comparison of secondary structures from linear and circular RNA while also considering pseudoknot and G4 motifs. We show that the super-n-motifs representation effectively captures the most important structural features of secondary structures, as compared to other representations such as ordered tree, arc-annotated and string representations. Finally, we demonstrate the time efficiency of our model, which is alignment free and capable of performing large-scale comparisons of 10 000 secondary structures with an efficiency up to 4 orders of magnitude faster than existing approaches. AVAILABILITY AND IMPLEMENTATION: The super-n-motifs model was implemented in C ++. Source code and Linux binary are freely available at http://jpsglouzon.github.io/supernmotifs/ . CONTACT: Shengrui.Wang@Usherbrooke.ca. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics o nline.

Small RNA Derived from the Virulence Modulating Region of the Potato spindle tuber viroid Silences callose synthase Genes of Tomato Plants.

The tomato (Solanum lycopersicum) callose synthase genes CalS11-like and CalS12-like encode proteins that are essential for the formation of callose, a major component of pollen mother cell walls; these enzymes also function in callose formation during pathogen infection. This article describes the targeting of these callose synthase mRNAs by a small RNA derived from the virulence modulating region of two Potato spindle tuber viroid variants. More specifically, viroid infection of tomato plants resulted in the suppression of the target mRNAs up to 1.5-fold, depending on the viroid variant used and the gene targeted. The targeting of these mRNAs by RNA silencing was validated by artificial microRNA experiments in a transient expression system and by RNA ligase-mediated rapid amplification of cDNA ends. Viroid mutants incapable of targeting callose synthase mRNAs failed to induce typical infection phenotypes, whereas a chimeric viroid obtained by swapping the virulence modulating regions of a mild and a severe variant of Potato spindle tuber viroid greatly affected the accumulation of viroids and the severity of disease symptoms. These data provide evidence of the silencing of multiple genes by a single small RNA derived from a viroid.

Alterations of the viroid regions that interact with the host defense genes attenuate viroid infection in host plant.

Understanding in intimate details how the viroid interaction with host's defense genes is a cornerstone for developing viroid resistant plants. In this present study, small RNAs (sRNA) derived from Potato spindle tuber viroid (PSTVd) were studied in silico in order to detect any interactions with the serine threonine kinase receptor, a transmembrane protein that plays a role in disease resistance in plants. Using molecular biology techniques, it was determined that PSTVd infection negatively affects at least three serine threonine kinase receptors as well as with three other genes that are known to be involved in the overall development of the tomato plants. The transient expression of these putative PSTVd-sRNAs, using the microRNA sequence as a backbone, in tomato plants induced phenotypes similar to viroid infection. Mutants created by altering the sequence of PSTVd in these regions failed to infect the tomato plant. The data presented here illustrates the importance of these regions in viroid survival, and suggests a possible avenue of exploration for the development of viroid resistant plants.

G-Quadruplexes influence pri-microRNA processing.

RNA G-Quadruplexes (G4) have been shown to possess many biological functions, including the regulation of microRNA (miRNA) biogenesis and function. However, their impact on pri-miRNA processing remains unknown. We identified G4 located near the Drosha cleavage site in three distinct pri-miRNAs: pri-mir200c, pri-mir451a, and pri-mir497. The folding of the potential G4 motifs was determined in solution. Subsequently, mutations disrupting G4 folding led to important changes in the mature miRNAs levels in cells. Moreover, using small antisense oligonucleotides binding to the pri-miRNA, it was possible to modulate, either positively or negatively, the mature miRNA levels. Together, these data demonstrate that G4 motifs could contribute to the regulation of pri-mRNA processing, a novel role for G4. Considering that bio-informatics screening indicates that between 9% and 50% of all pri-miRNAs contain a putative G4, these structures possess interesting potential as future therapeutic targets.

Irregular G-quadruplexes Found in the Untranslated Regions of Human mRNAs Influence Translation.

G-quadruplex structures are composed of coplanar guanines and are found in both DNA and RNA. They are formed by the stacking of two or more G-quartets that are linked together by three loops. The current belief is that RNA G-quadruplexes include loops of l to 7 nucleotides in length, although recent evidence indicates that the central loop (loop 2) can be longer if loops 1 and 3 are limited to a single nucleotide each. With the objective of broadening the definition of irregular RNA G-quadruplexes, a bioinformatic search was performed to find potential G-quadruplexes located in the untranslated regions of human mRNAs (i.e. in the 5' and 3'-UTRs) that contain either a long loop 1 or 3 of up to 40 nucleotides in length. RNA molecules including the potential sequences were then synthesized and examined in vitro by in-line probing for the formation of G-quadruplex structures. The sequences that adopted a G-quadruplex structure were cloned into a luciferase dual vector and examined for their ability to modulate translation in cellulo Some irregular G-quadruplexes were observed to either promote or repress translation regardless of the position or the size of the long loop they possessed. Even if the composition of a RNA G-quadruplex is not quite completely understood, the results presented in this report clearly demonstrate that what defines a RNA G-quadruplex is much broader than what we previously believed.

Small antisense oligonucleotides against G-quadruplexes: specific mRNA translational switches.

G-quadruplexes (G4) are intricate RNA structures found throughout the transcriptome. Because they are associated with a variety of biological cellular mechanisms, these fascinating structural motifs are seen as potential therapeutic targets against many diseases. While screening of chemical compounds specific to G4 motifs has yielded interesting results, no single compound successfully discriminates between G4 motifs based on nucleotide sequences alone. This level of specificity is best attained using antisense oligonucleotides (ASO). Indeed, oligonucleotide-based strategies are already used to modulate DNA G4 folding in vitro. Here, we report that, in human cells, the use of short ASO to promote and inhibit RNA G4 folding affects the translation of specific mRNAs, including one from the 5'UTR of the H2AFY gene, a histone variant associated with cellular differentiation and cancer. These results suggest that the relatively high specificity of ASO-based strategies holds significant potential for applications aimed at modulating G4-motif folding.

Global profiling of alternative RNA splicing events provides insights into molecular differences between various types of hepatocellular carcinoma.

BACKGROUND: Dysregulations in alternative splicing (AS) patterns have been associated with many human diseases including cancer. In the present study, alterations to the global RNA splicing landscape of cellular genes were investigated in a large-scale screen from 377 liver tissue samples using high-throughput RNA sequencing data. RESULTS: Our study identifies modifications in the AS patterns of transcripts encoded by more than 2500 genes such as tumor suppressor genes, transcription factors, and kinases. These findings provide insights into the molecular differences between various types of hepatocellular carcinoma (HCC). Our analysis allowed the identification of 761 unique transcripts for which AS is misregulated in HBV-associated HCC, while 68 are unique to HCV-associated HCC, 54 to HBV&HCV-associated HCC, and 299 to virus-free HCC. Moreover, we demonstrate that the expression pattern of the RNA splicing factor hnRNPC in HCC tissues significantly correlates with patient survival. We also show that the expression of the HBx protein from HBV leads to modifications in the AS profiles of cellular genes. Finally, using RNA interference and a reverse transcription-PCR screening platform, we examined the implications of cellular proteins involved in the splicing of transcripts involved in apoptosis and demonstrate the potential contribution of these proteins in AS control. CONCLUSIONS: This study provides the first comprehensive portrait of global changes in the RNA splicing signatures that occur in hepatocellular carcinoma. Moreover, these data allowed us to identify unique signatures of genes for which AS is misregulated in the different types of HCC.

The folding of 5'-UTR human G-quadruplexes possessing a long central loop.

G-quadruplexes are widespread four-stranded structures that are adopted by G-rich regions of both DNA and RNA and are involved in essential biological processes such as mRNA translation. They are formed by the stacking of two or more G-quartets that are linked together by three loops. Although the maximal loop length is usually fixed to 7 nt in most G-quadruplex-predicting software, it has already been demonstrated that artificial DNA G-quadruplexes containing two distal loops that are limited to 1 nt each and a central loop up to 30 nt long are likely to form in vitro. This report demonstrates that such structures possessing a long central loop are actually found in the 5'-UTRs of human mRNAs. Firstly, 1453 potential G-quadruplex-forming sequences (PG4s) were identified through a bioinformatic survey that searched for sequences respecting the requirement for two 1-nt long distal loops and a long central loop of 2-90 nt in length. Secondly, in vitro in-line probing experiments confirmed and characterized the folding of eight candidates possessing central loops of 10-70 nt long. Finally, the biological effect of several G-quadruplexes with a long central loop on mRNA expression was studied in cellulo using a luciferase gene reporter assay. Clearly, the actual definition of G-quadruplex-forming sequences is too conservative and must be expanded to include the long central loop. This greatly expands the number of expected PG4s in the transcriptome. Consideration of these new candidates might aid in elucidating the potentially important biological implications of the G-quadruplex structure.

Tumor microenvironment-associated modifications of alternative splicing.

Pre-mRNA alternative splicing is modified in cancer, but the origin and specificity of these changes remain unclear. Here, we probed ovarian tumors to identify cancer-associated splicing isoforms and define the mechanism by which splicing is modified in cancer cells. Using high-throughput quantitative PCR, we monitored the expression of splice variants in laser-dissected tissues from ovarian tumors. Surprisingly, changes in alternative splicing were not limited to the tumor tissues but were also found in the tumor microenvironment. Changes in the tumor-associated splicing events were found to be regulated by splicing factors that are differentially expressed in cancer tissues. Overall, ∼20% of the alternative splicing events affected by the down-regulation of the splicing factors QKI and RBFOX2 were altered in the microenvironment of ovarian tumors. Together, our results indicate that the tumor microenvironment undergoes specific changes in alternative splicing orchestrated by a limited number of splicing factors.