RNA-Binding Proteins: Splicing Factors and Disease

Understanding the genetics of sex determination in insects and its relevance to genetic pest management

Sex determination pathways regulate male and female-specific development and differentiation and offer potential targets for genetic pest management methods. Insect sex determination pathways are comprised of primary signals, relay genes and terminal genes. Primary signals of coleopteran, dipteran, hymenopteran and lepidopteran species are highly diverse and regulate the sex-specific splicing of relay genes based on the primary signal dosage, amino acid composition or the interaction with paternally inherited genes. In coleopterans, hymenopterans and some dipterans, relay genes are Transformer orthologs from the serine-arginine protein family that regulate sex-specific splicing of the terminal genes. Alternative genes regulate the splicing of the terminal genes in dipterans that lack Transformer orthologs and lepidopterans. Doublesex and Fruitless orthologs are the terminal genes. Doublesex and Fruitless orthologs are highly conserved zinc-finger proteins that regulate the expression of downstream proteins influencing physical traits and courtship behaviours in a sex-specific manner. Genetic pest management methods can use different mechanisms to exploit or disrupt female-specific regions of different sex determination genes. Female-specific regions of sex determination genes can be exploited to produce a lethal gene only in females or disrupted to impede female development or fertility. Reducing the number of fertile females in pest populations creates a male-biased sex ratio and eventually leads to the local elimination of the pest population. Knowledge on the genetic basis of sex determination is important to enable these sex determination pathways to be exploited for genetic pest management.

Alternative splicing categorizes organ development by stage and reveals unique human splicing variants linked to neuromuscular disorders

Alternative splicing (AS) diversifies protein expression and contributes to species-specific differences in organ development. Here, we focused on stage-specific splicing variants and their correlation with disease in human compared to mouse during brain and heart development. Temporal transcriptomic analysis revealed that splicing factors (SFs) can accurately classify organ developmental stages, and 5 SFs were identified specifically upregulated in human during organogenesis. Additionally, inter-stage splicing variations were identified across analogous human and mouse developmental stages. Developmentally dynamic alternative splicing genes (devASGs) were enriched in various neurodevelopmental disorders in both species, with the most significant changes observed in human newborn brain and 16 weeks post-conception heart. Intriguingly, diseases specifically enriched in humans were primarily associated with neuro-muscular dysfunction, and human-specific neuromuscular devASGs were linked to mannose glycosylation and ciliary motility. These findings highlight the significance of SFs and AS events in organogenesis, and inform the selection of appropriate models for translational research.

LincRNA-ASAO promotes dental pulp repair through interacting with PTBP1 to increase ALPL alternative splicing

BACKGROUND

Alternative splicing not only expands the genetic encoding of genes but also determines cellular activities. This study aimed to elucidate the regulation mechanism and biological functions of lincRNA-ASAO in the process of odontogenesis-related genes alternative splicing mediated odontogenic differentiation of hDPSCs.

METHODS

RACE, RNA-seq, FISH and bioinformatics techniques were used to identify novel lincRNA-ASAO. ALP staining, alizarin red staining, qRT-PCR and western blot were used to identify the role of lincRNA-ASAO in regulating the odontoblast differentiation of hDPSCs. The binding protein PTBP1 of lincRNA-ASAO was screened by RNA-Pulldown, protein profiling and bioinformatics. The target gene ALPL of lincRNA-ASAO/PTBP1 was identified by RNA-seq, bioinformatics technology and DNA agarose gel electrophoresis. FISH, IF, PAR-CLIP and bioinformatics techniques were used to determine the roles of lincRNA-ASAO, PTBP1 and ALPL pre-mRNA in the odontoblast differentiation of hDPSCs.

RESULTS

We identified a novel lincRNA-ASAO that could promote the odontogenic differentiation of human Dental Pulp Stem Cells (hDPSCs). And, the interaction between lincRNA-ASAO and alternative splicing factor PTBP1 promoted the odontoblast differentiation of hDPSCs. In addition, lincRNA-ASAO forms duplexes with ALPL pre-mRNA, targeting PTBP1 to exonic splicing silencer (ESS) of ALPL and regulating exon 2 skipping. Notably, lincRNA-ASAO/PTBP1 regulated ALPL production to increase the type 2 splice variant, which promoted the odontoblast differentiation of hDPSCs.

CONCLUSIONS

We have identified the novel lincRNA-ASAO, which can promote the odontoblast differentiation of hDPSCs. The mechanism study found that lincRNA-ASAO/PTBP1 mediated the exon 2 skipping of ALPL pre-mRNA, resulting in the type 2 splice variant of ALPL. Our results enrich the understanding of lncRNAs and alternative splicing in regulating the odontoblast differentiation of hDPSCs, and provide clues to improve the clinical therapeutic potential of hDPSCs for dental pulp restoration.

A compendium of Amplification-Related Gain Of Sensitivity genes in human cancer

While the effect of amplification-induced oncogene expression in cancer is known, the impact of copy-number gains on "bystander" genes is less understood. We create a comprehensive map of dosage compensation in cancer by integrating expression and copy number profiles from over 8000 tumors in The Cancer Genome Atlas and cell lines from the Cancer Cell Line Encyclopedia. Additionally, we analyze 17 cancer open reading frame screens to identify genes toxic to cancer cells when overexpressed. Combining these approaches, we propose a class of 'Amplification-Related Gain Of Sensitivity' (ARGOS) genes located in commonly amplified regions, yet expressed at lower levels than expected by their copy number, and toxic when overexpressed. We validate RBM14 as an ARGOS gene in lung and breast cancer cells, and suggest a toxicity mechanism involving altered DNA damage response and STING signaling. We additionally observe increased patient survival in a radiation-treated cancer cohort with RBM14 amplification.

OpenASO: RNA Rescue - designing splice-modulating antisense oligonucleotides through community science

Splice-modulating antisense oligonucleotides (ASOs) are precision RNA-based drugs that are becoming an established modality to treat human disease. Previously, we reported the discovery of ASOs that target a novel, putative intronic RNA structure to rescue splicing of multiple pathogenic variants of F8 exon 16 that cause hemophilia A. However, the conventional approach to discovering splice-modulating ASOs is both laborious and expensive. Here, we describe an alternative paradigm that integrates data-driven RNA structure prediction and community science to discover splice-modulating ASOs. Using a splicing-deficient pathogenic variant of F8 exon 16 as a model, we show that 25% of the top-scoring molecules designed in the Eterna OpenASO challenge have a statistically significant impact on enhancing exon 16 splicing. Additionally, we show that a distinct combination of ASOs designed by Eterna players can additively enhance the inclusion of the splicing-deficient exon 16 variant. Together, our data suggests that crowdsourcing designs from a community of citizen scientists may accelerate the discovery of splice-modulating ASOs with potential to treat human disease.

Global RNA Interaction and Transcriptome Profiles Demonstrate the Potential Anti-Oncogenic Targets and Pathways of RBM6 in HeLa Cells

BACKGROUND

The fate and functions of RNAs are coordinately regulated by RNA-binding proteins (RBPs), which are often dysregulated in various cancers. Known as a splicing regulator, RNA-binding motif protein 6 (RBM6) harbors tumor-suppressor activity in many cancers; however, there is a lack of research on the molecular targets and regulatory mechanisms of RBM6.

METHODS

In this study, we constructed an RBM6 knock-down (shRBM6) model in the HeLa cell line to investigate its functions and molecular targets. Then we applied improved RNA immunoprecipitation coupled with sequencing (iRIP-seq) and whole transcriptome sequencing approaches to investigate the potential role and RNA targets of RBM6.

RESULTS

Using The Cancer Genome Atlas dataset, we found that higher expression of RBM6 is associated with a better prognosis in many cancer types. In addition, we found that RBM6 knockdown promoted cell proliferation and inhibited apoptosis, demonstrating that RBM6 may act as an anti-oncogenic protein in cancer cells. RBM6 can regulate the alternative splicing (AS) of genes involved in DNA damage response, proliferation, and apoptosis-associated pathways. Meanwhile, RBM6 knockdown activated type I interferon signaling pathways and inhibited the expression of genes involved in the cell cycle, cellular responses to DNA damage, and DNA repair pathways. The differentially expressed genes (DEGs) by shRBM6 and their involved pathways were likely regulated by the transcription factors undergoing aberrant AS by RBM6 knockdown. For iRIP-seq analysis, we found that RBM6 could interact with a large number of mRNAs, with a tendency for binding motifs GGCGAUG and CUCU. RBM6 bound to the mRNA of cell proliferation- and apoptosis-associated genes with dysregulated AS after RBM6 knockdown.

CONCLUSIONS

In summary, our study highlights the important role of RBM6, as well as the downstream targets and regulated pathways, suggesting the potential regulatory mechanisms of RBM6 in the development of cancer.

SexAnnoDB, a knowledgebase of sex-specific regulations from multi-omics data of human cancers

