CLASH Analyst: A Web Server to Identify In Vivo RNA–RNA Interactions from CLASH Data

Non-coding RNAs, such as miRNAs and piRNAs, play critical roles in gene regulation through base-pairing interactions with their target molecules. The recent development of the crosslinking, ligation, and sequencing of hybrids (CLASH) method has allowed scientists to map transcriptome-wide RNA–RNA interactions by identifying chimeric reads consisting of fragments from regulatory RNAs and their targets. However, analyzing CLASH data requires scientists to use advanced bioinformatics, and currently available tools are limited for users with little bioinformatic experience. In addition, many published CLASH studies do not show the full scope of RNA–RNA interactions that were captured, highlighting the importance of reanalyzing published data. Here, we present CLASH Analyst, a web server that can analyze raw CLASH data within a fully customizable and easy-to-use interface. CLASH Analyst accepts raw CLASH data as input and identifies the RNA chimeras containing the regulatory and target RNAs according to the user’s interest. Detailed annotation of the captured RNA–RNA interactions is then presented for the user to visualize within the server or download for further analysis. We demonstrate that CLASH Analyst can identify miRNA- and piRNA-targeting sites reported from published CLASH data and should be applicable to analyze other RNA–RNA interactions. CLASH Analyst is freely available for academic use.

Hierarchical length and sequence preferences establish a single major piRNA 3'-end

PIWI-interacting RNAs (piRNAs) guard germline genomes against the deleterious action of retroviruses and other mobile genetic elements. How piRNAs faithfully discriminate between self and non-self to restrict all mobile elements while sparing essential genes remains a key outstanding question in genome biology. PiRNAs use extensive base-pairing to recognize their targets and variable 3'ends could change the specificity and efficacy of piRNA silencing. Here, we identify conserved rules that ensure the generation of a single major piRNA 3'end in flies and mice. Our data suggest that the PIWI proteins initially define a short interval on pre-piRNAs that grants access to the ZUC-processor complex. Within this Goldilocks zone, the preference of the ZUC-processor to cut in front of a Uridine determines the ultimate processing site. We observe a mouse-specific roadblock that relocates the Goldilocks zone and generates an opportunity for consecutive trimming by PNLDC1. Our data reveal a conserved hierarchy between length and sequence preferences that controls the piRNA sequence space. The unanticipated precision of 3'end formation bolsters the emerging understanding that the functional piRNA sequence space is tightly controlled to ensure effective defense.

Uveal Melanoma: a miR-16 disease?

Uveal melanoma (UM), the most common primary intraocular tumor in adults, has been extensively characterized by omics technologies during the last 5 years. Despite the discovery of gene signatures, the molecular actors driving the cancer aggressiveness are not fully understood and UM is still associated to a dismal overall survival at metastatic stage. Here, we showed that microRNA-16 (miR-16) is involved in uveal melanoma by an unexpected mechanism. By defining the miR-16-interactome, we revealed that miR-16 mainly interacts via non-canonical base-pairing to a subset of RNAs, promoting their expression levels (sponge RNAs). Consequently, the canonical miR-16 activity, involved in the RNA decay of oncogenes such as cyclin D1 and D3, is impaired. This miR-16 non-canonical base-pairing to sponge RNAs can explain both the derepression of miR-16 targets and the promotion of oncogenes expression observed for patients with poor overall survival in two cohorts. miR-16 activity assessment using our sponge-signature discriminates the patients overall survival as efficiently as the current method based on copy number variations and driver mutations detection. To conclude, miRNA loss of function due to miRNA sequestration seems to promote cancer burden by two combined events; loss of brake and an acceleration. Our results highlight the oncogenic role of the non-canonical base-pairing between miRNAs/mRNAs in uveal melanoma.

Widespread microRNA degradation elements in target mRNAs can assist the encoded proteins

Binding of microRNAs (miRNAs) to mRNAs normally results in post-transcriptional repression of gene expression. However, extensive base-pairing between miRNAs and target RNAs can trigger miRNA degradation, a phenomenon called target RNA-directed miRNA degradation (TDMD). Here, we systematically analyzed Argonaute-CLASH (cross-linking, ligation, and sequencing of miRNA–target RNA hybrids) data and identified numerous candidate TDMD triggers, focusing on their ability to induce nontemplated nucleotide addition at the miRNA 3′ end. When exogenously expressed in various cell lines, eight triggers induce degradation of corresponding miRNAs. Both the TDMD base-pairing and surrounding sequences are essential for TDMD. CRISPR knockout of endogenous trigger or ZSWIM8, a ubiquitin ligase essential for TDMD, reduced miRNA degradation. Furthermore, degradation of miR-221 and miR-222 by a trigger in BCL2L11, which encodes a proapoptotic protein, enhances apoptosis. Therefore, we uncovered widespread TDMD triggers in target RNAs and demonstrated an example that could functionally cooperate with the encoded protein.

Constraint of Base Pairing on HDV Genome Evolution

