Alternative splicing in disease and therapy

TREM2 variants that cause early dementia and increase Alzheimer's disease risk affect gene splicing

Loss-of-function variants in the triggering receptor expressed on myeloid cells 2 (TREM2) are responsible for a spectrum of neurodegenerative disorders. In the homozygous state, they cause severe pathologies with early onset dementia, such as Nasu-Hakola disease and behavioural variants of frontotemporal dementia (FTD), whereas heterozygous variants increase the risk of late-onset Alzheimer's disease (AD) and FTD. For over half of TREM2 variants found in families with recessive early onset dementia, the defect occurs at the transcript level via premature termination codons or aberrant splicing. The remaining variants are missense alterations thought to affect the protein; however, the underlying pathogenic mechanism is less clear. In this work, we tested whether these disease-associated TREM2 variants contribute to the pathology via altered splicing. Variants scored by SpliceAI algorithm were tested by a full-size TREM2 splicing reporter assay in different cell lines. The effect of variants was quantified by qRT-/RT-PCR and western blots. Nanostring nCounter was used to measure TREM2 RNA in the brains of NHD patients who carried spliceogenic variants. Exon skipping events were analysed from brain RNA-Seq datasets available through the Accelerating Medicines Partnership for Alzheimer's Disease Consortium. We found that for some Nasu-Hakola disease and early onset FTD-causing variants, splicing defects were the primary cause (D134G) or likely contributor to pathogenicity (V126G and K186N). Similar but milder effects on splicing of exons 2 and 3 were demonstrated for A130V, L133L and R136W enriched in patients with dementia. Moreover, the two most frequent missense variants associated with AD/FTD risk in European and African ancestries (R62H, 1% in Caucasians and T96K, 12% in Africans) had splicing defects via excessive skipping of exon 2 and overproduction of a potentially antagonistic TREM2 protein isoform. The effect of R62H on exon 2 skipping was confirmed in three independent brain RNA-Seq datasets. Our findings revealed an unanticipated complexity of pathogenic variation in TREM2, in which effects on post-transcriptional gene regulation and protein function often coexist. This necessitates the inclusion of computational and experimental analyses of splicing and mRNA processing for a better understanding of genetic variation in disease.

SPLICER: A Highly Efficient Base Editing Toolbox That Enables In Vivo Therapeutic Exon Skipping

Exon skipping technologies enable exclusion of targeted exons from mature mRNA transcripts, which has broad applications in molecular biology, medicine, and biotechnology. Existing exon skipping techniques include antisense oligonucleotides, targetable nucleases, and base editors, which, while effective for specific applications at some target exons, remain hindered by shortcomings, including transient effects for oligonucleotides, genotoxicity for nucleases and inconsistent exon skipping for base editors. To overcome these limitations, we created SPLICER, a toolbox of next-generation base editors consisting of near-PAMless Cas9 nickase variants fused to adenosine or cytosine deaminases for the simultaneous editing of splice acceptor (SA) and splice donor (SD) sequences. Synchronized SA and SD editing with SPLICER improves exon skipping, reduces aberrant outcomes, including cryptic splicing and intron retention, and enables skipping of exons refractory to single splice-site editing. To demonstrate the therapeutic potential of SPLICER, we targeted APP exon 17, which encodes the amino acid residues that are cleaved to form the Aβ plaques in Alzheimer's disease. SPLICER reduced the formation of Aβ42 peptides in vitro and enabled efficient exon skipping in a mouse model of Alzheimer's disease. Overall, SPLICER is a widely applicable and efficient toolbox for exon skipping with broad therapeutic applications.

Comprehensive splicing analysis of the alternatively spliced CHEK2 exons 8 and 10 reveals three enhancer/silencer-rich regions and 38 spliceogenic variants

Splicing is controlled by a large set of regulatory elements (SREs) including splicing enhancers and silencers, which are involved in exon recognition. Variants at these motifs may dysregulate splicing and trigger loss-of-function transcripts associated with disease. Our goal here was to study the alternatively spliced exons 8 and 10 of the breast cancer susceptibility gene CHEK2. For this purpose, we used a previously published minigene with exons 6-10 that produced the expected minigene full-length transcript and replicated the naturally occurring events of exon 8 [Δ(E8)] and exon 10 [Δ(E10)] skipping. We then introduced 12 internal microdeletions of exons 8 and 10 by mutagenesis in order to map SRE-rich intervals by splicing assays in MCF-7 cells. We identified three minimal (10-, 11-, 15-nt) regions essential for exon recognition: c.863_877del [ex8, Δ(E8): 75%] and c.1073_1083del and c.1083_1092del [ex10, Δ(E10): 97% and 62%, respectively]. Then 87 variants found within these intervals were introduced into the wild-type minigene and tested functionally. Thirty-eight of them (44%) impaired splicing, four of which (c.883G>A, c.883G>T, c.884A>T, and c.1080G>T) induced negligible amounts (<5%) of the minigene full-length transcript. Another six variants (c.886G>A, c.886G>T, c.1075G>A, c.1075G>T, c.1076A>T, and c.1078G>T) showed significantly strong impacts (20-50% of the minigene full-length transcript). Thirty-three of the 38 spliceogenic variants were annotated as missense, three as nonsense, and two as synonymous, underlying the fact that any exonic change is capable of disrupting splicing. Moreover, c.883G>A, c.883G>T, and c.884A>T were classified as pathogenic/likely pathogenic variants according to ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based criteria. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Systematic characterization of protein structural features of alternative splicing isoforms using AlphaFold 2

Alternative splicing is an important cellular process in eukaryotes, altering pre-mRNA to yield multiple protein isoforms from a single gene. However, our understanding of the impact of alternative splicing events on protein structures is currently constrained by a lack of sufficient protein structural data. To address this limitation, we employed AlphaFold 2, a cutting-edge protein structure prediction tool, to conduct a comprehensive analysis of alternative splicing for approximately 3,000 human genes, providing valuable insights into its impact on the protein structural. Our investigation employed state of the art high-performance computing infrastructure to systematically characterize structural features in alternatively spliced regions and identified changes in protein structure following alternative splicing events. Notably, we found that alternative splicing tends to alter the structure of residues primarily located in coils and beta-sheets. Our research highlighted a significant enrichment of loops and highly exposed residues within human alternatively spliced regions. Specifically, our examination of the Septin-9 protein revealed potential associations between loops and alternative splicing, providing insights into its evolutionary role. Furthermore, our analysis uncovered two missense mutations in the Tau protein that could influence alternative splicing, potentially contributing to the pathogenesis of Alzheimer's disease. In summary, our work, through a thorough statistical analysis of extensive protein structural data, sheds new light on the intricate relationship between alternative splicing, evolution, and human disease.

Oligomeric, phosphorylated, and truncated tau and spliceosome pathology within the entorhinal-hippocampal connectome across stages of Alzheimer's disease

Neurofibrillary tangles (NFTs) contain abnormally phosphorylated tau proteins, which spread within components of the medial temporal lobe (MTL) memory circuit in Alzheimer's disease (AD). Here, we used quantitative immunohistochemistry to determine the density of posttranslational oligomeric (TOC1 and TNT1), phosphorylated (AT8), and late truncated (TauC3) tau epitopes within the MTL subfields including entorhinal cortex (EC) layer II, subiculum, Cornu Ammonis (CA) subfields, and dentate gyrus (DG) in subjects who died with a clinical diagnosis of no cognitive impairment (NCI), mild cognitive impairment (MCI), and AD. We also examined whether alterations of the nuclear alternative splicing protein, SRSF2, are associated with tau pathology. Although a significant increase in TOC1, TNT1, and AT8 neuron density occurred in the EC in MCI and AD, subicular, DG granule cell, and CA1 and CA3 densities were only significantly higher in AD. TauC3 counts were not different between connectome regions and clinical groups. SRSF2 intensity in AT8-positive cells decreased significantly in all regions independent of the clinical groups examined. CA1 and subicular AT8, TauC3, and oligomeric densities correlated across clinical groups. EC AT8 counts correlated with CA subfields and subicular and DG values across clinical groups. Oligomeric and AT8 CA1, EC, and subicular density correlated with Braak stage. Decreased nuclear SRSF2 in the presence of cytoplasmic phosphorylated tau suggests a dual-hit process in NFT formation within the entorhinal hippocampal connectome during the onset of AD. Although oligomeric and phosphorylated tau follow a stereotypical pattern, clinical disease stage determined density of tau deposition and not anatomic location within the entorhinal-hippocampal connectome.

