Alternative splicing: multiple control mechanisms and involvement in human disease

Predicting splicing patterns from the transcription factor binding sites in the promoter with deep learning

BACKGROUND

Alternative splicing is a pivotal mechanism of post-transcriptional modification that contributes to the transcriptome plasticity and proteome diversity in metazoan cells. Although many splicing regulations around the exon/intron regions are known, the relationship between promoter-bound transcription factors and the downstream alternative splicing largely remains unexplored.

RESULTS

In this study, we present computational approaches to unravel the regulatory relationship between promoter-bound transcription factor binding sites (TFBSs) and the splicing patterns. We curated a fine dataset that includes DNase I hypersensitive site sequencing and transcriptomes across fifteen human tissues from ENCODE. Specifically, we proposed different representations of TF binding context and splicing patterns to examine the associations between the promoter and downstream splicing events. While machine learning models demonstrated potential in predicting splicing patterns based on TFBS occupancies, the limitations in the generalization of predicting the splicing forms of singleton genes across diverse tissues was observed with carefully examination using different cross-validation methods. We further investigated the association between alterations in individual TFBS at promoters and shifts in exon splicing efficiency. Our results demonstrate that the convolutional neural network (CNN) models, trained on TF binding changes in the promoters, can predict the changes in splicing patterns. Furthermore, a systemic in silico substitutions analysis on the CNN models highlighted several potential splicing regulators. Notably, using empirical validation using K562 CTCFL shRNA knock-down data, we showed the significant role of CTCFL in splicing regulation.

CONCLUSION

In conclusion, our finding highlights the potential role of promoter-bound TFBSs in influencing the regulation of downstream splicing patterns and provides insights for discovering alternative splicing regulations.

Insights into the ANKRD11 variants and short-stature phenotype through literature review and ClinVar database search

Ankyrin repeat domain containing-protein 11 (ANKRD11), a transcriptional factor predominantly localized in the cell nucleus, plays a crucial role in the expression regulation of key genes by recruiting chromatin remodelers and interacting with specific transcriptional repressors or activators during numerous biological processes. Its pathogenic variants are strongly linked to the pathogenesis and progression of multisystem disorder known as KBG syndrome. With the widespread application of high-throughput DNA sequencing technologies in clinical medicine, numerous pathogenic variants in the ANKRD11 gene have been reported. Patients with KBG syndrome usually exhibit a broad phenotypic spectrum with a variable degree of severity, even if having identical variants. In addition to distinctive dental, craniofacial and neurodevelopmental abnormalities, patients often present with skeletal anomalies, particularly postnatal short stature. The relationship between ANKRD11 variants and short stature is not well-understood, with limited knowledge regarding its occurrence rate or underlying biological mechanism involved. This review aims to provide an updated analysis of the molecular spectrum associated with ANKRD11 variants, investigate the prevalence of the short stature among patients harboring these variants, evaluate the efficacy of recombinant human growth hormone in treating children with short stature and ANKRD11 variants, and explore the biological mechanisms underlying short stature from both scientific and clinical perspectives. Our investigation indicated that frameshift and nonsense were the most frequent types in 583 pathogenic or likely pathogenic variants identified in the ANKRD11 gene. Among the 245 KBGS patients with height data, approximately 50% displayed short stature. Most patients showed a positive response to rhGH therapy, although the number of patients receiving treatment was limited. ANKRD11 deficiency potentially disrupts longitudinal bone growth by affecting the orderly differentiation of growth plate chondrocytes. Our review offers crucial insights into the association between ANKRD11 variants and short stature and provides valuable guidance for precise clinical diagnosis and treatment of patients with KBG syndrome.

Characterization of pre-mRNA Splicing Defects Caused by CLCN5 and OCRL Mutations and Identification of Novel Variants Associated with Dent Disease

Dent disease (DD) is an X-linked renal tubulopathy characterized by low-molecular-weight proteinuria, hypercalciuria, nephrocalcinosis, nephrolithiasis and progressive renal failure. Two-thirds of cases are associated with inactivating variants in the CLCN5 gene (Dent disease 1, DD1) and a few present variants in the OCRL gene (Dent disease 2, DD2). The aim of the present study was to test the effect on the pre-mRNA splicing process of DD variants, described here or in the literature, and describe the clinical and genotypic features of thirteen unrelated patients with suspected DD. All patients presented tubular proteinuria, ten presented hypercalciuria and five had nephrolithiasis or nephrocalcinosis. CLCN5 and OCRL genes were analyzed by Sanger sequencing. Nine patients showed variants in CLCN5 and four in OCRL; eight of these were new. Bioinformatics tools were used to select fifteen variants with a potential effect on pre-mRNA splicing from our patients' group and from the literature, and were experimentally tested using minigene assays. Results showed that three exonic missense mutations and two intronic variants affect the mRNA splicing process. Our findings widen the genotypic spectrum of DD and provide insight into the impact of variants causing DD.

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.

A patient-derived organoid-based study identified an ASO targeting SNORD14E for endometrial cancer through reducing aberrant FOXM1 Expression and β-catenin nuclear accumulation

BACKGROUND

Most of the endometrial cancer (EC) patients are diagnosis in early stage with a good prognosis while the patients with locally advanced recurrent or metastatic result in a poor prognosis. Adjuvant therapy could benefit the prognosis of patients with high-risk factors. Unfortunately, the molecular classification of great prognostic value has not yet reached an agreement and need to be further refined. The present study aims to identify new targets that have prognostic value in EC based on the method of EC patient-derived organ-like organs (PDOs), and further investigate their efficacy and mechanism.

METHODS

The Cancer Genome Atlas (TCGA) database was used to determine SNORD14E expression. The effects of SNORD14E were investigated using CCK8, Transwell, wound-healing assays, and a xenograft model experiment; apoptosis was measured by flow cytometry. Antisense oligonucleotide (ASO) targeting SNORD14E was designed and patient-derived organoids (PDO) models in EC patients was established. A xenograft mouse and PDO model were employed to evaluate the effects of ASO targeting SNORD14E. RNA-seq, Nm-seq, and RNA immunoprecipitation (RIP) experiments were employed to confirm the alternative splicing (AS) and modification induced by SNORD14E. A minigene reporter gene assay was conducted to confirm AS and splicing factors on a variable exon. Actinomycin-d (Act-D) and Reverse Transcription at Low deoxy-ribonucleoside triphosphate concentrations followed by PCR (RTL-P) were utilized to confirm the effects of 2'-O methylation modification on FOXM1.

RESULTS

