Post-transcriptional regulatory patterns revealed by protein-RNA interactions

RNA methylations in hepatic fibrosis, a gradually emerging new treatment strategy

BACKGROUND

Hepatic fibrosis (HF) is a pathological process caused by excessive accumulation of extracellular matrix caused by a series of causes, leading to the formation of fiber scar. RNA methylation is a newly discovered epigenetic modification that exists widely in eukaryotes and prokaryotes and plays a crucial role in the pathogenesis of many diseases.

RESULTS

The occurrence and development of HF are regulated by many factors, including excessive deposition of extracellular matrix, activation of hepatic stellate cells, inflammation, and oxidative stress. RNA methylations of different species have become a crucial regulatory mode of transcript expression, And participate in the pathogenesis of tumors, nervous system diseases, autoimmune diseases, and other diseases. In addition, there are five common types of RNA methylation, but only m6A plays a crucial regulatory role in HF. The pathophysiological regulation of m6A on HF is achieved by the combination of the methylated transferase, demethylated enzyme, and methylated reading protein.

CONCLUSIONS

RNA methylated methyltransferase, demethylase, and reading protein extensively affect the pathological mechanism of HF, which may be a new therapeutic and diagnostic target, representing a new class of therapeutic strategies.

A proteome-scale map of the SARS-CoV-2-human contactome

Understanding the mechanisms of coronavirus disease 2019 (COVID-19) disease severity to efficiently design therapies for emerging virus variants remains an urgent challenge of the ongoing pandemic. Infection and immune reactions are mediated by direct contacts between viral molecules and the host proteome, and the vast majority of these virus-host contacts (the 'contactome') have not been identified. Here, we present a systematic contactome map of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with the human host encompassing more than 200 binary virus-host and intraviral protein-protein interactions. We find that host proteins genetically associated with comorbidities of severe illness and long COVID are enriched in SARS-CoV-2 targeted network communities. Evaluating contactome-derived hypotheses, we demonstrate that viral NSP14 activates nuclear factor κB (NF-κB)-dependent transcription, even in the presence of cytokine signaling. Moreover, for several tested host proteins, genetic knock-down substantially reduces viral replication. Additionally, we show for USP25 that this effect is phenocopied by the small-molecule inhibitor AZ1. Our results connect viral proteins to human genetic architecture for COVID-19 severity and offer potential therapeutic targets.

Pan-cancer analysis of mRNA stability for decoding tumour post-transcriptional programs

Measuring mRNA decay in tumours is a prohibitive challenge, limiting our ability to map the post-transcriptional programs of cancer. Here, using a statistical framework to decouple transcriptional and post-transcriptional effects in RNA-seq data, we uncover the mRNA stability changes that accompany tumour development and progression. Analysis of 7760 samples across 18 cancer types suggests that mRNA stability changes are ~30% as frequent as transcriptional events, highlighting their widespread role in shaping the tumour transcriptome. Dysregulation of programs associated with >80 RNA-binding proteins (RBPs) and microRNAs (miRNAs) drive these changes, including multi-cancer inactivation of RBFOX and miR-29 families. Phenotypic activation or inhibition of RBFOX1 highlights its role in calcium signaling dysregulation, while modulation of miR-29 shows its impact on extracellular matrix organization and stemness genes. Overall, our study underlines the integral role of mRNA stability in shaping the cancer transcriptome, and provides a resource for systematic interrogation of cancer-associated stability pathways.

The role of mRNA stability in shaping the cancer transcriptome is revealed using a statistical analysis of transcriptomic data.

Differential Effects on the Translation of Immune-Related Alternatively Polyadenylated mRNAs in Melanoma and T Cells by eIF4A Inhibition

