Search

Circular RNAs (circRNAs) are a class of single-stranded closed non-coding RNA molecules (ncRNAs), which are formed as a result of reverse splicing of mRNAs. Despite their relative abundance, an interest in understanding their regulatory importance is rather recent. High stability, abundance and evolutionary conservation among species underline some of their important traits. CircRNAs perform a variety of cellular functions ranging from miRNA and proteins sponges to transcriptional modulation and splicing. Additionally, most circRNAs are expressed aberrantly in pathological conditions suggesting their possible exploitation as diagnostic biomarkers. Their covalent closed cyclic structure resulting in resistance to RNases further makes them suitable as cancer biomarkers. Studies involving human tumors have verified differences in the expression profiles of circRNAs, indicating a regulatory role in cancer pathogenesis and metastasis. As endogenous competitive RNA, circRNAs can regulate tumor proliferation and invasion. Further, some circRNAs located in the nucleus can regulate transcription of genes by binding to RNA polymerase II. In this review, we elaborate the characteristics, functions and mechanisms of action of circRNAs in cancer. We also discuss the possibility of using circRNAs as potential therapeutic targets and biomarkers for cancer.

Emerging roles of circular RNAs (circRNAs) in various biological processes have advanced our knowledge of transcriptional and posttranscriptional gene regulation. To date, no research has been conducted to explore their roles in the rice–Xanthomonas oryzae pv. oryzae interaction. Therefore, we identified 3,517 circRNAs from rice leaves infected with the highly virulent X. oryzae pv. oryzae strain PXO99A by using rRNA depleted RNA sequencing technique coupled with the CIRI2 and CIRCexplorer2 pipeline. Characterization analyses showed that these circRNAs were distributed across the whole genome of rice, and most circRNAs arose from exons (85.13%), ranged from 200 to 1,000 bp, and were with a noncanonical GT/AG (including CT/AC equivalent) splicing signal. Functional annotation and enrichment analysis of the host genes that produced the differentially expressed circRNAs (DEcircRNAs) suggested that these identified circRNAs might play an important role in reprogramming rice responses to PXO99A invasion, mainly by mediating photorespiration and chloroplast, peroxisome, and diterpenoid biosynthesis. Moreover, 31 DEcircRNAs were predicted to act as microRNA decoys in rice. The expression profile of four DEcircRNAs were validated by quantitative real-time PCR with divergent primers, and the back-splicing sites of seven DEcircRNAs were verified by PCR analysis and Sanger sequencing. Collectively, these results inferred a potential functional role of circRNAs in the regulation of rice immunity and provide novel clues about the molecular mechanisms of rice–PXO99A interaction.

CTGrepeat expansion (CTGexp) is associated with aberrant alternate splicing that contributes to cardiac dysfunction in myotonic dystrophy type 1 (DM1). Excision of thisCTGexprepeat using CRISPR-Cas resulted in the disappearance of punctate ribonuclear foci in cardiomyocyte-like cells derived from DM1-induced pluripotent stem cells (iPSCs). This was associated with correction of the underlying spliceopathy as determined byRNA sequencingand alternate splicing analysis. Certain genes were of particular interest due to their role in cardiac development, maturation, and function (TPM4,CYP2J2,DMD,MBNL3,CACNA1H,ROCK2,ACTB) or their association with splicing (SMN2,GCFC2,MBNL3). Moreover, while comparing isogenic CRISPR-Cas9-corrected versus non-corrected DM1 cardiomyocytes, a prominent difference in the splicing pattern for a number of candidate genes was apparent pertaining to genes that are associated with cardiac function (TNNT,TNNT2,TTN,TPM1,SYNE1,CACNA1A,MTMR1,NEBL,TPM1),cellular signaling(NCOR2,CLIP1,LRRFIP2,CLASP1,CAMK2G), and other DM1-related genes (i.e.,NUMA1,MBNL2,LDB3) in addition to the disease-causingDMPKgene itself. Subsequent validation using a selected gene subset, includingMBNL1,MBNL2,INSR,ADD3, andCRTC2, further confirmed correction of the spliceopathy followingCTGexprepeat excision. To our knowledge, the present study provides the first comprehensive unbiased transcriptome-wide analysis of the differential splicing landscape in DM1 patient-derived cardiac cells after excision of theCTGexprepeat using CRISPR-Cas9, showing reversal of the abnormal cardiac spliceopathy in DM1.

Genome-wide transcriptome analysis provides systems-level insights into plant biology. Due to the limited depth of quantitative proteomics our understanding of gene-protein-complex stoichiometry is largely unknown in plants. Recently, the complexity of the proteome and its cell-/tissue-specific distribution have boosted the research community to the integration of transcriptomics and proteomics landscapes in a proteogenomic approach. Herein, we generated a quantitative proteome and transcriptome abundance atlas of 15 major sweet cherry (Prunus avium L., cv 'Tragana Edessis') tissues represented by 29 247 genes and 7584 proteins. Additionally, 199 984 alternative splicing events, particularly exon skipping and alternative 3' splicing, were identified in 23 383 transcribed regions of the analyzed tissues. Common signatures as well as differences between mRNA and protein quantities, including genes encoding transcription factors and allergens, within and across the different tissues are reported. Using our integrated dataset, we identified key putative regulators of fruit development, notably genes involved in the biosynthesis of anthocyanins and flavonoids. We also provide proteogenomic-based evidence for the involvement of ethylene signaling and pectin degradation in cherry fruit ripening. Moreover, clusters of genes and proteins with similar and different expression and suppression trends across diverse tissues and developmental stages revealed a relatively low RNA abundance-to-protein correlation. The present proteogenomic analysis allows us to identify 17 novel sweet cherry proteins without prior protein-level annotation evidenced in the currently available databases. To facilitate use by the community, we also developed the Sweet Cherry Atlas Database (https://grcherrydb.com/) for viewing and data mining these resources. This work provides new insights into the proteogenomics workflow in plants and a rich knowledge resource for future investigation of gene and protein functions in Prunus species.

Circular RNAs with covalently linked ends are generated from many eukaryotic protein-coding genes when the pre-mRNA splicing machinery backsplices. These mature transcripts are resistant to digestion by exonucleases and typically have much longer half-lives than their associated linear mRNAs. Circular RNAs thus have great promise as sensitive biomarkers, including for detection of transcriptional activity. Here, we show that circular RNAs can serve as markers of readthrough transcription events in Drosophila and human cells, thereby revealing mechanistic insights into RNA polymerase II transcription termination as well as pre-mRNA 3' end processing. We describe methods that take advantage of plasmids that generate a circular RNA when an upstream polyadenylation signal fails to be used and/or RNA polymerase II fails to terminate. As a proof-of-principle, we show that RNAi-mediated depletion of well-established transcription termination factors, including the RNA endonuclease Cpsf73, results in increased circular RNA output from these plasmids in Drosophila and human cells. This method is generalizable as a circular RNA can be easily encoded downstream of any genomic region of interest. Circular RNA biomarkers, therefore, have great promise for identifying novel cellular factors and conditions that impact transcription termination processes.