The hepatitis delta virus is a single-stranded circular RNA virus, which is characterized by high self-complementarity. About 70% of the genome sequences can form base-pairs with internal nucleotides. There are many studies on the evolution of the hepatitis delta virus. However, the secondary structure has not been taken into account in these studies. In this study, we developed a method to examine the effect of base pairing as a constraint on the nucleotide substitutions during the evolution of the hepatitis delta virus. The method revealed that the base pairing can reduce the evolutionary rate in the non-coding region of the virus. In addition, it is suggested that the non-coding nucleotides without base pairing may be under some constraint, and that the intensity of the constraint is weaker than that by the base pairing but stronger than that on the synonymous site.

Complementary Nucleobase Interactions Drive Co-Assembly of Drugs and Nanocarriers for Selective Cancer Chemotherapy

A new concept in cooperative adenine–uracil (A–U) hydrogen bonding interactions between anticancer drugs and nanocarrier complexes was successfully demonstrated by invoking the co-assembly of water soluble, uracil end-capped polyethylene glycol polymer (BU-PEG) upon association with the hydrophobic drug adenine-modified rhodamine (A-R6G). This concept holds promise as a smart and versatile drug delivery system for the achievement of targeted, more efficient cancer chemotherapy. Due to A–U base pairing between BU-PEG and A-R6G, BU-PEG has high tendency to interact with A-R6G, which leads to the formation of self-assembled A-R6G/BU-PEG nanogels in aqueous solution. The resulting nanogels exhibit a number of unique physical properties, including extremely high A-R6G-loading capacity, well-controlled, pH-triggered A-R6G release behavior, and excellent structural stability in biological media. Importantly, a series of in vitro cellular experiments clearly demonstrated that A-R6G/BU-PEG nanogels improved the selective uptake of A-R6G by cancer cells via endocytosis and promoted the intracellular release of A-R6G to subsequently induce apoptotic cell death, while control rhodamine/BU-PEG nanogels did not exert selective toxicity in cancer or normal cell lines. Overall, these results indicate that cooperative A–U base pairing within nanogels is a critical factor that improves selective drug uptake and effectively promotes apoptotic programmed cell death in cancer cells.

Trans-Activator Binding Site Context in RCNMV Modulates Subgenomic mRNA Transcription

Many positive-sense RNA viruses transcribe subgenomic (sg) mRNAs during infections that template the translation of a subset of viral proteins. Red clover necrotic mosaic virus (RCNMV) expresses its capsid protein through the transcription of a sg mRNA from RNA1 genome segment. This transcription event is activated by an RNA structure formed by base pairing between a trans-activator (TA) in RNA2 and a trans-activator binding site (TABS) in RNA1. In this study, the impact of the structural context of the TABS in RNA1 on the TA–TABS interaction and sg mRNA transcription was investigated using in vitro and in vivo approaches. The results (i) generated RNA secondary structure models for the TA and TABS, (ii) revealed that the TABS is partially base paired with proximal upstream sequences, which limits TA access, (iii) demonstrated that the aforementioned intra-RNA1 base pairing involving the TABS modulates the TA–TABS interaction in vitro and sg mRNA levels during infections, and (iv) revealed that the TABS in RNA1 can be modified to mediate sg mRNA transcription in a TA-independent manner. These findings advance our understanding of transcriptional regulation in RCNMV and provide novel insights into the origin of the TA–TABS interaction.

Evidence for a long-r ange RNA-RNA interaction between ORF8 and the downstream region of the Spike polybasic insertion of SARS-CoV-2

SARS-CoV-2 has affected people all over the world as the causative agent of COVID-19. The virus is related to the highly lethal SARS-CoV responsible for the 2002-2003 SARS outbreak in Asia. Intense research is ongoing to understand why both viruses have different spreading capacities and mortality rates. Similar to other betacoronaviruses, long-range RNA-RNA interactions occur between different parts of the viral genomic RNA, resulting in discontinuous transcription and production of various sub-genomic RNAs. These sub-genomic RNAs are then translated into different viral proteins. An important difference between both viruses is a polybasic insertion in the Spike region of SARS-CoV-2, absent in SARS-CoV. Here we show that a 26-base-pair long-range RNA-RNA interaction occurs between the genomic region downstream of the Spike insertion and ORF8 in SARS-CoV-2. Predictions suggest that the corresponding ORF8 region forms the most energetically favorable interaction with that of Spike region from amongst all possible candidate regions within SARS-CoV-2 genomic RNA. We also found signs of sequence covariation in the predicted interaction using a large dataset with 27,592 full-length SARS-CoV-2 genomes. In particular, a synonymous mutation in ORF8 accommodated for base pairing with Spike [G23675 C28045U], and a non-synonymous mutation in Spike accommodated for base pairing with ORF8 [C23679U G28042] both of which were in close proximity of one another. The predicted interactions can potentially be related to regulation of sub-genomic RNA production rates.