Targeting the translation initiation complex eIF4F, which binds the 5' cap of mRNAs, is a promising anti-cancer approach. Silvestrol, a small molecule inhibitor of eIF4A, the RNA helicase component of eIF4F, inhibits the translation of the mRNA encoding the signal transducer and activator of transcription 1 (STAT1) transcription factor, which, in turn, reduces the transcription of the gene encoding one of the major immune checkpoint proteins, i.e., programmed death ligand-1 (PD-L1) in melanoma cells. A large proportion of human genes produce multiple mRNAs differing in their 3'-ends through the use of alternative polyadenylation (APA) sites, which, when located in alternative last exons, can generate protein isoforms, as in the STAT1 gene. Here, we provide evidence that the STAT1α, but not STAT1β protein isoform generated by APA, is required for silvestrol-dependent inhibition of PD-L1 expression in interferon-γ-treated melanoma cells. Using polysome profiling in activated T cells we find that, beyond STAT1, eIF4A inhibition downregulates the translation of some important immune-related mRNAs, such as the ones encoding TIM-3, LAG-3, IDO1, CD27 or CD137, but with little effect on the ones for BTLA and ADAR-1 and no effect on the ones encoding CTLA-4, PD-1 and CD40-L. We next apply RT-qPCR and 3'-seq (RNA-seq focused on mRNA 3' ends) on polysomal RNAs to analyze in a high throughput manner the effect of eIF4A inhibition on the translation of APA isoforms. We identify about 150 genes, including TIM-3, LAG-3, AHNAK and SEMA4D, for which silvestrol differentially inhibits the translation of APA isoforms in T cells. It is therefore crucial to consider 3'-end mRNA heterogeneity in the understanding of the anti-tumor activities of eIF4A inhibitors.

Desmosomes as Signaling Hubs in the Regulation of Cell Behavior

Desmosomes are intercellular junctions, which preserve tissue integrity during homeostatic and stress conditions. These functions rely on their unique structural properties, which enable them to respond to context-dependent signals and transmit them to change cell behavior. Desmosome composition and size vary depending on tissue specific expression and differentiation state. Their constituent proteins are highly regulated by posttranslational modifications that control their function in the desmosome itself and in addition regulate a multitude of desmosome-independent functions. This review will summarize our current knowledge how signaling pathways that control epithelial shape, polarity and function regulate desmosomes and how desmosomal proteins transduce these signals to modulate cell behavior.

Generation of 3'UTR knockout cell lines by CRISPR/Cas9-mediated genome editing

In addition to the protein code, messenger RNAs (mRNAs) also contain untranslated regions (UTRs). 3'UTRs span the region between the translational stop codon and the poly(A) tail. Sequence elements located in 3'UTRs are essential for pre-mRNA processing. 3'UTRs also contain elements that can regulate protein abundance, localization, and function. At least half of all human genes use alternative cleavage and polyadenylation (APA) to further diversify the regulatory potential of protein functions. Traditional gene editing approaches are designed to disrupt the production of functional proteins. Here, we describe a method that allows investigators to manipulate 3'UTR sequences of endogenous genes for both single- 3'UTR and multi-3'UTR genes. As 3'UTRs can regulate individual functions of proteins, techniques to manipulate 3'UTRs at endogenous gene loci will help to disentangle multi-functionality of proteins. Furthermore, the ability to directly examine the impact of gene regulatory elements in 3'UTRs will provide further insights into their functional significance.

Endogenous p53 expression in human and mouse is not regulated by its 3'UTR

The TP53 gene encodes the tumor suppressor p53 which is functionally inactivated in many human cancers. Numerous studies suggested that 3′UTR-mediated p53 expression regulation plays a role in tumorigenesis and could be exploited for therapeutic purposes. However, these studies did not investigate post-transcriptional regulation of the native TP53 gene. Here, we used CRISPR/Cas9 to delete the human and mouse TP53/Trp53 3′UTRs while preserving endogenous mRNA processing. This revealed that the endogenous 3′UTR is not involved in regulating p53 mRNA or protein expression neither in steady state nor after genotoxic stress. Using reporter assays, we confirmed the previously observed repressive effects of the isolated 3′UTR. However, addition of the TP53 coding region to the reporter had a dominant negative impact on expression as its repressive effect was stronger and abrogated the contribution of the 3′UTR. Our data highlight the importance of genetic models in the validation of post-transcriptional gene regulatory effects.

A mark of disease: how mRNA modifications shape genetic and acquired pathologies