Heterogeneous splicing patterns resulting from KIF5A variants associated with amyotrophic lateral sclerosis

Single-nucleotide variants (SNVs) in the gene encoding Kinesin Family Member 5A (KIF5A), a neuronal motor protein involved in anterograde transport along microtubules, have been associated with amyotrophic lateral sclerosis (ALS). ALS is a rapidly progressive and fatal neurodegenerative disease that primarily affects the motor neurons. Numerous ALS-associated KIF5A SNVs are clustered near the splice-site junctions of the penultimate exon 27 and are predicted to alter the carboxy-terminal (C-term) cargo-binding domain of KIF5A. Mis-splicing of exon 27, resulting in exon exclusion, is proposed to be the mechanism by which these SNVs cause ALS. Whether all SNVs proximal to exon 27 result in exon exclusion is unclear. To address this question, we designed an in vitro minigene splicing assay in human embryonic kidney 293 cells, which revealed heterogeneous site-specific effects on splicing: only 5' splice-site (5'ss) SNVs resulted in exon skipping. We also quantified splicing in select clustered, regularly interspaced, short palindromic repeats-edited human stem cells, differentiated to motor neurons, and in neuronal tissues from a 5'ss SNV knock-in mouse, which showed the same result. Moreover, the survival of representative 3' splice site, 5'ss, and truncated C-term variant KIF5A (v-KIF5A) motor neurons was severely reduced compared with wild-type motor neurons, and overt morphological changes were apparent. While the total KIF5A mRNA levels were comparable across the cell lines, the total KIF5A protein levels were decreased for v-KIF5A lines, suggesting an impairment of protein synthesis or stability. Thus, despite the heterogeneous effect on ribonucleic acid splicing, KIF5A SNVs similarly reduce the availability of the KIF5A protein, leading to axonal transport defects and motor neuron pathology.

Identification and functional regulation of three alternative splicing isoforms of the fthl27 gene in miiuy croaker, Miichthys miiuy

Alternative splicing is an important basic mechanism for eukaryotes to control gene expression. Different forms of alternative splicing may lead to the production of protein subtypes with different functions, leading to the expansion of protein diversity in organisms, affecting cell production and metabolism, and is even related to the occurrence of many diseases. Many studies have shown that ferritin is usually associated with inflammation, vascular proliferation, and tumors, which is the focus of immunological research. It not only plays a role in iron metabolism and storage in the body, but also plays an important regulatory role in pathways related to immune and inflammatory regulation. However, there are few studies on alternative splicing events of the ferritin gene nowadays. Therefore, this study identified three different splicing isoforms in its ferritin gene fthl27 of Miichthys miiuy through Sanger sequencing, qRT-PCR, and other experimental techniques, and we found that three different splicing isoforms of the ferritin gene fthl27 in M. Miiuy cells showed an upregulation trend after being stimulated by Lipopolysaccharide (LPS) and poly (I: C). The experiment also found that the three isoforms may have different regulatory effects on the expression of inflammatory factors and antiviral immune factors, playing an important role in the innate immune response of fish.

Genome and transcriptome engineering by compact and versatile CRISPR-Cas systems

Comparative genomics has enabled the discovery of tiny clustered regularly interspaced short palindromic repeat (CRISPR) bacterial immune system effectors with enormous potential for manipulating eukaryotic genomes. Recently, smaller Cas proteins, including miniature Cas9, Cas12, and Cas13 proteins, have been identified and validated as efficient genome editing and base editing tools in human cells. The compact size of these novel CRISPR effectors is highly desirable for generating CRISPR-based therapeutic approaches, mainly to overcome in vivo delivery constraints, providing a promising opportunity for editing pathogenic mutations of clinical relevance and knocking down RNAs in human cells without inducing chromosomal insertions or genome alterations. Thus, these tiny CRISPR-Cas systems represent new and highly programmable, specific, and efficient platforms, which expand the CRISPR toolkit for potential therapeutic opportunities.

Identification and verification of plasma protein biomarkers that accurately identify an ectopic pregnancy

BACKGROUND

Differentiating between a normal intrauterine pregnancy (IUP) and abnormal conditions including early pregnancy loss (EPL) or ectopic pregnancy (EP) is a major clinical challenge in early pregnancy. Currently, serial β-human chorionic gonadotropin (β-hCG) and progesterone are the most commonly used plasma biomarkers for evaluating pregnancy prognosis when ultrasound is inconclusive. However, neither biomarker can predict an EP with sufficient and reproducible accuracy. Hence, identification of new plasma biomarkers that can accurately diagnose EP would have great clinical value.

METHODS

Plasma was collected from a discovery cohort of 48 consenting women having an IUP, EPL, or EP. Samples were analyzed by liquid chromatography-tandem mass spectrometry (LC-MS/MS) followed by a label-free proteomics analysis to identify significant changes between pregnancy outcomes. A panel of 14 candidate biomarkers were then verified in an independent cohort of 74 women using absolute quantitation by targeted parallel reaction monitoring mass spectrometry (PRM-MS) which provided the capacity to distinguish between closely related protein isoforms. Logistic regression and Lasso feature selection were used to evaluate the performance of individual biomarkers and panels of multiple biomarkers to predict EP.

RESULTS

A total of 1391 proteins were identified in an unbiased plasma proteome discovery. A number of significant changes (FDR ≤ 5%) were identified when comparing EP vs. non-EP (IUP + EPL). Next, 14 candidate biomarkers (ADAM12, CGA, CGB, ISM2, NOTUM, PAEP, PAPPA, PSG1, PSG2, PSG3, PSG9, PSG11, PSG6/9, and PSG8/1) were verified as being significantly different between EP and non-EP in an independent cohort (FDR ≤ 5%). Using logistic regression models, a risk score for EP was calculated for each subject, and four multiple biomarker logistic models were identified that performed similarly and had higher AUCs than models with single predictors.

CONCLUSIONS

Overall, four multivariable logistic models were identified that had significantly better prediction of having EP than those logistic models with single biomarkers. Model 4 (NOTUM, PAEP, PAPPA, ADAM12) had the highest AUC (0.987) and accuracy (96%). However, because the models are statistically similar, all markers in the four models and other highly correlated markers should be considered in further validation studies.

A brain specific alternatively spliced isoform of nonmuscle myosin IIA lacks its mechanoenzymatic activities

Recent genomic studies reported that 90 to 95% of human genes can undergo alternative splicing, by which multiple isoforms of proteins are synthesized. However, the functional consequences of most of the isoforms are largely unknown. Here, we report a novel alternatively spliced isoform of nonmuscle myosin IIA (NM IIA), called NM IIA2, which is generated by the inclusion of 21 amino acids near the actin-binding region (loop 2) of the head domain of heavy chains. Expression of NM IIA2 is found exclusively in the brain tissue, where it reaches a maximum level at 24 h during the circadian rhythm. The actin-dependent Mg2+-ATPase activity and in vitro motility assays reveal that NM IIA2 lacks its motor activities but localizes with actin filaments in cells. Interestingly, NM IIA2 can also make heterofilaments with NM IIA0 (noninserted isoform of NM IIA) and can retard the in vitro motility of NM IIA, when the two are mixed. Altogether, our findings provide the functional importance of a previously unknown alternatively spliced isoform, NM IIA2, and its potential physiological role in regulating NM IIA activity in the brain.