BACKGROUND

Sexual differences across molecular levels profoundly impact cancer biology and outcomes. Patient gender significantly influences drug responses, with divergent reactions between men and women to the same drugs. Despite databases on sex differences in human tissues, understanding regulations of sex disparities in cancer is limited. These resources lack detailed mechanistic studies on sex-biased molecules.

METHODS

In this study, we conducted a comprehensive examination of molecular distinctions and regulatory networks across 27 cancer types, delving into sex-biased effects. Our analyses encompassed sex-biased competitive endogenous RNA networks, regulatory networks involving sex-biased RNA binding protein-exon skipping events, sex-biased transcription factor-gene regulatory networks, as well as sex-biased expression quantitative trait loci, sex-biased expression quantitative trait methylation, sex-biased splicing quantitative trait loci, and the identification of sex-biased cancer therapeutic drug target genes. All findings from these analyses are accessible on SexAnnoDB ( https://ccsm.uth.edu/SexAnnoDB/ ).

RESULTS

From these analyses, we defined 126 cancer therapeutic target sex-associated genes. Among them, 9 genes showed sex-biased at both the mRNA and protein levels. Specifically, S100A9 was the target of five drugs, of which calcium has been approved by the FDA for the treatment of colon and rectal cancers. Transcription factor (TF)-gene regulatory network analysis suggested that four TFs in the SARC male group targeted S100A9 and upregulated the expression of S100A9 in these patients. Promoter region methylation status was only associated with S100A9 expression in KIRP female patients. Hypermethylation inhibited S100A9 expression and was responsible for the downregulation of S100A9 in these female patients.

CONCLUSIONS

Comprehensive network and association analyses indicated that the sex differences at the transcriptome level were partially the result of corresponding sex-biased epigenetic and genetic molecules. Overall, SexAnnoDB offers a discipline-specific search platform that could potentially assist basic experimental researchers or physicians in developing personalized treatment plans.

Alteration of RNA m6A methylation mediates aberrant RNA binding protein expression and alternative splicing in condyloma acuminatum

Background

Condyloma acuminatum (CA) is caused by low-risk human papillomavirus, and is characterized by high recurrence after treatment. The RNA modification N6-methyladenosine (m6A) plays an important role during diverse viral infections, including high-risk HPV infection in cervical cancer. However, it is unclear whether low-risk HPV infection changes the RNA m6A methylation in CA.

Methods

High-throughputm6A-sequencing was performed to profile the transcriptome-wide mRNA modifications of CA tissues infected by LR-HPVs and the paired normal tissues from CA patients. We further investigated the regulation of alternative splicing by RNA binding proteins (RBPs) with altered m6A modification and constructed a regulatory network among these RBPs, regulated alternative splicing events (RASEs) and regulated alternative splicing genes (RASGs) in CA.

Results

The results show that the m6A level in CA tissues differed from that in the paired controls. Furthermore, cell cycle- and cell adhesion- associated genes with m6A modification were differentially expressed in CA tissues compared to the paired controls. In particular, seven RNA binding protein genes with specific m6A methylated sites, showed a higher or lower expression at the mRNA level in CA tissues than in the paired normal tissues. In addition, these differentially expressed RNA binding protein genes would regulate the alternative splicing pattern of apoptotic process genes in CA tissue.

Conclusions

Our study reveals a sophisticated m6A modification profile in CA tissue that affects the response of host cells to HPV infection, and provides cues for the further exploration of the roles of m6A and the development of a novel treatment strategy for CA.

MAPP unravels frequent co-regulation of splicing and polyadenylation by RNA-binding proteins and their dysregulation in cancer

Maturation of eukaryotic pre-mRNAs via splicing and polyadenylation is modulated across cell types and conditions by a variety of RNA-binding proteins (RBPs). Although there exist over 1,500 RBPs in human cells, their binding motifs and functions still remain to be elucidated, especially in the complex environment of tissues and in the context of diseases. To overcome the lack of methods for the systematic and automated detection of sequence motif-guided pre-mRNA processing regulation from RNA sequencing (RNA-Seq) data we have developed MAPP (Motif Activity on Pre-mRNA Processing). Applying MAPP to RBP knock-down experiments reveals that many RBPs regulate both splicing and polyadenylation of nascent transcripts by acting on similar sequence motifs. MAPP not only infers these sequence motifs, but also unravels the position-dependent impact of the RBPs on pre-mRNA processing. Interestingly, all investigated RBPs that act on both splicing and 3' end processing exhibit a consistently repressive or activating effect on both processes, providing a first glimpse on the underlying mechanism. Applying MAPP to normal and malignant brain tissue samples unveils that the motifs bound by the PTBP1 and RBFOX RBPs coordinately drive the oncogenic splicing program active in glioblastomas demonstrating that MAPP paves the way for characterizing pre-mRNA processing regulators under physiological and pathological conditions.

Emerging roles of RNA-binding proteins in fatty liver disease

A rampant and urgent global health issue of the 21st century is the emergence and progression of fatty liver disease (FLD), including alcoholic fatty liver disease and the more heterogenous metabolism-associated (or non-alcoholic) fatty liver disease (MAFLD/NAFLD) phenotypes. These conditions manifest as disease spectra, progressing from benign hepatic steatosis to symptomatic steatohepatitis, cirrhosis, and, ultimately, hepatocellular carcinoma. With numerous intricately regulated molecular pathways implicated in its pathophysiology, recent data have emphasized the critical roles of RNA-binding proteins (RBPs) in the onset and development of FLD. They regulate gene transcription and post-transcriptional processes, including pre-mRNA splicing, capping, and polyadenylation, as well as mature mRNA transport, stability, and translation. RBP dysfunction at every point along the mRNA life cycle has been associated with altered lipid metabolism and cellular stress response, resulting in hepatic inflammation and fibrosis. Here, we discuss the current understanding of the role of RBPs in the post-transcriptional processes associated with FLD and highlight the possible and emerging therapeutic strategies leveraging RBP function for FLD treatment. This article is categorized under: RNA in Disease and Development > RNA in Disease.

RNA-Binding Proteins: A Role in Neurotoxicity?

Despite sustained efforts to treat neurodegenerative diseases, little is known at the molecular level to understand and generate novel therapeutic approaches for these malignancies. Therefore, it is not surprising that neurogenerative diseases are among the leading causes of death in the aged population. Neurons require sophisticated cellular mechanisms to maintain proper protein homeostasis. These cells are generally sensitive to loss of gene expression control at the post-transcriptional level. Post-translational control responds to signals that can arise from intracellular processes or environmental factors that can be regulated through RNA-binding proteins. These proteins recognize RNA through one or more RNA-binding domains and form ribonucleoproteins that are critically involved in the regulation of post-transcriptional processes from splicing to the regulation of association of the translation machinery allowing a relatively rapid and precise modulation of the transcriptome. Neurotoxicity is the result of the biological, chemical, or physical interaction of agents with an adverse effect on the structure and function of the central nervous system. The disruption of the proper levels or function of RBPs in neurons and glial cells triggers neurotoxic events that are linked to neurodegenerative diseases such as spinal muscular atrophy (SMA), amyotrophic lateral sclerosis (ALS), fragile X syndrome (FXS), and frontotemporal dementia (FTD) among many others. The connection between RBPs and neurodegenerative diseases opens a new landscape for potentially novel therapeutic targets for the intervention of these neurodegenerative pathologies. In this contribution, a summary of the recent findings of the molecular mechanisms involved in the plausible role of RBPs in RNA processing in neurodegenerative disease is discussed.

In silico identification, characterization, and expression analysis of RNA recognition motif (RRM) containing RNA-binding proteins in Aedes aegypti

RNA-binding proteins (RBPs) are the proteins that bind RNAs and regulate their functioning. RBPs in mosquitoes are gaining attention due to their ability to bind flaviviruses and regulate their replication and transmission. Despite their relevance, RBPs in mosquitoes are not explored much. In this study, we screened the whole genome of Aedes aegypti, the primary vector of several pathogenic viruses, and identified the proteins containing RNA recognition motif (RRM), the most abundant protein domain in eukaryotes. Using several in silico strategies, a total of 135 RRM-containing RBPs were identified in Ae. aegypti. The proteins were characterized based on their available annotations and the sequence similarity with Drosophila melanogaster. Ae. aegypti RRM-containing RBPs included serine/arginine-rich (SR) proteins, polyadenylate-binding proteins (PABP), heteronuclear ribonucleoproteins (hnRNP), small nuclear ribonucleoproteins (snRNP), splicing factors, eukaryotic initiation factors, transformers, and nucleolysins. Phylogenetic analysis revealed that the proteins and the domain organization are conserved among Ae. aegypti, Bombyx mori, and Drosophila melanogaster. However, the gene length and the intron-exon organization varied across the insect species. Expression analysis of the genes encoding RBPs using publicly available RNA sequencing data for different developmental time points of the mosquito life cycle starting from the ovary and eggs up to the adults revealed stage-specific expression with several genes preferentially expressed in early embryonic stages and blood-fed female ovaries. This is the first database for the Ae. aegypti RBPs that can serve as the reference base for future investigations. Stage-specific genes can be further explored to determine their role in mosquito growth and development with a focus on developing novel mosquito control strategies.

Nucleocytoplasmic mRNA redistribution accompanies RNA binding protein mislocalization in ALS motor neurons and is restored by VCP ATPase inhibition

Amyotrophic lateral sclerosis (ALS) is characterized by nucleocytoplasmic mislocalization of the RNA-binding protein (RBP) TDP-43. However, emerging evidence suggests more widespread mRNA and protein mislocalization. Here, we employed nucleocytoplasmic fractionation, RNA sequencing, and mass spectrometry to investigate the localization of mRNA and protein in induced pluripotent stem cell-derived motor neurons (iPSMNs) from ALS patients with TARDBP and VCP mutations. ALS mutant iPSMNs exhibited extensive nucleocytoplasmic mRNA redistribution, RBP mislocalization, and splicing alterations. Mislocalized proteins exhibited a greater affinity for redistributed transcripts, suggesting a link between RBP mislocalization and mRNA redistribution. Notably, treatment with ML240, a VCP ATPase inhibitor, partially restored mRNA and protein localization in ALS mutant iPSMNs. ML240 induced changes in the VCP interactome and lysosomal localization and reduced oxidative stress and DNA damage. These findings emphasize the link between RBP mislocalization and mRNA redistribution in ALS motor neurons and highlight the therapeutic potential of VCP inhibition.

Development of Engineered-U1 snRNA Therapies: Current Status

Splicing of pre-mRNA is a crucial regulatory stage in the pathway of gene expression. The majority of human genes that encode proteins undergo alternative pre-mRNA splicing and mutations that affect splicing are more prevalent than previously thought. Targeting aberrant RNA(s) may thus provide an opportunity to correct faulty splicing and potentially treat numerous genetic disorders. To that purpose, the use of engineered U1 snRNA (either modified U1 snRNAs or exon-specific U1s-ExSpeU1s) has been applied as a potentially therapeutic strategy to correct splicing mutations, particularly those affecting the 5' splice-site (5'ss). Here we review and summarize a vast panoply of studies that used either modified U1 snRNAs or ExSpeU1s to mediate gene therapeutic correction of splicing defects underlying a considerable number of genetic diseases. We also focus on the pre-clinical validation of these therapeutic approaches both in vitro and in vivo, and summarize the main obstacles that need to be overcome to allow for their successful translation to clinic practice in the future.