Alternative splicing of pre-mRNAs is an important gene regulatory mechanism shaping the transcriptome. AtMC1 is an Arabidopsis thaliana type I metacaspase that positively regulates the hypersensitive response. Here, we found that AtMC1 is involved in the regulation of plant immunity to the bacterial pathogen Pseudomonas syringae pv. tomato DC3000 and is physically associated with Sm-like4 (LSM4), which is involved in pre-mRNA splicing. AtMC1 and LSM4 protein levels both increased with their coexpression as compared with their separate expression in vivo. Like AtMC1, LSM4 negatively regulates plant immunity to P. syringae pv. tomato DC3000 infection. By RNA sequencing, AtMC1 was shown to modulate the splicing of many pre-mRNAs, including 4CL3, which is a negative regulator of plant immunity. Thus, AtMC1 plays a regulatory role in pre-mRNA splicing, which might contribute to AtMC1-mediated plant immunity. [Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

Autosomal dominant optic atrophy (ADOA) is frequently caused by mutations in the optic atrophy 1 (OPA1) gene, with haploinsufficiency being the major genetic pathomechanism. Almost 30% of the OPA1-associated cases suffer from splice defects. We identified a novel OPA1 mutation, c.1065+5G>A, in patients with ADOA. In patient-derived fibroblasts, the mutation led to skipping of OPA1 exon 10, reducing the OPA1 protein expression by approximately 50%. We developed a molecular treatment to correct the splice defect in OPA1 using engineered U1 splice factors retargeted to different locations in OPA1 exon 10 or intron 10. The strongest therapeutic effect was detected when U1 binding was engineered to bind to intron 10 at position +18, a position predicted by bioinformatics to be a promising binding site. We were able to significantly silence the effect of the mutation (skipping of exon 10) and simultaneously increase the expression level of normal transcripts. Retargeting U1 to the canonical splice donor site did not lead to a detectable splice correction. This proof-of-concept study indicates for the first time the feasibility of splice mutation correction as a treatment option for ADOA. Increasing the amount of correctly spliced OPA1 transcripts may suffice to overcome the haploinsufficiency.

Covalently closed circular RNAs (circRNAs) produced by back-splicing of exon(s) are co-expressed with their cognate linear RNAs from the same gene loci. Most circRNAs are fully overlapped with their cognate linear RNAs in sequences except the back-spliced junction (BSJ) site, thus challenging the computational detection, experimental validation and hence functional evaluation of circRNAs. Nevertheless, specific bioinformatic pipelines were developed to identify fragments mapped to circRNA-featured BSJ sites, and circRNAs were pervasively identified from non-polyadenylated RNA-seq datasets in different cell lines/tissues and across species. Precise identification and quantification of circRNAs provide a basis to further understand their functions. Here, we describe detailed computational steps to annotate and quantify circRNAs using a series of CIRCexplorer pipelines.

Split intein-mediated protein trans-splicing expands AAV transfer capacity, thus overcoming the limited AAV cargo. However, non-mammalian inteins persist as trans-splicing by-products, and this could raise safety concerns for AAV intein clinical applications. In this study, we tested the ability of several degrons to selectively decrease levels of inteins after protein trans-splicing and found that a version of E. coli dihydrofolate reductase, which we have shortened to better fit into the AAV vector, is the most effective. We show that subretinal administration of AAV intein armed with this short degron is both safe and effective in a mouse model of Stargardt disease (STGD1), which is the most common form of inherited macular degeneration in humans. This supports the use of optimized AAV intein for gene therapy of both STGD1 and other conditions that require transfer of large genes., Graphical abstract, Intein-mediated protein trans-splicing (PTS) expands AAV gene transfer capacity. However, non-mammalian inteins persist as trans-splicing by-products. Here, we show that ecDHFR selectively degrades inteins after PTS and that AAV intein vectors armed with this degron are both safe and effective in the retina of a mouse model of genetic blindness.

Stable Intronic Sequence RNA (sisRNA) is a relatively new class of non-coding RNA. Found in many organisms, these sisRNA produced from their host genes are generally involved in regulatory roles, controlling gene expression at multiple levels through active involvement in regulatory feedback loops. Large scale identification of sisRNA via genome-wide RNA sequencing has been difficult, largely in part due to its low abundance. Done on its own, RNA sequencing often yields a large mass of information that is ironically uninformative; the potential sisRNA reads being masked by other highly abundant RNA species like ribosomal RNA and messenger RNA. In this review, we present a practical workflow for the enrichment of circular sisRNA through the use of transcriptionally quiescent systems, rRNA-depletion, and RNase R treatment prior to deep sequencing. This workflow allows circular sisRNA to be reliably detected. We also present various methods to experimentally validate the circularity and stability of the circular sisRNA identified, as well as a few methods for further functional characterisation.

The honeybee is a haplodiploid organism in which sexual development is determined by the complementary sex determiner (csd) gene and realized by sex-specific splicing processes involving the feminizer (fem) gene. We used high throughput transcriptome sequencing (RNA-Seq) to characterize the transcriptional differences between the sexes caused by the fertilization and sex determination processes in honeybee (Apis mellifera) embryos. We identified 758, 372 and 43 differentially expressed genes (DEGs) and 58, 176 and 233 differentially spliced genes (DSGs) in 10-15-h-old, 25-40-h-old and 55-70-h-old female and male embryos, respectively. The early difference in male and female embryos in response to the fertilization and nonfertilization processes resulted mainly in differential expression of genes (758 DEGs versus 58 DSGs). In the latest sampled embryonic stage, the transcriptional differences between the sexes were dominated by alternative splicing of transcripts (43 DEGs versus 233 DSGs). Interestingly, differentially spliced transcripts that encode RNA-binding properties were overrepresented in 55-70-h-old embryos, indicating a more diverse regulation via alternative splicing than previous work on the sex determination pathway suggested. These stage- and sex-specific transcriptome data from honeybee embryos provide a comprehensive resource for examining the roles of fertilization and sex determination in developmental programming in a haplodiploid system. This article is protected by copyright. All rights reserved.

Signal regulatory protein SIRPγ (CD172G) is expressed on the surface of lymphocytes where it acts by engaging its ligand, CD47. SIRPG, which encodes SIRPγ, contains a non-synonymous coding variant, rs6043409, which is significantly associated with risk for type 1 diabetes. SIRPG produces multiple transcript isoforms via alternative splicing, all encoding potentially functional proteins. We show that rs6043409 alters a predicted exonic splicing enhancer, resulting in significant shifts in the distribution of SIRPG transcript isoforms. All of these transcript isoforms produced protein upon transient expression in vitro. However, CRISPR targeting of one of the alternatively spliced exons in SIRPG eliminated all SIRPγ expression in Jurkat T cells. These targeted cells formed fewer cell-cell conjugates with each other than with wild type Jurkat cells, expressed reduced levels of genes associated with CD47 signaling and had significantly increased levels of cell surface CD47. In primary CD4+ and CD8+ T cells cell surface SIRPγ levels in response to anti-CD3 stimulation varied quantitatively by rs6043409 genotype. Our results suggest that SIRPG is the most likely causative gene for type 1 diabetes risk in the 20p13 region and highlight the role of alternative splicing in lymphocytes in mediating the genetic risk for autoimmunity.

Lettuce (Lactuca sativa) is one of the most important vegetable crops in the world. As a leafy vegetable, the polymorphism of lettuce leaves from dark to pale green is an important trait. However, the genetic and molecular mechanisms underlying such variations remain poorly understood. In this study, one major locus controlling the polymorphism of dark- and pale-green leaves in lettuce was identified using genome-wide association studies (GWAS). This locus was then fine mapped to an interval of 5375 bp on chromosome 4 using a segregating population containing 2480 progeny. Only one gene, homologous to the GLK genes in Arabidopsis and other plants, is present in the candidate region. A complementation test confirmed that the candidate gene, LsGLK, contributes to the variation of dark- and pale-green leaves. Sequence analysis showed that a CACTA transposon of 7434 bp was inserted 10 bp downstream of the stop codon of LsGLK, followed by a duplication of a 1826-bp fragment covering exons 3-6 of the LsGLK gene. The transposon insertion did not change the expression level of the LsGLK gene. However, because of alternative splicing, only 6% of the transcripts produced from the transposon insertion were wild-type transcripts, which led to the production of pale-green leaves. An evolutionary analysis revealed that the insertion of the CACTA transposon occurred in cultivated lettuce and might have been selected in particular cultivars to satisfy the diverse demands of consumers. In this study, we demonstrated that a transposon insertion near a gene may affect its splicing and consequently generate phenotypic variations.

Alternative splicing (AS) is a key process in which precursor RNAs produce different mature RNAs, and the disorder of AS is a key factor in promoting cancer development. Compared with coding RNA, studies on the functions of long non-coding RNAs (lncRNAs) are far from enough. In fact, lncRNA is an important participant and regulator in the process of AS. On the one hand, lncRNAs regulate cancer progression as AS products of precursor messenger RNA (mRNA), but on the other hand, precursor lncRNA generates cancer-related abnormal splicing variants through AS. In addition, lncRNAs directly or indirectly regulate the AS events of downstream target genes, thus affecting the occurrence and development of cancer. Here, we reviewed how lncRNAs regulate AS and influence oncogenesis in different ways.

Neuron-restrictive silencing factor (NRSF), also known as RE-1 silencing transcription factor (REST), has pivotal functions in many neuron-specific genes. Previous studies revealed that neuron-specific alternative splicing (AS) of REST produces divergent forms of REST variants and provides regulatory complexity in the nervous system. However, the biological significance of these variants in the regulation of neuronal activities remains to be clarified. Here, we revealed that Charlatan (Chn), a Drosophila REST-like molecule, is also regulated by neuron-specific AS. Neuron-specific AS produced six divergent variants of Chn proteins, one of which preferentially localized to axons. A small sequence of this variant was especially important for the axonal localization. Our data suggest that some variants have roles beyond the transcriptional regulation of neuronal activities.

MALAT1, an abundant lncRNA specifically localized to nuclear speckles, regulates alternative-splicing (AS). The molecular basis of its role in AS remains poorly understood. Here, we report three conserved, thermodynamically stable, parallel RNA-G-quadruplexes (rG4s) present in the 3′ region of MALAT1 which regulates this function. Using rG4 domain-specific RNA-pull-down followed by mass-spectrometry, RNA-immuno-precipitation, and imaging, we demonstrate the rG4 dependent localization of Nucleolin (NCL) and Nucleophosmin (NPM) to nuclear speckles. Specific G-to-A mutations that abolish rG4 structures, result in the localization loss of both the proteins from speckles. Functionally, disruption of rG4 in MALAT1 phenocopies NCL knockdown resulting in altered pre-mRNA splicing of endogenous genes. These results reveal a central role of rG4s within the 3′ region of MALAT1 orchestrating AS.

RNA sequencing (RNAseq) can reveal gene fusions, splicing variants, mutations/indels in addition to differential gene expression, thus providing a more complete genetic picture than DNA sequencing. This most widely used technology in genomics tool box has evolved from classic bulk RNA sequencing (RNAseq), popular single cell RNA sequencing (scRNAseq) to newly emerged spatial RNA sequencing (spRNAseq). Bulk RNAseq studies average global gene expression, scRNAseq investigates single cell RNA biology up to 20,000 individual cells simultaneously, while spRNAseq has ability to dissect RNA activities spatially, representing next generation of RNA sequencing. This article highlights these technologies, characteristic features and suitable applications in precision oncology.

Circular RNAs (circRNAs) are a complicated class of non-coding RNAs that have a covalently closed loop structure and are very stable and cautious. Multiple biological processes of malignancy, including tumorigenesis, development, invasion, metastasis, apoptosis, and vascularization, are disrupted by an increased number of circRNAs. Recent research has showed that circRNAs, functioning as microRNA (miRNA) sponges or protein scaffolds, interacting with RNA-binding proteins (RBPs), and autophagy regulators, affect the transcription and splicing regulation. Many circRNAs have tissue-specific expression patterns and are heavily conserved. CircRNA levels in neurons are dynamically modulated. Growing evidence suggests that circRNAs are highly abundant in neural tissues, perhaps owing to the proliferation of particular genes that promote circularization, implying that circRNA dysregulation is linked to nervous system disorders including glioma. The most widespread and deadly primary malignant brain tumor is glioma. CircRNA has a close connection to glioma, according to reported research. Here, the current knowledge about the properties of circRNAs is introduced and the biological and molecular functions of circRNAs are described. Then, the clinical association of circRNAs with glioma/glioblastoma and their level of expression and their regulatory mechanisms in tumorigenesis are discussed. Moreover, the potential of circRNAs as diagnostic biomarkers and predictors of brain cancer risk and possible therapeutic targets in medicine is examined.

Uterine carcinosarcoma (UCS) is a highly invasive malignant tumor that originated from the uterine epithelium. Many studies suggested that the abnormal changes of alternative splicing (AS) of pre-mRNA are related to the occurrence and metastasis of the tumor. This study investigates the mechanism of alternative splicing events (ASEs) in the tumorigenesis and metastasis of UCS. RNA-seq of UCS samples and alternative splicing event (ASE) data of UCS samples were downloaded from The Cancer Genome Atlas (TCGA) and TCGASpliceSeq databases, several times. Firstly, we performed the Cox regression analysis to identify the overall survival-related alternative splicing events (OSRASEs). Secondly, a multivariate model was applied to approach the prognostic values of the risk score. Afterwards, a coexpressed network between splicing factors (SFs) and OSRASEs was constructed. In order to explore the relationship between the potential prognostic signaling pathways and OSRASEs, we fabricated a network between these pathways and OSRASEs. Finally, validations from multidimension platforms were used to explain the results unambiguously. 1,040 OSRASEs were identified by Cox regression. Then, 6 OSRASEs were incorporated in a multivariable model by Lasso regression. The area under the curve (AUC) of the receiver operator characteristic (ROC) curve was 0.957. The risk score rendered from the multivariate model was corroborated to be an independent prognostic factor ( P < 0.001 ). In the network of SFs and ASEs, junction plakoglobin (JUP) noteworthily regulated RALGPS1-87608-AT ( P < 0.001 , R = 0.455 ). Additionally, RALGPS1-87608-AT ( P = 0.006 ) showed a prominent relationship with distant metastasis. KEGG pathways related to prognosis of UCS were selected by gene set variation analysis (GSVA). The pyrimidine metabolism ( P < 0.001 , R = − 0.470 ) was the key pathway coexpressed with RALGPS1. We considered that aberrant JUP significantly regulated RALGPS1-87608-AT and the pyrimidine metabolism pathway might play a significant part in the metastasis and prognosis of UCS.

High-risk human papillomaviruses (HPV), exemplified by HPV16/18, are causally linked to human cancers of the anogenital tract, skin, and upper aerodigestive tract. Previously, we identified Ecdysoneless (ECD) protein, the human homolog of the Drosophila ecdysoneless gene, as a novel HPV16 E6–interacting protein. Here, we show that ECD, through its C-terminal region, selectively binds to high-risk but not to low-risk HPV E6 proteins. We demonstrate that ECD is overexpressed in cervical and head and neck squamous cell carcinoma (HNSCC) cell lines as well as in tumor tissues. Using The Cancer Genome Atlas dataset, we show that ECD mRNA overexpression predicts shorter survival in patients with cervical and HNSCC. We demonstrate that ECD knockdown in cervical cancer cell lines led to impaired oncogenic behavior, and ECD co-overexpression with E7 immortalized primary human keratinocytes. RNA-sequencing analyses of SiHa cells upon ECD knockdown showed to aberrations in E6/E7 RNA splicing, as well as RNA splicing of several HPV oncogenesis–linked cellular genes, including splicing of components of mRNA splicing machinery itself. Taken together, our results support a novel role of ECD in viral and cellular mRNA splicing to support HPV-driven oncogenesis. Implications: This study links ECD overexpression to poor prognosis and shorter survival in HNSCC and cervical cancers and identifies a critical role of ECD in cervical oncogenesis through regulation of viral and cellular mRNA splicing.

© 2022 Prudêncio et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society. This article, published in RNA, is available under a Creative Commons License (Attribution 4.0 International), as described at http://creativecommons.org/licenses/by/4.0/., Widespread cotranscriptional splicing has been demonstrated from yeast to human. However, most studies to date addressing the kinetics of splicing relative to transcription used either Saccharomyces cerevisiae or metazoan cultured cell lines. Here, we adapted native elongating transcript sequencing technology (NET-seq) to measure cotranscriptional splicing dynamics during the early developmental stages of Drosophila melanogaster embryos. Our results reveal the position of RNA polymerase II (Pol II) when both canonical and recursive splicing occur. We found heterogeneity in splicing dynamics, with some RNAs spliced immediately after intron transcription, whereas for other transcripts no splicing was observed over the first 100 nt of the downstream exon. Introns that show splicing completion before Pol II has reached the end of the downstream exon are necessarily intron-defined. We studied the splicing dynamics of both nascent pre-mRNAs transcribed in the early embryo, which have few and short introns, as well as pre-mRNAs transcribed later in embryonic development, which contain multiple long introns. As expected, we found a relationship between the proportion of spliced reads and intron size. However, intron definition was observed at all intron sizes. We further observed that genes transcribed in the early embryo tend to be isolated in the genome whereas genes transcribed later are often overlapped by a neighboring convergent gene. In isolated genes, transcription termination occurred soon after the polyadenylation site, while in overlapped genes, Pol II persisted associated with the DNA template after cleavage and polyadenylation of the nascent transcript. Taken together, our data unravel novel dynamic features of Pol II transcription and splicing in the developing Drosophila embryo., This work was supported by funding to M.C.-F. (Fundação para a Ciência e Tecnologia, FCT/Ministério da Ciência, Tecnologia e Ensino Superior - Fundos do Orçamento de Estado [UIDB/50005/2020], and FCT/ FEDER/POR Lisboa 2020, Programa Operacional Regional de Lisboa PORTUGAL 2020, grant LISBOA-01-0145-FEDER-016394) and to R.G.M (FCT grant PTDC/BIA-BID/28441/2017). P.P. was a recipient of an FCT fellowship (SFRH/BD/109689/2015). R.S. was a recipient of an EMBO Long-Term Fellowship (EMBO ALTF 101-2019). This project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreements no. 842695 (Marie Skłodowska- Curie Actions) and no. 857119 (RiboMed)

Background. Breast cancer (BRCA) is one of the most common cancers and the leading cause of cancer-related death in women. RNA-binding proteins (RBPs) play an important role in the emergence and pathogenesis of tumors. The target RNAs of RBPs are very diverse; in addition to binding to mRNA, RBPs also bind to noncoding RNA. Noncoding RNA can cause secondary structures that can bind to RBPs and regulate multiple processes such as splicing, RNA modification, protein localization, and chromosomes remodeling, which can lead to tumor initiation, progression, and invasion. Methods. (1) BRCA data were downloaded from The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) databases and were used as training and testing datasets, respectively. (2) The prognostic RBPs-related genes were screened according to the overlapping differentially expressed genes (DEGs) from the TCGA database. (3) Univariate Cox proportional hazard regression was performed to identify the genes with significant prognostic value. (4) Further, we used the LASSO regression to construct a prognostic signature and validated the signature in the TCGA and ICGC cohort. (5) Besides, we also performed prognostic analysis, expression level verification, immune cell correlation analysis, and drug correlation analysis of the genes in the model. Results. Four genes (MRPL13, IGF2BP1, BRCA1, and MAEL) were identified as prognostic gene signatures. The prognostic model has been validated in the TCGA and ICGC cohorts. The risk score calculated with four genes signatures could largely predict overall survival for 1, 3, and 5 years in patients with BRCA. The calibration plot demonstrated outstanding consistency between the prediction and actual observation. The findings of online database verification revealed that these four genes were significantly highly expressed in tumors. Also, we observed their significant correlations with some immune cells and also potential correlations with some drugs. Conclusion. We constructed a 4-RBPs-based prognostic signature to predict the prognosis of BRCA patients, and it has the potential for treating and diagnosing BRCA.