RNA modifications have recently emerged as a widespread and complex facet of gene expression regulation. Counting more than 170 distinct chemical modifications with far-reaching implications for RNA fate, they are collectively referred to as the epitranscriptome. These modifications can occur in all RNA species, including messenger RNAs (mRNAs) and noncoding RNAs (ncRNAs). In mRNAs the deposition, removal, and recognition of chemical marks by writers, erasers and readers influence their structure, localization, stability, and translation. In turn, this modulates key molecular and cellular processes such as RNA metabolism, cell cycle, apoptosis, and others. Unsurprisingly, given their relevance for cellular and organismal functions, alterations of epitranscriptomic marks have been observed in a broad range of human diseases, including cancer, neurological and metabolic disorders. Here, we will review the major types of mRNA modifications and editing processes in conjunction with the enzymes involved in their metabolism and describe their impact on human diseases. We present the current knowledge in an updated catalog. We will also discuss the emerging evidence on the crosstalk of epitranscriptomic marks and what this interplay could imply for the dynamics of mRNA modifications. Understanding how this complex regulatory layer can affect the course of human pathologies will ultimately lead to its exploitation toward novel epitranscriptomic therapeutic strategies.

Transcriptome-wide high-throughput mapping of protein-RNA occupancy profiles using POP-seq

Interaction between proteins and RNA is critical for post-transcriptional regulatory processes. Existing high throughput methods based on crosslinking of the protein–RNA complexes and poly-A pull down are reported to contribute to biases and are not readily amenable for identifying interaction sites on non poly-A RNAs. We present Protein Occupancy Profile-Sequencing (POP-seq), a phase separation based method in three versions, one of which does not require crosslinking, thus providing unbiased protein occupancy profiles on whole cell transcriptome without the requirement of poly-A pulldown. Our study demonstrates that ~ 68% of the total POP-seq peaks exhibited an overlap with publicly available protein–RNA interaction profiles of 97 RNA binding proteins (RBPs) in K562 cells. We show that POP-seq variants consistently capture protein–RNA interaction sites across a broad range of genes including on transcripts encoding for transcription factors (TFs), RNA-Binding Proteins (RBPs) and long non-coding RNAs (lncRNAs). POP-seq identified peaks exhibited a significant enrichment (p value < 2.2e−16) for GWAS SNPs, phenotypic, clinically relevant germline as well as somatic variants reported in cancer genomes, suggesting the prevalence of uncharacterized genomic variation in protein occupied sites on RNA. We demonstrate that the abundance of POP-seq peaks increases with an increase in expression of lncRNAs, suggesting that highly expressed lncRNA are likely to act as sponges for RBPs, contributing to the rewiring of protein–RNA interaction network in cancer cells. Overall, our data supports POP-seq as a robust and cost-effective method that could be applied to primary tissues for mapping global protein occupancies.

Role of SARS-CoV-2 in Altering the RNA-Binding Protein and miRNA-Directed Post-Transcriptional Regulatory Networks in Humans

The outbreak of a novel coronavirus SARS-CoV-2 responsible for the COVID-19 pandemic has caused a worldwide public health emergency. Due to the constantly evolving nature of the coronaviruses, SARS-CoV-2-mediated alterations on post-transcriptional gene regulations across human tissues remain elusive. In this study, we analyzed publicly available genomic datasets to systematically dissect the crosstalk and dysregulation of the human post-transcriptional regulatory networks governed by RNA-binding proteins (RBPs) and micro-RNAs (miRs) due to SARS-CoV-2 infection. We uncovered that 13 out of 29 SARS-CoV-2-encoded proteins directly interacted with 51 human RBPs, of which the majority of them were abundantly expressed in gonadal tissues and immune cells. We further performed a functional analysis of differentially expressed genes in mock-treated versus SARS-CoV-2-infected lung cells that revealed enrichment for the immune response, cytokine-mediated signaling, and metabolism-associated genes. This study also characterized the alternative splicing events in SARS-CoV-2-infected cells compared to the control, demonstrating that skipped exons and mutually exclusive exons were the most abundant events that potentially contributed to differential outcomes in response to the viral infection. A motif enrichment analysis on the RNA genomic sequence of SARS-CoV-2 clearly revealed the enrichment for RBPs such as SRSFs, PCBPs, ELAVs, and HNRNPs, suggesting the sponging of RBPs by the SARS-CoV-2 genome. A similar analysis to study the interactions of miRs with SARS-CoV-2 revealed functionally important miRs that were highly expressed in immune cells, suggesting that these interactions may contribute to the progression of the viral infection and modulate the host immune response across other human tissues. Given the need to understand the interactions of SARS-CoV-2 with key post-transcriptional regulators in the human genome, this study provided a systematic computational analysis to dissect the role of dysregulated post-transcriptional regulatory networks controlled by RBPs and miRs across tissue types during a SARS-CoV-2 infection.

