Disentangling the splicing factor programs underlying complex molecular phenotypes.

The regulation of exon inclusion through alternative splicing tunes the cell's behavior by increasing the functional diversity of the transcriptome and the proteome. Splicing factors work in concert to generate gene isoform pools that contribute to cell phenotypes yet their activity is controlled by multiple regulatory and signaling layers. This hinders identification of functional, phenotype-specific splicing factors using traditional single-omic measurements, such as their mutational state or expression. To address this challenge, we propose repurposing the virtual inference of protein activity by enriched regulon analysis (VIPER) to measure splicing factor activity solely from their downstream exon transcriptomic inclusion signatures. This approach is effective in assessing the effect of co-occurring splicing factor perturbations, as well as their post-translational regulation. As proof of concept, we dissect recurrent splicing factor programs underlying tumorigenesis including aberrantly activated factors acting as oncogenes and inactivated ones acting as tumor suppressors, which are undetectable by more conventional methodologies. Activation and inactivation of these cancer splicing programs effectively stratifies overall survival, as well as cancer hallmarks such as proliferation and immune evasion. Altogether, repurposing network-based inference of protein activity for splicing factor networks distills common, functionally relevant splicing programs in otherwise heterogeneous molecular contexts.

Oncogenic NOVA1 expression dysregulates alternative splicing in breast cancer.

Neuro-Oncological Ventral Antigen 1 ( NOVA1 ) is best known for its role in mediating an alternative splicing (AS) program in neurons, yet was first discovered as an antigen expressed in breast tumors, causing rare autoimmune reactions and paraneoplastic neurological disorders (PNDs). The PND model suggests a plausible role of the tumor antigen expression in tumor suppression, whereas it has emerged that NOVA may function as an oncogene in a variety of cancers. In addition, whether NOVA mediates AS in breast cancer remains unanswered. Here we examine the AS profiles of breast invasive carcinoma (BRCA) tumor samples and demonstrate that ectopic NOVA1 expression led to the activation of neuron-like splicing patterns in many genes, including exons targeted by NOVA in the brain. The splicing dysregulation is especially prevalent in cell periphery and cytoskeleton genes related to cell-cell communication, actin-based movement, and neuronal functions. We find that NOVA1-mediated AS is most prominent in Luminal A tumors and high NOVA1 expression in this subtype is associated with poorer prognosis. Our results suggest that ectopic NOVA1 in tumors has regulatory activity affecting pathways with high relevance to tumor progression and that this might be a more general mechanism for PND antigens.

A clinicopathological-imaging nomogram for the prediction of pathological complete response in breast cancer cases administered neoadjuvant therapy.

OBJECTIVE: To construct a user-friendly nomogram with MRI and clinicopathological parameters for the prediction of pathological complete response (pCR) after neoadjuvant therapy (NAT) in patients with breast cancer (BC). METHODS: We retrospectively enrolled consecutive female patients pathologically confirmed with breast cancer who received NAT followed by surgery between January 2018 and December 2022 as the development cohort. Additionally, we prospectively collected eligible candidates between January 2023 and December 2023 as an external validation group at our institution. Pretreatment MRI features and clinicopathological variables were collected, and the pre- and post-treatment background parenchymal enhancement (BPE) and the changes in BPE on two MRIs were compared between patients who achieved pCR and those who did not. Multivariable logistic regression analysis was used to identify independent variables associated with pCR in the development cohort. These independent variables were combined into a predictive nomogram for which performance was assessed using the area under the receiver operating characteristic curve (AUC), calibration plot, decision curve analysis, and external validation. RESULTS: In the development cohort, there were a total of 276 female patients with a mean age of 48.3 ± 8.7 years, while in the validation cohort, there were 87 female patients with a mean age of 49.0 ± 9.5 years. Independent prognostic factors of pCR included small tumor size, HER2(+), high Ki-67 index,high signal enhancement ratio (SER), low minimum value of apparent diffusion coefficient (ADCmin), and significantly decreased BPE after NAT(change of BPE). The nomogram, which incorporates the above parameters, demonstrated excellent predictive performance in both the development and external validation cohorts, with AUC values of 0.900 and 0.850, respectively. Additionally, the nomogram showed excellent calibration capacities, as indicated by Hosmer-Lemeshow test p values of 0.508 and 0.423 in the two cohorts. Furthermore, the nomogram provided greater net benefits compared to the default simple schemes in both cohorts. CONCLUSION: A nomogram constructed using tumor size, HER2 status, Ki-67 index, SER, ADCmin, and changes in pre- and post-NAT BPE demonstrated strong predictive performance, calibration ability, and greater net benefits for predicting pCR in patients with BC after NAT. This suggests that the user-friendly nomogram could be a valuable imaging biomarker for identifying suitable candidates for NAT.