The p53 tumor suppressor integrates upstream signals such as DNA damage and active oncogenes to initiate cell cycle arrest or apoptosis. This response is critical to halting inappropriate growth signals. As such, p53 activity is lost in cancer. In melanoma, however, the p53 gene is intact in a reported 94% of human cases. Rather than direct mutation, p53 is held inactive through interaction with inhibitory proteins. Here, we examine the expression of the two primary inhibitors of p53, MDM2 and MDM4, in genomic databases and biopsy specimens. We find that MDM4 is frequently overexpressed. Moreover, changes in splicing of MDM4 occur frequently and early in melanomagenesis. These changes in splicing must be considered in the design of therapeutic inhibitors of the MDM2/4 proteins for melanoma.

Thymic epithelial cells (TEC) control the selection of a T cell repertoire reactive to pathogens but tolerant of self. This process is known to involve the promiscuous expression of virtually the entire protein-coding gene repertoire, but the extent to which TEC recapitulate peripheral isoforms, and the mechanisms by which they do so, remain largely unknown. We performed the first assembly-based transcriptomic census of transcript structures and splicing factor (SF) expression in mouse medullary TEC (mTEC) and 21 peripheral tissues. Mature mTEC expressed 60.1% of all protein-coding transcripts, more than was detected in any of the peripheral tissues. However, for genes with tissue-restricted expression, mTEC produced fewer isoforms than did the relevant peripheral tissues. Analysis of exon inclusion revealed an absence of brain-specific microexons in mTEC. We did not find unusual numbers of novel transcripts in TEC, and we show that Aire, the facilitator of promiscuous gene expression, promotes the generation of long “classical” transcripts (with 5′ and 3′ UTRs) but has only a limited impact on alternative splicing in mTEC. Comprehensive assessment of SF expression in mTEC identified a small set of nonpromiscuously expressed SF genes, among which we confirmed RBFOX to be present with AIRE in mTEC nuclei. Using a conditional loss-of-function approach, we show that Rbfox2 promotes mTEC development and regulates the alternative splicing of promiscuously expressed genes. These data indicate that TEC recommission a small number of peripheral SFs, including members of the RBFOX family, to generate a broad but selective representation of the peripheral splice isoform repertoire., Genome Research, 31 (11), ISSN:1088-9051, ISSN:1549-5469