Towards in silico CLIP-seq: predicting protein-RNA interaction via sequence-to-signal learning

We present RBPNet, a novel deep learning method, which predicts CLIP-seq crosslink count distribution from RNA sequence at single-nucleotide resolution. By training on up to a million regions, RBPNet achieves high generalization on eCLIP, iCLIP and miCLIP assays, outperforming state-of-the-art classifiers. RBPNet performs bias correction by modeling the raw signal as a mixture of the protein-specific and background signal. Through model interrogation via Integrated Gradients, RBPNet identifies predictive sub-sequences that correspond to known and novel binding motifs and enables variant-impact scoring via in silico mutagenesis. Together, RBPNet improves imputation of protein-RNA interactions, as well as mechanistic interpretation of predictions.

Multi-omics data analysis reveals the biological implications of alternative splicing events in lung adenocarcinoma

Cancer is characterized by the dysregulation of alternative splicing (AS). However, the comprehensive regulatory mechanisms of AS in lung adenocarcinoma (LUAD) are poorly understood. Here, we displayed the AS landscape in LUAD based on the integrated analyses of LUAD's multi-omics data. We identified 13,995 AS events in 6309 genes as differentially expressed alternative splicing events (DEASEs) mainly covering protein-coding genes. These DEASEs were strongly linked to "cancer hallmarks", such as apoptosis, DNA repair, cell cycle, cell proliferation, angiogenesis, immune response, generation of precursor metabolites and energy, p53 signaling pathway and PI3K-AKT signaling pathway. We further built a regulatory network connecting splicing factors (SFs) and DEASEs. In addition, RNA-binding protein (RBP) mutations that can affect DEASEs were investigated to find some potential cancer drivers. Further association analysis demonstrated that DNA methylation levels were highly correlated with DEASEs. In summary, our results can bring new insight into understanding the mechanism of AS and provide novel biomarkers for personalized medicine of LUAD.

Widespread dysregulation of mRNA splicing implicates RNA processing in the development and progression of Huntington's disease

BACKGROUND

In Huntington's disease (HD), a CAG repeat expansion mutation in the Huntingtin (HTT) gene drives a gain-of-function toxicity that disrupts mRNA processing. Although dysregulation of gene splicing has been shown in human HD post-mortem brain tissue, post-mortem analyses are likely confounded by cell type composition changes in late-stage HD, limiting the ability to identify dysregulation related to early pathogenesis.

METHODS

To investigate gene splicing changes in early HD, we performed alternative splicing analyses coupled with a proteogenomics approach to identify early CAG length-associated splicing changes in an established isogenic HD cell model.

FINDINGS

We report widespread neuronal differentiation stage- and CAG length-dependent splicing changes, and find an enrichment of RNA processing, neuronal function, and epigenetic modification-related genes with mutant HTT-associated splicing. When integrated with a proteomics dataset, we identified several of these differential splicing events at the protein level. By comparing with human post-mortem and mouse model data, we identified common patterns of altered splicing from embryonic stem cells through to post-mortem striatal tissue.

INTERPRETATION

We show that widespread splicing dysregulation in HD occurs in an early cell model of neuronal development. Importantly, we observe HD-associated splicing changes in our HD cell model that were also identified in human HD striatum and mouse model HD striatum, suggesting that splicing-associated pathogenesis possibly occurs early in neuronal development and persists to later stages of disease. Together, our results highlight splicing dysregulation in HD which may lead to disrupted neuronal function and neuropathology.

FUNDING

This research is supported by the Lee Kong Chian School of Medicine, Nanyang Technological University Singapore Nanyang Assistant Professorship Start-Up Grant, the Singapore Ministry of Education under its Singapore Ministry of Education Academic Research Fund Tier 1 (RG23/22), the BC Children's Hospital Research Institute Investigator Grant Award (IGAP), and a Scholar Award from the Michael Smith Health Research BC.

Contribution of A-to-I RNA editing, M6A RNA Methylation, and Alternative Splicing to physiological brain aging and neurodegenerative diseases

Aging is a physiological and progressive phenomenon in all organisms' life cycle, characterized by the accumulation of degenerative processes triggered by several alterations within molecular pathways. These changes compromise cell fate, resulting in the loss of functions in tissues throughout the body, including the brain. Physiological brain aging has been linked to structural and functional alterations, as well as to an increased risk of neurodegenerative diseases. Post-transcriptional RNA modifications modulate mRNA coding properties, stability, translatability, expanding the coding capacity of the genome, and are involved in all cellular processes. Among mRNA post-transcriptional modifications, the A-to-I RNA editing, m6A RNA Methylation and Alternative Splicing play a critical role in all the phases of a neuronal cell life cycle and alterations in their mechanisms of action significantly contribute to aging and neurodegeneration. Here we review our current understanding of the contribution of A-to-I RNA editing, m6A RNA Methylation, and Alternative Splicing to physiological brain aging process and neurodegenerative diseases.

Overexpression of the QKI Gene Promotes Differentiation of Goat Myoblasts into Myotubes

The QKI genes encode RNA-binding proteins regulating cell proliferation, differentiation, and apoptosis. The Goat QKI has six isoforms, but their roles in myogenesis are unclear. In this study, the six isoforms of the QKI gene were overexpressed in goat myoblast. Immunofluorescence, qPCR and Western blot were used to evaluate the effect of QKI on the differentiation of goat myoblast. An RNA-Seq was performed on the cells with the gain of the function from the major isoforms to screen differentially expressed genes (DEGs). The results show that six isoforms had different degrees of deletion in exons 6 and 7, and caused the appearance of different types of encoded amino acids. The expression levels of the QKI-1 and QKI-5 groups were upregulated in the biceps femoris and latissimus dorsi muscle tissues compared with those of the QKI-4, QKI-7, QKI-3 and QKI-6 groups. After 6 d of myoblast differentiation, QKI-5 and the myogenic differentiators MyoG, MyoD, and MyHC were upregulated. Compared to the negative control group, QKI promoted myotube differentiation and the myoblasts overexpressing QKI-5 formed large, abundant myotubes. In summary, we identified that the overexpression of the QKI gene promotes goat-myoblast differentiation and that QKI-5 is the major isoform, with a key role. The RNA-Seq screened 76 upregulated and 123 downregulated DEGs between the negative control and the QKI-5-overexpressing goat myoblasts after d 6 of differentiation. The GO and KEGG analyses associated the downregulated DEGs with muscle-related biological functions. Only the pathways related to muscle growth and development were enriched. This study provides a theoretical basis for further exploring the regulatory mechanism of QKI in skeletal-muscle development in goats.