Role of SARS-CoV-2 in altering the RNA binding protein and miRNA directed post-transcriptional regulatory networks in humans

The outbreak of a novel coronavirus SARS-CoV-2 responsible for COVID-19 pandemic has caused worldwide public health emergency. Due to the constantly evolving nature of the coronaviruses, SARS-CoV-2 mediated alteration on post-transcriptional gene regulation across human tissues remains elusive. In this study, we analyze publicly available genomic datasets to systematically dissect the crosstalk and dysregulation of human post-transcriptional regulatory networks governed by RNA binding proteins (RBPs) and micro-RNAs (miRs), due to SARS-CoV-2 infection. We uncovered that 13 out of 29 SARS-CoV-2 encoded proteins directly interact with 51 human RBPs of which majority of them were abundantly expressed in gonadal tissues and immune cells. We further performed a functional analysis of differentially expressed genes in mock-treated versus SARS-CoV-2 infected lung cells that revealed enrichment for immune response, cytokine-mediated signaling, and metabolism associated genes. This study also characterized the alternative splicing events in SARS-CoV-2 infected cells compared to control demonstrating that skipped exons and mutually exclusive exons were the most abundant events that potentially contributed to differential outcomes in response to viral infection. Motif enrichment analysis on the RNA genomic sequence of SARS-CoV-2 clearly revealed the enrichment for RBPs such as SRSFs, PCBPs, ELAVs, and HNRNPs suggesting the sponging of RBPs by SARS-CoV-2 genome. A similar analysis to study the interactions of miRs with SARS-CoV-2 revealed functionally important miRs that were highly expressed in immune cells, suggesting that these interactions may contribute to the progression of the viral infection and modulate host immune response across other human tissues. Given the need to understand the interactions of SARS-CoV-2 with key post-transcriptional regulators in the human genome, this study provides a systematic computational analysis to dissect the role of dysregulated post-transcriptional regulatory networks controlled by RBPs and miRs, across tissues types during SARS-CoV-2 infection.

The role of 3'UTR-protein complexes in the regulation of protein multifunctionality and subcellular localization

Multifunctional proteins often perform their different functions when localized in different subcellular compartments. However, the mechanisms leading to their localization are largely unknown. Recently, 3'UTRs were found to regulate the cellular localization of newly synthesized proteins through the formation of 3'UTR-protein complexes. Here, we investigate the formation of 3'UTR-protein complexes involving multifunctional proteins by exploiting large-scale protein-protein and protein-RNA interaction networks. Focusing on 238 human 'extreme multifunctional' (EMF) proteins, we predicted 1411 3'UTR-protein complexes involving 54% of those proteins and evaluated their role in regulating protein cellular localization and multifunctionality. We find that EMF proteins lacking localization addressing signals, yet present at both the nucleus and cell surface, often form 3'UTR-protein complexes, and that the formation of these complexes could provide EMF proteins with the diversity of interaction partners necessary to their multifunctionality. Our findings are reinforced by archetypal moonlighting proteins predicted to form 3'UTR-protein complexes. Finally, the formation of 3'UTR-protein complexes that involves up to 17% of the proteins in the human protein-protein interaction network, may be a common and yet underestimated protein trafficking mechanism, particularly suited to regulate the localization of multifunctional proteins.