Alternative splicing in neurodegenerative disease and the promise of RNA therapies

Alternative splicing generates a myriad of RNA products and protein isoforms of different functions from a single gene. Dysregulated alternative splicing has emerged as a new mechanism broadly implicated in the pathogenesis of neurodegenerative diseases such as Alzheimer disease, amyotrophic lateral sclerosis, frontotemporal dementia, Parkinson disease and repeat expansion diseases. Understanding the mechanisms and functional outcomes of abnormal splicing in neurological disorders is vital in developing effective therapies to treat mis-splicing pathology. In this Review, we discuss emerging research and evidence of the roles of alternative splicing defects in major neurodegenerative diseases and summarize the latest advances in RNA-based therapeutic strategies to target these disorders.

Investigating Open Reading Frames in Known and Novel Transcripts using ORFanage

ORFanage is a system designed to assign open reading frames (ORFs) to known and novel gene transcripts while maximizing similarity to annotated proteins. The primary intended use of ORFanage is the identification of ORFs in the assembled results of RNA sequencing experiments, a capability that most transcriptome assembly methods do not have. Our experiments demonstrate how ORFanage can be used to find novel protein variants in RNA-seq datasets, and to improve the annotations of ORFs in tens of thousands of transcript models in the human annotation databases. Through its implementation of a highly accurate and efficient pseudo-alignment algorithm, ORFanage is substantially faster than other ORF annotation methods, enabling its application to very large datasets. When used to analyze transcriptome assemblies, ORFanage can aid in the separation of signal from transcriptional noise and the identification of likely functional transcript variants, ultimately advancing our understanding of biology and medicine.

Mapping splice QTLs reveals distinct transcriptional and post-transcriptional regulatory variation of gene expression and identifies putative alternative splicing variation mediating complex trait variation in pigs

BACKGROUND

Alternative splicing is an important step in gene expression, generating multiple isoforms for the same genes and greatly expanding the diversity of proteomes. Genetic variation in alternative splicing contributes to phenotypic diversity in natural populations. However, the genetic basis of variation in alternative splicing in livestock including pigs remains poorly understood.

RESULTS