The basic-helix-loop-helix (bHLH) transcription factors superfamily is present in eukaryotes and has been widely characterized in several species. The study of these proteins has focused on processes involved in plant growth and development. However, it has also been shown that bHLH plays a relevant role in processes of response to abiotic stress. In this study, the bHLH family in Eucalyptus was systematically characterized for the first time, using the Eucalyptus grandis genome. 153 bHLH genes in the E. grandis genome were identified, which were characterized phylogenetically and structurally. Then, a cold-tolerant hybrid of Eucalyptus (Eucalyptus nitens × Eucalyptus globulus) to generate a reference assembly was used to study the alternative splicing events produced in the cold response. 44 bHLH candidate genes to generate alternative splicing events in the hybrid were identified. These genes had exons differentially expressed in cold acclimation conditions. The splicing event prediction was further validated by qRT-PCR in three bHLH genes. Through this experiment was possible to determine that two of them were differentially expressed, so that new isoforms of these genes were being produced under cold acclimation conditions.

Background Abnormal alternative splicing is frequently associated with carcinogenesis. In B-cell acute lymphoblastic leukemia (B-ALL), double homeobox 4 fused with immunoglobulin heavy chain (DUX4/IGH) can lead to the aberrant production of E-26 transformation-specific family related gene abnormal transcript (ERGalt ) and other splicing variants. However, the molecular mechanism underpinning this process remains elusive. Here, we aimed to know how DUX4/IGH triggers abnormal splicing in leukemia. Methods The differential intron retention analysis was conducted to identify novel DUX4/IGH-driven splicing in B-ALL patients. X-ray crystallography, small angle X-ray scattering (SAXS), and analytical ultracentrifugation were used to investigate how DUX4/IGH recognize double DUX4 responsive element (DRE)-DRE sites. The ERGalt biogenesis and B-cell differentiation assays were performed to characterize the DUX4/IGH crosslinking activity. To check whether recombination-activating gene 1/2 (RAG1/2) was required for DUX4/IGH-driven splicing, the proximity ligation assay, co-immunoprecipitation, mammalian two hybrid characterizations, in vitro RAG1/2 cleavage, and shRNA knock-down assays were performed. Results We reported previously unrecognized intron retention events in C-type lectin domain family 12, member A abnormal transcript (CLEC12Aalt ) and chromosome 6 open reading frame 89 abnormal transcript (C6orf89alt ), where also harbored repetitive DRE-DRE sites. Supportively, X-ray crystallography and SAXS characterization revealed that DUX4 homeobox domain (HD)1-HD2 might dimerize into a dumbbell-shape trans configuration to crosslink two adjacent DRE sites. Impaired DUX4/IGH-mediated crosslinking abolishes ERGalt , CLEC12Aalt , and C6orf89alt biogenesis, resulting in marked alleviation of its inhibitory effect on B-cell differentiation. Furthermore, we also observed a rare RAG1/2-mediated recombination signal sequence-like DNA edition in DUX4/IGH target genes. Supportively, shRNA knock-down of RAG1/2 in leukemic Reh cells consistently impaired the biogenesis of ERGalt , CLEC12Aalt , and C6orf89alt . Conclusions All these results suggest that DUX4/IGH-driven DNA crosslinking is required for RAG1/2 recruitment onto the double tandem DRE-DRE sites, catalyzing V(D)J-like recombination and oncogenic splicing in acute lymphoblastic leukemia.

