RBPmap: A Tool for Mapping and Predicting the Binding Sites of RNA-Binding Proteins Considering the Motif Environment

Demonstration of the pathogenicity of a common non-exomic mutation in ABCA4 using iPSC-derived retinal organoids and retrospective clinical data

Mutations in ABCA4 are the most common cause of Mendelian retinal disease. Clinical evaluation of this gene is challenging because of its extreme allelic diversity, the large fraction of non-exomic mutations, and the wide range of associated disease. We used patient-derived retinal organoids as well as DNA samples and clinical data from a large cohort of patients with ABCA4-associated retinal disease to investigate the pathogenicity of a variant in ABCA4 (IVS30 + 1321 A>G) that occurs heterozygously in 2% of Europeans. We found that this variant causes mis-splicing of the gene in photoreceptor cells such that the resulting protein contains 36 incorrect amino acids followed by a premature stop. We also investigated the phenotype of 10 patients with compound genotypes that included this mutation. Their median age of first vision loss was 39 years, which is in the mildest quintile of a large cohort of patients with ABCA4 disease. We conclude that the IVS30 + 1321 A>G variant can cause disease when paired with a sufficiently deleterious opposing allele in a sufficiently permissive genetic background.

Transcriptomic analysis reveals regulation of adipogenesis via long non-coding RNA, alternative splicing, and alternative polyadenylation

Obesity is characterized by dysregulated adipogenesis that leads to increased number and/or size of adipocytes. Understanding the molecular mechanisms governing adipogenesis is therefore key to designing therapeutic interventions against obesity. In our study, we analyzed 3'-end sequencing data that we generated from human preadipocytes and adipocytes, as well as previously published RNA-seq datasets, to elucidate mechanisms of regulation via long non-coding RNA (lncRNA), alternative splicing (AS) and alternative polyadenylation (APA). We discovered lncRNAs that have not been previously characterized but may be key regulators of white adipogenesis. We also detected 100 AS events and, using motif enrichment analysis, identified RNA binding proteins (RBPs) that could mediate exon skipping-the most prevalent AS event. In addition, we show that usage of alternative poly(A) sites in introns or 3'-UTRs of key adipogenesis genes leads to isoform diversity, which can have significant biological consequences on differentiation efficiency. We also identified RBPs that may modulate APA and defined how 3'-UTR APA can regulate gene expression through gain or loss of specific microRNA binding sites. Taken together, our bioinformatics-based analysis reveals potential therapeutic avenues for obesity through manipulation of lncRNA levels and the profile of mRNA isoforms via alternative splicing and polyadenylation.

Mechanistic insights into the basis of widespread RNA localization

The importance of subcellular mRNA localization is well established, but the underlying mechanisms mostly remain an enigma. Early studies suggested that specific mRNA sequences recruit RNA-binding proteins (RBPs) to regulate mRNA localization. However, despite the observation of thousands of localized mRNAs, only a handful of these sequences and RBPs have been identified. This suggests the existence of alternative, and possibly predominant, mechanisms for mRNA localization. Here I re-examine currently described mRNA localization mechanisms and explore alternative models that could account for its widespread occurrence.

The leader RNA of SARS-CoV-2 sequesters polypyrimidine tract binding protein (PTBP1) and influences pre-mRNA splicing in infected cells

The large amount of viral RNA produced during infections has the potential to interact with and effectively sequester cellular RNA binding proteins, thereby influencing aspects of post-transcriptional gene regulation in the infected cell. Here we demonstrate that the abundant 5' leader RNA region of SARS-CoV-2 viral RNAs can interact with the cellular polypyrimidine tract binding protein (PTBP1). Interestingly, the effect of a knockdown of PTBP1 protein on cellular gene expression is also mimicked during SARS-CoV-2 infection, suggesting that this protein may be functionally sequestered by viral RNAs. Consistent with this model, the alternative splicing of mRNAs that is normally controlled by PTBP1 is dysregulated during SARS-CoV-2 infection. Collectively, these data suggest that the SARS-CoV-2 leader RNA sequesters the cellular PTBP1 protein during infection, resulting in significant impacts on the RNA biology of the host cell. These alterations in post-transcriptional gene regulation may play a role in SARS-CoV-2 mediated molecular pathogenesis.

