Alternative splicing perturbation landscape identifies RNA binding proteins as potential therapeutic targets in cancer

Single-cell RNA-binding protein pattern-mediated molecular subtypes depict the hallmarks of the tumor microenvironment in bladder urothelial carcinoma

Objectives

Bladder carcinoma (BC) is a common malignancy of the urinary tract. As a new hallmark of cancer for drug therapy, RNA-binding proteins (RBPs) are key regulatory factors in alternative splicing events. This work is to uncover the relationship between BC and RBP in order to find drug targets in BC.

Methods

In this work, data from single-cell RNA-seq GSE1355337, PRJNA662018, and the TCGA-Bladder urothelial carcinoma (BLCA) cohorts are integrated to identify their relationships. A scoring system is constructed according to RBPs gene expression and patients’ survival. A network is constructed to analyze the alternative splicing events and RBP genes.

Results

A scoring system identified 321 RBPs significantly associated with the prognosis of patients. Subsequent typing of these RBP genes in two single-cell datasets demonstrated that most of the RBP genes had variable copy numbers. Three RBP clusters were identified. Using RBP genes as a signature in BC epithelial cells allows for differentiation between different grades of BC samples. The novel RBP genes-based subtype system reflects BC clinical staging. Notably, CellChat analysis revealed that the RBP genes-associated cell subtypes of T cells had extensive interactions with epithelial cells. Further analysis showed that the ligand-receptor pair MIF-CXCR4 mediated the communication between RBP-associated subtypes of BC epithelial cells and T cells.

Conclusions

Taken together, RBP genes are associated with BC progress and offer new indicators for precision medicine in BC.

Differential alternative splicing landscape identifies potentially functional RNA binding proteins in early embryonic development in mammals

Alternative splicing (AS) as one of the important post-transcriptional regulatory mechanisms has been poorly studied during embryogenesis. In this study, we comprehensively collected and analyzed the transcriptome data of early embryos from human and mouse. We found that AS plays an important role in this process and predicted candidate RNA binding protein (RBP) regulators that are associated with reproductive development. The predicted RBPs such as EIF4A3, MAK16, SRSF2, and UTP23 were found to be associated with reproductive disorders. By Smart-seq2 sequencing analysis, we identified 5445 aberrant alternative splicing events in Eif4a3-knockdown embryos. These events were preferentially associated with RNA processing. In conclusion, our work on the landscape and potential function of alternative splicing events will boost further investigation of detailed mechanisms and key factors regulating mammalian early embryo development and promote the inspiration of pharmaceutical approaches for disorders in this crucial biology process.

Splicing modulators impair DNA damage response and induce killing of cohesin-mutant MDS and AML

Splicing modulation is a promising treatment strategy pursued to date only in splicing factor-mutant cancers; however, its therapeutic potential is poorly understood outside of this context. Like splicing factors, genes encoding components of the cohesin complex are frequently mutated in cancer, including myelodysplastic syndromes (MDS) and secondary acute myeloid leukemia (AML), where they are associated with poor outcomes. Here, we showed that cohesin mutations are biomarkers of sensitivity to drugs targeting the splicing factor 3B subunit 1 (SF3B1) H3B-8800 and E-7107. We identified drug-induced alterations in splicing, and corresponding reduced gene expression, of a number of DNA repair genes, including BRCA1 and BRCA2, as the mechanism underlying this sensitivity in cell line models, primary patient samples and patient-derived xenograft (PDX) models of AML. We found that DNA damage repair genes are particularly sensitive to exon skipping induced by SF3B1 modulators due to their long length and large number of exons per transcript. Furthermore, we demonstrated that treatment of cohesin-mutant cells with SF3B1 modulators not only resulted in impaired DNA damage response and accumulation of DNA damage, but it sensitized cells to subsequent killing by poly(ADP-ribose) polymerase (PARP) inhibitors and chemotherapy and led to improved overall survival of PDX models of cohesin-mutant AML in vivo. Our findings expand the potential therapeutic benefits of SF3B1 splicing modulators to include cohesin-mutant MDS and AML.

Human papillomavirus integration transforms chromatin to drive oncogenesis

BACKGROUND

Human papillomavirus (HPV) drives almost all cervical cancers and up to 70% of head and neck cancers. Frequent integration into the host genome occurs predominantly in tumorigenic types of HPV. We hypothesize that changes in chromatin state at the location of integration can result in changes in gene expression that contribute to the tumorigenicity of HPV.

RESULTS

We find that viral integration events often occur along with changes in chromatin state and expression of genes near the integration site. We investigate whether introduction of new transcription factor binding sites due to HPV integration could invoke these changes. Some regions within the HPV genome, particularly the position of a conserved CTCF binding site, show enriched chromatin accessibility signal. ChIP-seq reveals that the conserved CTCF binding site within the HPV genome binds CTCF in 4 HPV^+ cancer cell lines. Significant changes in CTCF binding pattern and increases in chromatin accessibility occur exclusively within 100 kbp of HPV integration sites. The chromatin changes co-occur with out-sized changes in transcription and alternative splicing of local genes. Analysis of The Cancer Genome Atlas (TCGA) HPV⁺ tumors indicates that HPV integration upregulates genes which have significantly higher essentiality scores compared to randomly selected upregulated genes from the same tumors.