D-arabinose induces cell cycle arrest by promoting autophagy via p38 MAPK signaling pathway in breast cancer.

Breast cancer patients often have a poor prognosis largely due to lack of effective targeted therapy. It is now well established that monosaccharide enhances growth retardation and chemotherapy sensitivity in tumor cells. We investigated whether D-arabinose has capability to restrict the proliferation of tumor cells and its mechanism. Here, we report that D-arabinose induced cytotoxicity is modulated by autophagy and p38 MAPK signaling pathway in breast cancer cell lines. The proliferation of cells was evaluated by CCK-8 and Colony formation assay. The distribution of cells in cell cycle phases was analyzed by flow cytometry. Cell cycle, autophagy and MAPK signaling related proteins were detected by western blotting. Mouse xenograft model was used to evaluate the efficacy of D-arabinose in vivo. The proliferation of cells was dramatically inhibited by D-arabinose exposure in a dose-dependent manner, which was relevant to cell cycle arrest, as demonstrated by G2/M cell cycle restriction and ectopic expression of cell cycle related proteins. Mechanistically, we further identified that D-arabinose is positively associated with autophagy and the activation of the p38 MAPK signaling in breast cancer. In contrast, 3-Ma or SB203580, the inhibitor of autophagy or p38 MAPK, reversed the efficacy of D-arabinose. Additionally, D-arabinose in vivo treatment could significantly inhibit xenograft growth of breast cancer cells. Our findings were the first to reveal that D-arabinose triggered cell cycle arrest by inducing autophagy through the activation of p38 MAPK signaling pathway in breast cancer cells.

The features associated with mammography-occult MRI-detected newly diagnosed breast cancer analysed by comparing machine learning models with a logistic regression model.

PURPOSE: To compare machine learning (ML) models with logistic regression model in order to identify the optimal factors associated with mammography-occult (i.e. false-negative mammographic findings) magnetic resonance imaging (MRI)-detected newly diagnosed breast cancer (BC). MATERIAL AND METHODS: The present single-centre retrospective study included consecutive women with BC who underwent mammography and MRI (no more than 45 days apart) for breast cancer between January 2018 and May 2023. Various ML algorithms and binary logistic regression analysis were utilized to extract features linked to mammography-occult BC. These features were subsequently employed to create different models. The predictive value of these models was assessed using receiver operating characteristic curve analysis. RESULTS: This study included 1957 malignant lesions from 1914 patients, with an average age of 51.64 ± 9.92 years and a range of 20-86 years. Among these lesions, there were 485 mammography-occult BCs. The optimal features of mammography-occult BC included calcification status, tumour size, mammographic density, age, lesion enhancement type on MRI, and histological type. Among the different ML models (ANN, L1-LR, RF, and SVM) and the LR-based combined model, the ANN model with RF features was found to be the optimal model. It demonstrated the best discriminative performance in predicting mammography false- negative findings, with an AUC of 0.912, an accuracy of 86.90%, a sensitivity of 85.85%, and a specificity of 84.18%. CONCLUSION: Mammography-occult MRI-detected breast cancers have features that should be considered when performing breast MRI to improve the detection rate for breast cancer and aid in clinician management.