In this study, using a Duroc x Pietrain F2 pig population, we performed genome-wide analysis of alternative splicing estimated from stranded RNA-Seq data in skeletal muscle. We characterized the genetic architecture of alternative splicing and compared its basic features with those of overall gene expression. We detected a large number of novel alternative splicing events that were not previously annotated. We found heritability of quantitative alternative splicing scores (percent spliced in or PSI) to be lower than that of overall gene expression. In addition, heritabilities showed little correlation between alternative splicing and overall gene expression. We mapped expression QTLs (eQTLs) and splice QTLs (sQTLs) and found them to be largely non-overlapping. Finally, we integrated sQTL mapping with phenotype QTL (pQTL mapping to identify potential mediator of pQTL effect by alternative splicing.

CONCLUSIONS

Our results suggest that regulatory variation exists at multiple levels and that their genetic controls are distinct, offering opportunities for genetic improvement.

Novel thrombospondin-1 transcript exhibits distinctive expression and activity in thyroid tumorigenesis

Thrombospondin 1 (TSP1) is known for its cell-specific functions in cancer progression, such as proliferation and migration. It contains 22 exons that may potentially produce several different transcripts. Here, we identified TSP1V as a novel TSP1-splicing variant produced by intron retention (IR) in human thyroid cancer cells and tissues. We observed that TSP1V functionally inhibited tumorigenesis contrary to TSP1 wild-type, as identified in vivo and in vitro. These activities of TSP1V are caused by inhibiting phospho-Smad and phospho-focal adhesion kinase. Reverse transcription polymerase chain reaction and minigene experiments revealed that some phytochemicals/non-steroidal anti-inflammatory drugs enhanced IR. We further found that RNA-binding motif protein 5 (RBM5) suppressed IR induced by sulindac sulfide treatment. Additionally, sulindac sulfide reduced phospho-RBM5 levels in a time-dependent manner. Furthermore, trans-chalcone demethylated TSP1V, thereby preventing methyl-CpG-binding protein 2 binding to TSP1V gene. In addition, TSP1V levels were significantly lower in patients with differentiated thyroid carcinoma than in those with benign thyroid nodule, indicating its potential application as a diagnostic biomarker in tumor progression.

Investigating Open Reading Frames in Known and Novel Transcripts using ORFanage

ORFanage is a system designed to assign open reading frames (ORFs) to both known and novel gene transcripts while maximizing similarity to annotated proteins. The primary intended use of ORFanage is the identification of ORFs in the assembled results of RNA sequencing (RNA-seq) experiments, a capability that most transcriptome assembly methods do not have. Our experiments demonstrate how ORFanage can be used to find novel protein variants in RNA-seq datasets, and to improve the annotations of ORFs in tens of thousands of transcript models in the RefSeq and GENCODE human annotation databases. Through its implementation of a highly accurate and efficient pseudo-alignment algorithm, ORFanage is substantially faster than other ORF annotation methods, enabling its application to very large datasets. When used to analyze transcriptome assemblies, ORFanage can aid in the separation of signal from transcriptional noise and the identification of likely functional transcript variants, ultimately advancing our understanding of biology and medicine.

Single nucleotide polymorphism (SNP) in the doublesex ( dsx) gene splice sites and relevance for its alternative splicing in the malaria vector Anopheles gambiae

Background: Malaria burden continues to be significant in tropical regions, and conventional vector control methods are faced with challenges such as insecticide resistance. To overcome these challenges, additional vector control interventions are vital and include modern genetic approaches as well as classical methods like the sterile insect technique (SIT). In the major human malaria vector Anopheles gambiae, a candidate gene favourable for sterility induction is the doublesex ( dsx) gene, involved in mosquitos' somatic sexually dimorphic traits determination. However, the pathways that trigger the signal of dsx gene exon skipping alternative splicing mechanism in anopheline mosquitoes are not well characterized. This study aims to screen the An. gambiae dsx gene splice site sequences for single-nucleotide polymorphisms (SNPs) that could be critical to its alternative splicing. Methods: Variant annotation data from Ag1000G project phase 2 was analysed, in order to identify splice-relevant SNPs within acceptor and donor splice sites of the An. gambiae dsx gene ( Agdsx). Results: SNPs were found in both donor and acceptor sites of the Agdsx. No splice-relevant SNPs were identified in the female-specific intron 4 acceptor site and the corresponding region in males. Two SNPs (rs48712947, rs48712962) were found in the female-specific donor site of exon 5. They were not specific to either males or females as the rs48712947 was found in female mosquitoes from Cameroon, and in both males and females from Burkina Faso. In the other splice sites, the intron 3 acceptor site carried the greatest abundance of SNPs.   Conclusions: There were no gender association between the identified SNPs and the random distribution of these SNPs in mosquito populations. The SNPs in Agdsx splice sites are not critical for the alternative splicing. Other molecular mechanisms should be considered and investigated.

Single nucleotide polymorphism (SNP) in the doublesex (dsx) gene splice sites and relevance for its alternative splicing in the malaria vector Anopheles gambiae

Background: Malaria burden continues to be significant in tropical regions, and conventional vector control methods are faced with challenges such as insecticide resistance. To overcome these challenges, additional vector control interventions are vital and include modern genetic approaches as well as classical methods like the sterile insect technique (SIT). In the major human malaria vector Anopheles gambiae, a candidate gene favourable for sterility induction is the doublesex (dsx) gene, involved in mosquitos’ somatic sexually dimorphic traits determination. However, the pathways that trigger the signal of dsx gene exon skipping alternative splicing mechanism in anopheline mosquitoes are not well characterized. This study aims to screen the An. gambiae dsx gene splice site sequences for single-nucleotide polymorphisms (SNPs) that could be critical to its alternative splicing. Methods: Variant annotation data from Ag1000G project phase 2 was analysed, in order to identify splice-relevant SNPs within acceptor and donor splice sites of the An. gambiae dsx gene (Agdsx). Results: SNPs were found in both donor and acceptor sites of the Agdsx. No splice-relevant SNPs were identified in the female-specific intron 4 acceptor site and the corresponding region in males. Two SNPs (rs48712947, rs48712962) were found in the female-specific donor site of exon 5. They were not specific to either males or females as the rs48712947 was found in female mosquitoes from Cameroon, and in both males and females from Burkina Faso. In the other splice sites, the intron 3 acceptor site carried the greatest abundance of SNPs.   Conclusions: There were no gender association between the identified SNPs and the random distribution of these SNPs in mosquito populations. The SNPs in Agdsx splice sites are not critical for the alternative splicing. Other molecular mechanisms should be considered and investigated.

A Chemical Biology Perspective to Therapeutic Regulation of RNA Splicing in Spinal Muscular Atrophy (SMA)

Manipulation of RNA splicing machinery has emerged as a drug modality. Here, we illustrate the potential of this novel paradigm to correct aberrant splicing events focused on the recent therapeutic advances in spinal muscular atrophy (SMA). SMA is an incurable neuromuscular disorder and at present the primary genetic cause of early infant death. This Review summarizes the exciting journey from the first reported SMA cases to the currently approved splicing-switching treatments, i.e., antisense oligonucleotides and small-molecule modifiers. We emphasize both chemical structures and molecular bases for recognition. We briefly discuss the advantages and disadvantages of these treatments and include the remaining challenges and future directions. Finally, we also predict that these success stories will contribute to further therapies for human diseases by RNA-splicing control.

Gene expression adjustment of inflammatory mechanisms in dairy cow mammary gland parenchyma during host defense against staphylococci

The aim of the study was to identify differences in the expression of splice variants of the PRMT2, LTF and C4A genes in the mammary glands of healthy dairy cows and those infected with staphylococci. An expression study was conducted on 38 Polish Holstein-Friesian dairy cows who were removed from the herd owing to subclinical and chronic mastitic or reproductive issues. Two days before slaughter, milk samples were taken for microbiological analysis and examined for the presence of bacteria. The mammary gland parenchyma samples with a predominance of secretory tissue were taken; these were divided into three groups according to the health status of the mammary gland: H (without pathogenic bacteria in milk), CoNS (with coagulase-negative staphylococci in milk), and CoPS (with coagulase-positive staphylococci in milk). Two of the investigated genes, LTF and C4A, demonstrated variants unequivocally expressed in infected tissue. Two LTF gene variants were found to be associated with cow health status, and with the type of bacteria causing mastitis (CoPS or CoNS). In addition, the expression of C4A isoforms differed with regard to mastitis etiology groups. The comprehensive evaluation of PRMT2 transcript suggested that the gene may also be involved in course of mastitis: two of four PRMT2 transcripts showed increased expression in the mammary gland of the CoPS group compared to controls. The obtained results are important for the knowledge on the etiology of bovine mastitis. The effects of the identified mastitis-relevant splice variants need to be further explored on the protein level to verify the suitability of splice variants and recognize their contribution towards the disease phenotypes and course.

Electroacupuncture at PC6 (Neiguan) Attenuates Angina Pectoris in Rats with Myocardial Ischemia-Reperfusion Injury Through Regulating the Alternative Splicing of the Major Inhibitory Neurotransmitter Receptor GABRG2

Angina pectoris is the most common manifestation of coronary heart disease, causing suffering in patients. Electroacupuncture at PC6 can effectively alleviate angina by regulating the expression of genes, whether the alternative splicing (AS) of genes is affected by acupuncture is still unknown. We established a rat model of myocardial ischemia-reperfusion by coronary artery ligation and confirmed electroacupuncture alleviated the abnormal discharge caused by angina pectoris measured in EMG electromyograms. Analysis of the GSE61840 dataset established that AS events were altered after I/R and regulated by electroacupuncture. I/R decreased the expression of splicing factor Nova1 while electroacupuncture rescued it. Further experiments in dorsal root ganglion cells showed Nova1 regulated the AS of the GABRG2, specifically on its exon 9 where an important phosphorylation site is present. In vivo, results also showed that electroacupuncture can restore AS of GABRG2. Our results proved that electroacupuncture alleviates angina results by regulating alternative splicing.

Nuclear magnetic resonance reveals a two hairpin equilibrium near the 3'-splice site of influenza A segment 7 mRNA that can be shifted by oligonucleotides

Influenza A kills hundreds of thousands of people globally every year and has the potential to generate more severe pandemics. Influenza A's RNA genome and transcriptome provide many potential therapeutic targets. Here, nuclear magnetic resonance (NMR) experiments suggest that one such target could be a hairpin loop of 8 nucleotides in a pseudoknot that sequesters a 3' splice site in canonical pairs until a conformational change releases it into a dynamic 2 × 2-nt internal loop. NMR experiments reveal that the hairpin loop is dynamic and able to bind oligonucleotides as short as pentamers. A 3D NMR structure of the complex contains 4 and likely 5 bp between pentamer and loop. Moreover, a hairpin sequence was discovered that mimics the equilibrium of the influenza hairpin between its structure in the pseudoknot and upon release of the splice site. Oligonucleotide binding shifts the equilibrium completely to the hairpin secondary structure required for pseudoknot folding. The results suggest this hairpin can be used to screen for compounds that stabilize the pseudoknot and potentially reduce splicing.

Beyond Genes: Inclusion of Alternative Splicing and Alternative Polyadenylation to Assess the Genetic Architecture of Predisposition to Voluntary Alcohol Consumption in Brain of the HXB/BXH Recombinant Inbred Rat Panel

Post transcriptional modifications of RNA are powerful mechanisms by which eukaryotes expand their genetic diversity. For instance, researchers estimate that most transcripts in humans undergo alternative splicing and alternative polyadenylation. These splicing events produce distinct RNA molecules, which in turn yield distinct protein isoforms and/or influence RNA stability, translation, nuclear export, and RNA/protein cellular localization. Due to their pervasiveness and impact, we hypothesized that alternative splicing and alternative polyadenylation in brain can contribute to a predisposition for voluntary alcohol consumption. Using the HXB/BXH recombinant inbred rat panel (a subset of the Hybrid Rat Diversity Panel), we generated over one terabyte of brain RNA sequencing data (total RNA) and identified novel splice variants (via StringTie) and alternative polyadenylation sites (via aptardi) to determine the transcriptional landscape in the brains of these animals. After establishing an analysis pipeline to ascertain high quality transcripts, we quantitated transcripts and integrated genotype data to identify candidate transcript coexpression networks and individual candidate transcripts associated with predisposition to voluntary alcohol consumption in the two-bottle choice paradigm. For genes that were previously associated with this trait (e.g., Lrap, Ift81, and P2rx4) (Saba et al., Febs. J., 282, 3556–3578, Saba et al., Genes. Brain. Behav., 20, e12698), we were able to distinguish between transcript variants to provide further information about the specific isoforms related to the trait. We also identified additional candidate transcripts associated with the trait of voluntary alcohol consumption (i.e., isoforms of Mapkapk5, Aldh1a7, and Map3k7). Consistent with our previous work, our results indicate that transcripts and networks related to inflammation and the immune system in brain can be linked to voluntary alcohol consumption. Overall, we have established a pipeline for including the quantitation of alternative splicing and alternative polyadenylation variants in the transcriptome in the analysis of the relationship between the transcriptome and complex traits.

The Splicing of the Mitochondrial Calcium Uniporter Genuine Activator MICU1 Is Driven by RBFOX2 Splicing Factor during Myogenic Differentiation

Alternative splicing, the process by which exons within a pre-mRNA transcript are differentially joined or skipped, is crucial in skeletal muscle since it is required both during myogenesis and in post-natal life to reprogram the transcripts of contractile proteins, metabolic enzymes, and transcription factors in functionally distinct muscle fiber types. The importance of such events is underlined by the numerosity of pathological conditions caused by alternative splicing aberrations. Importantly, many skeletal muscle Ca²⁺ homeostasis genes are also regulated by alternative splicing mechanisms, among which is the Mitochondrial Ca²⁺ Uniporter (MCU) genuine activator MICU1 which regulates MCU opening upon cell stimulation. We have previously shown that murine skeletal muscle MICU1 is subjected to alternative splicing, thereby generating a splice variant-which was named MICU1.1-that confers unique properties to the mitochondrial Ca²⁺ uptake and ensuring sufficient ATP production for muscle contraction. Here we extended the analysis of MICU1 alternative splicing to human tissues, finding two additional splicing variants that were characterized by their ability to regulate mitochondrial Ca²⁺ uptake. Furthermore, we found that MICU1 alternative splicing is induced during myogenesis by the splicing factor RBFOX2. These results highlight the complexity of the alternative splicing mechanisms in skeletal muscle and the regulation of mitochondrial Ca²⁺ among tissues.

Single nucleotide polymorphism (SNP) in the doublesex (dsx) gene splice sites and relevance for its alternative splicing in the malaria vector Anopheles gambiae

Background: The malaria burden continues to be significant in tropical regions, and conventional vector control methods are faced with challenges such as insecticide resistance. To overcome these challenges, additional vector control interventions are vital and include modern genetic approaches as well as classical methods like the sterile insect technique (SIT). In the major human malaria vector Anopheles gambiae, a candidate gene favourable for sterility induction is the doublesex (dsx) gene, encoding somatic sexually dimorphic traits in mosquitoes. However, the mechanism that regulates the expression of this gene in anopheline mosquitoes is poorly understood. This study aimed to screen the An. gambiae dsx gene splice site sequences for single nucleotide polymorphisms (SNPs) that could be critical to its alternative splicing. Methods: Variant annotation data from Ag1000G project phase 2 was analysed, in order to identify splice-relevant SNPs within acceptor and donor splice sites of the An. gambiae dsx gene (Agdsx). Results: SNPs were found in both donor and acceptor sites of the Agdsx. No splice-relevant SNPs were identified in the female-specific intron 4 acceptor site and the corresponding region in males. Two SNPs (rs48712947, rs48712962) were found in the female-specific donor site of exon 5. They were not specific to either males or females as the rs48712947 was found in female mosquitoes from Cameroon, and in both males and females from Burkina Faso. In the other splice sites, the intron 3 acceptor site carried the greatest abundance of SNPs.   Conclusions: There were no gender association between the identified SNPs and the random distribution of these SNPs in mosquito populations. The SNPs in Agdsx splice sites are not critical for the alternative splicing. Other molecular mechanisms should be considered and investigated.

Evaluation of FRET X for single-molecule protein fingerprinting

Single-molecule protein identification is an unrealized concept with potentially ground-breaking applications in biological research. We propose a method called FRET X (Förster Resonance Energy Transfer via DNA eXchange) fingerprinting, in which the FRET efficiency is read out between exchangeable dyes on protein-bound DNA docking strands and accumulated FRET efficiencies constitute the fingerprint for a protein. To evaluate the feasibility of this approach, we simulated fingerprints for hundreds of proteins using a coarse-grained lattice model and experimentally demonstrated FRET X fingerprinting on model peptides. Measured fingerprints are in agreement with our simulations, corroborating the validity of our modeling approach. In a simulated complex mixture of >300 human proteins of which only cysteines, lysines, and arginines were labeled, a support vector machine was able to identify constituents with 95% accuracy. We anticipate that our FRET X fingerprinting approach will form the basis of an analysis tool for targeted proteomics.

•

We propose a FRET-based single-molecule protein identification method

•

Peptides are experimentally distinguishable by their fingerprints

•

Our approach can classify the constituents of complex samples with 95% accuracy

In silico analysis of GATA4 variants demonstrates main contribution to congenital heart disease

Introduction: Congenital heart disease (CHD) is the most common congenital abnormality and the main cause of infant mortality worldwide. Some of the mutations that occur in the GATA4 gene region may result in different types of CHD. Here, we report our in silico analysis of gene variants to determine the effects of the GATA4 gene on the development of CHD.

Methods: Online 1000 Genomes Project, ExAC, gnomAD, GO-ESP, TOPMed, Iranome, GME, ClinVar, and HGMD databases were drawn upon to collect information on all the reported GATA4 variations.The functional importance of the genetic variants was assessed by using SIFT, MutationTaster, CADD,PolyPhen-2, PROVEAN, and GERP prediction tools. Thereafter, network analysis of the GATA4protein via STRING, normal/mutant protein structure prediction via HOPE and I-TASSER, and phylogenetic assessment of the GATA4 sequence alignment via ClustalW were performed.

Results: The most frequent variant was c.874T>C (45.58%), which was reported in Germany.Ventricular septal defect was the most frequent type of CHD. Out of all the reported variants of GATA4,38 variants were pathogenic. A high level of pathogenicity was shown for p.Gly221Arg (CADD score=31), which was further analyzed.

Conclusion: The GATA4 gene plays a significant role in CHD; we, therefore, suggest that it be accorded priority in CHD genetic screening.

AI applications in functional genomics

We review the current applications of artificial intelligence (AI) in functional genomics. The recent explosion of AI follows the remarkable achievements made possible by "deep learning", along with a burst of "big data" that can meet its hunger. Biology is about to overthrow astronomy as the paradigmatic representative of big data producer. This has been made possible by huge advancements in the field of high throughput technologies, applied to determine how the individual components of a biological system work together to accomplish different processes. The disciplines contributing to this bulk of data are collectively known as functional genomics. They consist in studies of: i) the information contained in the DNA (genomics); ii) the modifications that DNA can reversibly undergo (epigenomics); iii) the RNA transcripts originated by a genome (transcriptomics); iv) the ensemble of chemical modifications decorating different types of RNA transcripts (epitranscriptomics); v) the products of protein-coding transcripts (proteomics); and vi) the small molecules produced from cell metabolism (metabolomics) present in an organism or system at a given time, in physiological or pathological conditions. After reviewing main applications of AI in functional genomics, we discuss important accompanying issues, including ethical, legal and economic issues and the importance of explainability.