BackgroundAlternative splicing (AS) produces transcript variants playing potential roles in proteome diversification and gene expression regulation. AS modulation is thus essential to respond to developmental and environmental stimuli. In grapevine, a better understanding of berry development is crucial for implementing breeding and viticultural strategies allowing adaptation to climate changes. Although profound changes in gene transcription have been shown to occur in the course of berry ripening, no detailed study on splicing modifications during this period has been published so far. We report here on the regulation of gene AS in developing berries of two grapevine (Vitis vinifera L.) varieties, Gewurztraminer (Gw) and Riesling (Ri), showing distinctive phenotypic characteristics. Using the software rMATS, the transcriptomes of berries at four developmental steps, from the green stage to mid-ripening, were analysed in pairwise comparisons between stages and varieties.ResultsA total of 305 differential AS (DAS) events, affecting 258 genes, were identified. Interestingly, 22% of these AS events had not been reported before. Among the 80 genes that underwent the most significant variations during ripening, 22 showed a similar splicing profile in Gw and Ri, which suggests their involvement in berry development. Conversely, 23 genes were subjected to splicing regulation in only one variety. In addition, the ratios of alternative isoforms were different in Gw and Ri for 35 other genes, without any change during ripening. This last result indicates substantial AS differences between the two varieties. Remarkably, 8 AS events were specific to one variety, due to the lack of a splice site in the other variety. Furthermore, the transcription rates of the genes affected by stage-dependent splicing regulation were mostly unchanged, identifying AS modulation as an independent way of shaping the transcriptome.ConclusionsThe analysis of AS profiles in grapevine varieties with contrasting phenotypes revealed some similarity in the regulation of several genes with developmental functions, suggesting their involvement in berry ripening. Additionally, many splicing differences were discovered between the two varieties, that could be linked to phenotypic specificities and distinct adaptive capacities. Together, these findings open perspectives for a better understanding of berry development and for the selection of grapevine genotypes adapted to climate change.

Rett syndrome (RTT) is a severe neurological disorder and a leading cause of intellectual disability in young females. RTT is mainly caused by mutations found in the X-linked gene encoding methyl-CpG binding protein 2 (MeCP2). Despite extensive studies, the molecular mechanism underlying RTT pathogenesis is still poorly understood. Here, we report MeCP2 as a key subunit of a higher-order multiunit protein complex Rbfox/LASR. Defective MeCP2 in RTT mouse models disrupts the assembly of the MeCP2/Rbfox/LASR complex, leading to reduced binding of Rbfox proteins to target pre-mRNAs and aberrant splicing of Nrxns and Nlgn1 critical for synaptic plasticity. We further show that MeCP2 disease mutants display defective condensate properties and fail to promote phase-separated condensates with Rbfox proteins in vitro and in cultured cells. These data link an impaired function of MeCP2 with disease mutation in splicing control to its defective properties in mediating the higher-order assembly of the MeCP2/Rbfox/LASR complex., MeCP2 mutations can cause Rett syndrome, a severe childhood neurological disorder. Here the authors show that MeCP2 mediates the higher-order assembly of a large splicing complex Rbfox/LASR, which is disrupted in the mouse models of Rett syndrome.

Identifying the molecular mechanisms contributing to phenotypic variation in natural populations is a major goal of molecular ecology. However, the multiple regulatory steps between genotype and phenotype mean that many potential mechanisms can lead to trait divergence. To date, the role of transcriptional regulation in local adaptation has received much focus, as we can readily measure mRNA quantity and have a reasonable grasp of how variation in the expression of many protein-coding genes can influence phenotype. Thus, studying the evolution of protein-coding gene mRNA abundance in candidate tissues has led to successes in detecting the molecular mechanisms underlying local adaptation (reviewed by Hill et al., 2021). However, the contribution of differential splicing of precursor mRNA (pre-mRNA) to adaptive differentiation, as well as the loci controlling this variation, remains largely unexplored in wild populations. In their "From the Cover'" article in this issue of Molecular Ecology, Jacobs and Elmer (2021) reanalyse muscle RNA sequencing (RNA-seq) data to quantify the relative contributions of variation in mRNA quantity (differentially expressed "DE" genes) and splice variant identity (differentially spliced "DS" genes) to parallel divergence of wild "benthic" and "pelagic" ecotypes of a salmonid fish, the Arctic charr (Salvelinus alpinus). They found little overlap in the identity and biological functions of DE and DS genes, suggesting that these two regulatory mechanisms act on different cellular traits to complementarily alter organismal phenotype. Furthermore, many DE and DS genes could be mapped to cis-acting QTL, arguing that some of this regulatory divergence is genetically based. DE and DS genes were also more likely to be "hub genes" than their nondivergent counterparts, hinting that this regulatory variation may have a variety of phenotypic effects. The comparison of three independently evolved pairs of benthic and pelagic charr uncovered greater than expected parallelism in both expression and splicing between ecotypes across different lakes, supporting a role for these molecular phenotypes in adaptive divergence. Overall, the findings of Jacobs and Elmer (2021) highlight the importance of alternative splicing as a potential mechanism underlying local adaptation and provide a framework for others hoping to make the most of their RNA-seq data.