Modular architecture and functional annotation of human RNA-binding proteins containing RNA recognition motif

RNA-binding proteins (RBPs) are structurally and functionally diverse macromolecules with significant involvement in several post-transcriptional gene regulatory processes and human diseases. RNA recognition motif (RRM) is one of the most abundant RNA-binding domains in human RBPs. The unique modular architecture of each RBP containing RRM is crucial for its diverse target recognition and function. Genome-wide study of these structurally conserved and functionally diverse domains can enhance our understanding of their functional implications. In this study, modular architecture of RRM containing RBPs in human proteome is identified and systematically analysed. We observe that 30% of human RBPs with RNA-binding function contain RRM in single or multiple repeats or with other domains with maximum of six repeats. Zinc-fingers are the most frequently co-occurring domain partner of RRMs. Human RRM containing RBPs mostly belong to RNA metabolism class of proteins and are significantly enriched in two functional pathways including spliceosome and mRNA surveillance. Various human diseases are associated with 18% of the RRM containing RBPs. Single RRM containing RBPs are highly enriched in disorder regions. Gene ontology (GO) molecular functions including poly(A), poly(U) and miRNA binding are highly depleted in RBPs with single RRM, indicating the significance of modular nature of RRMs in specific function. The current study reports all the possible domain architectures of RRM containing human RBPs and their functional enrichment. The idea of domain architecture, and how they confer specificity and new functionalities to RBPs, can help in re-designing of modular RRM containing RBPs with re-engineered function.

Hippocampal ceRNA networks from chronic intermittent ethanol vapor-exposed male mice and functional analysis of top-ranked lncRNA genes for ethanol drinking phenotypes

The molecular mechanisms regulating the development and progression of alcohol use disorder (AUD) are largely unknown. While noncoding RNAs have previously been implicated as playing key roles in AUD, long-noncoding RNA (lncRNA) remains understudied in relation to AUD. In this study, we first identified ethanol-responsive lncRNAs in the mouse hippocampus that are transcriptional network hub genes. Microarray analysis of lncRNA, miRNA, circular RNA, and protein coding gene expression in the hippocampus from chronic intermittent ethanol vapor- or air- (control) exposed mice was used to identify ethanol-responsive competing endogenous RNA (ceRNA) networks. Highly interconnected lncRNAs (genes that had the strongest overall correlation to all other dysregulated genes identified) were ranked. The top four lncRNAs were novel, previously uncharacterized genes named Gm42575, 4930413E15Rik, Gm15767, and Gm33447, hereafter referred to as Pitt1, Pitt2, Pitt3, and Pitt4, respectively. We subsequently tested the hypothesis that CRISPR/Cas9 mutagenesis of the putative promoter and first exon of these lncRNAs in C57BL/6J mice would alter ethanol drinking behavior. The Drinking in the Dark (DID) assay was used to examine binge-like drinking behavior, and the Every-Other-Day Two-Bottle Choice (EOD-2BC) assay was used to examine intermittent ethanol consumption and preference. No significant differences between control and mutant mice were observed in the DID assay. Female-specific reductions in ethanol consumption were observed in the EOD-2BC assay for Pitt1, Pitt3, and Pitt4 mutant mice compared to controls. Male-specific alterations in ethanol preference were observed for Pitt1 and Pitt2. Female-specific increases in ethanol preference were observed for Pitt3 and Pitt4. Total fluid consumption was reduced in Pitt1 and Pitt2 mutants at 15% v/v ethanol and in Pitt3 and Pitt4 at 20% v/v ethanol in females only. We conclude that all lncRNAs targeted altered ethanol drinking behavior, and that lncRNAs Pitt1, Pitt3, and Pitt4 influenced ethanol consumption in a sex-specific manner. Further research is necessary to elucidate the biological mechanisms for these effects. These findings add to the literature implicating noncoding RNAs in AUD and suggest lncRNAs also play an important regulatory role in the disease.

Analysis of RPGR gene mutations in 41 Chinese families affected by X-linked inherited retinal dystrophy

Background: This study analyzed the phenotypes and genotypes of 41 Chinese families with inherited retinal dystrophy (IRD) and RPGR gene mutations.

Methods: This retrospective analysis evaluated a cohort of 41 patients who were subjected to a specific Hereditary Eye Disease Enrichment Panel (HEDEP) analysis. All (likely) pathogenic variants were determined by Sanger sequencing, and co-segregation analyses were performed on the available family members. All cases were subjected to Sanger sequencing for RPGR open reading frame 15 (ORF15) mutations.

Results: A total of 41 probands from different families with a clinical diagnosis of retinitis pigmentosa (RP; 34 cases) and cone-rod dystrophy (CORD; 7 cases) were included in this cohort. According to clinical information, 2, 18, and 21 cases were first assigned as autosomal dominant (AD), sporadic, and X-linked (XL) inheritance, respectively. Several cases of affected females who presented with a male phenotype have been described, posing challenges at diagnosis related to the apparent family history of AD. Mutations were located in RPGR exons or introns 1–14 and in ORF15 of 12 of 41 (29.3%) and 29 of 41 (70.7%) subjects, respectively. Thirty-four (likely) pathogenic mutations were identified. Frameshifts were the most frequently observed variants, followed by nonsense, splice, and missense mutations. Herein, a detailed description of four RP patients carrying RPGR intronic mutations is reported, and in vitro splice assays were performed to confirm the pathogenicity of these intronic mutations.

Conclusion: Our findings provide useful insights for the genetic and clinical counseling of patients with XL IRD, which will be useful for ongoing and future gene therapy trials.

RNA binding proteins in osteoarthritis

Osteoarthritis (OA) is a common chronic degenerative joint disease worldwide. The pathological features of OA are the erosion of articular cartilage, subchondral bone sclerosis, synovitis, and metabolic disorder. Its progression is characterized by aberrant expression of genes involved in inflammation, proliferation, and metabolism of chondrocytes. Effective therapeutic strategies are limited, as mechanisms underlying OA pathophysiology remain unclear. Significant research efforts are ongoing to elucidate the complex molecular mechanisms underlying OA focused on gene transcription. However, posttranscriptional alterations also play significant function in inflammation and metabolic changes related diseases. RNA binding proteins (RBPs) have been recognized as important regulators in posttranscriptional regulation. RBPs regulate RNA subcellular localization, stability, and translational efficiency by binding to their target mRNAs, thereby controlling their protein expression. However, their role in OA is less clear. Identifying RBPs in OA is of great importance to better understand OA pathophysiology and to figure out potential targets for OA treatment. Hence, in this manuscript, we summarize the recent knowledge on the role of dysregulated RBPs in OA and hope it will provide new insight for OA study and targeted treatment.

RNA-seq reveals insights into molecular mechanisms of metabolic restoration via tryptophan supplementation in low birth weight piglet model

Low birth weight (LBW) is associated with metabolic disorders in early life. While dietary l-tryptophan (Trp) can ameliorate postprandial plasma triglycerides (TG) disposal in LBW piglets, the genetic and biological basis underlying Trp-caused alterations in lipid metabolism is poorly understood. In this study, we collected 24 liver samples from 1-mo-old LBW and normal birth weight (NBW) piglets supplemented with different concentrations of dietary Trp (NBW with 0% Trp, N0; LBW with 0% Trp, L0; LBW with 0.4% Trp, L4; LBW with 0.8% Trp, L8; N = 6 in each group.) and conducted systematic, transcriptome-wide analysis using RNA sequencing (RNA-seq). We identified 39 differentially expressed genes (DEG) between N0 and L0, and genes within "increased dose effect" clusters based on dose-series expression profile analysis, enriched in fatty acid response of gene ontology (GO) biological process (BP). We then identified RNA-binding proteins including SRSF1, DAZAP1, PUM2, PCBP3, IGF2BP2, and IGF2BP3 significantly (P < 0.05) enriched in alternative splicing events (ASE) in comparison with L0 as control. There were significant positive and negative relationships between candidate genes from co-expression networks (including PID1, ANKRD44, RUSC1, and CYP2J34) and postprandial plasma TG concentration. Further, we determined whether these candidate hub genes were also significantly associated with metabolic and cardiovascular traits in humans via human phenome-wide association study (Phe-WAS), and analysis of mammalian orthologs suggests a functional conservation between human and pig. Our work demonstrates that transcriptomic changes during dietary Trp supplementation in LBW piglets. We detected candidate genes and related BP that may play roles on lipid metabolism restoration. These findings will help to better understand the amino acid support in LBW metabolic complications.

Non-canonical splice junction processing increases the diversity of RBFOX2 splicing isoforms