Maternal cecal microbiota transfer rescues early-life antibiotic-induced enhancement of type 1 diabetes in mice

Early-life antibiotic exposure perturbs the intestinal microbiota and accelerates type 1 diabetes (T1D) development in the NOD mouse model. Here, we found that maternal cecal microbiota transfer (CMT) to NOD mice after early-life antibiotic perturbation largely rescued the induced T1D enhancement. Restoration of the intestinal microbiome was significant and persistent, remediating the antibiotic-depleted diversity, relative abundance of particular taxa, and metabolic pathways. CMT also protected against perturbed metabolites and normalized innate and adaptive immune effectors. CMT restored major patterns of ileal microRNA and histone regulation of gene expression. Further experiments suggest a gut-microbiota-regulated T1D protection mechanism centered on Reg3γ, in an innate intestinal immune network involving CD44, TLR2, and Reg3γ. This regulation affects downstream immunological tone, which may lead to protection against tissue-specific T1D injury.

Programmable System of Cas13-Mediated RNA Modification and Its Biological and Biomedical Applications

Clustered regularly interspaced short palindromic repeats (CRISPR)-Cas13 has drawn broad interest to control gene expression and cell fate at the RNA level in general. Apart from RNA interference mediated by its endonuclease activity, the nuclease-deactivated form of Cas13 further provides a versatile RNA-guided RNA-targeting platform for manipulating kinds of RNA modifications post-transcriptionally. Chemical modifications modulate various aspects of RNA fate, including translation efficiency, alternative splicing, RNA–protein affinity, RNA–RNA interaction, RNA stability and RNA translocation, which ultimately orchestrate cellular biologic activities. This review summarizes the history of the CRISPR-Cas13 system, fundamental components of RNA modifications and the related physiological and pathological functions. We focus on the development of epi-transcriptional editing toolkits based on catalytically inactive Cas13, including RNA Editing for Programmable A to I Replacement (REPAIR) and xABE (adenosine base editor) for adenosine deamination, RNA Editing for Specific C-to-U Exchange (RESCUE) and xCBE (cytidine base editor) for cytidine deamination and dm⁶ACRISPR, as well as the targeted RNA methylation (TRM) and photoactivatable RNA m⁶A editing system using CRISPR-dCas13 (PAMEC) for m⁶A editing. We further highlight the emerging applications of these useful toolkits in cell biology, disease and imaging. Finally, we discuss the potential limitations, such as off-target editing, low editing efficiency and limitation for AAV delivery, and provide possible optimization strategies.