CONCLUSIONS

Our results suggest that introduction of a new CTCF binding site due to HPV integration reorganizes chromatin state and upregulates genes essential for tumor viability in some HPV⁺ tumors. These findings emphasize a newly recognized role of HPV integration in oncogenesis.

The internal ribosome entry site determines the neurotropic potential of enterovirus A71

The mechanisms underlying tissue-specific replication of enteroviruses are currently unclear. Although enterovirus A71 (EV-A71) and coxsackievirus A16 (CV-A16) are both common pathogens that cause hand-foot-mouth disease, they display quite different neurotropic properties. Herein, we characterized the role of the internal ribosomal entry site (IRES) in determining neurovirulence using an oral inoculation model of EV-A71. The receptor transgenic (hSCARB2-Tg) mice developed neurological symptoms after being infected with a mouse-adapted EV-A71 strain (MP4) via different administrative routes. Intragastric administration of the MP4 strain caused pathological changes such as neuronal loss and neuropil vacuolation, whereas replacing EV-A71 IRES with CV-A16 abolished the neuropathological phenotypes. The attenuated neurotropic potential of IRES-swapped EV-A71 was observed even in mice that received intraperitoneal and intracerebral inoculations. Fewer chimeric MP4 viral RNAs and proteins were detected in the mouse tissues, regardless of the inoculation route. Tissue-specific replication can be reflected in cell-based characterization. While chimeric MP4 virus replicated poorly in human intestinal C2BBe1 and neuroblastoma SH-SY5Y cells, its replication in susceptible rhabdomyosarcoma cells was not affected. Overall, our results demonstrated that the IRES determined the neurotropic potential of EV-A71 and CV-A16, emphasizing the importance of the IRES in tissue tropism, along with the host receptors.

Single-cell RNA binding protein regulatory network analyses reveal oncogenic HNRNPK-MYC signalling pathway in cancer

RNA-binding proteins (RBPs) are key players of gene expression and perturbations of RBP-RNA regulatory network have been observed in various cancer types. Here, we propose a computational method, RBPreg, to identify the RBP regulators by integration of single cell RNA-Seq (N = 233,591) and RBP binding data. Pan-cancer analyses suggest that RBP regulators exhibit cancer and cell specificity and perturbations of RBP regulatory network are involved in cancer hallmark-related functions. We prioritize an oncogenic RBP-HNRNPK, which is highly expressed in tumors and associated with poor prognosis of patients. Functional assays performed in cancer cells reveal that HNRNPK promotes cancer cell proliferation, migration, and invasion in vitro and in vivo. Mechanistic investigations further demonstrate that HNRNPK promotes tumorigenesis and progression by directly binding to MYC and perturbed the MYC targets pathway in lung cancer. Our results provide a valuable resource for characterizing RBP regulatory networks in cancer, yielding potential biomarkers for precision medicine.

Data analysis guidelines for single-cell RNA-seq in biomedical studies and clinical applications

The application of single-cell RNA sequencing (scRNA-seq) in biomedical research has advanced our understanding of the pathogenesis of disease and provided valuable insights into new diagnostic and therapeutic strategies. With the expansion of capacity for high-throughput scRNA-seq, including clinical samples, the analysis of these huge volumes of data has become a daunting prospect for researchers entering this field. Here, we review the workflow for typical scRNA-seq data analysis, covering raw data processing and quality control, basic data analysis applicable for almost all scRNA-seq data sets, and advanced data analysis that should be tailored to specific scientific questions. While summarizing the current methods for each analysis step, we also provide an online repository of software and wrapped-up scripts to support the implementation. Recommendations and caveats are pointed out for some specific analysis tasks and approaches. We hope this resource will be helpful to researchers engaging with scRNA-seq, in particular for emerging clinical applications.

Alternative splicing signatures of congenital heart disease and induced pluripotent stem cell-derived cardiomyocytes from congenital heart disease patients

Congenital heart disease (CHD) is the most serious congenital defect in newborns with higher mortality. Alternative splicing (AS) plays an essential role in numerous heart diseases. However, our understanding of the link between mRNA splicing and CHD in humans is limited. Here, we try to investigate the genome-wide AS events in CHD using bioinformatics methods. We collected available RNA-seq datasets of CHD-induced pluripotent stem cell-cardiomyocytes (iPSC-CMs) (including single ventricle disease [SVD] and tetralogy of Fallot [TOF]) and non-CHD from the Gene Expression Omnibus database. Then, we unprecedentedly performed AS profiles in CHD-iPSC-CMs and non-CHD-iPSC-CMs. The rMAPS was used to generate RNA-maps for the analysis of RNA-binding proteins' (RBPs) binding sites. We used StringTie to identify and quantify the transcripts from aligned RNA-Seq reads. A quantification matrix was generated with respect to different groups by extracting the transcripts per million values from StringTie outputs. Then, this matrix was used for correlation analysis between the expression level of RBP and AS level. Finally, we validated our AS results using RNA-seq data from CHD and non-CHD patient tissue samples. We identified CHD-related AS events using CHD-iPSC-CMs and CHD samples from patients. The results showed that functional enrichment of abnormal AS in SVD and TOF was transcription factor-related. Using rMAPS, RNA-binding proteins which regulated these AS were also determined, and RBP-AS regulatory network was constructed. Overall, we identified abnormal AS in CHD-iPSC-CMs and CHD samples from patients. We predicted AS regulators in SVD and TOF, respectively. At last, we concluded that AS played a key role in the pathogenesis of CHD.