We found that SNORD14E was overexpressed in EC tissues and patients with high expressed SNORD14E were distributed in the TCGA biomolecular classification subgroups without difference. Further, SNORD14E could reduce disease-free survival (DFS) and recurrence free survival (RFS) of EC patients. SNORD14E promoted proliferation, migration, and invasion and inhibited the apoptosis of EC cells in vitro. ASOs targeting SNORD14E inhibited cell proliferation, migration, invasion while promoted cell apoptosis. ASOs targeting SNORD14E inhibited tumor growth in the xenograft mouse model. TCGA-UCEC database showed that the proportion of patients with high expression of SNORD14E in middle-high risk and high-risk patients recommended by EMSO-ESGO-ESTRO guidelines for adjuvant therapy is more than 50%. Next, we enrolled 8 cases of high-risk and high-risk EC patients according to EMSO-ESGO-ESTRO guidelines and successfully constructed EC-PDOs. ASOs targeting SNORD14E inhibited the EC-PDO growth. Mechanistically, SNORD14E could recognize the mRNA of FOXM1 and recruit SRSF1 to promote the shearing of the variable exon VIIa of FOXM1, resulting in the overexpression of the FOXM1 malignant subtypes FOXM1b and FOXM1c. In addition, SNORD14E modified FOXM1 mRNA with 2`-O-methylation, which prolonged the half-life of FOXM1 mRNA. The nucleus accumulation of β-catenin caused by aberrant FOXM1 expression led to EC progression.

CONCLUSIONS

ASO targeting SNORD14E can be an effective treatment for EC.

Alternative splicing regulates adaptor protein binding, trafficking, and activity of the Vps10p domain receptor SorCS2 in neuronal development

The Vps10p domain receptor SorCS2 is crucial for the development and function of the nervous system and essential for brain-derived neurotrophic factor (BDNF)-induced changes in neuronal morphology and plasticity. SorCS2 regulates the subcellular trafficking of the BDNF signaling receptor TrkB as well as selected neurotransmitter receptors in a manner that is dependent on the SorCS2 intracellular domain (ICD). However, the cellular machinery and adaptor protein (AP) interactions that regulate receptor trafficking via the SorCS2 ICD are unknown. We here identify four splice variants of human SorCS2 differing in the insertion of an acidic cluster motif and/or a serine residue within the ICD. We show that each variant undergoes posttranslational proteolytic processing into a one- or two-chain receptor, giving rise to eight protein isoforms, the expression of which differs between neuronal and nonneuronal tissues and is affected by cellular stressors. We found that the only variants without the serine were able to rescue BDNF-induced branching of SorCS2 knockout hippocampal neurons, while variants without the acidic cluster showed increased interactions with clathrin-associated APs AP-1, AP-2, and AP-3. Using yeast two-hybrid screens, we further discovered that all variants bound dynein light chain Tctex-type 3; however, only variants with an acidic cluster motif bound kinesin light chain 1. Accordingly, splice variants showed markedly different trafficking properties and localized to different subcellular compartments. Taken together, our findings demonstrate the existence of eight functional SorCS2 isoforms with differential capacity for interactions with cytosolic ligands dynein light chain Tctex-type 3 and kinesin light chain 1, which potentially allows cell-type specific SorCS2 trafficking and BDNF signaling.

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.

The role of alternative splicing in lung cancer

Aberrant alternative splicing (AS) events are frequently observed in lung cancer, which can be attributed to aberrant gene AS, alterations in splicing regulatory factors, or changes in splicing regulatory mechanisms. Consequently, the dysregulation of alternative RNA splicing is the fundamental cause of lung cancer. In this review, we have summarized the pivotal role of AS in the development, progression, invasion, metastasis, angiogenesis, and drug resistance of lung cancer. Ultimately, this review emphasizes the potential of AS as biomarkers in lung cancer prognosis and diagnosis, and introduces some applications of AS isoform in the treatment of lung cancer. The comprehension of the AS may provide a glimmer of hope for the eradication of lung cancer.

The identification of protein and RNA interactors of the splicing factor Caper in the adult Drosophila nervous system

Post-transcriptional gene regulation is a fundamental mechanism that helps regulate the development and healthy aging of the nervous system. Mutations that disrupt the function of RNA-binding proteins (RBPs), which regulate post-transcriptional gene regulation, have increasingly been implicated in neurological disorders including amyotrophic lateral sclerosis, Fragile X Syndrome, and spinal muscular atrophy. Interestingly, although the majority of RBPs are expressed widely within diverse tissue types, the nervous system is often particularly sensitive to their dysfunction. It is therefore critical to elucidate how aberrant RNA regulation that results from the dysfunction of ubiquitously expressed RBPs leads to tissue specific pathologies that underlie neurological diseases. The highly conserved RBP and alternative splicing factor Caper is widely expressed throughout development and is required for the development of Drosophila sensory and motor neurons. Furthermore, caper dysfunction results in larval and adult locomotor deficits. Nonetheless, little is known about which proteins interact with Caper, and which RNAs are regulated by Caper. Here we identify proteins that interact with Caper in both neural and muscle tissue, along with neural specific Caper target RNAs. Furthermore, we show that a subset of these Caper-interacting proteins and RNAs genetically interact with caper to regulate Drosophila gravitaxis behavior.

Resurrecting the alternative splicing landscape of archaic hominins using machine learning

Alternative splicing contributes to adaptation and divergence in many species. However, it has not been possible to directly compare splicing between modern and archaic hominins. Here, we unmask the recent evolution of this previously unobservable regulatory mechanism by applying SpliceAI, a machine-learning algorithm that identifies splice-altering variants (SAVs), to high-coverage genomes from three Neanderthals and a Denisovan. We discover 5,950 putative archaic SAVs, of which 2,186 are archaic-specific and 3,607 also occur in modern humans via introgression (244) or shared ancestry (3,520). Archaic-specific SAVs are enriched in genes that contribute to traits potentially relevant to hominin phenotypic divergence, such as the epidermis, respiration and spinal rigidity. Compared to shared SAVs, archaic-specific SAVs occur in sites under weaker selection and are more common in genes with tissue-specific expression. Further underscoring the importance of negative selection on SAVs, Neanderthal lineages with low effective population sizes are enriched for SAVs compared to Denisovan and shared SAVs. Finally, we find that nearly all introgressed SAVs in humans were shared across the three Neanderthals, suggesting that older SAVs were more tolerated in human genomes. Our results reveal the splicing landscape of archaic hominins and identify potential contributions of splicing to phenotypic differences among hominins.

Aberrant Cyclin D1 splicing in cancer: from molecular mechanism to therapeutic modulation

Cyclin D1 (CCND1), a crucial mediator of cell cycle progression, possesses many mutation types with different mutation frequencies in human cancers. The G870A mutation is the most common mutation in CCND1, which produces two isoforms: full-length CCND1a and divergent C-terminal CCND1b. The dysregulation of the CCND1 isoforms is associated with multiple human cancers. Exploring the molecular mechanism of CCND1 isoforms has offer new insight for cancer treatment. On this basis, the alterations of CCND1 gene are described, including amplification, overexpression, and mutation, especially the G870A mutation. Subsequently, we review the characteristics of CCND1 isoforms caused by G870A mutation. Additionally, we summarize cis-regulatory elements, trans-acting factors, and the splice mutation involved in splicing regulation of CCND1. Furthermore, we highlight the function of CCND1 isoforms in cell cycle, invasion, and metastasis in cancers. Importantly, the clinical role of CCND1 isoforms is also discussed, particularly concerning prognosis, chemotherapy, and radiotherapy. Last, emphasis is given to the corrective strategies that modulate the cancerous CCND1 isoforms. Thus, it is highlighting significance of aberrant isoforms of CCND1 as targets for cancer therapy.

The implications of alternative pre-mRNA splicing in cell signal transduction

Cells produce multiple mRNAs through alternative splicing, which ensures proteome diversity. Because most human genes undergo alternative splicing, key components of signal transduction pathways are no exception. Cells regulate various signal transduction pathways, including those associated with cell proliferation, development, differentiation, migration, and apoptosis. Since proteins produced through alternative splicing can exhibit diverse biological functions, splicing regulatory mechanisms affect all signal transduction pathways. Studies have demonstrated that proteins generated by the selective combination of exons encoding important domains can enhance or attenuate signal transduction and can stably and precisely regulate various signal transduction pathways. However, aberrant splicing regulation via genetic mutation or abnormal expression of splicing factors negatively affects signal transduction pathways and is associated with the onset and progression of various diseases, including cancer. In this review, we describe the effects of alternative splicing regulation on major signal transduction pathways and highlight the significance of alternative splicing.

Splicing accuracy varies across human introns, tissues and age

Alternative splicing impacts most multi-exonic human genes. Inaccuracies during this process may have an important role in ageing and disease. Here, we investigated mis-splicing using RNA-sequencing data from ~14K control samples and 42 human body sites, focusing on split reads partially mapping to known transcripts in annotation. We show that mis-splicing occurs at different rates across introns and tissues and that these splicing inaccuracies are primarily affected by the abundance of core components of the spliceosome assembly and its regulators. Using publicly available data on short-hairpin RNA-knockdowns of numerous spliceosomal components and related regulators, we found support for the importance of RNA-binding proteins in mis-splicing. We also demonstrated that age is positively correlated with mis-splicing, and it affects genes implicated in neurodegenerative diseases. This in-depth characterisation of mis-splicing can have important implications for our understanding of the role of splicing inaccuracies in human disease and the interpretation of long-read RNA-sequencing data.

Serine-arginine splicing factor 2 promotes oesophageal cancer progression by regulating alternative splicing of interferon regulatory factor 3

OBJECTIVE

Often, alternative splicing is used by cancer cells to produce or increase proteins that promote growth and survival through alternative splicing. Although RNA-binding proteins are known to regulate alternative splicing events associated with tumorigenesis, their role in oesophageal cancer (EC) has rarely been explored.

METHODS

We analysed the expression pattern of several relatively well characterized splicing regulators on 183 samples from TCGA cohort of oesophageal cancer; the effectiveness of the knockdown of SRSF2 was subsequently verified by immunoblotting; we measured the ability of cells treated with lenti-sh-SRSF2/lenti-sh2-SRSF2 to invade through an extracellular matrix coating by transwell invasion assay; using RNA-seq data to identify its potential target genes; we performed qRT-PCR to detect the changes of exon 2 usage in lenti-sh-SRSF2 transduced KYSE30 cells to determine the possible effect of SRSF2 on splicing regulation of IRF3; RNA Electrophoretic mobility shift assay (RNA-EMSA) was performed by the incubation of purified SRSF2 protein and biotinylated RNA probes; we performed luciferase assay to confirm the effect of SRSF2 on IFN1 promoter activity.

RESULTS

We found upregulation of SRSF2 is correlated with the development of EC; Knock-down of SRSF2 inhibits EC cell proliferation, migration, and invasion; SRSF2 regulates the splicing pattern of IRF3 in EC cells; SRSF2 interacts with exon 2 of IRF3 to regulate its exclusion; SRSF2 inhibits the transcription of IFN1 in EC cells.

CONCLUSION

This study identified a novel regulatory axis involved in EC from the various aspects of splicing regulation.

Identification of Candidate mRNA Isoforms for Prostate Cancer-Risk SNPs Utilizing Iso-eQTL and sQTL Methods

Single nucleotide polymorphisms (SNPs) impacting the alternative splicing (AS) process (sQTLs) or isoform expression (iso-eQTL) are implicated as important cancer regulatory elements. To find the sQTL and iso-eQTL, we retrieved prostate cancer (PrCa) tissue RNA-seq and genotype data originating from 385 PrCa European patients from The Cancer Genome Atlas. We conducted RNA-seq analysis with isoform-based and splice event-based approaches. The MatrixEQTL was used to identify PrCa-associated sQTLs and iso-eQTLs. The overlap between sQTL and iso-eQTL with GWAS loci and those that are differentially expressed between cancer and normal tissue were identified. The cis-acting associations (FDR < 0.05) for PrCa-risk SNPs identified 42, 123, and 90 PrCa-associated cassette exons, intron retention, and mRNA isoforms belonging to 25, 95, and 83 genes, respectively; while assessment of trans-acting association (FDR < 0.05) yielded 59, 65, and 196 PrCa-associated cassette exons, intron retention and mRNA isoforms belonging to 35, 55, and 181 genes, respectively. The results suggest that functional PrCa-associated SNPs can play a role in PrCa genesis by making an important contribution to the dysregulation of AS and, consequently, impacting the expression of the mRNA isoforms.

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.

Transcriptomic analysis to elucidate the effects of high stocking density on grass carp (Ctenopharyngodon idella)

Background

Grass carp (Ctenopharyngodon idella) is one of the most widely cultivated fishes in China. High stocking density can reportedly affect fish growth and immunity. Herein we performed PacBio long-read single-molecule real-time (SMRT) sequencing and Illumina RNA sequencing to evaluate the effects of high stocking density on grass carp transcriptome.

Results

SMRT sequencing led to the identification of 33,773 genes (14,946 known and 18,827 new genes). From the structure analysis, 8,009 genes were detected with alternative splicing events, 10,219 genes showed alternative polyadenylation sites and 15,521 long noncoding RNAs. Further, 1,235, 962, and 213 differentially expressed genes (DEGs) were identified in the intestine, muscle, and brain tissues, respectively. We performed functional enrichment analyses of DEGs, and they were identified to be significantly enriched in nutrient metabolism and immune function. The expression levels of several genes encoding apolipoproteins and activities of enzymes involved in carbohydrate enzymolysis were found to be upregulated in the high stocking density group, indicating that lipid metabolism and carbohydrate decomposition were accelerated. Besides, four isoforms of grass carp major histocompatibility complex class II antigen alpha and beta chains in the aforementioned three tissue was showed at least a 4-fold decrease.

Conclusions

The results suggesting that fish farmed at high stocking densities face issues associated with the metabolism and immune system. To conclude, our results emphasize the importance of maintaining reasonable density in grass carp aquaculture.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07924-4.

Transcription/Replication Conflicts in Tumorigenesis and Their Potential Role as Novel Therapeutic Targets in Multiple Myeloma

Simple Summary

Multiple myeloma is a hematologic cancer characterized by the accumulation of malignant plasma cells in the bone marrow. It remains a mostly incurable disease due to the inability to overcome refractory disease and drug-resistant relapse. Oncogenic transformation of PC in multiple myeloma is thought to occur within the secondary lymphoid organs. However, the precise molecular events leading to myelomagenesis remain obscure. Here, we identified genes involved in the prevention and the resolution of conflicts between the replication and transcription significantly overexpressed during the plasma cell differentiation process and in multiple myeloma cells. We discussed the potential role of these factors in myelomagenesis and myeloma biology. The specific targeting of these factors might constitute a new therapeutic strategy in multiple myeloma.

Abstract

Plasma cells (PCs) have an essential role in humoral immune response by secretion of antibodies, and represent the final stage of B lymphocytes differentiation. During this differentiation, the pre-plasmablastic stage is characterized by highly proliferative cells that start to secrete immunoglobulins (Igs). Thus, replication and transcription must be tightly regulated in these cells to avoid transcription/replication conflicts (TRCs), which could increase replication stress and lead to genomic instability. In this review, we analyzed expression of genes involved in TRCs resolution during B to PC differentiation and identified 41 genes significantly overexpressed in the pre-plasmablastic stage. This illustrates the importance of mechanisms required for adequate processing of TRCs during PCs differentiation. Furthermore, we identified that several of these factors were also found overexpressed in purified PCs from patients with multiple myeloma (MM) compared to normal PCs. Malignant PCs produce high levels of Igs concomitantly with cell cycle deregulation. Therefore, increasing the TRCs occurring in MM cells could represent a potent therapeutic strategy for MM patients. Here, we describe the potential roles of TRCs resolution factors in myelomagenesis and discuss the therapeutic interest of targeting the TRCs resolution machinery in MM.

The upstream 5' splice site remains associated to the transcription machinery during intron synthesis

In the earliest step of spliceosome assembly, the two splice sites flanking an intron are brought into proximity by U1 snRNP and U2AF along with other proteins. The mechanism that facilitates this intron looping is poorly understood. Using a CRISPR interference-based approach to halt RNA polymerase II transcription in the middle of introns in human cells, we discovered that the nascent 5′ splice site base pairs with a U1 snRNA that is tethered to RNA polymerase II during intron synthesis. This association functionally corresponds with splicing outcome, involves bona fide 5′ splice sites and cryptic intronic sites, and occurs transcriptome-wide. Overall, our findings reveal that the upstream 5′ splice sites remain attached to the transcriptional machinery during intron synthesis and are thus brought into proximity of the 3′ splice sites; potentially mediating the rapid splicing of long introns.

We know that most splicing reactions take place co-transcriptionally, but how the transcription machinery facilitate splicing of introns is unknown. Here the authors show that the 5′ splice site remains associated with the transcription machinery during intron synthesis through U1 snRNP, providing a basis for the rapid splicing reaction of introns.

RAD51D Aberrant Splicing in Breast Cancer: Identification of Splicing Regulatory Elements and Minigene-Based Evaluation of 53 DNA Variants

RAD51D loss-of-function variants increase lifetime risk of breast and ovarian cancer. Splicing disruption is a frequent pathogenic mechanism associated with variants in susceptibility genes. Herein, we have assessed the splicing and clinical impact of splice-site and exonic splicing enhancer (ESE) variants identified through the study of ~113,000 women of the BRIDGES cohort. A RAD51D minigene with exons 2-9 was constructed in splicing vector pSAD. Eleven BRIDGES splice-site variants (selected by MaxEntScan) were introduced into the minigene by site-directed mutagenesis and tested in MCF-7 cells. The 11 variants disrupted splicing, collectively generating 25 different aberrant transcripts. All variants but one produced negligible levels (<3.4%) of the full-length (FL) transcript. In addition, ESE elements of the alternative exon 3 were mapped by testing four overlapping exonic microdeletions (≥30-bp), revealing an ESE-rich interval (c.202_235del) with critical sequences for exon 3 recognition that might have been affected by germline variants. Next, 26 BRIDGES variants and 16 artificial exon 3 single-nucleotide substitutions were also assayed. Thirty variants impaired splicing with variable amounts (0-65.1%) of the FL transcript, although only c.202G>A demonstrated a complete aberrant splicing pattern without the FL transcript. On the other hand, c.214T>C increased efficiency of exon 3 recognition, so only the FL transcript was detected (100%). In conclusion, 41 RAD51D spliceogenic variants (28 of which were from the BRIDGES cohort) were identified by minigene assays. We show that minigene-based mapping of ESEs is a powerful approach for identifying ESE hotspots and ESE-disrupting variants. Finally, we have classified nine variants as likely pathogenic according to ACMG/AMP-based guidelines, highlighting the complex relationship between splicing alterations and variant interpretation.

Plant responses to high temperature: a view from pre-mRNA alternative splicing

This review focused on the recent breakthroughs in plant high temperature responses from an alternative splicing angle. With the inevitable global warming, high temperature triggers plants to change their growth and developmental programs for adapting temperature increase. In the past decades, the signaling mechanisms from plant thermo-sensing to downstream transcriptional cascades have been extensively studied. Plenty of elegant review papers have summarized these breakthroughs from signal transduction to cross-talk within plant hormones and environmental cues. Precursor messenger RNA (pre-mRNA) splicing enables plants to produce a series of functional un-related proteins and thus enhances the regulation flexibility. Plants take advantage of this strategy to modulate their proteome diversity under high ambient temperature and elicit developmental plasticity. In this review, we particularly focus on pre-mRNA splicing regulation underlying plant high temperature responses, and will shed new light on the understanding of post-transcriptional regulation on plant growth and development.

Cyclin-Dependent Kinases (CDK) and Their Role in Diseases Development-Review

Cyclin-dependent kinases (CDKs) are involved in many crucial processes, such as cell cycle and transcription, as well as communication, metabolism, and apoptosis. The kinases are organized in a pathway to ensure that, during cell division, each cell accurately replicates its DNA, and ensure its segregation equally between the two daughter cells. Deregulation of any of the stages of the cell cycle or transcription leads to apoptosis but, if uncorrected, can result in a series of diseases, such as cancer, neurodegenerative diseases (Alzheimer’s or Parkinson’s disease), and stroke. This review presents the current state of knowledge about the characteristics of cyclin-dependent kinases as potential pharmacological targets.

Identification of a novel and potent small molecule inhibitor of SRPK1: mechanism of dual inhibition of SRPK1 for the inhibition of cancer progression

Protein kinases are the family of attractive enzyme targets for drug design with relevance to cancer biology. Serine arginine protein kinase 1 (SRPK1) is responsible for the phosphorylation of serine/arginine (SR)-rich proteins. Alternative Splicing Factor/Splicing Factor 2 (ASF/SF2) involved in mRNA editing. ASF/SF2 is over expressed in many cancers and plays crucial roles in the cell survival. Phosphorylation of ASF/SF2 is decisive for its functions in cancer. In search of potential anticancer therapeutic agents for attenuating phosphorylation of ASF/SF2, we have explored specific and potential inhibitors of SRPK1 from natural and drug like compounds databases using in-silico methods. Compound ZINC02154892 (C02) was found to be the most potent inhibitor for SRPK1. In-vitro molecular and cell biology studies have shown C02 as a potent and specific inhibitor of phosphorylation of ASF/SF2 and cell survival in leukemic cell line. Structural analysis of SRPK1 with compound C02 revealed a unique pattern of binding targeting ATP binding site along with inhibiting recruitment of ASF/SF2 by SRPK1. The possibilities of compound C02 to be used as a lead compound paving way for the development of potent and specific inhibitors of SRPK1 for designing of novel potential anticancer inhibitor is inferred from the current studies.

Gene Expression and RNA Splicing Imputation Identifies Novel Candidate Genes Associated with Osteoporosis

CONTEXT

Though genome-wide association studies (GWASs) have identified hundreds of genetic variants associated with osteoporosis related traits, such as bone mineral density (BMD) and fracture, it remains a challenge to interpret their biological functions and underlying biological mechanisms.

OBJECTIVE

Integrate diverse expression quantitative trait loci and splicing quantitative trait loci data with several powerful GWAS datasets to identify novel candidate genes associated with osteoporosis.

DESIGN, SETTING, AND PARTICIPANTS

Here, we conducted a transcriptome-wide association study (TWAS) for total body BMD (TB-BMD) (n = 66 628 for discovery and 7697 for validation) and fracture (53 184 fracture cases and 373 611 controls for discovery and 37 857 cases and 227 116 controls for validation), respectively. We also conducted multi-SNP-based summarized mendelian randomization analysis to further validate our findings.

RESULTS

In total, we detected 88 genes significantly associated with TB-BMD or fracture through expression or ribonucleic acid splicing. Summarized mendelian randomization analysis revealed that 78 of the significant genes may have potential causal effects on TB-BMD or fracture in at least 1 specific tissue. Among them, 64 genes have been reported in previous GWASs or TWASs for osteoporosis, such as ING3, CPED1, and WNT16, as well as 14 novel genes, such as DBF4B, GRN, TMUB2, and UNC93B1.

CONCLUSIONS

Overall, our findings provide novel insights into the pathogenesis mechanisms of osteoporosis and highlight the power of a TWAS to identify and prioritize potential causal genes.

Pan-cancer analysis of alternative splicing regulator heterogeneous nuclear ribonucleoproteins (hnRNPs) family and their prognostic potential

As the most critical alternative splicing regulator, heterogeneous nuclear ribonucleoproteins (hnRNPs) have been reported to be implicated in various aspects of cancer. However, the comprehensive understanding of hnRNPs in cancer is still lacking. The molecular alterations and clinical relevance of hnRNP genes were systematically analysed in 33 cancer types based on next-generation sequence data. The expression, mutation, copy number variation, functional pathways, immune cell correlations and prognostic value of hnRNPs were investigated across different cancer types. HNRNPA1 and HNRNPAB were highly expressed in most tumours. HNRNPM, HNRNPUL1, and HNRNPL showed high mutation frequencies, and most hnRNP genes were frequently mutated in uterine corpus endometrial carcinoma (UCEC). HNRNPA2B1 showed widespread copy number amplification across various cancer types. HNRNPs participated in cancer-related pathways including protein secretion, mitotic spindle, G2/M checkpoint, DNA repair, IL6/JAK/STAT3 signal and coagulation, of which hnRNP genes of HNRNPF, HNRNPH2, HNRNPU and HNRNPUL1 are more likely to be implicated. Significant correlation of hnRNP genes with T help cells, NK cells, CD8 positive T cells and neutrophils was identified. Most hnRNPs were associated with worse survival of adrenocortical carcinoma (ACC), liver hepatocellular carcinoma (LIHC) and lung adenocarcinoma (LUAD), whereas hnRNPs predicted better prognosis in kidney renal clear cell carcinoma (KIRC) and thymoma (THYM). The prognosis analysis of KIRC suggested that hnRNPs gene cluster was significantly associated with overall survival (HR = 0.5, 95% CI = 0.35-0.73, P = 0.003). These findings provide novel evidence for further investigation of hnRNPs in the development and therapy of cancer in the future.

A novel splicing mutation in F8 causes various aberrant transcripts in a hemophilia A patient and identifies a new transcript from healthy individuals

: Hemophilia A is an X-linked hemorrhagic disorder caused by deficiency or dysfunction of the coagulation factor VIII (FVIII), and a great variety of mutations in the factor VIII gene (F8) are identified. We aimed to identify the genetic defects of the F8 gene in a Chinese patient with moderate hemophilia A. We have identified a novel intronic variant in the hemophilia A patient by DNA sequence analysis, cDNA sequencing, and TA clone sequencing. An intronic variant, c.5816-1G>A, was identified and the cDNA sequencing confirmed the pathogenicity of the transition. TA clone sequencing showed that the splicing mutation produced two aberrant premRNA skipping exons (18 and exon 18 + 19, respectively). These aberrant mRNA forms maintain the reading frame and are predicted to code for deleted FVIII isoforms and the shorter abnormal transcript accounted for one-eighth of the total mRNA. There was a new unreported transcript with E22 spliced out in healthy individuals and our patient, whose specific functions need to be determined in further studies. Our study widens the mutation spectrum of the F8 gene. In addition, the study findings could provide the opportunity to reveal alternative splicing patterns.

The Crosstalk Between Epigenetic Mechanisms and Alternative RNA Processing Regulation

As a co-transcriptional process, RNA processing, including alternative splicing and alternative polyadenylation, is crucial for the generation of multiple mRNA isoforms. RNA processing mechanisms are widespread across all higher eukaryotes and play critical roles in cell differentiation, organ development and disease response. Recently, significant progresses have been made in understanding the mechanism of RNA processing. RNA processing is regulated by trans-acting factors such as splicing factors, RNA-binding proteins and cis-sequences in pre-mRNA, and increasing evidence suggests that epigenetic mechanisms, which are important for the dynamic regulation and state of specific chromatic regions, are also involved in co-transcriptional RNA processing. In contrast, recent studies also suggest that alternative RNA processing also has a feedback regulation on epigenetic mechanisms. In this review, we discuss recent studies and summarize the current knowledge on the epigenetic regulation of alternative RNA processing. In addition, a feedback regulation of RNA processing on epigenetic regulators is also discussed.

Aberrant spliceosome expression and altered alternative splicing events correlate with maturation deficiency in human oocytes

Different strategies of ovarian stimulation are widely used in IVF to retrieve mature metaphase II (MII) oocytes for fertilization. On average, approximately 70% of recovered oocytes are mature, while personalized administration of hCG and/or GnRH agonist trigger and in vitro maturation (IVM) management can further improve the maturation rate. However, even under such conditions, a complete absence of oocyte maturation is still observed sporadically. The probable causes for such maturation-deficient (MD) oocytes - which arrest abnormally at metaphase I (MI) stage - are still under investigation. In the present study, using single-cell transcriptomic RNA sequencing (RNA-seq) and differential expression analysis, we showed that gene expression profiles were aberrant, and alternative splicing (AS) patterns were changed in MD oocytes when compared with normally mature (MN) oocytes. Gene ontology (GO) enrichment demonstrated that the differently expressed genes (DEGs) were mostly correlated with pre-mRNA splicing, RNA transportation, RNA processing, and mRNA regulation. Subsequently, analysis of AS events revealed that genes with altered AS patterns were primarily associated with metabolism and cell cycle. With these findings, we have demonstrated aberrant gene expression in complete maturation-deficient oocytes, and we propose that alterations in post-transcriptional regulation constitute a potential underlying mechanism governing oocyte maturation.

Splicing Anomalies in Myeloproliferative Neoplasms: Paving the Way for New Therapeutic Venues

Since the discovery of spliceosome mutations in myeloid malignancies, abnormal pre-mRNA splicing, which has been well studied in various cancers, has attracted novel interest in hematology. However, despite the common occurrence of spliceosome mutations in myelo-proliferative neoplasms (MPN), not much is known regarding the characterization and mechanisms of splicing anomalies in MPN. In this article, we review the current scientific literature regarding "splicing and myeloproliferative neoplasms". We first analyse the clinical series reporting spliceosome mutations in MPN and their clinical correlates. We then present the current knowledge about molecular mechanisms by which these mutations participate in the pathogenesis of MPN or other myeloid malignancies. Beside spliceosome mutations, splicing anomalies have been described in myeloproliferative neoplasms, as well as in acute myeloid leukemias, a dreadful complication of these chronic diseases. Based on splicing anomalies reported in chronic myelogenous leukemia as well as in acute leukemia, and the mechanisms presiding splicing deregulation, we propose that abnormal splicing plays a major role in the evolution of myeloproliferative neoplasms and may be the target of specific therapeutic strategies.

SURF: integrative analysis of a compendium of RNA-seq and CLIP-seq datasets highlights complex governing of alternative transcriptional regulation by RNA-binding proteins

Advances in high-throughput profiling of RNA-binding proteins (RBPs) have resulted inCLIP-seq datasets coupled with transcriptome profiling by RNA-seq. However, analysis methods that integrate both types of data are lacking. We describe SURF, Statistical Utility for RBP Functions, for integrative analysis of large collections of CLIP-seq and RNA-seq data. We demonstrate SURF's ability to accurately detect differential alternative transcriptional regulation events and associate them to local protein-RNA interactions. We apply SURF to ENCODE RBP compendium and carry out downstream analysis with additional reference datasets. The results of this application are browsable at http://www.statlab.wisc.edu/shiny/surf/.

Heterogeneous nuclear ribonucleoprotein A3 controls mitotic progression of neural progenitors via interaction with cohesin

Cortex development is controlled by temporal patterning of neural progenitor (NP) competence with sequential generation of deep and superficial layer neurons, but underlying mechanisms remain elusive. Here, we report a role for heterogeneous nuclear ribonucleoprotein A3 (HNRNPA3) in regulating the division of early cortical NPs that mainly give rise to deep-layer neurons via direct neurogenesis. HNRNPA3 is expressed at high levels in NPs of mouse and human cortex at early stages, with a unique peri-chromosome pattern. Intriguingly, downregulation of HNRNPA3 caused chromosome disarrangement, which hindered normal separation of chromosomes during NP division, leading to mitotic delay. Furthermore, HNRNPA3 is associated with the cohesin-core subunit SMC1A and controls its association with chromosomes, implicating a mechanism for the role of HNRNPA3 in regulating chromosome segregation in dividing NPs. Hnrnpa3-deficient mice exhibited reduced cortical thickness, especially of deep layers. Moreover, downregulation of HNRNPA3 in cultured human cerebral organoids led to marked reduction in NPs and deep-layer neurons. Thus, this study has identified a crucial role for HNRNPA3 in NP division and highlighted the relationship between mitosis progression and early neurogenesis.

A novel splice variant of Gαq-coupled Bombyx CAPA-PVK receptor 1 functions as a specific Gαi/o-linked receptor for CAPA-PK

Alternative splicing enables G protein-coupled receptor (GPCR) genes to greatly increase the number of structurally and functionally distinct receptor isoforms. However, the functional role and relevance of the individual GPCR splice variants in regulating physiological processes are still to be assessed. A naturally occurring alternative splice variant of Bombyx CAPA-PVK receptor, BomCAPA-PVK-R1-Δ341, has been shown to act as a dominant-negative protein to regulate cell surface expression and function of the canonical CAPA-PVK receptor. Herein, using functional assays, we identify the splice variant Δ341 as a specific receptor for neuropeptide CAPA-PK, and upon activation, Δ341 signals to ERK1/2 pathway. Further characterization demonstrates that Δ341 couples to Gαi/o, distinct from the Gαq-coupled canonical CAPA-PVK receptor, triggering ERK1/2 phosphorylation through Gβγ-PI3K-PKCζ signaling cascade. Moreover, our ELISA data show that the ligand-dependent internalization of the splice variant Δ341 is significantly impaired due to lack of GRKs-mediated phosphorylation sites. Our findings highlight the potential of this knowledge for molecular, pharmacological and physiological studies on GPCR splice variants in the future.

Significance of LncRNA CASC8 genetic polymorphisms on the tuberculosis susceptibility in Chinese population

Background

Tuberculosis remains an important disease threatening the security of public health, and no effective targets have been found for the immunological diagnosis or therapy of tuberculosis. The aim of this study was to explore the associations between lncRNA CASC8 genetic polymorphism and tuberculosis risk.

Method

A total of 900 tuberculosis patients and 1534 healthy individuals in the Western Chinese Han population were recruited for our study. Candidate SNPs of CASC8 were initially filtered by importing the 1000 genomes database into Haploview, and subsequently genotyped using modified multiplex ligation detection reactions.

Results

The lncRNA CASC8 genetic variant rs7836840 was associated with an increased tuberculosis risk with a P‐value of .034, but .134 after Bonferroni correction. Using subtype analysis, the C allele in rs7836840 showed a significant association with tuberculosis susceptibility (OR = 1.196, 95% CI = 1.05‐1.362, P = .02739 after Bonferroni correction). Patients carrying genotype AG and GG of rs7825118 and rs9297758 exhibited lower Hb concentrations (P = .006) and neutrophil counts (P = .015), respectively, while genotype AG and AA in rs6981424 demonstrated higher levels of ALT (P = .005) and AST (P = .033) in a dominant model, which were consistent with a tendency toward increased TB risk.

Conclusions

This study was the first to explore the association between lncRNA CASC8 polymorphisms and TB infection risk and clinical manifestations. Our results provide evidence that CASC8 may act as a biomarker for the progression of clinical tuberculosis.

Poison-Exon Inclusion in DHX9 Reduces Its Expression and Sensitizes Ewing Sarcoma Cells to Chemotherapeutic Treatment

Alternative splicing is a combinatorial mechanism by which exons are joined to produce multiple mRNA variants, thus expanding the coding potential and plasticity of eukaryotic genomes. Defects in alternative splicing regulation are associated with several human diseases, including cancer. Ewing sarcoma is an aggressive tumor of bone and soft tissue, mainly affecting adolescents and young adults. DHX9 is a key player in Ewing sarcoma malignancy, and its expression correlates with worse prognosis in patients. In this study, by screening a library of siRNAs, we have identified splicing factors that regulate the alternative inclusion of a poison exon in DHX9 mRNA, leading to its downregulation. In particular, we found that hnRNPM and SRSF3 bind in vivo to this poison exon and suppress its inclusion. Notably, DHX9 expression correlates with that of SRSF3 and hnRNPM in Ewing sarcoma patients. Furthermore, downregulation of SRSF3 or hnRNPM inhibited DHX9 expression and Ewing sarcoma cell proliferation, while sensitizing cells to chemotherapeutic treatment. Hence, our study suggests that inhibition of hnRNPM and SRSF3 expression or activity could be exploited as a therapeutic tool to enhance the efficacy of chemotherapy in Ewing sarcoma.

An Alu Element Insertion in Intron 1 Results in Aberrant Alternative Splicing of APOBEC3G Pre-mRNA in Northern Pig-Tailed Macaques (Macaca leonina) That May Reduce APOBEC3G-Mediated Hypermutation Pressure on HIV-1

APOBEC3 family members, particularly APOBEC3F and APOBEC3G, inhibit the replication and spread of various retroviruses by inducing hypermutation in newly synthesized viral DNA. Viral hypermutation by APOBEC3 is associated with viral evolution, viral transmission, and disease progression. In recent years, increasing attention has been paid to targeting APOBEC3G for AIDS therapy. Thus, a controllable model system using species such as macaques, which provide a relatively ideal in vivo system, is needed for the study of APOBEC3-related issues. To appropriately utilize this animal model for biomedical research, important differences between human and macaque APOBEC3s must be considered. In this study, we found that the ratio of APOBEC3G-mediated/APOBEC3-mediated HIV-1 hypermutation footprints was much lower in peripheral blood mononuclear cells (PBMCs) from northern pig-tailed macaques than in PBMCs from humans. Next, we identified a novel and conserved APOBEC3G pre-mRNA alternative splicing pattern in macaques, which differed from that in humans and resulted from an Alu element insertion into macaque APOBEC3G gene intron 1. This alternative splicing pattern generating an aberrant APOBEC3G mRNA isoform may significantly dilute full-length APOBEC3G and reduce APOBEC3G-mediated hypermutation pressure on HIV-1 in northern pig-tailed macaques, which was supported by the elimination of other possibilities accounting for this hypermutation difference between the two hosts.IMPORTANCE APOBEC3 family members, particularly APOBEC3F and APOBEC3G, are important cellular antiviral factors. Recently, more attention has been paid to targeting APOBEC3G for AIDS therapy. To appropriately utilize macaque animal models for the study of APOBEC3-related issues, it is important that the differences between human and macaque APOBEC3s are clarified. In this study, we identified a novel and conserved APOBEC3G pre-mRNA alternative splicing pattern in macaques, which differed from that in humans and which may reduce the APOBEC3G-mediated hypermutation pressure on HIV-1 in northern pig-tailed macaques (NPMs). Our work provides important information for the proper application of macaque animal models for APOBEC3-related issues in AIDS research and a better understanding of the biological functions of APOBEC3 proteins.

Splicing Busts a Move: Isoform Switching Regulates Migration

Cell migration is essential for normal development, neural patterning, pathogen eradication, and cancer metastasis. Pre-mRNA processing events such as alternative splicing and alternative polyadenylation result in greater transcript and protein diversity as well as function and activity. A critical role for alternative pre-mRNA processing in cell migration has emerged in axon outgrowth during neuronal development, immune cell migration, and cancer metastasis. These findings suggest that migratory signals result in expression changes of post-translational modifications of splicing or polyadenylation factors, leading to splicing events that generate promigratory isoforms. We summarize this recent progress and suggest emerging technologies that may facilitate a deeper understanding of the role of alternative splicing and polyadenylation in cell migration.

Systematic profiling of alternative splicing signature reveals prognostic predictor for cervical cancer

Aim

Cervical cancer is a common malignant carcinoma of the gynecological tract with high morbidity and mortality. Therefore, it is crucial to elucidate the pathogenesis, prevention, diagnosis and prognosis of cervical cancer by searching for the involved key genes.

Method

In this study, the alternative splicing (AS) events of 253 patients with cervical cancer were analyzed, and 41,766 AS events were detected in 9961 genes. Univariate analysis was performed to screen prognostic AS events. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was used to identify the pathways in which these AS events were involved.

Results

We found that exon skip (ES) is the main AS event in patients with cervical cancer. There was pronounced consistency between the genes involved in overall survival and those involved in recurrence. At the same time, we found that a gene may exhibit several different types of AS events, and these different AS events may be related to prognosis. Four characteristic genes, HSPA14, SDHAF2, CAMKK2 and TM9SF1, that can be used as prognostic markers for cervical cancer were selected. Conclusion: The importance of AS events in the development of cervical cancer and prediction of prognosis was revealed by a large amount of data at the whole genome level, which may provide a potential target for cervical cancer treatment. We also provide a new method for exploring the pathogenesis of cervical cancer to determine clinical treatment and prognosis more accurately.

Thermosensitive alternative splicing senses and mediates temperature adaptation in Drosophila

Circadian rhythms are generated by the cyclic transcription, translation, and degradation of clock gene products, including timeless (tim), but how the circadian clock senses and adapts to temperature changes is not completely understood. Here, we show that temperature dramatically changes the splicing pattern of tim in Drosophila. We found that at 18°C, TIM levels are low because of the induction of two cold-specific isoforms: tim-cold and tim-short and cold. At 29°C, another isoform, tim-medium, is upregulated. Isoform switching regulates the levels and activity of TIM as each isoform has a specific function. We found that tim-short and cold encodes a protein that rescues the behavioral defects of tim⁰¹ mutants, and that flies in which tim-short and cold is abrogated have abnormal locomotor activity. In addition, miRNA-mediated control limits the expression of some of these isoforms. Finally, data that we obtained using minigenes suggest that tim alternative splicing might act as a thermometer for the circadian clock.

Polymorphisms within the Boule Gene Detected by Tetra-Primer Amplification Refractory Mutation System PCR (T-ARMS-PCR) are Significantly Associated with Goat Litter Size

As a gene contributing to spermatogenesis, the Boule gene (also called Boll), whose mutations result in azoospermia and sterility of flies and mice, was conserved in reductional maturation divisions. However, in goats, the polymorphisms of Boule, especially with regard to their fundamental roles in female reproduction traits, are still unknown. Therefore, the aims of this study were to detect a potential mutation (rs661484476: g.7254T>C) located in intron 2 of the Boule gene by tetra-primer amplification refractory mutation system PCR (T-ARMS-PCR) and to explore its potential association with the litter size of Shaanbei White-Cashmere goats (SBWGs). In this study, g.7254T>C was firstly detected. The TT genotype was the dominant genotype in the single-lamb group, and T was also the dominant allele in all tested groups. Additionally, the detected locus displayed moderate polymorphism with polymorphism information content (PIC) values among all studied goats ranging from 0.303 to 0.344. Notably, according to the χ² test, the distribution differences for the genotypic frequencies between the single- and multi-lamb groups was significant (p = 0.014). Furthermore, the polymorphisms of the goat Boule gene were significantly associated with the goat litter size in SBWGs (p < 0.05), which indicated that g.7254T>C could be a potential marker in the marker-assisted selection process for potential litter size in goats. These results also indicated that the Boule gene might exercise important functions in female goat reproduction, which provided new insight for female goat breeding and could accelerate the process of goat breeding.

Expression of ZNF695 Transcript Variants in Childhood B-Cell Acute Lymphoblastic Leukemia

B-cell acute lymphoblastic leukemia is the most commonly diagnosed childhood malignancy worldwide; more than 50% of these cases are diagnosed in Mexico. Although the five-year survival rate is >80%, 30% of patients experience relapse with poor prognosis. Cancer-associated gene expression profiles have been identified in several malignancies, and some transcripts have been used to predict disease prognosis. The human transcriptome is incompletely elucidated; moreover, more than 80% of transcripts can be processed via alternative splicing (AS), which increases transcript and protein diversity. The human transcriptome is divided; coding RNA accounts for 2%, and the remaining 98% is noncoding RNA. Noncoding RNA can undergo AS, promoting the diversity of noncoding transcripts. We designed specific primers to amplify previously reported alternative transcript variants of ZNF695 and showed that six ZNF695 transcript variants are co-expressed in cancer cell lines. The amplicons were sequenced and identified. Additionally, we analyzed the expression of these six transcript variants in bone marrow from B-cell acute lymphoblastic leukemia patients and observed that ZNF695 transcript variants one and three were the predominant variants expressed in leukemia. Moreover, our results showed the co-expression of coding and long noncoding RNA. Finally, we observed that long noncoding RNA ZNF695 expression predicted survival rates.

Bridging the gap between reference and real transcriptomes

Genetic, transcriptional, and post-transcriptional variations shape the transcriptome of individual cells, rendering establishing an exhaustive set of reference RNAs a complicated matter. Current reference transcriptomes, which are based on carefully curated transcripts, are lagging behind the extensive RNA variation revealed by massively parallel sequencing. Much may be missed by ignoring this unreferenced RNA diversity. There is plentiful evidence for non-reference transcripts with important phenotypic effects. Although reference transcriptomes are inestimable for gene expression analysis, they may turn limiting in important medical applications. We discuss computational strategies for retrieving hidden transcript diversity.

Identification, evolution and alternative splicing profile analysis of the splicing factor 30 (SPF30) in plant species

The work offers a comprehensive evaluation on the phylogenetics and conservation of splicing patterns of the plant SPF30 splicing factor gene family. In eukaryotes, one pre-mRNA can generate multiple mRNA transcripts by alternative splicing (AS), which expands transcriptome and proteome diversity. Splicing factor 30 (SPF30), also known as survival motor neuron domain containing protein 1 (SMNDC1), is a spliceosomal protein that plays an essential role in spliceosomal assembly. Although SPF30 genes have been well characterised in human and yeast, little is known about their homologues in plants. Here, we report the genome-wide identification and phylogenetic analysis of SPF30 genes in the plant kingdom. In total, 82 SPF30 genes were found in 64 plant species from algae to land plants. Alternative transcripts were found in many SPF30 genes and splicing profile analysis revealed that the second intron in SPF30 genome is frequently associated with AS events and contributed to the birth of novel exons in a few SPF30 members. In addition, different conserved sequences were observed at these putative splice sites among moss, monocots and dicots, respectively. Our findings will facilitate further functional characterization of plant SPF30 genes as putative splicing factors.

Potential epigenomic co-management in rare diseases and epigenetic therapy

The purpose of this review is to highlight the impact of the alternative splicing process on human disease. Epigenetic regulation determines not only what parts of the genome are expressed but also how they are spliced. The recent progress in the field of epigenetics has important implications for the study of rare diseases. The role of epigenetics in rare diseases is a key issue in molecular physiology and medicine because not only rare diseases can benefit from epigenetic research, but can also provide useful principles for other common and complex disorders such as cancer, cardiovascular, type 2 diabetes, obesity, and neurological diseases. Predominantly, epigenetic modifications include DNA methylation, histone modification, and RNA-associated silencing. These modifications in the genome regulate numerous cellular activities. Disruption of epigenetic regulation process can contribute to the etiology of numerous diseases during both prenatal and postnatal life. Here, I discuss current knowledge about this matter including some current epigenetic therapies and future directions in the field by emphasizing on the RNA-based therapy via antisense oligonucleotides to correct splicing defects.

Novel alternative splice variants of the human protein arginine methyltransferase 1 (PRMT1) gene, discovered using next-generation sequencing

Next-generation sequencing (NGS) technology is highly expected to help researchers disclose the complexity of alternative splicing and understand its association with carcinogenesis. Alternative splicing alterations are firmly associated with multiple malignancies, in terms of functional roles in malignant transformation, motility, and/or metastasis of cancer cells. One perfect example illustrating the connection between alternative splicing and cancer is the human protein arginine methyltransferase 1 (PRMT1) gene, previously cloned from members of our research group and involved in a variety of processes including transcription, DNA repair, and signal transduction. Two splice variants of PRMT1 (variants v.1 and v.2) are downregulated in breast cancer. In addition, PRMT1 v.2 promotes the survival and invasiveness of breast cancer cells, while it could serve as a biomarker of unfavorable prognosis in colon cancer patients. The aim of this study was the molecular cloning of novel alternative splice variants of PRMT1 with the use of 3' RACE coupled with NGS technology. Extensive bioinformatics and computational analysis revealed a significant number of 19 novel alternative splicing events between annotated exons of PRMT1 as well as one novel exon, resulting in the discovery of multiple PRMT1 transcripts. In order to validate the full sequence of the novel transcripts, RT-PCR was carried out with the use of variant-specific primers. As a result, 58 novel PRMT1 transcripts were identified, 34 of which are mRNAs encoding new protein isoforms, whereas the rest 24 transcripts are candidates for nonsense-mediated mRNA decay (NMD).

Functional Study of Carboxylesterase 1 Protein Isoforms

Carboxylesterase 1 (CES1) is a primary human hepatic hydrolase involved in hydrolytic biotransformation of numerous medications. Considerable interindividual variability in CES1 expression and activity has been consistently reported. Four isoforms of the CES1 protein are produced by alternative splicing (AS). In the current study, the activity and expression of each CES1 isoform are examined using transfected cell lines, and CES1 isoform composition and its impact on CES1 activity in human livers are determined. In transfected cells, isoforms 3 and 4 show mRNA and protein expressions comparable to isoforms 1 and 2, but have significantly impaired activity when hydrolyzing enalapril and clopidogrel. In individual human liver samples, isoforms 1 and 2 are the major forms, contributing 73-90% of the total CES1 protein expression. In addition, the protein expression ratios of isoforms 1 and 2 to isoforms 3 and 4 are positively associated with CES1 activity in the liver, suggesting that CES1 isoform composition is a factor contributing to the variability in hepatic CES1 function. Further investigations of the regulation of CES1 AS would improve the understanding of CES1 variability and help develop a strategy to optimize the pharmacotherapy of many CES1 substrate medications.

An exonic splicing enhancer mutation in DUOX2 causes aberrant alternative splicing and severe congenital hypothyroidism in Bama pigs

Pigs share many similarities with humans in terms of anatomy, physiology and genetics, and have long been recognized as important experimental animals in biomedical research. Using an N-ethyl-N-nitrosourea (ENU) mutagenesis screen, we previously identified a large number of pig mutants, which could be further established as human disease models. However, the identification of causative mutations in large animals with great heterogeneity remains a challenging endeavor. Here, we select one pig mutant, showing congenital nude skin and thyroid deficiency in a recessive inheritance pattern. We were able to efficiently map the causative mutation using family-based genome-wide association studies combined with whole-exome sequencing and a small sample size. A loss-of-function variant (c.1226 A>G) that resulted in a highly conserved amino acid substitution (D409G) was identified in the DUOX2 gene. This mutation, located within an exonic splicing enhancer motif, caused aberrant splicing of DUOX2 transcripts and resulted in lower H₂O₂ production, which might cause a severe defect in thyroid hormone production. Our findings suggest that exome sequencing is an efficient way to map causative mutations and that DUOX2^D409G/D409G mutant pigs could be a potential large animal model for human congenital hypothyroidism.

Summary: Here, we show that an exonic splicing enhancer variant in DUOX2 (c.1226 A>G) causes aberrant splicing of DUOX2 transcripts, resulting in lower H₂O₂ production, to cause severe congenital hypothyroidism in Bama pigs.

View from an mRNP: The Roles of SR Proteins in Assembly, Maturation and Turnover

Serine- and arginine-rich proteins (SR proteins) are a family of multitasking RNA-binding proteins (RBPs) that are key determinants of messenger ribonucleoprotein (mRNP) formation, identity and fate. Apart from their essential functions in pre-mRNA splicing, SR proteins display additional pre- and post-splicing activities and connect nuclear and cytoplasmic gene expression machineries. Through changes in their post-translational modifications (PTMs) and their subcellular localization, they provide functional specificity and adjustability to mRNPs. Transcriptome-wide UV crosslinking and immunoprecipitation (CLIP-Seq) studies revealed that individual SR proteins are present in distinct mRNPs and act in specific pairs to regulate different gene expression programmes. Adopting an mRNP-centric viewpoint, we discuss the roles of SR proteins in the assembly, maturation, quality control and turnover of mRNPs and describe the mechanisms by which they integrate external signals, coordinate their multiple tasks and couple subsequent mRNA processing steps.

Targeting VEGF pathway to normalize the vasculature: an emerging insight in cancer therapy

Vascular normalization is a new concept of targeting angiogenesis to restore vessel structure and function and to increase blood perfusion and delivery of drugs. It has been confirmed that vascular normalization can decrease relapse and benefit other cancer therapy, including chemotherapy, radiotherapy, and immune cell therapy. The key point of this therapy is to inhibit pro-angiogenic factors and make it be balanced with anti-angiogenic factors, resulting in a mature and normal vessel characteristic. Vascular endothelial growth factor (VEGF) is a key player in the process of tumor angiogenesis, and inhibiting VEGF is a primary approach to tumor vessel normalization. Herein, we review newly uncovered mechanisms governing angiogenesis and vascular normalization of cancer and place emphasis on targeting VEGF pathway to normalize the vasculature. Also, important methods to depress VEGF pathway and make tumor vascular are discussed.

Circular RNAs as Therapeutic Agents and Targets

It has recently been reported that thousands of covalently linked circular RNAs (circRNAs) are expressed from human genomes. circRNAs emerge during RNA splicing. circRNAs are circularized in a reaction termed "backsplicing," whereby the spliceosome fuses a splice donor site in a downstream exon to a splice acceptor site in an upstream exon. Although a young field of research, first studies indicate that backsplicing is not an erroneous reaction of the spliceosome. Instead, circRNAs are produced in cells with high cell-type specificity and can exert biologically meaningful and specific functions. These observations and the finding that circRNAs are stable against exonucleolytic decay are raising the question whether circRNAs may be relevant as therapeutic agents and targets. In this review, we start out with a short introduction into classification, biogenesis and general molecular mechanisms of circRNAs. We then describe reports, where manipulating circRNA abundance has been shown to have therapeutic value in animal disease models in vivo, with a focus on cardiovascular disease (CVD). Starting from existing approaches, we outline particular challenges and opportunities for future circRNA-based therapeutic approaches that exploit stability and molecular effector functions of native circRNAs. We end with considerations which designer functions could be engineered into artificial therapeutic circular RNAs.