The Functional Impact of Alternative Splicing and Single Nucleotide Polymorphisms in Rheumatoid Arthritis

Advances in genomics and proteomics aid the identification of genes associated with various diseases. Genome-Wide Association Studies (GWAS) have identified multiple loci as risk alleles for susceptibility to Rheumatoid Arthritis (RA). A bisection of RA risk can be attributed to genetic factors. Over 100 associated genetic loci that encompass immune regulatory factors have been found to be linked with RA. Aberrant Single Nucleotide Polymorphisms (SNPs) and alternative splicing mechanisms in such loci induce RA. These aberrations are viewed as potential therapeutic targets due to their association with a multitude of diseases. This review presents a few imperious genes whose alterations can cause severe bone deformities culminating in RA.

Antisense oligonucleotide-based drug development for Cystic Fibrosis patients carrying the 3849+10 kb C-to-T splicing mutation

Background:

Antisense oligonucleotide (ASO)-based drugs for splicing modulation were recently approved for various genetic diseases with unmet need. Here we aimed to develop an ASO-based splicing modulation therapy for Cystic Fibrosis (CF) patients carrying the 3849 + 10 kb C-to-T splicing mutation in the CFTR gene.

Methods:

We have screened, in FRT cells expressing the 3849 + 10 kb C-to-T splicing mutation, ~30 2ʹ-O-Methyl-modified phosphorothioate ASOs, targeted to prevent the recognition and inclusion of a cryptic exon generated due to the mutation. The effect of highly potent ASO candidates on the splicing pattern, protein maturation and CFTR function was further analyzed in well differentiated primary human nasal and bronchial epithelial cells, derived from patients carrying at least one 3849 + 10 kb C-to-T allele.

Results:

A highly potent lead ASO, efficiently delivered by free uptake, was able to significantly increase the level of correctly spliced mRNA and completely restore the CFTR function to wild type levels in cells from a homozygote patient. This ASO led to CFTR function with an average of 43% of wild type levels in cells from various heterozygote patients. Optimized efficiency of the lead ASO was further obtained with 2ʹ-Methoxy Ethyl modification (2ʹMOE).

Conclusion:

The highly efficient splicing modulation and functional correction, achieved by free uptake of the selected lead ASO in various patients, demonstrate the ASO therapeutic potential benefit for CF patients carrying splicing mutations and is aimed to serve as the basis for our current clinical development.

Tubulin-folding cofactor E deficiency promotes vascular dysfunction by increased endoplasmic reticulum stress

AIMS

Assessment of endothelial function in humans by measuring flow-mediated dilation (FMD) risk-stratifies individuals with established cardiovascular disease, whereas its predictive value is limited in primary prevention. We therefore aimed to establish and evaluate novel markers of FMD at the population level.

METHODS AND RESULTS

In order to identify novel targets that were negatively correlated with FMD and investigate their contribution to vascular function, we performed a genome-wide association study (GWAS) of 4175 participants of the population based Gutenberg Health Study. Subsequently, conditional knockout mouse models deleting the gene of interest were generated and characterized. GWAS analysis revealed that single-nucleotide polymorphisms (SNPs) in the tubulin-folding cofactor E (TBCE) gene were negatively correlated with endothelial function and TBCE expression. Vascular smooth muscle cell (VSMC)-targeted TBCE deficiency was associated with endothelial dysfunction, aortic wall hypertrophy, and endoplasmic reticulum (ER) stress-mediated VSMC hyperproliferation in mice, paralleled by calnexin up-regulation and exacerbated by the blood pressure hormone angiotensin II. Treating SMMHC-ERT2-Cre+/-TBCEfl/fl mice with the ER stress modulator tauroursodeoxycholic acid amplified Raptor/Beclin-1-dependent autophagy and reversed vascular dysfunction.

CONCLUSION

TBCE and tubulin homeostasis seem to be novel predictors of vascular function and offer a new drug target to ameliorate ER stress-dependent vascular dysfunction.

Role of promoters in regulating alternative splicing

Alternative splicing (AS) plays a critical role in enhancing proteome complexity in higher eukaryotes. Almost all the multi intron-containing genes undergo AS in humans. Splicing mainly occurs co-transcriptionally, where RNA polymerase II (RNA pol II) plays a crucial role in coordinating transcription and pre-mRNA splicing. Aberrant AS leads to non-functional proteins causative in various pathophysiological conditions such as cancers, neurodegenerative diseases, and muscular dystrophies. Transcription and pre-mRNA splicing are deeply interconnected and can influence each other's functions. Several studies evinced that specific promoters employed by RNA pol II dictate the RNA processing decisions. Promoter-specific recruitment of certain transcriptional factors or transcriptional coactivators influences splicing, and the extent to which these factors affect splicing has not been discussed in detail. Here, in this review, various DNA-binding proteins and their influence on promoter-specific AS are extensively discussed. Besides, this review highlights how the promoter-specific epigenetic changes might regulate AS.

Isoform-resolved correlation analysis between mRNA abundance regulation and protein level degradation

Profiling of biological relationships between different molecular layers dissects regulatory mechanisms that ultimately determine cellular function. To thoroughly assess the role of protein post‐translational turnover, we devised a strategy combining pulse stable isotope‐labeled amino acids in cells (pSILAC), data‐independent acquisition mass spectrometry (DIA‐MS), and a novel data analysis framework that resolves protein degradation rate on the level of mRNA alternative splicing isoforms and isoform groups. We demonstrated our approach by the genome‐wide correlation analysis between mRNA amounts and protein degradation across different strains of HeLa cells that harbor a high grade of gene dosage variation. The dataset revealed that specific biological processes, cellular organelles, spatial compartments of organelles, and individual protein isoforms of the same genes could have distinctive degradation rate. The protein degradation diversity thus dissects the corresponding buffering or concerting protein turnover control across cancer cell lines. The data further indicate that specific mRNA splicing events such as intron retention significantly impact the protein abundance levels. Our findings support the tight association between transcriptome variability and proteostasis and provide a methodological foundation for studying functional protein degradation.

Emerging Methods and Resources for Biological Interrogation of Neuropsychiatric Polygenic Signal