Overexpression of long non-coding RNA ANRIL promotes post-ischaemic angiogenesis and improves cardiac functions by targeting Akt

Angiogenesis is critical for re‐establishing the blood supply to the surviving myocardium after myocardial infarction (MI). Long non‐coding RNA ANRIL (lncRNA‐ANRIL) has been reported to regulate endothelial functions in cardiovascular diseases. This study was to determine the role of lncRNA‐ANRIL in Akt regulation and cardiac functions after MI. Human umbilical vein endothelial cells (HUVECs) were exposed to oxygen‐glucose deprivation (OGD) to mimic in vivo ischaemia. The MI model in mice was induced by ligating left anterior descending coronary artery. OGD remarkably decreased lncRNA‐ANRIL expression level, reduced the phosphorylated levels of Akt and eNOS proteins, and inhibited NO release and cell viability, which were duplicated by shRNA‐mediated gene knockdown of lncRNA‐ANRIL. Conversely, all these effects induced by OGD were abolished by adenovirus‐mediated overexpression of lncRNA‐ANRIL in HUVECs. Further, OGD impaired cell migrations and tube formations in HUVECs, which were reversed by lncRNA‐ANRIL overexpression or Akt up‐regulation. RNA immunoprecipitation analysis indicated that the affinity of lncRNA‐ANRIL to Akt protein was increased in OGD‐treated cells. In animal studies, adenovirus‐mediated lncRNA‐ANRIL overexpression increased the phosphorylated levels of Akt and eNOS, promoted post‐ischaemic angiogenesis and improved heart functions in mice with MI surgery. LncRNA‐ANRIL regulates Akt phosphorylation to improve endothelial functions, which promotes angiogenesis and improves cardiac functions in mice following MI. In this perspective, targeting lncRNA‐ANRIL/Akt may be considered to develop a drug to treat angiogenesis‐related diseases.

The Interplay Between the DNA Damage Response, RNA Processing and Extracellular Vesicles

RNA processing was recently found to affect DNA damage response. The RNA processing factors THRAP3 and BCLAF1 play critical role in keeping DNA genomic stability by regulating the transcription, mRNA splicing and export of DNA repair proteins BRCA2, PALB2, Rad51, FANCD2, and FANCL in response to DNA damage. RNA processing factors THRAP3 and BCLAF1 play critical roles in maintaining DNA genomic stability. These factors regulate transcription, mRNA splicing and nuclear RNA export of DNA repair proteins BRCA2, PALB2, Rad51, FANCD2, and FANCL in response to DNA damage. Splicing factors SRSF10 and Sam68 were found to control the DNA damage agent-induced mRNA splicing of transcripts including BCLAF1, BRCA1, BCL2L1, CASP8, CHK2, and RBBP8 to regulate apoptosis, cell-cycle transition and DNA repair. Splicing factors and RNA binding proteins (RBPs) were also found to play a critical role in DNA/RNA hybrids (R-loops) formed during transcription and RNA processing to prevent RNA-induced genome instability. At the same time, DNA repair proteins FANCI and FANCD2 were found to regulate the nuclear localization of splicing factors SF3B1 in the DNA damage response. In addition, tumor-derived extracellular vesicles (Evs) enhanced by chemotherapeutic agents in cancer were found to promote cancer metastasis and drug resistance. Inhibiting Evs from cancer cells significantly reduced cancer metastasis and drug resistance. Furthermore, cross-talk between the DNA damage response and the immune response was observed including the enhancement of the efficacy of immune checkpoint blockade by PARP inhibitors and the effect of PD-L1 on mRNA stability of various mRNAs involved in DNA damage response by acting as a novel RNA binding protein to increase drug resistance in cancer cells. This review will introduce recent progress on the interplay of the DNA damage response, the RNA processing and the extracellular vesicles mediated metastasis.

Shaping the regulation of the p53 mRNA tumour suppressor: the co-evolution of genetic signatures