Kabuki syndrome (KS) is an autosomal dominant genetic disorder in which most cases are caused by de novo mutations. KS type 1 is caused by mutations in KMT2D (OMIM: #147920) and is more common. KS type 2 is caused by mutations in KDM6A (OMIM: #300867). Both genes encode proteins that modify histones and are involved in epigenetic regulation. The enzyme histone-lysine N-methyltransferase 2D, the product of KMT2D, is expressed in most adult tissues and is essential for early embryonic development. The main clinical manifestations of KS include dysmorphic facial features, such as elongated palpebral fissures, eversion of the lateral third of the lower eyelids, and short nasal columella with a broad and depressed nasal tip. Additionally, patients also present with skeletal abnormalities, dermatoglyphic features, mild-to-moderate intellectual disability, hearing loss, and postnatal growth deficiency. We describe an 11-year-old girl from Colombia, who presented with characteristic clinical signs of KS. Genetic studies showed a KMT2D intronic variant (KMT2D NM_003482.3: c.511‐2A> T) as a cause of KS.

Exonic (i.e. coding) variants in genes associated with disease can exert pathogenic effects both at the protein and mRNA level, either by altering the amino acid sequence or by affecting pre-mRNA splicing. The latter is often neglected due to the lack of RNA analyses in genetic diagnostic testing. In this study we considered both pathomechanisms and performed a comprehensive analysis of nine exonic nucleotide changes in OPA1, which is the major gene underlying autosomal dominant optic atrophy (DOA) and is characterized by pronounced allelic heterogeneity. We focused on the GTPase-encoding domain of OPA1, which harbors most of the missense variants associated with DOA. Given that the consensus splice sites extend into the exons, we chose a split codon, namely codon 438, for our analyses. Variants at this codon are the second most common cause of disease in our large cohort of DOA patients harboring disease-causing variants in OPA1. In silico splice predictions, heterologous splice assays, analysis of patient’s RNA when available, and protein modeling revealed different molecular outcomes for variants at codon 438. The wildtype aspartate residue at amino acid position 438 is directly involved in the dimerization of OPA1 monomers. We found that six amino acid substitutions at codon 438 (i.e. all substitutions of the first and second nucleotide of the codon) destabilized dimerization while only substitutions of the first nucleotide of the codon caused exon skipping. Our study highlights the value of combining RNA analysis and protein modeling approaches to accurately assign patients to future precision therapies.

Background: A mutation/deletion involving donor or acceptor sites for exon 14 results in splicing out of exon 14 of the mesenchymal epithelial transition (MET) gene and is known as “MET exon 14 skipping” (ΔMET14). The two recent approvals with substantial objective responses and improved progression-free survival to MET inhibitors namely capmatinib and tepotinib necessitate the identification of this alteration upfront. We herein describe our experience of ΔMET14 detection by an mRNA-based assay using polymerase chain reaction followed by fragment sizing. Methods: This is a home brew assay which was developed with the concept that the transcripts from true ΔMET14 will be shorter by ~140 bases than their wild type counterparts. The cases which were called MET exon 14 skipping positive on next-generation sequencing (NGS) were subjected to this assay, along with 13 healthy controls in order to establish the validity for true negatives. Results: Thirteen cases of ΔMET14 mutation were detected on NGS using RNA-based sequencing. Considering NGS as a gold standard, the sizing assay using both gel and capillary electrophoresis that showed 100% specificity for both with concordance rates of 84.6% and 88.2% with NGS, respectively, were obtained. Conclusions: Owing to the cost-effective nature and easy to use procedures, this assay will prove beneficial for small- and medium-sized laboratories where skilled technical personnel and NGS platforms are unavailable.

TP53 is one of the most ubiquitously altered genes in human cancer. The biological impact of rare variants, particularly those located within noncoding regions, remains poorly understood. From interrogation of clinical massively parallel sequencing data from over 55,000 tumors, which included 23,330 tumors with known TP53 mutations, TP53 intron 4 nucleotide substitutions at position c.375+5G were identified in 45 tumors (0.2% of TP53-mutated cancers), comprising cancers of different organ sites. Loss-of-heterozygosity or a second-hit somatic TP53 mutation was observed in 34 of 40 (85%) informative cases. RT-PCR analysis showed the c.375+5G>T variant to be associated with aberrantly spliced TP53 mRNA transcripts with concomitant loss of the normal transcript. Immunohistochemical staining for p53 was performed on a representative subset of tumors with TP53 c.375+5G variants (n = 14), all of which showed loss of protein expression (100%; n = 13 complete loss, n = 1 subclonal loss). Our data are consistent with classification of TP53 c.375+5G variants as deleterious intronic mutations that interfere with proper mRNA splicing, ultimately resulting in loss of expression of functional p53 protein. The clinical scenario of a tumor with loss of p53 immunohistochemical staining, yet lacking a detectable TP53 exonic mutation, should therefore prompt consideration of splice-altering intronic variants.

Intronic polyadenylation (IpA) usually leads to changes in the coding region of an mRNA, and its implication in diseases has been recognized, although at its very beginning status. Conveniently and accurately identifying IpA is of great importance for further evaluating its biological significance. Here, we developed IPAFinder, a bioinformatic method for the de novo identification of intronic poly(A) sites and their dynamic changes from standard RNA-seq data. Applying IPAFinder to 256 pan-cancer tumor/normal pairs across six tumor types, we discovered 490 recurrent dynamically changed IpA events, some of which are novel and derived from cancer-associated genes such as TSC1, SPERD2, and CCND2. Furthermore, IPAFinder revealed that IpA could be regulated by factors related to splicing and m6A modification. In summary, IPAFinder enables the global discovery and characterization of biologically regulated IpA with standard RNA-seq data and should reveal the biological significance of IpA in various processes.

Many gene expression quantitative trait locus (eQTL) studies have published their summary statistics, which can be used to gain insight into complex human traits by downstream analyses, such as fine mapping and co-localization. However, technical differences between these datasets are a barrier to their widespread use. Consequently, target genes for most genome-wide association study (GWAS) signals have still not been identified. In the present study, we present the eQTL Catalogue (https://www.ebi.ac.uk/eqtl), a resource of quality-controlled, uniformly re-computed gene expression and splicing QTLs from 21 studies. We find that, for matching cell types and tissues, the eQTL effect sizes are highly reproducible between studies. Although most QTLs were shared between most bulk tissues, we identified a greater diversity of cell-type-specific QTLs from purified cell types, a subset of which also manifested as new disease co-localizations. Our summary statistics are freely available to enable the systematic interpretation of human GWAS associations across many cell types and tissues., The eQTL Catalogue provides uniform processing of 21 eQTL studies, allowing the identification of highly reproducible eQTLs affecting whole gene and transcript expression levels.

Somatic single nucleotide variants (SNVs) in cancer genome affect gene expression through various mechanisms depending on their genomic location. While somatic SNVs near canonical splice sites have been reported to cause abnormal splicing of cancer-related genes, whether these SNVs can affect gene expression through other mechanisms remains an open question. Here, we analyzed RNA sequencing and exome data from 4,998 cancer patients covering ten cancer types and identified 152 somatic SNVs near splice sites that were associated with abnormal intronic polyadenylation (IPA). IPA-associated somatic variants favored the localization near the donor splice sites compared to the acceptor splice sites. A proportion of SNV-associated IPA events overlapped with premature cleavage and polyadenylation events triggered by U1 small nuclear ribonucleoproteins (snRNP) inhibition. GC content, intron length and polyadenylation signal were three genomic features that differentiated between SNV-associated IPA and intron retention. Notably, IPA-associated SNVs were enriched in tumor suppressor genes (TSGs), including the well-known TSGs such as PTEN and CDH1 with recurrent SNV-associated IPA events. Minigene assay confirmed that SNVs from PTEN, CDH1, VEGFA, GRHL2, CUL3 and WWC2 could lead to IPA. This work reveals that IPA acts as a novel mechanism explaining the functional consequence of somatic SNVs in human cancer.

Post-transcriptional regulation of RNA processing (RNAP), including splicing and alternative polyadenylation (APA), controls eukaryotic gene function. Conservative estimates based on bulk tissue studies conclude that at least 50% of mammalian genes undergo APA. Single-cell RNA sequencing (scRNA-seq) could enable a near complete estimate of the extent, function, and regulation of these and other forms of RNA processing. Yet, statistical methods to detect regulated RNAP are limited in their detection power because they suffer from reliance on (a) incomplete annotations of 39 untranslated regions (39 UTRs), (b) peak calling heuristics, (c) analysis based on measurements collapsed over all cells in a cell type (pseudobulking), or (d) APA-specific detection. Here, we introduce ReadZS, a computationally-efficient, and annotation-free statistical approach to identify regulated RNAP, including but not limited to APA, in single cells. ReadZS rediscovers and substantially extends the scope of known cell type-specific RNAP in the human lung and during human spermatogenesis. The unique single-cell resolution and statistical properties of ReadZS enable discovery of new evolutionarily conserved, developmentally regulated RNAP and subpopulations of lung-resident macrophages, homogenous by gene expression alone.

Quaking (QKI) controls RNA metabolism in many biological processes including innate immunity, where its roles remain incompletely understood. To illuminate these roles, we performed genome scale transcriptome profiling in QKI knockout cells with or without poly(I:C) transfection, a double-stranded RNA analog that mimics viral infection. Analysis of RNA-sequencing data shows that QKI knockout upregulates genes induced by interferons, suggesting that QKI is an immune suppressor. Furthermore, differential splicing analysis shows that QKI primarily controls cassette exons, and among these events, we noted that QKI silences splicing of the extra domain A (EDA) exon in fibronectin (FN1) transcripts. QKI knockout results in elevated production and secretion of FN1-EDA protein, which is a known activator of interferons. Consistent with an upregulation of the interferon response in QKI knockout cells, our results show reduced production of dengue virus-2 and Japanese encephalitis virus in these cells. In conclusion, we demonstrate that QKI downregulates the interferon system and attenuates the antiviral state.

Clinically used opioids, such as morphine, activate the mu opioid receptor (MOR) encoded by Opioid Receptor Mu 1 (OPRM1) gene. Examination of the opioid receptor genes showed that the human OPRM1 pre-mRNA undergoes extensive alternative splicing events and capable of expressing 21 isoforms. However, characterization of OPRM1 signaling is generalized, and only one isoform (MOR-1) has been extensively studied. Compounding this issue is the increasing significance of intravenous drug abuse in HIV neuropathogenesis. Here, we investigated the molecular impact of morphine and HIV-1 on regulation of OPRM1 pre-mRNA splicing in in vitro and in vivo models. Our results suggested that morphine treatment specifically induces the alternative splicing of MOR-1X isoform among the other isoforms analyzed in neuronal cells. Interestingly, alternative splicing and expression of MOR-1X isoform was also induced in postmortem brain tissues obtained from people with HIV (PWH). Additionally, treatment of control rats with morphine induced alternative splicing of MOR-1X in the brain regions involved in the reward pathways. More interestingly, HIV-1 transgenic (HIV-1Tg) rats, showed an additive induction of MOR-1X isoform with the exposure to morphine. To further assess the possible role of HIV secretory proteins in alternative splicing of OPRM1 gene, we analyzed the impact of HIV-1 Tat, gp120 and Nef proteins on alternative splicing of MOR-1X isoform. While the Tat and gp120 had no visible effects, treatment of neurons with Nef induced MOR-1X alternative splicing that was comparable to treatment with morphine. Altogether, our results suggest that HIV-1 may alter MOR isoform expression with Nef protein by amplifying the rate of MOR-1X alternative splicing induced by morphine.

In eukaryotes, alternative pre-mRNA splicing allows a single gene to encode different protein isoforms that function in many biological processes, and they are used as biomarkers or therapeutic targets for diseases. Although protein isoforms in the human genome are well annotated, we speculate that some low-abundance protein isoforms may still be under-annotated because most genes have a primary coding product and alternative protein isoforms tend to be under-expressed. A peptide coencoded by a novel exon and an annotated exon separated by an intron is known as a novel junction peptide. In the absence of known transcripts and homologous proteins, traditional whole-genome six-frame translation-based proteogenomics cannot identify novel junction peptides, and it cannot capture novel alternative splice sites. In this article, we first propose a strategy and tool for identifying novel junction peptides, called CJunction, which we then integrate into a proteogenomics process specifically designed for novel protein isoform discovery and apply to the analysis of a deep-coverage HeLa mass spectrometry data set with identifier PXD004452 in ProteomeXchange. We succeeded in identifying and validating three novel protein isoforms of two functionally important genes, NHSL1 (causative gene of Nance-Horan syndrome) and EEF1B2 (translation elongation factor), which validate our hypothesis. These novel protein isoforms have significant sequence differences from the annotated gene-coding products introduced by the novel N-terminal, suggesting that they may play importantly different functions.

Pentatricopeptide repeat (PPR) proteins play important roles in the post-transcriptional modification of organellar RNAs in plants. However, the function of most PPR proteins remains unknown. Here, we characterized the rice (Oryza sativa L.) chlorophyll deficient 4 (cde4) mutant which exhibits an albino phenotype during early leaf development, with decreased chlorophyll contents and abnormal chloroplasts at low-temperature (20°C). Positional cloning revealed that CDE4 encodes a P-type PPR protein localized in chloroplasts. In the cde4 mutant, plastid-encoded polymerase (PEP)-dependent transcript levels were significantly reduced, but transcript levels of nuclear-encoded genes were increased compared to wild-type plants at 20°C. CDE4 directly binds to the transcripts of the chloroplast genes rpl2, ndhA, and ndhB. Intron splicing of these transcripts was defective in the cde4 mutant at 20°C, but was normal at 32°C. Moreover, CDE4 interacts with the guanylate kinase VIRESCENT 2 (V2); overexpression of V2 enhanced CDE4 protein stability, thereby rescuing the cde4 phenotype at 20°C. Our results suggest that CDE4 participates in plastid RNA splicing and plays an important role in rice chloroplast development under low-temperature conditions.

Introduction The prevalence of gene mutations in hemophagocytic lymphohistiocytosis (HLH) varied between studies. Thus far, data on the genetic background of HLH in Vietnamese patients are limited. Methods We recruited 94 HLH patients and analyzed for the 4 genes using Sanger sequencing technology. Results Pathogenic variants were observed in 36 (38.29%) patients, including 27 in UNC13D, 5 in STXBP2, 3 in PRF1, and 2 in STX11 (one patient with digenic variants in both UNC13D and STX11). Monoallelic variants accounted for 77.8% of all cases with mutation. A total of 23 different types of pathogenic variants were documented in the 4 genes tested, including 15 in UNC13D, 3 in PRF1, 3 in STXBP2, and 2 in STX11. Interestingly, the novel splicing variant c.3151G>A in UNC13D was recurrently identified in 8 unrelated patients. Conclusion Vietnamese patients with HLH showed a distinct genetic variant spectrum, in which UNC13D is the predominant genetic lesion associated with HLH.

Aspergillus neoniger NCIM 1400 whose cell-free fraction was earlier established for transglycosylation activity conferred by α-glucosidase gene (agdA), was subjected to sequence analysis. Preliminary results revealed certain dynamics in the intron splicing mechanism, and to ascertain these molecular events, a detailed study was carried. The electrophoresis results from the cDNA portion (B-fragment) of agdA showed multiple bands, indicating the amplification of one or more fragments. The sequence results of cDNA cloned vector revealed the retention type of alternative splicing in the agdA. The splicing mechanism of agdA in NCIM 1400 was compared to different A. niger strains, which harbours agdA orthologues, using PCR. It was observed that effective intron splicing leads to higher α-glucosidase activity from these selected Aspergillus spp. To explore the dynamics of intron retention in A. neoniger NCIM 1400, time-course analysis of intron retention, enzyme activity, and sugar consumption were carried over a period of 168 h of fungal growth. RT-qPCR results revealed that introns retention was not detected during the initial growth phase when the maltose and its hydrolysed product, glucose were consumed. Here we demonstrate that exhaustion of maltose causes increase in retention of introns in the mRNA transcripts of agdA gene, and this could be the possible mode of regulating this gene.

Ab repertoire diversity plays a critical role in the host’s ability to fight pathogens. CDR3 is partially responsible for Ab–Ag binding and is a significant source of diversity in the repertoire. CDR3 diversity is generated during VDJ rearrangement because of gene segment selection, gene segment trimming and splicing, and the addition of nucleotides. We analyzed the Ab repertoire diversity across multiple experiments examining the effects of spaceflight on the Ab repertoire after vaccination. Five datasets from four experiments were analyzed using rank-abundance curves and Shannon indices as measures of diversity. We discovered a trend toward lower diversity as a result of spaceflight but did not find the same decrease in our physiological model of microgravity in either the spleen or bone marrow. However, the bone marrow repertoire showed a reduction in diversity after vaccination. We also detected differences in Shannon indices between experiments and tissues. We did not detect a pattern of CDR3 usage across the experiments. Overall, we were able to find differences in the Ab repertoire diversity across experimental groups and tissues.

Production of biopharmaceuticals relies on the expression of mammalian cDNAs in host organisms. Here we show that the expression of a human cDNA in the moss Physcomitrium patens generates the expected full-length and four additional transcripts due to unexpected splicing. This mRNA splicing results in non-functional protein isoforms, cellular misallocation of the proteins and low product yields. We integrated these results together with the results of our analysis of all 32,926 protein-encoding Physcomitrella genes and their 87,533 annotated transcripts in a web application, physCO, for automatized optimization. A thus optimized cDNA results in about twelve times more protein, which correctly localizes to the ER. An analysis of codon preferences of different production hosts suggests that similar effects occur also in non-plant hosts. We anticipate that the use of our methodology will prevent so far undetected mRNA heterosplicing resulting in maximized functional protein amounts for basic biology and biotechnology., Top et al. report the expression of a human blood-clotting factor IX-encoding cDNA in the moss Physcomitrium patens and demonstrate that in addition to the unspliced full-length transcript, the splicing machinery produced several different transcripts, and thus different protein isoforms, which the authors designate as heterosplicing. The authors present a web application, physCO, to automatize the process of P. patens-specific codon optimization, a tool that can assist in maximized functional protein amounts for biotechnology.

Chloroplast genomes in land plants include approximately 20 intron-containing genes. Most of the introns are similar to the group II introns found in fungi, algae and some bacteria, but no self-splicing has been reported. To analyze splicing reactions in chloroplasts, we developed a tobacco (Nicotiana tabacum) chloroplast-based in vitro system. We optimized the splicing reaction using atpF precursor messenger RNA (pre-mRNA). Our system requires a high ATP concentration, whereas ATP is not necessary for self-splicing group II introns. Self-splicing group II introns possess two exon-binding sites (EBS1 and 2) complementary to two intron-binding sites (IBS1 and 2) in the 3′ end of 5′ exons, which are involved in 5′ splice-site selection. Using our in vitro system and atpF pre-mRNA, we analyzed short sequences corresponding to the above EBSs and IBSs. Mutation analyses revealed that EBS1–IBS1 pairing is essential, while EBS2–IBS2 pairing is important but not crucial for splicing. The first 3′ exon nucleotide determines the 3′ splice sites of self-splicing introns. However, mutations to this nucleotide in atpF pre-mRNA did not affect splicing. This result suggests that the mechanism underlying chloroplast pre-mRNA splicing differs partly from that mediating the self-splicing of group II introns.

Intron splicing is an essential event in post-transcriptional RNA processing in plant mitochondria, which requires the participation of diverse nuclear-encoded splicing factors. However, it is presently unclear how these proteins cooperatively take part in the splicing of specific introns. In this study, we characterized a nuclear-encoded mitochondrial P-type pentatricopeptide repeat (PPR) protein named EMP603. This protein is essential for splicing of intron 2 in the Nad1 gene and interacts with the mitochondria-localized DEAD-box RNA helicase PMH2-5140, the RAD52-like proteins ODB1-0814 and ODB1-5061, and the CRM domain-containing protein Zm-mCSF1. Further study revealed that the N-terminal region of EMP603 interacts with the DEAD-box of PMH2-5140, the CRM domain of Zm-mCSF1, and OBD1-5061, but not with OBD1-0814, whereas the PPR domain of EMP603 can interact with ODB1-0814, ODB1-5061, and PMH2-5140, but not with Zm-mCSF1. Defects in EMP603 severely disrupt the assembly and activity of mitochondrial complex I, leading to impaired mitochondrial function, and delayed seed development. The interactions revealed between EMP603 and PMH2-5140, ODB1-0814, ODB1-5061, and Zm-mCSF1 indicate a possible involvement of a dynamic ‘spliceosome-like’ complex in intron splicing, and may accelerate the elucidation of the intron splicing mechanism in plant mitochondria.

Comprehensive characterization of differentially spliced RNA transcripts with nanopore sequencing is limited by bioinformatics tools that are reliant on existing annotations. We have developed FLAME, a bioinformatics pipeline for alternative splicing analysis of gene-specific or transcriptome-wide long-read sequencing data. FLAME is a Python-based tool aimed at providing comprehensible quantification of full-length splice variants, reliable de novo recognition of splice sites and exons, and representation of consecutive exon connectivity in the form of a weighted adjacency matrix. Notably, this workflow circumvents issues related to inadequate reference annotations and allows for incorporation of short-read sequencing data to improve the confidence of nanopore sequencing reads. In this study, the Epstein-Barr virus long noncoding RNA RPMS1 was used to demonstrate the utility of the pipeline. RPMS1 is ubiquitously expressed in Epstein-Barr virus associated cancer and known to undergo ample differential splicing. To fully resolve the RPMS1 spliceome, we combined gene-specific nanopore sequencing reads from a primary gastric adenocarcinoma and a nasopharyngeal carcinoma cell line with matched publicly available short-read sequencing data sets. All previously reported splice variants, including putative ORFs, were detected using FLAME. In addition, 32 novel exons, including two intron retentions and a cassette exon, were discovered within the RPMS1 gene.

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.

Nonsense-mediated RNA decay (NMD) is recognized as an RNA surveillance pathway that targets aberrant mRNAs with premature translation termination codons (PTC) for degradation, however, its molecular mechanisms and roles in health and disease remain incompletely understood. In this study, we developed a novel reporter system to accurately measure NMD activity in individual cells. A genome-wide CRISPR-Cas9 knockout screen using this reporter system identified novel NMD-promoting factors, including multiple components of the SF3B complex and other U2 spliceosome factors. Interestingly, cells with mutations in the spliceosome genes SF3B1 and U2AF1, which are commonly found in myelodysplastic syndrome (MDS) and cancers, have overall attenuated NMD activity. Compared with wild-type (WT) cells, SF3B1- and U2AF1-mutant cells were more sensitive to NMD inhibition, a phenotype that is accompanied by elevated DNA replication obstruction, DNA damage, and chromosomal instability. Remarkably, the sensitivity of spliceosome mutant cells to NMD inhibition was rescued by overexpression of RNase H1, which removes R-loops in the genome. Together, these findings shed new light on the functional interplay between NMD and RNA splicing and suggest a novel synthetic lethal strategy for the treatment of MDS and cancers with spliceosome mutations. Significance: This study has developed a novel NMD reporter system and identified a potential therapeutic approach of targeting the NMD pathway to treat cancer with spliceosome gene mutations.