Most neuropsychiatric disorders are highly polygenic, implicating hundreds to thousands of causal genetic variants that span much of the genome. This widespread polygenicity complicates biological understanding because no single variant can explain disease etiology. A strategy to advance biological insight is to seek convergent functions among the large set of variants and map them to a smaller set of disease-relevant genes and pathways. Accordingly, functional genomic resources that provide data on intermediate molecular phenotypes, such as gene-expression and methylation status, can be leveraged to functionally annotate variants and map them to genes. Such molecular quantitative trait locus mappings can be integrated with genome-wide association studies to make sense of the polygenic signal that underlies complex disease. Other resources that provide data on the 3-dimensional structure of chromatin and functional importance of specific genomic regions can be integrated similarly. In addition, mapped genes can then be tested for convergence in biological function, tissue, cell type, or developmental stage. In this review, we provide an overview of functional genomic resources and methods that can be used to interpret results from genome-wide association studies, and we discuss current challenges for biological understanding and future requirements to overcome them.

Alternative splicing: Human disease and quantitative analysis from high-throughput sequencing

Alternative splicing contributes to the majority of protein diversity in higher eukaryotes by allowing one gene to generate multiple distinct protein isoforms. It adds another regulation layer of gene expression. Up to 95% of human multi-exon genes undergo alternative splicing to encode proteins with different functions. Moreover, around 15% of human hereditary diseases and cancers are associated with alternative splicing. Regulation of alternative splicing is attributed to a set of delicate machineries interacting with each other in aid of important biological processes such as cell development and differentiation. Given the importance of alternative splicing events, their accurate mapping and quantification are paramount for downstream analysis, especially for associating disease with alternative splicing. However, deriving accurate isoform expression from high-throughput RNA-seq data remains a challenging task. In this mini-review, we aim to illustrate I) mechanisms and regulation of alternative splicing, II) alternative splicing associated human disease, III) computational tools for the quantification of isoforms and alternative splicing from RNA-seq.

RNA Splicing: Basic Aspects Underlie Antitumor Targeting

BACKGROUND

RNA splicing, a fundamental step in gene expression, is aimed at intron removal and ordering of exons to form the protein's reading frame.

OBJECTIVE

This review is focused on the role of RNA splicing in cancer biology; the splicing abnormalities that lead to tumor progression emerge as targets for therapeutic intervention.

METHODS

We discuss the role of aberrant mRNA splicing in carcinogenesis and drug response.

RESULTS AND CONCLUSION

Pharmacological modulation of RNA splicing sets the stage for treatment approaches in situations where mRNA splicing is a clinically meaningful mechanism of the disease.

Identification and characterization of a novel isoform of heparin cofactor II in human liver

Heparin cofactor II (HCII) is predominantly expressed in the liver and inhibits thrombin in blood plasma to influence the blood coagulation cascade. Its deficiency is associated with arterial thrombosis. Its cleavage by neutrophil elastase produces fragment that helps in neutrophil chemotaxis in the acute inflammatory response in human. In the present study, we have identified a novel alternatively spliced transcript of the HCII gene in human liver. This novel transcript includes an additional novel region in continuation with exon 3 called exon 3b. Exon 3b acts like an alternate last exon, and hence its inclusion in the transcript due to alternative splicing removes exon 4 and encodes for a different C-terminal region to give a novel protein, HCII-N. MD simulations of HCII-N and three-dimensional structure showed a unique 51 amino acid sequence at the C-terminal having unique RCL-like structure. The HCII-N protein purified from bacterial culture showed a protein migrating at lower molecular weight (MW 55 kDa) as compared to native HCII (MW 66 kDa). A fluorescence-based analysis revealed a more compact structure of HCII-N that was in a more hydrophilic environment. The HCII-N protein, however, showed no inhibitory activity against thrombin. Due to large conformational variation observed in comparison with native HCII, HCII-N may have alternate protease specificity or a non-inhibitory role. Western blot of HCII purified from large plasma volume showed the presence of a low MW 59 kDa band with no thrombin activity. This study provides the first evidence of alternatively spliced novel isoform of the HCII gene.

SRSF3: Newly discovered functions and roles in human health and diseases

The serine/arginine rich proteins (SR proteins) are members of a family of RNA binding proteins involved in regulating various features of RNA metabolism, including pre-mRNA constitutive and alternative splicing. In humans, a total of 12 SR splicing factors (SRSFs) namely SRSF1-SRSF12 have been reported. SRSF3, the smallest member of the SR family and the focus of this review, regulates critical steps in mRNA metabolism and has been shown to have mRNA-independent functions as well. Recent studies on SRSF3 have uncovered its role in a wide array of complex biological processes. We have also reviewed the involvement of SRSF3 in disease conditions like cancer, ageing, neurological and cardiac disorders. Finally, we have discussed in detail the autoregulation of SRSF3 and its implications in cancer and commented on the potential of SRSF3 as a therapeutic target, especially in the context of cancer.

Transcriptome profiling reveals the high incidence of hnRNPA1 exon 8 inclusion in chronic myeloid leukemia

Chronic myeloid leukemia (CML) is a malignancy that evolves through a multi-step process. Alternative splicing of several genes has been linked to the progression of the disease, but involvement of alternations in splicing profiles has not been reported. RNA-seq of peripheral blood mononuclear cell (PBMC) samples characterized the differentially expressed and spliced transcripts in five CML chronic phase (CP) and five blast phase (BP) patients, and five healthy controls. Global splicing alteration analysis detected 6474 altered splicing events altered between CML and healthy samples, including many of the previously reported splicing variants and showing a more profound altered splicing deregulation in BP samples. Functional clustering of differentially spliced genes in CP revealed a preferred enrichment relating to cell signaling, while the spliceosome pathway was most overrepresented in BP samples. One differentially spliced spliceosome gene hnRNPA1 showed two splice isoforms; the longer isoform contained exon 8 was preferentially expressed in the BP patients, and the short one excluding exon 8 was specific to healthy controls. Our findings suggested that alternative splicing deregulation played a central role during the progression of CML from CP to BP, and the longer isoform of hnRNPA1 might represent a diagnostic marker and therapeutic target for CML.

RNA modifications in brain tumorigenesis

RNA modifications are emerging as critical regulators in cancer biology, thanks to their ability to influence gene expression and the predominant protein isoforms expressed during cell proliferation, migration, and other pro-oncogenic properties. The reversibility and dynamic nature of post-transcriptional RNA modifications allow cells to quickly adapt to microenvironmental changes. Recent literature has revealed that the deregulation of RNA modifications can promote a plethora of developmental diseases, including tumorigenesis. In this review, we will focus on four key post-transcriptional RNA modifications which have been identified as contributors to the pathogenesis of brain tumors: m6A, alternative polyadenylation, alternative splicing and adenosine to inosine modifications. In addition to the role of RNA modifications in brain tumor progression, we will also discuss potential opportunities to target these processes to improve the dismal prognosis for brain tumors.

Upstream analysis of alternative splicing: a review of computational approaches to predict context-dependent splicing factors

Alternative splicing (AS) has shown to play a pivotal role in the development of diseases, including cancer. Specifically, all the hallmarks of cancer (angiogenesis, cell immortality, avoiding immune system response, etc.) are found to have a counterpart in aberrant splicing of key genes. Identifying the context-specific regulators of splicing provides valuable information to find new biomarkers, as well as to define alternative therapeutic strategies. The computational models to identify these regulators are not trivial and require three conceptual steps: the detection of AS events, the identification of splicing factors that potentially regulate these events and the contextualization of these pieces of information for a specific experiment. In this work, we review the different algorithmic methodologies developed for each of these tasks. Main weaknesses and strengths of the different steps of the pipeline are discussed. Finally, a case study is detailed to help the reader be aware of the potential and limitations of this computational approach.

Interactions of the GABRG2 polymorphisms and childhood trauma on suicide attempt and related traits in depressed patients

BACKGROUND

Previously, we reported that the longest variant of the GABA A receptor γ2 subunit (GABRG2) was associated with suicidal behavior. The present study therefore aimed to determine whether polymorphisms near the alternatively spliced exon of GABRG2 are associated with suicide attempt (SA) and its related traits, and how these variants might interact with reported childhood trauma (CT) in their association with suicidal behavior.