Structured RNA regulatory motifs exist from the prebiotic stages of the RNA world to the more complex eukaryotic systems. In cases where a functional RNA structure is within the coding sequence a selective pressure drives a parallel co-evolution of the RNA structure and the encoded peptide domain. The p53-MDM2 axis, describing the interactions between the p53 tumor suppressor and the MDM2 E3 ubiquitin ligase, serves as particularly useful model revealing how secondary RNA structures have co-evolved along with corresponding interacting protein motifs, thus having an impact on protein - RNA and protein - protein interactions; and how such structures developed signal-dependent regulation in mammalian systems. The p53(BOX-I) RNA sequence binds the C-terminus of MDM2 and controls p53 synthesis while the encoded peptide domain binds MDM2 and controls p53 degradation. The BOX-I peptide domain is also located within p53 transcription activation domain. The folding of the p53 mRNA structure has evolved from temperature-regulated in pre-vertebrates to an ATM kinase signal-dependent pathway in mammalian cells. The protein - protein interaction evolved in vertebrates and became regulated by the same signaling pathway. At the same time the protein - RNA and protein - protein interactions evolved, the p53 trans-activation domain progressed to become integrated into a range of cellular pathways. We discuss how a single synonymous mutation in the BOX-1, the p53(L22 L), observed in a chronic lymphocyte leukaemia patient, prevents the activation of p53 following DNA damage. The concepts analysed and discussed in this review may serve as a conceptual mechanistic paradigm of the co-evolution and function of molecules having roles in cellular regulation, or the aetiology of genetic diseases and how synonymous mutations can affect the encoded protein.

AKR1B10 (Aldo-keto reductase family 1 B10) promotes brain metastasis of lung cancer cells in a multi-organ microfluidic chip model

Brain metastasis (BM) is a leading cause of mortality in patients with non-small cell lung cancer (NSCLC). However, the molecular mechanisms underlying BM of NSCLC remain largely unknown because of the lack of models to accurately investigate such a dynamic and complex process. Here we developed a multi-organ microfluidic chip as a new methodological platform to study BM. The chip consisted of two bionic organ units - an upstream "lung" and a downstream "brain" characterized by a functional "blood-brain barrier (BBB)" structure, allowing real-time visual monitoring of the entire BM process, from the growth of primary tumor to its breaking through the BBB, and finally reaching the brain parenchyma. The chip was verified by lung cancer cell lines with differing metastatic abilities and then applied for the BM research where we first demonstrated that the protein expression of Aldo-keto reductase family 1 B10 (AKR1B10) was significantly elevated in lung cancer BM. Silencing AKR1B10 in brain metastatic tumor cells suppressed their extravasation through the BBB in the in vitro Transwell model, in our ex vivo microfluidic chip, as well as the in vivo model of brain metastasis in nude mice. Moreover, AKR1B10 downregulated the expression of matrix metalloproteinase (MMP)-2 and MMP-9 via MEK/ERK signaling in metastatic lung cancers. These data suggest that our multi-organ microfluidic chip is a practical alternative to study BM pathogenesis, and AKR1B10 is a diagnostic biomarker and a prospective therapeutic target for NSCLC BM. STATEMENT OF SIGNIFICANCE: Brain metastasis (BM) of non-small cell lung cancer (NSCLC) is a complex cascade, and in particular, the process of lung cancer cells penetrating the blood-brain barrier (BBB) is very unique. However, due to the lack of reliable models that can faithfully mimic the dynamic process of BBB breaking, its molecular mechanisms have not well elucidated so far. In addition, although Aldo-keto reductase family 1 B10 (AKR1B10) has been implicated to the tumor development of liver cancer and many other cancers, little is known on its roles in the BM. Here, we established a multi-organ microfluidic bionic chip platform to recapitulate the entire BM process, and applied it to the BM pathology research, especially BBB extravasation. By using the chip and traditional models synergistically, we first demonstrated that AKR1B10 was significantly elevated in lung cancer BM, and defined the value of AKR1B10 as a diagnostic serum biomarker for lung cancer patients suffering from BM. Further, we investigated the role and mechanisms of AKR1B10 in BM that it promotes the extravasation of cancer cells through the BBB.