CDK11 requires a critical activator SAP30BP to regulate pre-mRNA splicing

CDK11 is an emerging druggable target for cancer therapy due to its prevalent roles in phosphorylating critical transcription and splicing factors and in facilitating cell cycle progression in cancer cells. Like other cyclin-dependent kinases, CDK11 requires its cognate cyclin, cyclin L1 or cyclin L2, for activation. However, little is known about how CDK11 activities might be modulated by other regulators. In this study, we show that CDK11 forms a tight complex with cyclins L1/L2 and SAP30BP, the latter of which is a poorly characterized factor. Acute degradation of SAP30BP mirrors that of CDK11 in causing widespread and strong defects in pre-mRNA splicing. Furthermore, we demonstrate that SAP30BP facilitates CDK11 kinase activities in vitro and in vivo, through ensuring the stabilities and the assembly of cyclins L1/L2 with CDK11. Together, these findings uncover SAP30BP as a critical CDK11 activator that regulates global pre-mRNA splicing.

Long non-coding RNA generated from CDKN1A gene by alternative polyadenylation regulates p21 expression during DNA damage response

Alternative Polyadenylation (APA) is an emerging mechanism for dynamic changes in gene expression. Previously, we described widespread APA occurrence in introns during the DNA damage response (DDR). Here, we show that a DDR-activated APA event occurs in the first intron of CDKN1A, inducing an alternate last exon-containing lncRNA. We named this lncRNA SPUD (Selective Polyadenylation Upon DNA Damage). SPUD localizes to polysomes in the cytoplasm and is detectable as multiple isoforms in available high-throughput studies. SPUD has low abundance compared to the CDKN1A full-length isoform under non-stress conditions, and SPUD is induced in cancer and normal cells under a variety of DNA damaging conditions in part through p53. The RNA binding protein HuR binds to and promotes the stability of SPUD precursor RNA. SPUD induction increases p21 protein, but not mRNA levels, affecting p21 functions in cell-cycle, CDK2 expression and cell growth. Like CDKN1A full-length isoform, SPUD can bind two competitive p21 translational regulators, the inhibitor calreticulin and the activator CUGBP1; SPUD alters their association with CDKN1A full-length in a DDR-dependent manner, promoting CDKN1A translation. Together, these results show a new regulatory mechanism by which a lncRNA controls p21 expression post-transcriptionally, highlighting lncRNA relevance in DDR progression and cell-cycle.

Long-read sequencing reveals the landscape of aberrant alternative splicing and novel therapeutic target in colorectal cancer

BACKGROUND

Alternative splicing complexity plays a vital role in carcinogenesis and cancer progression. Improved understanding of novel splicing events and the underlying regulatory mechanisms may contribute new insights into developing new therapeutic strategies for colorectal cancer (CRC).

METHODS

Here, we combined long-read sequencing technology with short-read RNA-seq methods to investigate the transcriptome complexity in CRC. By using experiment assays, we explored the function of newly identified splicing isoform TIMP1 Δ4-5. Moreover, a CRISPR/dCasRx-based strategy to induce the TIMP1 exon 4-5 exclusion was introduced to inhibit neoplasm growth.

RESULTS

A total of 90,703 transcripts were identified, of which > 62% were novel compared with current transcriptome annotations. These novel transcripts were more likely to be sample specific, expressed at relatively lower levels with more exons, and oncogenes displayed a characteristic to generate more transcripts in CRC. Clinical outcome data analysis showed that 1472 differentially expressed alternative splicing events (DEAS) were tightly associated with CRC patients' prognosis, and many novel isoforms were likely to be important determinants for patient survival. Among these, newly identified splicing isoform TIMP1 Δ4-5 was significantly downregulated in CRC. Further in vitro and in vivo assays demonstrated that ectopic expression of TIMP1 Δ4-5 significantly suppresses tumor cell growth and metastasis. Serine/arginine-rich splicing factor 1 (SRSF1) acts as a onco-splicing regulator through sustaining the inclusion of TIMP1 exon 4-5. Furthermore, CRISPR/dCasRx-based strategies designed to induce TIMP1 exon 4-5 exclusion have the potential to restrain the CRC growth.