METHODS

We examined 5 single nucleotide polymorphisms (SNPs) of GABRG2. Subjects were suicide Attempters (N = 94), non-suicide attempters (N = 168) with MDD or Bipolar depression, and healthy volunteers (N = 100). Data on demographics, depression severity and suicide attempts were collected. Participants also completed a set of instruments assessing CT, and lifetime aggression and impulsivity.. GABRG2 polymorphisms were genotyped using Sanger sequencing.

RESULTS

Allele A of rs211034 was a protective factor for SA (OR = 0.50 (0.32, 0.80), p = 0.003), and had an interaction effect with emotional neglect (OR = 0.89 (0.82, 0.97), p = 0.006) on depression. One haploblock (consisting of rs211035 and rs211034) was identified within these SNPs, and subjects with haplotype GA (frequency = 7.3%), had lower rate of SA (OR=0.26(0.10, 0.67), p = 0.006). Cognitive impulsivity (OR=1.38)1.24,1.55), p < 0.001), non-planning impulsivity (OR = 1.18 (1.10,1.25), p < 0.001), anger (OR = 1.13 (1.07,1.19), p < 0.001), impulsivity total score (OR = 1.10(1.06,1.15), p < 0.001), hostility (OR = 1.10 (1.04, 1.15), p < 0.001), aggression total score (OR = 1.05 (1.03,1.07), p < 0.001) were associated with depression, meanwhile, hopelessness (OR = 2.18 (1.56, 3.04), p < 0.001) and impulsivity total score (OR = 1.05 (1.02,1.08), p < 0.001) were associated with the risk of SA, adjusted by age and gender. There was no mediation effect in the relationship among CT, gene polymorphisms and SA or depression through increased impulsivity or aggression.

LIMITATIONS

The main limitation of this study is its modest sample size. More genetic variants as well as epigenetic markers should be examined in future studies.

CONCLUSIONS

These findings add to evidence for the involvement of GABRG2 and impulsivity and hopelessness in SA independent from their association with depression. More research is needed on possible mediators of the relationship between GABA-related gene and SA.

A functional network of gastric-cancer-associated splicing events controlled by dysregulated splicing factors

Comprehensive understanding of aberrant splicing in gastric cancer is lacking. We RNA-sequenced 19 gastric tumor–normal pairs and identified 118 high-confidence tumor-associated (TA) alternative splicing events (ASEs) based on high-coverage sequencing and stringent filtering, and also identified 8 differentially expressed splicing factors (SFs). The TA ASEs occurred in genes primarily involved in cytoskeletal organization. We constructed a correlative network between TA ASE splicing ratios and SF expression, replicated it in independent gastric cancer data from The Cancer Genome Atlas and experimentally validated it by knockdown of the nodal SFs (PTBP1, ESRP2 and MBNL1). Each SF knockdown drove splicing alterations in several corresponding TA ASEs and led to alterations in cellular migration consistent with the role of TA ASEs in cytoskeletal organization. We have therefore established a robust network of dysregulated splicing associated with tumor invasion in gastric cancer. Our work is a resource for identifying oncogenic splice forms, SFs and splicing-generated tumor antigens as biomarkers and therapeutic targets.

Real-time detection of mRNA splicing variants with specifically designed reverse-transcription loop-mediated isothermal amplification

Alternative splicing is a ubiquitous and crucial process in cellular processes and has a specific linkage with diseases. To date, developing cost-effective methods with high sensitivity and specificity for detection of splicing variants has been needed. Herein, we report a novel splicing variant assay based on specifically designed reverse-transcription loop-mediated isothermal amplification. After reverse transcribing the splicing variant into cDNA, four DNA primers are specifically designed to recognize six distinct regions. The four DNA primers can hybridize with corresponding sequences for extension and strand displacement DNA synthesis to form stem-loop DNA and then LAMP amplification is started. The proposed method can detect as low as 100 aM splicing variants in real-time fashion with high specificity, showing great potential in biological function and clinical studies.

A novel splicing variant assay is developed based on specifically designed reverse-transcription (RT) loop-mediated isothermal amplification.

Deoxynivalenol globally affects the selection of 3' splice sites in human cells by suppressing the splicing factors, U2AF1 and SF1

Deoxynivalenol (DON) is one of the most abundant mycotoxins and has adverse effects on several biological processes, posing risks of protein synthesis-disrupting effects and ribotoxic response. Therefore, chronic exposure to DON would fundamentally reshape the global expression pattern. Whether DON causes toxic effects on mRNA splicing, a fundamental biological process, remains unclear. In this study, we found that administration of the relative low dosage of DON dramatically changed the alternative splicing of pre-mRNA in HepG2 cells. The overall number of transcripts with aberrant selection of 3' splice sites was significantly increased in DON-exposed HepG2 cells. This effect was further confirmed in two other human cell lines, HEK293 and Caco-2, suggesting that this DON-induced alteration in splicing patterns was universal in human cells. Among these DON-induced changes in alternative splicing, the expression levels of two related splicing factors, SF1 and U2AF1, which are essential for 3' splice site recognitions, were strongly suppressed. Overexpression of either of the two splicing factors strongly alleviated the DON-induced aberrant selection of 3' splice sites. Moreover, SF1 was required for human cell proliferation in DON exposure, and the restoration of SF1 expression partially reinstated the proliferation potential for DON-treated cells. In conclusion, our study suggests that DON, even at a low dosage, has great potential to change gene expression globally by affecting not only protein synthesis but also mRNA processing in human cells.

Identification of a novel splice variant isoform of interferon regulatory factor 10, IRF10, in orange spotted grouper Epinephelus coioides

Interferon regulatory factors (IRFs) are a family of transcription factors essential to the control of antiviral immune response, cell growth, differentiation and apoptosis. IRF10 was originally found in chicken, which was induced by the v-Rel oncoprotein in lymphoid cell lines and involved in the upregulation of major histocompatibility complex (MHC) class I and guanylate-binding protein. In fish, IRF10 plays negative roles in regulation of the interferon (IFN) response. Here, we identified a splice variant of IRF10, named as EcIRF10-SF in orange spotted grouper, which shares the first three exons with the long form (EcIRF10-LF) and retains part of intron 3, creating a premature termination codon. Furthermore, we observed that the EcIRF10-SF exhibits similar expression pattern compared to its native counterparts. Functional studies demonstrate that the two EcIRF10 isoforms repress DrIFNϕ1 and DrIFNϕ3 promoter activity and negatively regulate fish antiviral gene expression. Subcellular localization analysis shows that the amino acids from 57 to 86 within DBD are required for IRF10 nuclear import. Overall, our description of transcript diversification of IRF10 in the grouper provides a coherent framework to further dissect its roles in immune response.

Functional significance of U2AF1 S34F mutations in lung adenocarcinomas

The functional role of U2AF1 mutations in lung adenocarcinomas (LUADs) remains incompletely understood. Here, we report a significant co-occurrence of U2AF1 S34F mutations with ROS1 translocations in LUADs. To characterize this interaction, we profiled effects of S34F on the transcriptome-wide distribution of RNA binding and alternative splicing in cells harboring the ROS1 translocation. Compared to its wild-type counterpart, U2AF1 S34F preferentially binds and modulates splicing of introns containing CAG trinucleotides at their 3' splice junctions. The presence of S34F caused a shift in cross-linking at 3' splice sites, which was significantly associated with alternative splicing of skipped exons. U2AF1 S34F induced expression of genes involved in the epithelial-mesenchymal transition (EMT) and increased tumor cell invasion. Finally, S34F increased splicing of the long over the short SLC34A2-ROS1 isoform, which was also associated with enhanced invasiveness. Taken together, our results suggest a mechanistic interaction between mutant U2AF1 and ROS1 in LUAD.