The underlying mechanisms of splicing regulation through non-canonical splice junction processing remain largely unknown. Here, we identified two RBFOX2 splicing isoforms by alternative 3' splice site selection of exon 9; the non-canonical splice junction processed RBFOX2 transcript (RBFOX2-N.C.) was expressed by the selection of the 3' splice GG acceptor sequence. The cytoplasmic localization of RBFOX2-C., a canonical splice junction-processed RBFOX2 transcript, was different from that of RBFOX2-N.C., which showed nuclear localization. In addition, we confirmed that RBFOX2-C. showed a significantly stronger localization into stress granules than RBFOX2-N.C. upon sodium arsenite treatment. Next, we investigated the importance of non-canonical 3' splice GG sequence selection of specific cis-regulatory elements using minigene constructs of the RBFOX2 gene. We found that the non-canonical 3' splice GG sequence and suboptimal branch point site adjacent region were critical for RBFOX2-N.C. expression through a non-canonical 3' splice selection. Our results suggest a regulatory mechanism for the non-canonical 3' splice selection in the RBFOX2 gene, providing a basis for studies related to the regulation of alternative pre-mRNA splicing through non-canonical splice junction processing.

RBM22, a Key Player of Pre-mRNA Splicing and Gene Expression Regulation, Is Altered in Cancer

RNA-Binding Proteins (RBP) are very diverse and cover a large number of functions in the cells. This review focuses on RBM22, a gene encoding an RBP and belonging to the RNA-Binding Motif (RBM) family of genes. RBM22 presents a Zinc Finger like and a Zinc Finger domain, an RNA-Recognition Motif (RRM), and a Proline-Rich domain with a general structure suggesting a fusion of two yeast genes during evolution: Cwc2 and Ecm2. RBM22 is mainly involved in pre-mRNA splicing, playing the essential role of maintaining the conformation of the catalytic core of the spliceosome and acting as a bridge between the catalytic core and other essential protein components of the spliceosome. RBM22 is also involved in gene regulation, and is able to bind DNA, acting as a bona fide transcription factor on a large number of target genes. Undoubtedly due to its wide scope in the regulation of gene expression, RBM22 has been associated with several pathologies and, notably, with the aggressiveness of cancer cells and with the phenotype of a myelodysplastic syndrome. Mutations, enforced expression level, and haploinsufficiency of RBM22 gene are observed in those diseases. RBM22 could represent a potential therapeutic target in specific diseases, and, notably, in cancer.

Case report: Altered pre-mRNA splicing caused by intronic variant c.1499 + 1G > A in the SLC4A4 gene

Proximal renal tubular acidosis (pRTA) with ocular abnormalities is an autosomal recessive disease caused by variants in the Solute Carrier Family 4 Member 4 (SLC4A4) gene. Patients present with metabolic acidosis and low plasma bicarbonate concentration (3∼17 mmol/L). In addition, they are often accompanied by ocular abnormalities, intellectual disability, and growth retardation. The patient underwent whole exome sequencing (WES) and bioinformatics analysis of variant pathogenicity in this study. Then, a minigene assay was conducted to analyze the splicing site variant further. Compound heterozygous variants in the SLC4A4 gene (NM_003759.3), c.145C > T (p.Arg49*) and c.1499 + 1G > A, were detected by WES. The minigene assay showed an mRNA splicing aberration caused by the c.1499 + 1G > A variant. Compared with the wild type, the mutant type caused 4-base insertion between exons 10 and 11 of SLC4A4 after expression in HEK293 cells. In conclusion, the c.1499 + 1G > A variant in the SLC4A4 gene may be one of the genetic causes in the patient. Moreover, our study provides the foundation for future gene therapy of such pathogenic variants.

Impact of Alternative Splicing Variants on Liver Cancer Biology

Simple Summary

Among the top ten deadly solid tumors are the two most frequent liver cancers, hepatocellular carcinoma, and intrahepatic cholangiocarcinoma, whose development and malignancy are favored by multifactorial conditions, which include aberrant maturation of pre-mRNA due to abnormalities in either the machinery involved in the splicing, i.e., the spliceosome and associated factors, or the nucleotide sequences of essential sites for the exon recognition process. As a consequence of cancer-associated aberrant splicing in hepatocytes- and cholangiocytes-derived cancer cells, abnormal proteins are synthesized. They contribute to the dysregulated proliferation and eventually transformation of these cells to phenotypes with enhanced invasiveness, migration, and multidrug resistance, which contributes to the poor prognosis that characterizes these liver cancers.

Abstract

The two most frequent primary cancers affecting the liver, whose incidence is growing worldwide, are hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA), which are among the five most lethal solid tumors with meager 5-year survival rates. The common difficulty in most cases to reach an early diagnosis, the aggressive invasiveness of both tumors, and the lack of favorable response to pharmacotherapy, either classical chemotherapy or modern targeted therapy, account for the poor outcome of these patients. Alternative splicing (AS) during pre-mRNA maturation results in changes that might affect proteins involved in different aspects of cancer biology, such as cell cycle dysregulation, cytoskeleton disorganization, migration, and adhesion, which favors carcinogenesis, tumor promotion, and progression, allowing cancer cells to escape from pharmacological treatments. Reasons accounting for cancer-associated aberrant splicing include mutations that create or disrupt splicing sites or splicing enhancers or silencers, abnormal expression of splicing factors, and impaired signaling pathways affecting the activity of the splicing machinery. Here we have reviewed the available information regarding the impact of AS on liver carcinogenesis and the development of malignant characteristics of HCC and iCCA, whose understanding is required to develop novel therapeutical approaches aimed at manipulating the phenotype of cancer cells.

HNRNPUL1 inhibits cisplatin sensitivity of esophageal squamous cell carcinoma through regulating the formation of circMAN1A2

Cisplatin (CDDP) is widely used for chemotherapy of esophageal squamous cell carcinoma (ESCC) but the drug resistance limits its therapeutic benefit. Heterogeneous nuclear ribonucleoprotein U-like 1 (HNRNPUL1) belongs to the family of RNA-binding proteins (RBPs) and is involved in DNA damage repair. To investigate whether and how HNRNPUL1 affects CDDP resistance of ESCC, we evaluated the expression of HNRNPUL1 and found that it was associated with recurrence in ESCC patients receiving postoperative platinum-based chemotherapy and was an independent prognostic factor for disease-free survival (DFS). Besides, we showed that the reduced expression of HNRNPUL1 enhanced the CDDP sensitivity of ESCC cells. Furthermore, RNA immunoprecipitation coupled with high-throughput sequencing (RIP-seq) were performed and a range of HNRNPUL1-binding RNAs influenced by CDDP treatment were identified followed by bioinformatics analysis. In terms of mechanism, we found that HNRNPUL1 inhibited CDDP sensitivity of ESCC cells by regulating the CDDP sensitivity-inhibited circular RNA (circRNA) MAN1A2 formation. Taken together, our results first demonstrated the role of HNRNPUL1 in CDDP resistance of ESCC and suggested that HNRNPUL1 may be a potential target of ESCC chemotherapy.

Inherited retinal diseases: Linking genes, disease-causing variants, and relevant therapeutic modalities

Inherited retinal diseases (IRDs) are a clinically complex and heterogenous group of visual impairment phenotypes caused by pathogenic variants in at least 277 nuclear and mitochondrial genes, affecting different retinal regions, and depleting the vision of affected individuals. Genes that cause IRDs when mutated are unique by possessing differing genotype-phenotype correlations, varying inheritance patterns, hypomorphic alleles, and modifier genes thus complicating genetic interpretation. Next-generation sequencing has greatly advanced the identification of novel IRD-related genes and pathogenic variants in the last decade. For this review, we performed an in-depth literature search which allowed for compilation of the Global Retinal Inherited Disease (GRID) dataset containing 4,798 discrete variants and 17,299 alleles published in 31 papers, showing a wide range of frequencies and complexities among the 194 genes reported in GRID, with 65% of pathogenic variants being unique to a single individual. A better understanding of IRD-related gene distribution, gene complexity, and variant types allow for improved genetic testing and therapies. Current genetic therapeutic methods are also quite diverse and rely on variant identification, and range from whole gene replacement to single nucleotide editing at the DNA or RNA levels. IRDs and their suitable therapies thus require a range of effective disease modelling in human cells, granting insight into disease mechanisms and testing of possible treatments. This review summarizes genetic and therapeutic modalities of IRDs, provides new analyses of IRD-related genes (GRID and complexity scores), and provides information to match genetic-based therapies such as gene-specific and variant-specific therapies to the appropriate individuals.

Impact of alternative splicing on mechanisms of resistance to anticancer drugs