CONCLUSIONS

This data provides a rich resource for deeper studies of gastrointestinal malignancies. Newly identified splicing isoform TIMP1 Δ4-5 plays an important role in mediating CRC progression and may be a potential therapy target in CRC.

Massively parallel identification of mRNA localization elements in primary cortical neurons

Cells adopt highly polarized shapes and form distinct subcellular compartments in many cases due to the localization of many mRNAs to specific areas, where they are translated into proteins with local functions. This mRNA localization is mediated by specific cis-regulatory elements in mRNAs, commonly called 'zipcodes'. Although there are hundreds of localized mRNAs, only a few zipcodes have been characterized. Here we describe a novel neuronal zipcode identification protocol (N-zip) that can identify zipcodes across hundreds of 3' untranslated regions. This approach combines a method of separating the principal subcellular compartments of neurons-cell bodies and neurites-with a massively parallel reporter assay. N-zip identifies the let-7 binding site and (AU)_n motif as de novo zipcodes in mouse primary cortical neurons. Our analysis also provides, to our knowledge, the first demonstration of an miRNA affecting mRNA localization and suggests a strategy for detecting many more zipcodes.

MAPKAPK2-centric transcriptome profiling reveals its major role in governing molecular crosstalk of IGFBP2, MUC4, and PRKAR2B during HNSCC pathogenesis

Transcriptome analysis of head and neck squamous cell carcinoma (HNSCC) has been pivotal to comprehending the convoluted biology of HNSCC tumors. MAPKAPK2 or MK2 is a critical modulator of the mRNA turnover of crucial genes involved in HNSCC progression. However, MK2-centric transcriptome profiles of tumors are not well known. This study delves into HNSCC progression with MK2 at the nexus to delineate the biological relevance and intricate crosstalk of MK2 in the tumor milieu. We performed next-generation sequencing-based transcriptome profiling of HNSCC cells and xenograft tumors to ascertain mRNA expression profiles in MK2-wild type and MK2-knockdown conditions. The findings were validated using gene expression assays, immunohistochemistry, and transcript turnover studies. Here, we identified a pool of crucial MK2-regulated candidate genes by annotation and differential gene expression analyses. Regulatory network and pathway enrichment revealed their significance and involvement in the HNSCC pathogenesis. Additionally, 3'-UTR-based filtering recognized important MK2-regulated downstream target genes and validated them by nCounter gene expression assays. Finally, immunohistochemistry and transcript stability studies revealed the putative role of MK2 in regulating the transcript turnover of IGFBP2, MUC4, and PRKAR2B in HNSCC. Conclusively, MK2-regulated candidate genes were identified in this study, and their plausible involvement in HNSCC pathogenesis was elucidated. These genes possess investigative values as targets for diagnosis and therapeutic interventions for HNSCC.

Long Non-Coding RNA Generated from CDKN1A Gene by Alternative Polyadenylation Regulates p21 Expression during DNA Damage Response

Alternative Polyadenylation (APA) is an emerging mechanism for dynamic changes in gene expression. Previously, we described widespread APA occurrence in introns during the DNA damage response (DDR). Here, we show that a DNA damage activated APA event occurs in the first intron of CDKN1A , inducing an alternate last exon (ALE)-containing lncRNA. We named this lncRNA SPUD (Selective Polyadenylation Upon Damage). SPUD localizes to polysomes in the cytoplasm and is detectable as multiple isoforms in available high throughput studies. SPUD has low abundance compared to the CDKN1A full-length isoform and is induced in cancer and normal cells under a variety of DNA damaging conditions in part through p53 transcriptional activation. RNA binding protein (RBP) HuR and the transcriptional repressor CTCF regulate SPUD levels. SPUD induction increases p21 protein, but not CDKN1A full-length levels, affecting p21 functions in cell-cycle, CDK2 expression, and cell viability. Like CDKN1A full-length isoform, SPUD can bind two competitive p21 translational regulators, the inhibitor calreticulin and the activator CUGBP1; SPUD can change their association with CDKN1A full-length in a DDR-dependent manner. Together, these results show a new regulatory mechanism by which a lncRNA controls p21 expression post-transcriptionally, highlighting lncRNA relevance in DDR progression and cellcycle.