Identification of Breast Cancer Immune-related Prognostic Characteristics in Tumor Microenvironment.

BACKGROUND: Accumulated evidence suggest that tumor microenvironment (TME) plays a crucial role in breast cancer (BRCA) progression and therapeutic effects. OBJECTIVE: This study aimed to characterize immune-related BRCA subtypes in TME, and identify genes with prognostic value. METHODS: RNA sequencing profiles with corresponding clinical data from The Cancer Genome Atlas (TCGA) database of BRCA patients were downloaded to evaluate immune infiltration using the single-sample gene set enrichment (ssGAEA) algorithm. Further, BRCA was clustered according to immune infiltration status by consensus clustering analysis. Using Venn analysis, differentially expressed genes (DEGs) were overlapped to obtain candidate genes. Kaplan-Meier (K-M) analysis was performed to identify prognostic genes, and the results were verified in the GEO and METABRIC datasets. RT-qPCR was conducted to detect the mRNA expression of prognostic genes. RESULTS: In the TCGA database, 3 immune-related BRCA subtypes were identified [cluster1 (C1), cluster2 (C2), and cluster3 (C2)]. The C2 subtype had better overall survival (OS) compared to the C1 subtype. Higher levels of immune markers and checkpoint protein were found in the C2 subtype than in others. By combining DEGs between BRCA and normal tissues, with the C1 and C2 subtypes associated with different OS, 25 BRCA candidate genes were identified. Among these, 8 genes were identified as prognostic genes for BRCA. RT-qPCR showed that the expressions of 2 genes were significantly elevated in BRCA tissues, while that of other genes were decreased. CONCLUSION: Three BRCA subtypes were identified with the immune index, which may help design advanced treatment of BRCA. The data code used for the analysis in this article was available on GitHub (https://github.com/tangzhn/BRCA1.git).

Multiomics analyses reveal DARS1-AS1/YBX1-controlled posttranscriptional circuits promoting glioblastoma tumorigenesis/radioresistance.

The glioblastoma (GBM) stem cell-like cells (GSCs) are critical for tumorigenesis/therapeutic resistance of GBM. Mounting evidence supports tumor-promoting function of long noncoding RNAs (lncRNAs), but their role in GSCs remains poorly understood. By combining CRISPRi screen with orthogonal multiomics approaches, we identified a lncRNA DARS1-AS1-controlled posttranscriptional circuitry that promoted the malignant properties of GBM cells/GSCs. Depleting DARS1-AS1 inhibited the proliferation of GBM cells/GSCs and self-renewal of GSCs, prolonging survival in orthotopic GBM models. DARS1-AS1 depletion also impaired the homologous recombination (HR)-mediated double-strand break (DSB) repair and enhanced the radiosensitivity of GBM cells/GSCs. Mechanistically, DARS1-AS1 interacted with YBX1 to promote target mRNA binding and stabilization, forming a mixed transcriptional/posttranscriptional feed-forward loop to up-regulate expression of the key regulators of G1-S transition, including E2F1 and CCND1. DARS1-AS1/YBX1 also stabilized the mRNA of FOXM1, a master transcription factor regulating GSC self-renewal and DSB repair. Our findings suggest DARS1-AS1/YBX1 axis as a potential therapeutic target for sensitizing GBM to radiation/HR deficiency-targeted therapy.

Systematic dissection of coordinated stromal remodeling identifies Sox10+ glial cells as a therapeutic target in myelofibrosis.

Remodeling of the tissue niche is often evident in diseases, yet, the stromal alterations and their contribution to pathogenesis are poorly characterized. Bone marrow fibrosis is a maladaptive feature of primary myelofibrosis (PMF). We performed lineage tracing and found that most collagen-expressing myofibroblasts were derived from leptin-receptor-positive (LepR+) mesenchymal cells, whereas a minority were from Gli1-lineage cells. Deletion of Gli1 did not impact PMF. Unbiased single-cell RNA sequencing (scRNA-seq) confirmed that virtually all myofibroblasts originated from LepR-lineage cells, with reduced expression of hematopoietic niche factors and increased expression of fibrogenic factors. Concurrently, endothelial cells upregulated arteriolar-signature genes. Pericytes and Sox10+ glial cells expanded drastically with heightened cell-cell signaling, suggesting important functional roles in PMF. Chemical or genetic ablation of bone marrow glial cells ameliorated fibrosis and improved other pathology in PMF. Thus, PMF involves complex remodeling of the bone marrow microenvironment, and glial cells represent a promising therapeutic target.

RANKL Promotes Chemotherapy Resistance in Breast Cancer Cells Through STAT3 Mediated Autophagy Induction.

BACKGROUND: This study was to investigate the functional role and mechanism of receptor activator of nuclear factor-kappa B ligand (RANKL) associated autophagy and chemoresistance in breast cancer. MATERIALS AND METHODS: Cell Counting Kit-8 (CCK-8) assay was used to detect the cell viability. Real-time polymerase chain reaction (PCR) was used for determining the relative mRNA levels of key genes and protein expression was assessed by Western blotting. Immunofluorescence was performed to evaluate the changes in the autophagy flux. Short hairpin (shRNA) was used to knockdown the expression of the target genes in breast cancer cells. Based on The Cancer Genome Atlas (TCGA) database, we explored the expression of receptor activator of nuclear factor-kappa B (RANK), autophagy and signal transducer and activator of transcription 3 (STAT3) signaling associated genes and analyzed their correlation with the prognosis of breast cancer patients. RESULTS: The findings showed that receptor activator of nuclear factor-kappa B ligand (RANKL), the ligand of RANK, could effectively enhance the chemoresistance potential of breast cancer cells. Our results showed that RANKL induced autophagy and enhanced the expression of autophagy associated genes in breast cancer cells. The knockdown of RANK suppressed RANKL mediated autophagy induction in these cells. Furthermore, the inhibition of autophagy suppressed RANKL mediated chemoresistance in breast cancer cells. We found STAT3 signaling pathway was involved in RANKL-induced autophagy. Analysis of the expression of RANK, and autophagy and STAT3 signaling associated genes in breast cancer tissues showed that the expression of autophagy and STAT3 signaling associated genes was correlated with the prognosis of breast cancer patients. CONCLUSION: The present study suggests that the RANKL/RANK axis may potentially mediate chemoresistance in breast cancer cells by inducing autophagy through the STAT3 signaling pathway.

Gain of toxic function by long-term AAV9-mediated SMN overexpression in the sensorimotor circuit.

The neurodegenerative disease spinal muscular atrophy (SMA) is caused by deficiency in the survival motor neuron (SMN) protein. Currently approved SMA treatments aim to restore SMN, but the potential for SMN expression beyond physiological levels is a unique feature of adeno-associated virus serotype 9 (AAV9)-SMN gene therapy. Here, we show that long-term AAV9-mediated SMN overexpression in mouse models induces dose-dependent, late-onset motor dysfunction associated with loss of proprioceptive synapses and neurodegeneration. Mechanistically, aggregation of overexpressed SMN in the cytoplasm of motor circuit neurons sequesters components of small nuclear ribonucleoproteins, leading to splicing dysregulation and widespread transcriptome abnormalities with prominent signatures of neuroinflammation and the innate immune response. Thus, long-term SMN overexpression interferes with RNA regulation and triggers SMA-like pathogenic events through toxic gain-of-function mechanisms. These unanticipated, SMN-dependent and neuron-specific liabilities warrant caution on the long-term safety of treating individuals with SMA with AAV9-SMN and the risks of uncontrolled protein expression by gene therapy.

Vascular-derived SPARC and SerpinE1 regulate interneuron tangential migration and accelerate functional maturation of human stem cell-derived interneurons.

Cortical interneurons establish inhibitory microcircuits throughout the neocortex and their dysfunction has been implicated in epilepsy and neuropsychiatric diseases. Developmentally, interneurons migrate from a distal progenitor domain in order to populate the neocortex - a process that occurs at a slower rate in humans than in mice. In this study, we sought to identify factors that regulate the rate of interneuron maturation across the two species. Using embryonic mouse development as a model system, we found that the process of initiating interneuron migration is regulated by blood vessels of the medial ganglionic eminence (MGE), an interneuron progenitor domain. We identified two endothelial cell-derived paracrine factors, SPARC and SerpinE1, that enhance interneuron migration in mouse MGE explants and organotypic cultures. Moreover, pre-treatment of human stem cell-derived interneurons (hSC-interneurons) with SPARC and SerpinE1 prior to transplantation into neonatal mouse cortex enhanced their migration and morphological elaboration in the host cortex. Further, SPARC and SerpinE1-treated hSC-interneurons also exhibited more mature electrophysiological characteristics compared to controls. Overall, our studies suggest a critical role for CNS vasculature in regulating interneuron developmental maturation in both mice and humans.

Hyaluronic acid functionalized biodegradable mesoporous silica nanocomposites for efficient photothermal and chemotherapy in breast cancer.

Chemotherapy is one of conventional treatment methods for breast cancer, but drug toxicity and side effects have severely limited its clinical applications. Photothermal therapy has emerged as a promising method that, upon combination with chemotherapy, can better treat breast cancer. In this context, a biodegradable mesoporous silica nanoparticle (bMSN NPs) system was developed for loading doxorubicin (DOX) and IR780, to be potentially applied in the treatment of breast cancer. IR780 is encapsulated in the pores of bMSN NPs by hydrophobic adsorption, while DOX is adsorbed on the surface of the bMSN NPs by hyaluronic acid electrostatically, to form the bMID NPs. Transmission electron microscopy, fluorescence spectrum and UV absorption spectrum are used to prove the successful encapsulation of IR780 and the loading of DOX. In vitro experiments have shown bMID NPs present an excellent therapeutic effect on breast cancer cells. In vivo fluorescence imaging results have indicated that bMID NPs can accumulate in tumor sites gradually and achieve in vivo long-term circulation and continuous drug release. Furthermore, bMID NPs have provided obvious antitumor effects in breast cancer mouse models, thus evolving as an efficient platform for breast cancer therapy.

Research on Temperature Variation during Coal and Gas Outbursts: Implications for Outburst Prediction in Coal Mines.

Coal and gas outbursts are among the most severe disasters threatening the safety of coal mines around the world. They are dynamic phenomena characterized by large quantities of coal and gas ejected from working faces within a short time. Numerous researchers have conducted studies on outburst prediction, and a variety of indices have been developed to this end. However, these indices are usually empirical or based on local experience, and the accurate prediction of outbursts is not feasible due to the complicated mechanisms of outbursts. This study conducts outburst experiments using large-scale multifunctional equipment developed in the laboratory to develop a more robust outburst prediction method. In this study, the coal temperature during the outburst process was monitored using temperature sensors. The results show that the coal temperature decreased rapidly as the outburst progressed. Meanwhile, the coal temperature in locations far from the outburst mouth increased. The coal broken in the stress concentration state is the main factor causing the abnormal temperature rise. The discovery of these phenomena lays a theoretical foundation and provides an experimental basis for an effective outburst prediction method. An outburst prediction method based on monitoring temperature was proposed, and has a simpler and faster operation process and is not easily disturbed by coal mining activities. What is more, the critical values of coal temperature rises or temperature gradients can be flexibly adjusted according to the actual situations of different coal mines to predict outbursts more effectively and accurately.

Modeling RNA-Binding Protein Specificity In Vivo by Precisely Registering Protein-RNA Crosslink Sites.

RNA-binding proteins (RBPs) regulate post-transcriptional gene expression by recognizing short and degenerate sequence motifs in their target transcripts, but precisely defining their binding specificity remains challenging. Crosslinking and immunoprecipitation (CLIP) allows for mapping of the exact protein-RNA crosslink sites, which frequently reside at specific positions in RBP motifs at single-nucleotide resolution. Here, we have developed a computational method, named mCross, to jointly model RBP binding specificity while precisely registering the crosslinking position in motif sites. We applied mCross to 112 RBPs using ENCODE eCLIP data and validated the reliability of the discovered motifs by genome-wide analysis of allelic binding sites. Our analyses revealed that the prototypical SR protein SRSF1 recognizes clusters of GGA half-sites in addition to its canonical GGAGGA motif. Therefore, SRSF1 regulates splicing of a much larger repertoire of transcripts than previously appreciated, including HNRNPD and HNRNPDL, which are involved in multivalent protein assemblies and phase separation.

Stage-specific requirement for Mettl3-dependent m6A mRNA methylation during haematopoietic stem cell differentiation.

Haematopoietic stem cells (HSCs) maintain balanced self-renewal and differentiation, but how these functions are precisely regulated is not fully understood. N6-methyladenosine (m6A) messenger RNA methylation has emerged as an important mode of epitranscriptional gene expression regulation affecting many biological processes. We show that deletion of the m6A methyltransferase Mettl3 from the adult haematopoietic system led to an accumulation of HSCs in the bone marrow and a marked reduction of reconstitution potential due to a blockage of HSC differentiation. Interestingly, deleting Mettl3 from myeloid cells using Lysm-cre did not impact myeloid cell number or function. RNA sequencing revealed 2,073 genes with significant m6A modifications in HSCs. Myc was identified as a direct target of m6A in HSCs. Mettl3-deficient HSCs failed to upregulate MYC expression following stimulation to differentiate and enforced expression of Myc rescued differentiation defects of Mettl3-deficient HSCs. Our results reveal a key role of m6A in governing HSC differentiation.

Cell type-specific CLIP reveals that NOVA regulates cytoskeleton interactions in motoneurons.

BACKGROUND: Alternative RNA processing plays an essential role in shaping cell identity and connectivity in the central nervous system. This is believed to involve differential regulation of RNA processing in various cell types. However, in vivo study of cell type-specific post-transcriptional regulation has been a challenge. Here, we describe a sensitive and stringent method combining genetics and CLIP (crosslinking and immunoprecipitation) to globally identify regulatory interactions between NOVA and RNA in the mouse spinal cord motoneurons. RESULTS: We developed a means of undertaking motoneuron-specific CLIP to explore motoneuron-specific protein-RNA interactions relative to studies of the whole spinal cord in mouse. This allowed us to pinpoint differential RNA regulation specific to motoneurons, revealing a major role for NOVA in regulating cytoskeleton interactions in motoneurons. In particular, NOVA specifically promotes the palmitoylated isoform of the cytoskeleton protein Septin 8 in motoneurons, which enhances dendritic arborization. CONCLUSIONS: Our study demonstrates that cell type-specific RNA regulation is important for fine tuning motoneuron physiology and highlights the value of defining RNA processing regulation at single cell type resolution.

LIN28 Selectively Modulates a Subclass of Let-7 MicroRNAs.

LIN28 is a bipartite RNA-binding protein that post-transcriptionally inhibits the biogenesis of let-7 microRNAs to regulate development and influence disease states. However, the mechanisms of let-7 suppression remain poorly understood because LIN28 recognition depends on coordinated targeting by both the zinc knuckle domain (ZKD), which binds a GGAG-like element in the precursor, and the cold shock domain (CSD), whose binding sites have not been systematically characterized. By leveraging single-nucleotide-resolution mapping of LIN28 binding sites in vivo, we determined that the CSD recognizes a (U)GAU motif. This motif partitions the let-7 microRNAs into two subclasses, precursors with both CSD and ZKD binding sites (CSD+) and precursors with ZKD but no CSD binding sites (CSD-). LIN28 in vivo recognition-and subsequent 3' uridylation and degradation-of CSD+ precursors is more efficient, leading to their stronger suppression in LIN28-activated cells and cancers. Thus, CSD binding sites amplify the regulatory effects of LIN28.

Systematic reconstruction of autism biology from massive genetic mutation profiles.

Autism spectrum disorder (ASD) affects 1% of world population and has become a pressing medical and social problem worldwide. As a paradigmatic complex genetic disease, ASD has been intensively studied and thousands of gene mutations have been reported. Because these mutations rarely recur, it is difficult to (i) pinpoint the fewer disease-causing versus majority random events and (ii) replicate or verify independent studies. A coherent and systematic understanding of autism biology has not been achieved. We analyzed 3392 and 4792 autism-related mutations from two large-scale whole-exome studies across multiple resolution levels, that is, variants (single-nucleotide), genes (protein-coding unit), and pathways (molecular module). These mutations do not recur or replicate at the variant level, but significantly and increasingly do so at gene and pathway levels. Genetic association reveals a novel gene + pathway dual-hit model, where the mutation burden becomes less relevant. In multiple independent analyses, hundreds of variants or genes repeatedly converge to several canonical pathways, either novel or literature-supported. These pathways define recurrent and systematic ASD biology, distinct from previously reported gene groups or networks. They also present a catalog of novel ASD risk factors including 118 variants and 72 genes. At a subpathway level, most variants disrupt the pathway-related gene functions, and in the same gene, they tend to hit residues extremely close to each other and in the same domain. Multiple interacting variants spotlight key modules, including the cAMP (adenosine 3',5'-monophosphate) second-messenger system and mGluR (metabotropic glutamate receptor) signaling regulation by GRKs (G protein-coupled receptor kinases). At a superpathway level, distinct pathways further interconnect and converge to three biology themes: synaptic function, morphology, and plasticity.

Prediction of clustered RNA-binding protein motif sites in the mammalian genome.

Sequence-specific interactions of RNA-binding proteins (RBPs) with their target transcripts are essential for post-transcriptional gene expression regulation in mammals. However, accurate prediction of RBP motif sites has been difficult because many RBPs recognize short and degenerate sequences. Here we describe a hidden Markov model (HMM)-based algorithm mCarts to predict clustered functional RBP-binding sites by effectively integrating the number and spacing of individual motif sites, their accessibility in local RNA secondary structures and cross-species conservation. This algorithm learns and quantifies rules of these features, taking advantage of a large number of in vivo RBP-binding sites obtained from cross-linking and immunoprecipitation data. We applied this algorithm to study two representative RBP families, Nova and Mbnl, which regulate tissue-specific alternative splicing through interacting with clustered YCAY and YGCY elements, respectively, and predicted their binding sites in the mouse transcriptome. Despite the low information content in individual motif elements, our algorithm made specific predictions for successful experimental validation. Analysis of predicted sites also revealed cases of extensive and distal RBP-binding sites important for splicing regulation. This algorithm can be readily applied to other RBPs to infer their RNA-regulatory networks. The software is freely available at http://zhanglab.c2b2.columbia.edu/index.php/MCarts.

NOVA-dependent regulation of cryptic NMD exons controls synaptic protein levels after seizure.

The neuronal RNA binding protein NOVA regulates splicing, shuttles to the cytoplasm, and co-localizes with target transcripts in dendrites, suggesting links between splicing and local translation. Here we identified >200 transcripts showing NOVA-dependent changes in abundance, but, surprisingly, HITS-CLIP revealed NOVA binds these RNAs in introns rather than 3' UTRs. This led us to discover NOVA-regulated splicing of cryptic exons within these introns. These exons triggered nonsense mediated decay (NMD), as UPF1 and protein synthesis were required for NOVA's effect on RNA levels. Their regulation was dynamic and physiologically relevant. The NMD exons were regulated by seizures, which also induced changes in Nova subcellular localization and mediated large changes in synaptic proteins, including proteins implicated in familial epilepsy. Moreover, Nova haploinsufficient mice had spontaneous epilepsy. The data reveal a hidden means of dynamic RNA regulation linking electrical activity to splicing and protein output, and of mediating homeostatic excitation/inhibition balance in neurons.DOI:http://dx.doi.org/10.7554/eLife.00178.001.