A shared characteristic of many tumors is the lack of response to anticancer drugs. Multiple mechanisms of pharmacoresistance (MPRs) are involved in permitting cancer cells to overcome the effect of these agents. Pharmacoresistance can be primary (intrinsic) or secondary (acquired), i.e., triggered or enhanced in response to the treatment. Moreover, MPRs usually result in the lack of sensitivity to several agents, which accounts for diverse multidrug-resistant (MDR) phenotypes. MPRs are based on the dynamic expression of more than one hundred genes, constituting the so-called resistome. Alternative splicing (AS) during pre-mRNA maturation results in changes affecting proteins involved in the resistome. The resulting splicing variants (SVs) reduce the efficacy of anticancer drugs by lowering the intracellular levels of active agents, altering molecular targets, enhancing both DNA repair ability and defensive mechanism of tumors, inducing changes in the balance between pro-survival and pro-apoptosis signals, modifying interactions with the tumor microenvironment, and favoring malignant phenotypic transitions. Reasons accounting for cancer-associated aberrant splicing include mutations that create or disrupt splicing sites or splicing enhancers or silencers, abnormal expression of splicing factors, and impaired signaling pathways affecting the activity of the splicing machinery. Here we have reviewed the impact of AS on MPR in cancer cells.

DRAK2 aggravates nonalcoholic fatty liver disease progression through SRSF6-associated RNA alternative splicing

Nonalcoholic steatohepatitis (NASH) is an advanced stage of nonalcoholic fatty liver disease (NAFLD) with serious consequences that currently lacks approved pharmacological therapies. Recent studies suggest the close relationship between the pathogenesis of NAFLD and the dysregulation of RNA splicing machinery. Here, we reveal death-associated protein kinase-related apoptosis-inducing kinase-2 (DRAK2) is markedly upregulated in the livers of both NAFLD/NASH patients and NAFLD/NASH diet-fed mice. Hepatic deletion of DRAK2 suppresses the progression of hepatic steatosis to NASH. Comprehensive analyses of the phosphoproteome and transcriptome indicated a crucial role of DRAK2 in RNA splicing and identified the splicing factor SRSF6 as a direct binding protein of DRAK2. Further studies demonstrated that binding to DRAK2 inhibits SRSF6 phosphorylation by the SRSF kinase SRPK1 and regulates alternative splicing of mitochondrial function-related genes. In conclusion, our findings reveal an indispensable role of DRAK2 in NAFLD/NASH and offer a potential therapeutic target for this disease.

A Candidate RNAi Screen Reveals Diverse RNA-Binding Protein Phenotypes in Drosophila Flight Muscle

The proper regulation of RNA processing is critical for muscle development and the fine-tuning of contractile ability among muscle fiber-types. RNA binding proteins (RBPs) regulate the diverse steps in RNA processing, including alternative splicing, which generates fiber-type specific isoforms of structural proteins that confer contractile sarcomeres with distinct biomechanical properties. Alternative splicing is disrupted in muscle diseases such as myotonic dystrophy and dilated cardiomyopathy and is altered after intense exercise as well as with aging. It is therefore important to understand splicing and RBP function, but currently, only a small fraction of the hundreds of annotated RBPs expressed in muscle have been characterized. Here, we demonstrate the utility of Drosophila as a genetic model system to investigate basic developmental mechanisms of RBP function in myogenesis. We find that RBPs exhibit dynamic temporal and fiber-type specific expression patterns in mRNA-Seq data and display muscle-specific phenotypes. We performed knockdown with 105 RNAi hairpins targeting 35 RBPs and report associated lethality, flight, myofiber and sarcomere defects, including flight muscle phenotypes for Doa, Rm62, mub, mbl, sbr, and clu. Knockdown phenotypes of spliceosome components, as highlighted by phenotypes for A-complex components SF1 and Hrb87F (hnRNPA1), revealed level- and temporal-dependent myofibril defects. We further show that splicing mediated by SF1 and Hrb87F is necessary for Z-disc stability and proper myofibril development, and strong knockdown of either gene results in impaired localization of kettin to the Z-disc. Our results expand the number of RBPs with a described phenotype in muscle and underscore the diversity in myofibril and transcriptomic phenotypes associated with splicing defects. Drosophila is thus a powerful model to gain disease-relevant insight into cellular and molecular phenotypes observed when expression levels of splicing factors, spliceosome components and splicing dynamics are altered.

Heterozygous loss of WBP11 function causes multiple congenital defects in humans and mice

The genetic causes of multiple congenital anomalies are incompletely understood. Here, we report novel heterozygous predicted loss-of-function (LoF) and predicted damaging missense variants in the WW domain binding protein 11 (WBP11) gene in seven unrelated families with a variety of overlapping congenital malformations, including cardiac, vertebral, tracheo-esophageal, renal and limb defects. WBP11 encodes a component of the spliceosome with the ability to activate pre-messenger RNA splicing. We generated a Wbp11 null allele in mouse using CRISPR-Cas9 targeting. Wbp11 homozygous null embryos die prior to E8.5, indicating that Wbp11 is essential for development. Fewer Wbp11 heterozygous null mice are found than expected due to embryonic and postnatal death. Importantly, Wbp11 heterozygous null mice are small and exhibit defects in axial skeleton, kidneys and esophagus, similar to the affected individuals, supporting the role of WBP11 haploinsufficiency in the development of congenital malformations in humans. LoF WBP11 variants should be considered as a possible cause of VACTERL association as well as isolated Klippel-Feil syndrome, renal agenesis or esophageal atresia.

Systematic Profiling of mRNA Splicing Reveals the Prognostic Predictor and Potential Therapeutic Target for Glioblastoma Multiforme

Despite many changes in alternative splicing events (ASEs) are frequently involved in various cancers, prognosis-related ASEs and drug treatment targets in glioblastoma multiforme (GBM) have not been well explored. ASEs participate in many biological behaviors in the initiation and progression of tumors, the aberrant ASE has been considered another hallmark of cancer, and the systematic study of alternative splicing may provide potential biomarkers for malignancies. In this study, we carried out a systematic analysis to characterize the ASE signatures in GBM cohort. Through comparing GBM tissues and nontumor tissues, a total of 48,191 differently expressed ASEs from 10,727 genes were obtained, and these aberrant ASEs play an important role in the oncogenic process. Then, we identified 514 ASEs independently associated with patient survival in GBM by univariate and multivariate Cox regression, including exon skip in CD3D, alternate acceptor site in POLD2, and exon skip in DCN. Those prognostic models built on ASEs of each splice type can accurately predict the outcome of GBM patients, and values for the area under curve were 0.97 in the predictive model based on alternate acceptor site. In addition, the splicing-regulatory network revealed an interesting correlation between survival-associated splicing factors and prognostic ASE corresponding genes. Moreover, these three hub splicing factors in splicing regulation network are the potential targets of some drugs. In conclusion, a systematic analysis of ASE signatures in GBM could serve as an indicator for identifying novel prognostic biomarkers and guiding clinical treatment.

Systematic analysis of the function and prognostic value of RNA binding protein in head and neck squamous cell carcinoma

OBJECTIVE

Dysregulation of RNA binding proteins (RBPs) plays an important role in controlling processes in cancer development. However, the function of RBPs has not been thoroughly and systematically documented in head and neck cancer. We aim to explore the role of RPB in the pathogenesis of HNSC.

METHODS

We obtained HNSC gene expression data and corresponding clinical information from The Cancer Genome Atlas (TCGA) and the GEO databases, and identified aberrantly expressed RBPs between tumors and normal tissues. Meanwhile, we performed a series of bioinformatics to explore the function and prognostic value of these RBPs.

RESULTS

A total of 249 abnormally expressed RBPs were identified, including 101 downregulated RBPs and 148 upregulated RBPs. Using least absolute shrinkage and selection operator (LASSO) and univariate Cox regression analysis, the 15 RPBs were identified as hub genes. With the 15 RPBS, the prognostic prediction model was constructed. Further analysis showed that the high-risk score of the patients expressed a better survival outcome. The prediction model was validated in another HNSC dataset, and similar findings were observed.

CONCLUSIONS

Our results provide novel insights into the pathogenesis of HNSC. The fifteen RBP gene signature exhibited the predictive value of moderate HNSC prognosis, and have potential application value in clinical decision-making and individualized treatment.

The pan-cancer lncRNA PLANE regulates an alternative splicing program to promote cancer pathogenesis

Genomic amplification of the distal portion of chromosome 3q, which encodes a number of oncogenic proteins, is one of the most frequent chromosomal abnormalities in malignancy. Here we functionally characterise a non-protein product of the 3q region, the long noncoding RNA (lncRNA) PLANE, which is upregulated in diverse cancer types through copy number gain as well as E2F1-mediated transcriptional activation. PLANE forms an RNA-RNA duplex with the nuclear receptor co-repressor 2 (NCOR2) pre-mRNA at intron 45, binds to heterogeneous ribonucleoprotein M (hnRNPM) and facilitates the association of hnRNPM with the intron, thus leading to repression of the alternative splicing (AS) event generating NCOR2-202, a major protein-coding NCOR2 AS variant. This is, at least in part, responsible for PLANE-mediated promotion of cancer cell proliferation and tumorigenicity. These results uncover the function and regulation of PLANE and suggest that PLANE may constitute a therapeutic target in the pan-cancer context.