Pathogenic Roles of RNA-Binding Proteins in Sarcomas

RNA-binding proteins (RBPs) are proteins that physically and functionally bind to RNA to regulate the RNA metabolism such as alternative splicing, polyadenylation, transport, maintenance of stability, localization, and translation. There is accumulating evidence that dysregulated RBPs play an essential role in the pathogenesis of malignant tumors including a variety of types of sarcomas. On the other hand, prognosis of patients with sarcoma, especially with sarcoma in advanced stages, is very poor, and almost no effective standard treatment has been established for most of types of sarcomas so far, highlighting the urgent need for identifying novel therapeutic targets based on the deep understanding of pathogenesis. Therefore, defining the network of interactions between RBPs and disease-related RNA targets will contribute to a better understanding of sarcomagenesis and identification of a novel therapeutic target for sarcomas.

TDP43 promotes stemness of breast cancer stem cells through CD44 variant splicing isoforms

Alternative splicing (AS) is a promising clinical target for cancer treatment at the post-transcriptional level. We previously identified a unique AS profile in triple-negative breast cancer (TNBC), which is regulated by the splicing regulator TAR DNA-binding protein-43 (TDP43), thus indicating the crucial role of TDP43 in heterogeneous TNBC. Cluster of differentiation 44 (CD44), a widely recognized marker for breast cancer stem cells (BCSCs), is extensively spliced into CD44 variant AS isoforms (CD44v) during the development of breast cancer. At present, however, the regulatory mechanism of CD44v is not fully understood. In the current study, we found that loss of TDP43 inhibits BCSC stemness by reducing the abundance of CD44v. In addition, serine-arginine-rich splicing factor 3 (SRSF3), another splicing factor and partner of TDP43, acts as an upstream regulator of TDP43 to maintain CD44v isoforms and thereafter BCSC stemness. Mechanistically, SRSF3 stabilizes the mRNA of TDP43 by inhibiting nonsense-mediated decay (NMD). These findings illustrate the important role of complicated regulatory networks formed by splicing factors in TNBC progression, thus providing potential therapeutic targets from an AS perspective.

Shedding light on the hidden human proteome expands immunopeptidome in cancer

Unrestrained cellular growth and immune escape of a tumor are associated with the incidental errors of the genome and transcriptome. Advances in next-generation sequencing have identified thousands of genomic and transcriptomic aberrations that generate variant peptides that assemble the hidden proteome, further expanding the immunopeptidome. Emerging next-generation sequencing technologies and a number of computational methods estimated the abundance of immune infiltration from bulk transcriptome have advanced our understanding of tumor microenvironments. Here, we will characterize several major types of tumor-specific antigens arising from single-nucleotide variants, insertions and deletions, gene fusion, alternative splicing, RNA editing and non-coding RNAs. Finally, we summarize the current state-of-the-art computational and experimental approaches or resources and provide an integrative pipeline for the identification of candidate tumor antigens. Together, the systematic investigation of the hidden proteome in cancer will help facilitate the development of effective and durable immunotherapy targets for cancer.

ImmReg: the regulon atlas of immune-related pathways across cancer types

Immune system gene regulation perturbation has been found to be a major cause of the development of various types of cancer. Numbers of mechanisms contribute to gene expression regulation, thus, systematically identification of potential regulons of immune-related pathways is critical to cancer immunotherapy. Here, we comprehensively chart the landscape of transcription factors, microRNAs, RNA binding proteins and long noncoding RNAs regulation in 17 immune-related pathways across 33 cancers. The potential immunology regulons are likely to exhibit higher expressions in immune cells, show expression perturbations in cancer, and are significantly correlated with immune cell infiltrations. We also identify a panel of clinically relevant immunology regulons across cancers. Moreover, the regulon atlas of immune-related pathways helps prioritizing cancer-related genes (i.e. ETV7, miR-146a-5p, ZFP36 and HCP5). We further identified two molecular subtypes of glioma (cold and hot tumour phenotypes), which were characterized by differences in immune cell infiltrations, expression of checkpoints, and prognosis. Finally, we developed a user-friendly resource, ImmReg (http://bio-bigdata.hrbmu.edu.cn/ImmReg/), with multiple modules to visualize, browse, and download immunology regulation. Our study provides a comprehensive landscape of immunology regulons, which will shed light on future development of RNA-based cancer immunotherapies.