MiR-141-3p overexpression suppresses the malignancy of osteosarcoma by targeting FUS to degrade LDHB

Osteosarcoma (OS) is a common malignant bone cancer. Lactate dehydrogenase B (LDHB) has been revealed to act as a tumor promoter in several cancers. It is also revealed to be correlated with poor prognosis in OS, but its molecular mechanism in OS remains veiled. Our work illustrated that LDHB was overexpressed in OS tissues and cells, and it could enhance cell proliferation, migration, and invasion in OS. Subsequently, it was confirmed that fused in sarcoma (FUS) could bind with LDHB to positively regulate the stability of LDHB messenger RNA (mRNA). Besides, FUS expression was revealed to be elevated in OS tissues and positively correlate with LDHB expression. Furthermore, miR-141-3p, down-regulated in OS cells, was identified as the upstream regulator of FUS in OS cells. Besides, miR-141-3p overexpression decreased mRNA and protein levels of FUS and LDHB. More importantly, overexpression of miR-141-3p could impair FUS overexpression-mediated promotion on LDHB mRNA stability and expression. Finally, rescue assays indicated that miR-141-3p regulated OS cells cellular process via regulating LDHB. In sum, miR-141-3p targets FUS to degrade LDHB, thereby attenuating the malignancy of OS cells.

The Role of the U5 snRNP in Genetic Disorders and Cancer

Pre-mRNA splicing is performed by the spliceosome, a dynamic macromolecular complex consisting of five small uridine-rich ribonucleoprotein complexes (the U1, U2, U4, U5, and U6 snRNPs) and numerous auxiliary splicing factors. A plethora of human disorders are caused by genetic variants affecting the function and/or expression of splicing factors, including the core snRNP proteins. Variants in the genes encoding proteins of the U5 snRNP cause two distinct and tissue-specific human disease phenotypes – variants in PRPF6, PRPF8, and SNRP200 are associated with retinitis pigmentosa (RP), while variants in EFTUD2 and TXNL4A cause the craniofacial disorders mandibulofacial dysostosis Guion-Almeida type (MFDGA) and Burn-McKeown syndrome (BMKS), respectively. Furthermore, recurrent somatic mutations or changes in the expression levels of a number of U5 snRNP proteins (PRPF6, PRPF8, EFTUD2, DDX23, and SNRNP40) have been associated with human cancers. How and why variants in ubiquitously expressed spliceosome proteins required for pre-mRNA splicing in all human cells result in tissue-restricted disease phenotypes is not clear. Additionally, why variants in different, yet interacting, proteins making up the same core spliceosome snRNP result in completely distinct disease outcomes – RP, craniofacial defects or cancer – is unclear. In this review, we define the roles of different U5 snRNP proteins in RP, craniofacial disorders and cancer, including how disease-associated genetic variants affect pre-mRNA splicing and the proposed disease mechanisms. We then propose potential hypotheses for how U5 snRNP variants cause tissue specificity resulting in the restricted and distinct human disorders.

Characterization of RNP Networks of PUM1 and PUM2 Post-Transcriptional Regulators in TCam-2 Cells, a Human Male Germ Cell Model

Mammalian Pumilio (PUM) proteins are sequence-specific, RNA-binding proteins (RBPs) with wide-ranging roles. They are involved in germ cell development, which has functional implications in development and fertility. Although human PUM1 and PUM2 are closely related to each other and recognize the same RNA binding motif, there is some evidence for functional diversity. To address that problem, first we used RIP-Seq and RNA-Seq approaches, and identified mRNA pools regulated by PUM1 and PUM2 proteins in the TCam-2 cell line, a human male germ cell model. Second, applying global mass spectrometry-based profiling, we identified distinct PUM1- and PUM2-interacting putative protein cofactors, most of them involved in RNA processing. Third, combinatorial analysis of RIP and RNA-Seq, mass spectrometry, and RNA motif enrichment analysis revealed that PUM1 and PUM2 form partially varied RNP-regulatory networks (RNA regulons), which indicate different roles in human reproduction and testicular tumorigenesis. Altogether, this work proposes that protein paralogues with very similar and evolutionary highly conserved functional domains may play divergent roles in the cell by combining with different sets of protein cofactors. Our findings highlight the versatility of PUM paralogue-based post-transcriptional regulation, offering insight into the mechanisms underlying their diverse biological roles and diseases resulting from their dysfunction.

A novel alternative splicing-based prediction model for uteri corpus endometrial carcinoma

Alternative splicing (AS) is crucial a mechanism by which the complexity of mammalian and viral proteom increased overwhelmingly. There lacks systematic and comprehensive analysis of the prognostic significance of AS profiling landscape for uteri corpus endometrial carcinoma (UCEC). In this study, univariate and multivariate Cox regression analyses were conducted to identify candidate survival-associated AS events curated from SpliceSeq for the construction of prognostic index (PI) models. A correlation network between splicing factor-related AS events and significant survival-associated AS events were constructed using Cytoscape 3.5. As consequence, 28281 AS events from 8137 genes were detected from 506 UCEC patients, including 2630 survival-associated AS events. Kaplan Meier survival analysis revealed that six of the seven PI models (AD, AP, AT, ME, RI and ALL) exhibited good performance in stratifying the prognosis of low risk and high risk group (P<0.05). Among the six PI models, PI-AT performed best with an area under curves (AUC) value of 0.758 from time-dependent receiver operating characteristic. Correlation network implicated potential regulatory mechanism of AS events in UCEC. PI models based on survival-associated AS events for UCEC in this study showed preferable prognosis-predicting ability and may be promising prognostic indicators for UCEC patients.

Summary: This is the first study to systematically investigate the prognostic value of AS in UCEC. Findings in the presents study supported the clinical potential of AS for UCEC and shed light on the potential AS-associated molecular basis of UCEC.

Dysregulated lncRNA-miRNA-mRNA Network Reveals Patient Survival-Associated Modules and RNA Binding Proteins in Invasive Breast Carcinoma

Breast cancer is the most common cancer in women, but few biomarkers are effective in clinic. Previous studies have shown the important roles of non-coding RNAs in diagnosis, prognosis, and therapy selection for breast cancer and have suggested the significance of integrating molecules at different levels to interpret the mechanism of breast cancer. Here, we collected transcriptome data including long non-coding RNA (lncRNA), microRNA (miRNA), and mRNA for ~1,200 samples, including 1079 invasive breast carcinoma samples and 104 normal samples, from The Cancer Genome Atlas (TCGA) project. We identified differentially expressed lncRNAs, miRNAs, and mRNAs that distinguished invasive carcinoma samples from normal samples. We further constructed an integrated dysregulated network consisting of differentially expressed lncRNAs, miRNAs, and mRNAs and found housekeeping and cancer-related functions. Moreover, 58 RNA binding proteins (RBPs) involved in biological processes that are essential to maintain cell survival were found in the dysregulated network, and 10 were correlated with overall survival. In addition, we identified two modules that stratify patients into high- and low-risk subgroups. The expression patterns of these two modules were significantly different in invasive carcinoma versus normal samples, and some molecules were high-confidence biomarkers of breast cancer. Together, these data demonstrated an important clinical application for improving outcome prediction for invasive breast cancers.

Rbm14 maintains the integrity of genomic DNA during early mouse embryogenesis via mediating alternative splicing

OBJECTIVE

In this study, we generated an Rbm14 knockout mouse model to explore its functions during early mouse embryogenesis.

MATERIALS AND METHODS

The Rbm14 knockout mouse model was generated by a combination of clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 and microinjection techniques. The developmental defects of the knockout embryos were characterized by histological analyses. The accumulation of DNA damage in mouse embryonic stem cells (ESCs) was detected by γH2AX staining and comet assay. The altered mRNA splicing of DNA damage response (DDR)-related genes was detected by RNA-Seq analysis and confirmed by semi-quantitative PCR. The interaction of RBM14 with alternative splicing-related genes was detected by immunoprecipitation-mass spectra (IP-MS) and confirmed by co-immunoprecipitation (Co-IP).

RESULTS

Rbm14 knockout in mice results in apoptosis and cell proliferation defects in early post-implantation epiblast cells, leading to gastrulation disruption and embryonic lethality. FACS and immunostaining demonstrate accumulation of DNA damage in Rbm14 knockout ES cells. We also identified altered splicing of DDR-related genes in the knockout mouse ESCs by RNA-Seq, indicating that RBM14-mediated alternative splicing is required for the maintenance of genome integrity during early mouse embryogenesis.

CONCLUSIONS

Our work reveals that Rbm14 plays an essential role in the maintenance of genome integrity during early mouse embryonic development by regulating alternative splicing of DDR-related genes.

Activation of Hypoxia-Inducible Factor Signaling Modulates the RNA Protein Interactome in Caenorhabditis elegans

The cellular response to hypoxia is crucial to organismal survival, and hypoxia-inducible factors (HIF) are the key mediators of this response. HIF-signaling is central to many human diseases and mediates longevity in the nematode. Despite the rapidly increasing knowledge on RNA-binding proteins (RBPs), little is known about their contribution to hypoxia-induced cellular adaptation. We used RNA interactome capture (RIC) in wild-type Caenorhabditis elegans and vhl-1 loss-of-function mutants to fill this gap. This approach identifies more than 1,300 nematode RBPs, 270 of which can be considered novel RBPs. Interestingly, loss of vhl-1 modulates the RBPome. This difference is not primarily explained by protein abundance suggesting differential RNA-binding. Taken together, our study provides a global view on the nematode RBPome and proteome as well as their modulation by HIF-signaling. The resulting RBP atlas is also provided as an interactive online data mining tool (http://shiny.cecad.uni-koeln.de:3838/celegans_rbpome).

Investigation of RNA Editing Sites within Bound Regions of RNA-Binding Proteins

Studies in epitranscriptomics indicate that RNA is modified by a variety of enzymes. Among these RNA modifications, adenosine to inosine (A-to-I) RNA editing occurs frequently in the mammalian transcriptome. These RNA editing sites can be detected directly from RNA sequencing (RNA-seq) data by examining nucleotide changes from adenosine (A) to guanine (G), which substitutes for inosine (I). However, a careful investigation of such nucleotide changes must be conducted to distinguish sequencing errors and genomic mutations from the genuine editing sites. Building upon our recent introduction of an easy-to-use bioinformatics tool, RNA Editor, to detect RNA editing events from RNA-seq data, we examined the extent by which RNA editing events affect the binding of RNA-binding proteins (RBP). Through employing bioinformatic techniques, we uncovered that RNA editing sites occur frequently in RBP-bound regions. Moreover, the presence of RNA editing sites are more frequent when RNA editing islands were examined, which are regions in which RNA editing sites are present in clusters. When the binding of one RBP, human antigen R [HuR; encoded by ELAV-like protein 1 (ELAV1)], was quantified experimentally, its binding was reduced upon silencing of the RNA editing enzyme adenosine deaminases acting on RNA (ADAR) compared to the control-suggesting that the presence of RNA editing islands influence HuR binding to its target regions. These data indicate RNA editing as an important mediator of RBP-RNA interactions-a mechanism which likely constitutes an additional mode of post-transcription gene regulation in biological systems.

Disease modeling of core pre-mRNA splicing factor haploinsufficiency

The craniofacial disorder mandibulofacial dysostosis Guion-Almeida type is caused by haploinsufficiency of the U5 snRNP gene EFTUD2/SNU114. However, it is unclear how reduced expression of this core pre-mRNA splicing factor leads to craniofacial defects. Here we use a CRISPR-Cas9 nickase strategy to generate a human EFTUD2-knockdown cell line and show that reduced expression of EFTUD2 leads to diminished proliferative ability of these cells, increased sensitivity to endoplasmic reticulum (ER) stress and the mis-expression of several genes involved in the ER stress response. RNA-Seq analysis of the EFTUD2-knockdown cell line revealed transcriptome-wide changes in gene expression, with an enrichment for genes associated with processes involved in craniofacial development. Additionally, our RNA-Seq data identified widespread mis-splicing in EFTUD2-knockdown cells. Analysis of the functional and physical characteristics of mis-spliced pre-mRNAs highlighted conserved properties, including length and splice site strengths, of retained introns and skipped exons in our disease model. We also identified enriched processes associated with the affected genes, including cell death, cell and organ morphology and embryonic development. Together, these data support a model in which EFTUD2 haploinsufficiency leads to the mis-splicing of a distinct subset of pre-mRNAs with a widespread effect on gene expression, including altering the expression of ER stress response genes and genes involved in the development of the craniofacial region. The increased burden of unfolded proteins in the ER resulting from mis-splicing would exceed the capacity of the defective ER stress response, inducing apoptosis in cranial neural crest cells that would result in craniofacial abnormalities during development.

Regulatory RNA binding proteins contribute to the transcriptome-wide splicing alterations in human cellular senescence

Dysregulation of mRNA splicing has been observed in certain cellular senescence process. However, the common splicing alterations on the whole transcriptome shared by various types of senescence are poorly understood. In order to systematically identify senescence-associated transcriptomic changes in genome-wide scale, we collected RNA sequencing datasets of different human cell types with a variety of senescence-inducing methods from public databases and performed meta-analysis. First, we discovered that a group of RNA binding proteins were consistently down-regulated in diverse senescent samples and identified 406 senescence-associated common differential splicing events. Then, eight differentially expressed RNA binding proteins were predicted to regulate these senescence-associated splicing alterations through an enrichment analysis of their RNA binding information, including motif scanning and enhanced cross-linking immunoprecipitation data. In addition, we constructed the splicing regulatory modules that might contribute to senescence-associated biological processes. Finally, it was confirmed that knockdown of the predicted senescence-associated potential splicing regulators through shRNAs in HepG2 cell line could result in senescence-like splicing changes. Taken together, our work demonstrated a broad range of common changes in mRNA splicing switches and detected their central regulatory RNA binding proteins during senescence. These findings would help to better understand the coordinating splicing alterations in cellular senescence.

Deciphering human ribonucleoprotein regulatory networks

RNA-binding proteins (RBPs) control and coordinate each stage in the life cycle of RNAs. Although in vivo binding sites of RBPs can now be determined genome-wide, most studies typically focused on individual RBPs. Here, we examined a large compendium of 114 high-quality transcriptome-wide in vivo RBP-RNA cross-linking interaction datasets generated by the same protocol in the same cell line and representing 64 distinct RBPs. Comparative analysis of categories of target RNA binding preference, sequence preference, and transcript region specificity was performed, and identified potential posttranscriptional regulatory modules, i.e. specific combinations of RBPs that bind to specific sets of RNAs and targeted regions. These regulatory modules represented functionally related proteins and exhibited distinct differences in RNA metabolism, expression variance, as well as subcellular localization. This integrative investigation of experimental RBP-RNA interaction evidence and RBP regulatory function in a human cell line will be a valuable resource for understanding the complexity of post-transcriptional regulation.

Dysregulation of the Splicing Machinery Is Associated to the Development of Nonalcoholic Fatty Liver Disease

CONTEXT

Nonalcoholic fatty liver disease (NAFLD) is a common obesity-associated pathology characterized by hepatic fat accumulation, which can progress to fibrosis, cirrhosis, and hepatocellular carcinoma. Obesity is associated with profound changes in gene-expression patterns of the liver, which could contribute to the onset of comorbidities.

OBJECTIVE

As these alterations might be linked to a dysregulation of the splicing process, we aimed to determine whether the dysregulation in the expression of splicing machinery components could be associated with NAFLD.

PARTICIPANTS

We collected 41 liver biopsies from nonalcoholic individuals with obesity, with or without hepatic steatosis, who underwent bariatric surgery.

INTERVENTIONS

The expression pattern of splicing machinery components was determined using a microfluidic quantitative PCR-based array. An in vitro approximation to determine lipid accumulation using HepG2 cells was also implemented.

RESULTS

The liver of patients with obesity and steatosis exhibited a severe dysregulation of certain splicing machinery components compared with patients with obesity without steatosis. Nonsupervised clustering analysis allowed the identification of three molecular phenotypes of NAFLD with a unique fingerprint of alterations in splicing machinery components, which also presented distinctive hepatic and clinical-metabolic alterations and a differential response to bariatric surgery after 1 year. In addition, in vitro silencing of certain splicing machinery components (i.e., PTBP1, RBM45, SND1) reduced fat accumulation and modulated the expression of key de novo lipogenesis enzymes, whereas conversely, fat accumulation did not alter spliceosome components expression.

CONCLUSION

There is a close relationship between splicing machinery dysregulation and NAFLD development, which should be further investigated to identify alternative therapeutic targets.

Pre-mRNA Splicing Modulation by Antisense Oligonucleotides

Pre-mRNA splicing, a dynamic process of intron removal and exon joining, is governed by a combinatorial control exerted by overlapping cis-elements that are unique to each exon and its flanking intronic sequences. Splicing cis-elements are usually 4-to-8-nucleotide-long linear motifs that provide binding sites for specific proteins. Pre-mRNA splicing is also influenced by secondary and higher order RNA structures that affect accessibility of splicing cis-elements. Antisense oligonucleotides (ASOs) that block splicing cis-elements and/or affect RNA structure have been shown to modulate splicing in vivo. Therefore, ASO-based strategies have emerged as a powerful tool for therapeutic manipulation of splicing in pathological conditions. Here we describe an ASO-based approach to increase the production of the full-length SMN2 mRNA in spinal muscular atrophy patient cells.