Genome-wide survey of human alternative pre-mRNA splicing with exon junction microarrays

Interdependent regulation of alternative splicing by SR and hnRNP proteins

Alternative pre-mRNA splicing is a combinatorial process involving SR and hnRNP splicing factors. These proteins can silence or enhance splicing based on their expression levels and binding positions. To better understand their combinatorial and interdependent regulation, computational analyses were performed using HepG2 and K562 cell knockdown and binding datasets from the ENCODE Project. Analyses of diMerential splicing for 6 SR proteins and 13 hnRNP knockdowns revealed statistically significant exon overlap among most RBP combinations, albeit at diMerent levels. Neither SR proteins nor hnRNPs showed strong preferences for collaborating with specific RBP classes in mediating exon inclusion. While SRSF1, hnRNPK, and hnRNPC stand out as major influencers of alternative splicing, they do so predominantly independent of other RBPs. Meanwhile, minor influencers of alternative splicing such as hnRNPAB and hnRNPA0 predominantly regulate exon inclusion in concert with other RBPs, indicating that inclusion can be mediated by both single and multiple RBPs. Interestingly, the higher the number of RBPs that regulate the inclusion of an exon, the more variable exon inclusion preferences become. Interdependently regulated exons are more modular and have diMerent physical characteristics such as reduced exon length compared to their independent counterparts. A comparison of RBP interdependence between HepG2 and K562 cells provides the framework that explains cell-type-specific alternative splicing. Our study highlights the importance of the interdependent regulation of alternative exons and identifies characteristics of interdependently regulated exons that diMer from independently regulated exons.

Assessing the causal relationship between plasma proteins and osteoporosis: novel insights into pathological mechanisms and therapeutic implications

Identifying dysregulated plasma proteins in osteoporosis (OP) progression offers insights into prevention and treatment. This study found 8 such proteins associated with OP, suggesting them as therapy targets. This discovery may cut drug development costs and improve personalized treatments.

PURPOSE

This study aims to identify potential therapeutic targets for OP using summary data-based Mendelian randomization (SMR) and colocalization analysis methods. Furthermore, we seek to explore the biological significance and pharmacological value of these drug targets.

METHODS

To identify potential therapeutic targets for OP, we conducted SMR and colocalization analysis. Plasma protein (pQTL, exposure) data were sourced from the study by Ferkingstad et al. (n = 35,559). Summary statistics for bone mineral density (BMD, outcome) were obtained from the GWAS Catalog (n = 56,284). Additionally, we utilized enrichment analysis, protein-protein interaction (PPI) network analysis, drug prediction, and molecular docking to further analyze the biological significance and pharmacological value of these drug targets.

RESULTS

In the SMR analysis, while 20 proteins showed significance, only 8 potential drug targets (GCKR, ERBB3, CFHR1, GPN1, SDF2, VTN, BET1L, and SERPING1) received support from colocalization (PP.H4 > 0.8). These proteins are closely associated with immune function in terms of biological significance. Molecular docking also demonstrated favorable binding of drugs to proteins, consistent with existing structural data, further substantiating the pharmacological value of these targets.

CONCLUSIONS

The study identified 8 potential drug targets for OP. These prospective targets are believed to have a higher chance of success in clinical trials, thus aiding in prioritizing OP drug development and reducing development costs.

Insulin receptor isoform B is required for efficient proinsulin processing in pancreatic β cells

The insulin receptor (INSR, IR) has two isoforms, IRA and IRB, through alternative splicing. However, their distinct functions in vivo remain unclear. Here we generated β cell-specific IRB knockout (KO) mice (βIRBKO). The KO mice displayed worsened hyperinsulinemia and hyperproinsulinemia in diet-induced obesity due to impaired proinsulin processing in β cells. Mechanistically, loss of IRB suppresses eukaryotic translation initiation factor 4G1 (eIF4G1) by stabilizing the transcriptional receptor sterol-regulatory element binding protein 1 (SREBP1). Moreover, excessive autocrine proinsulin in βIRBKO mice enhances the activity of extracellular signal-regulated kinase (ERK) through the remaining IRA to further stabilize nuclear SREBP1, forming a feedback loop. Collectively, our study paves the way to dissecting the isoform-specific function of IR in vivo and highlights the important roles of IRB in insulin processing and protecting β cells from lipotoxicity in obesity.

A reliable and quick method for screening alternative splicing variants for low-abundance genes

Alternative splicing (AS) is a universal phenomenon in eukaryotes, and it is still challenging to identify AS events. Several methods have been developed to identify AS events, such as expressed sequence tags (EST), microarrays and RNA-seq. However, EST has limitations in identifying low-abundance genes, while microarray and RNA-seq are high-throughput technologies, and PCR-based technology is needed for validation. To overcome the limitations of EST and shortcomings of high-throughput technologies, we established a method to identify AS events, especially for low-abundance genes, by reverse transcription (RT) PCR with gene-specific primers (GSPs) followed by nested PCR. This process includes two major steps: 1) the use of GSPs to amplify as long as the specific gene segment and 2) multiple rounds of nested PCR to screen the AS and confirm the unknown splicing variants. With this method, we successfully identified three new splicing variants, namely, GenBank Accession No. HM623886 for the bdnf gene (GenBank GeneID: 12064), GenBank Accession No. JF417977 for the trkc gene (GenBank GeneID: 18213) and GenBank Accession No. HM623888 for the glb-18 gene (GenBank GeneID: 172485). In addition to its reliability and simplicity, the method is also cost-effective and labor-intensive. In conclusion, we developed an RT-nested PCR method using gene-specific primers to efficiently identify known and novel AS variants. This approach overcomes the limitations of existing methods for detecting rare transcripts. By enabling the discovery of new isoforms, especially for low-abundance genes, this technique can aid research into aberrant splicing in disease. Future studies can apply this method to uncover AS variants involved in cancer, neurodegeneration, and other splicing-related disorders.

Landscape analysis of alternative splicing in kidney renal clear cell carcinoma and their clinical significance

A growing number of studies reveal that alternative splicing (AS) is associated with tumorigenesis, progression, and metastasis. Systematic analysis of alternative splicing signatures in renal cancer is lacking. In our study, we investigated the AS landscape of kidney renal clear cell carcinoma (KIRC) and identified AS predictive model to improve the prognostic prediction of KIRC. We obtained clinical data and gene expression profiles of KIRC patients from the TCGA database to evaluate AS events. The calculation results for seven types of AS events indicated that 46276 AS events from 10577 genes were identified. Next, we applied Cox regression analysis to identify 5864 prognostic-associated AS events. We used the Metascape database to verify the potential pathways of prognostic-associated AS. Moreover, we constructed KIRC prediction systems with prognostic-associated AS events by the LASSO Cox regression model. AUCs demonstrated that these prediction systems had excellent prognostic accuracy simultaneously. We identified 34 prognostic associated splicing factors (SFs) and constructed homologous regulatory networks. Furthermore, in vitro experiments were performed to validate the favorable effect of SFs FMR1 in KIRC. In conclusion, we overviewed AS events in KIRC and identified AS-based prognostic models to assist the survival prediction of KIRC patients. Our study may provide a novel predictive signature to improve the prognostic prediction of KIRC, which might facilitate KIRC patient counseling and individualized management.

Specific recognition and sensitive quantification of mRNA splice variants via one-pot ligation-dependent loop-mediated isothermal amplification

Specific recognition and sensitive quantification of mRNA alternative splice variants have been a necessity for exploring the regulatory mechanism of RNA splicing and revealing the association between pre-mRNA splicing and transcriptome function, as well as disease diagnosis. However, their wide abundance range and high sequence homology pose enormous challenges for high sensitivity and selectivity quantification of splice variants. Herein, taking advantage of the excellent specificity of ligation and the powerful nucleic acid replication feature of loop-mediated isothermal amplification (LAMP), we developed a one-pot method (termed one-pot ligation-LAMP) for specific recognition and sensitive quantification of mRNA splicing variants based on two splicing junction-specific stem-loop DNA probe ligation and the subsequently initiating LAMP. The one-pot ligation-LAMP can specifically detect as low as 100 aM mRNA splice variants without any nonspecific signals and quantify them with a wide dynamics range spanning at least six orders of magnitude. We have demonstrated that the one-pot ligation-LAMP is a versatile and practical strategy for accurately quantifying different splicing variants in complex biological samples with high sensitivity all in one tube within 90 min, thereby providing an attractive tool for mRNA splice variant-related studies.

Mutations Causing X-Linked Amelogenesis Imperfecta Alter miRNA Formation from Amelogenin Exon4

Amelogenin plays a crucial role in tooth enamel formation, and mutations on X-chromosomal amelogenin cause X-linked amelogenesis imperfecta (AI). Amelogenin pre-messenger RNA (mRNA) is highly alternatively spliced, and during alternative splicing, exon4 is mostly skipped, leading to the formation of a microRNA (miR-exon4) that has been suggested to function in enamel and bone formation. While delivering the functional variation of amelogenin proteins, alternative splicing of exon4 is the decisive first step to producing miR-exon4. However, the factors that regulate the splicing of exon4 are not well understood. This study aimed to investigate the association between known mutations in exon4 and exon5 of X chromosome amelogenin that causes X-linked AI, the splicing of exon4, and miR-exon4 formation. Our results showed mutations in exon4 and exon5 of the amelogenin gene, including c.120T>C, c.152C>T, c.155C>G, and c.155delC, significantly affected the splicing of exon4 and subsequent miR-exon4 production. Using an amelogenin minigene transfected in HEK-293 cells, we observed increased inclusion of exon4 in amelogenin mRNA and reduced miR-exon4 production with these mutations. In silico analysis predicted that Ser/Arg-rich RNA splicing factor (SRSF) 2 and SRSF5 were the regulatory factors for exon4 and exon5 splicing, respectively. Electrophoretic mobility shift assay confirmed that SRSF2 binds to exon4 and SRSF5 binds to exon5, and mutations in each exon can alter SRSF binding. Transfection of the amelogenin minigene to LS8 ameloblastic cells suppressed expression of the known miR-exon4 direct targets, Nfia and Prkch, related to multiple pathways. Given the mutations on the minigene, the expression of Prkch has been significantly upregulated with c.155C>G and c.155delC mutations. Together, we confirmed that exon4 splicing is critical for miR-exon4 production, and mutations causing X-linked AI in exon4 and exon5 significantly affect exon4 splicing and following miR-exon4 production. The change in miR-exon4 would be an additional etiology of enamel defects seen in some X-linked AI.

Splam: a deep-learning-based splice site predictor that improves spliced alignments

The process of splicing messenger RNA to remove introns plays a central role in creating genes and gene variants. Here we describe Splam, a novel method for predicting splice junctions in DNA based on deep residual convolutional neural networks. Unlike some previous models, Splam looks at a relatively limited window of 400 base pairs flanking each splice site, motivated by the observation that the biological process of splicing relies primarily on signals within this window. Additionally, Splam introduces the idea of training the network on donor and acceptor pairs together, based on the principle that the splicing machinery recognizes both ends of each intron at once. We compare Splam's accuracy to recent state-of-the-art splice site prediction methods, particularly SpliceAI, another method that uses deep neural networks. Our results show that Splam is consistently more accurate than SpliceAI, with an overall accuracy of 96% at predicting human splice junctions. Splam generalizes even to non-human species, including distant ones like the flowering plant Arabidopsis thaliana. Finally, we demonstrate the use of Splam on a novel application: processing the spliced alignments of RNA-seq data to identify and eliminate errors. We show that when used in this manner, Splam yields substantial improvements in the accuracy of downstream transcriptome analysis of both poly(A) and ribo-depleted RNA-seq libraries. Overall, Splam offers a faster and more accurate approach to detecting splice junctions, while also providing a reliable and efficient solution for cleaning up erroneous spliced alignments.

Alternative splicing of PSMD13 mediated by genetic variants is significantly associated with endometrial cancer risk

OBJECTIVE

Accumulating evidence has shown that aberrant alternative splicing events are closely associated with the onset and development of cancer. However, whether genetic variants-associated alternative splicing is linked to risk of endometrial cancer remains largely uncertain.

METHODS

We identified single nucleotide polymorphisms (SNPs) locates in the splicing number trait locus (sQTL) of endometrial cancer using the CancerSplicing QTL database. In parallel with bioinformatics analysis, we conducted a case-control study comprising 2,000 cases and 2,013 controls to assess the association between identified SNP which possesses mRNA splicing function and endometrial cancer susceptibility. Furthermore, we used the Kaplan-Meier Plotter, The Human Protein Atlas, SPNR, and Spliceman2 databases for sQTL and differential gene expression analyses to identify the genetic variant which most potentially influence the risk of endometrial cancer through alternative splicing to reveal the potential mechanism by which candidate SNPs regulate the risk of endometrial cancer.

RESULTS

The results indicated that SNP rs7128029 A

CONCLUSION

These findings suggest that SNP rs7128029-mediated alternative splicing events in PSMD13 are associated with endometrial cancer risk and may be a potential early screening biomarker for endometrial cancer-susceptible populations.

Collapse

Key Words

• 26S Proteasome Non-ATPase Regulatory Subunit 13
• Alternative Splicing
• Endometrial Cancer
• Single Nucleotide Polymorphism

Collapse

MESH Headings

• Female
• Humans
• Alternative Splicing/genetics
• Case-Control Studies
• Databases, Factual
• Endometrial Neoplasms/genetics
• Odds Ratio
• Risk
• Polymorphism, Single Nucleotide

Collapse

Grants

• (2022)626 Natural Science Foundation of Guangzhou City
• (2022)402 Zunyi City Science and Technology Plan Project
• 2018 F-913 Zunyi Medical University

Collapse

Affiliation(s)

Sisi He Department of Oncology, The Second Affiliated Hospital of Zunyi Medical University, Zunyi, China
Rong Cao The Second Clinical College, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Yan Mao Department of Oncology, The Second Affiliated Hospital of Zunyi Medical University, Zunyi, China
Na Li Department of Gynecology and Obstetrics, The Second Affiliated Hospital of Zunyi Medical University, Zunyi, China
Yanzhe Wang Department of Oncology, The Second Affiliated Hospital of Zunyi Medical University, Zunyi, China
Hu Ma Department of Oncology, The Second Affiliated Hospital of Zunyi Medical University, Zunyi, China
Kunming Tian Department of Gynecology and Obstetrics, The Second Affiliated Hospital of Zunyi Medical University, Zunyi, China Department of Preventive Medicine, School of Public Health, Zunyi Medical University, Zunyi, China.

Collapse

Spatiotemporal Monitoring of Subcellular mRNAs In Situ via Polyadenine-Mediated Dual-Color Sticky Flares

Spatiotemporal monitoring of multiple low-abundance messenger RNAs (mRNAs) is vitally important for the diagnosis and pathologic analysis of cancer. However, it remains a clinical challenge to monitor and track multiple mRNAs location simultaneously in situ at subcellular level with high efficiency. Herein, we proposed polyA-mediated dual-color sticky flares for simultaneous imaging of two kinds of intracellular mRNA biomarkers. Two kinds of fluorescent DNA specific for GalNac-T mRNA and c-Myc mRNA were functionalized onto gold nanoparticles (AuNPs) through efficient polyadenine (polyA) attachment. By tuning polyA length, the lateral spacing and densities of DNA on AuNPs could be precisely engineered. Compared to the traditional thio-DNA-modified nanoprobes, the uniformity, detection sensitivity, and response kinetics of sticky flares were greatly improved, which enables live-cell imaging of mRNAs with enhanced efficiency. With a sticky-end design, the fluorescent DNA could dynamically trace mRNAs after binding with target mRNAs, which realized spatiotemporal monitoring of subcellular mRNAs in situ. Compared to one target mRNA imaging mode, the multiple target imaging mode allows more accurate diagnosis of cancer. Furthermore, the proposed polyA-mediated dual-color sticky flares exhibit excellent cell entry efficiency and low cytotoxicity with a low-cost and simple assembling process, which provide a pivotal tool for multiple targets imaging in living cells.

Minigene splicing assays reveal new insights into exonic variants of the SLC12A3 gene in Gitelman syndrome

BACKGROUND

Gitelman syndrome (GS) is a type of salt-losing tubular disease, most of which is caused by SLC12A3 gene variants, and missense variants account for the majority. Recently, the phenomenon of exon skipping, in which variants disrupt normal pre-mRNA splicing, has been related to a variety of diseases. Therefore, we hypothesize that a certain proportion of SLC12A3 variants can result in disease via interfering with the normal splicing process.

METHODS

We analyzed 342 previously presumed SLC12A3 missense variants using bioinformatics programs and identified candidate variants that may alter the splicing of pre-mRNA through minigene assays.

RESULTS

Our study revealed that, among ten candidate variants, six variants (c.602G>A, c.602G>T, c.1667C>T, c.1925G>A, c.2548G>C, and c.2549G>C) led to complete or incomplete exon skipping by affecting exonic splicing regulatory elements and/or disturbing canonical splice sites.

CONCLUSION

It is worth mentioning that this is the largest study on pre-mRNA splicing of SLC12A3 exonic variants. In addition, our study emphasizes the importance of detecting splicing function at the mRNA level in GS and indicates that minigene analysis is a valuable tool for splicing functional assays of variants in vitro.

Cooperation and competition by RNA-binding proteins in cancer

Post-transcriptional regulation of gene expression plays a major role in determining the cellular proteome in health and disease. Post-transcriptional control mechanisms are disrupted in many cancers, contributing to multiple processes of tumorigenesis. RNA-binding proteins (RBPs), the main post-transcriptional regulators, often show altered expression and activity in cancer cells. Dysregulation of RBPs contributes to many cancer phenotypes, functioning in complex regulatory networks with other cellular players such as non-coding RNAs, signaling mediators and transcription factors to alter the expression of oncogenes and tumor suppressor genes. RBPs often function combinatorially, based on their binding to target sequences/structures on shared mRNA targets, to regulate the expression of cancer-related genes. This gives rise to cooperativity and competition between RBPs in mRNA binding and resultant functional outcomes in post-transcriptional processes such as mRNA splicing, stability, export and translation. Cooperation and competition is also observed in the case of interaction of RBPs and microRNAs with mRNA targets. RNA structural change is a common mechanism mediating the cooperative/competitive interplay between RBPs and between RBPs and microRNAs. RNA modifications, leading to changes in RNA structure, add a new dimension to cooperative/competitive binding of RBPs to mRNAs, further expanding the RBP regulatory landscape. Therefore, cooperative/competitive interplay between RBPs is a major determinant of the RBP interactome and post-transcriptional regulation of gene expression in cancer cells.

Transcriptome Sequencing Reveals Pathways Related to Proliferation and Differentiation of Shitou Goose Myoblasts

Chinese Shitou goose is a type of large goose with high meat yield. Understanding the genetic regulation of muscle development in Shitou goose would be beneficial to improve the meat production traits of geese. Muscle development is regulated by genes related to myoblast proliferation and differentiation. In this study, the RNA-seq method was used to construct the mRNA and lncRNA expression profiles of Shitou goose myoblasts and myotubes. A total of 1664 differentially expressed (DE) mRNAs and 244 DE-lncRNAs were identified. The alternative mRNA splicing in proliferation and differentiation stages was also analyzed. Notably, pathways enriched in DE-mRNAs, DE-splicing transcripts, and DE-lncRNAs all point to the Wnt signaling pathway, indicating that the Wnt signaling is a key regulatory pathway of muscle development in Shitou goose. We also constructed the interactive network of DE-lncRNAs and DE-mRNAs and revealed some key genes of lncRNAs regulating the proliferation and differentiation of myoblasts. These results provide new insights for the study of the muscle development of the Shitou goose.

F1 Male Sterility in Cattle-Yak Examined through Changes in Testis Tissue and Transcriptome Profiles

Simple Summary

Cattle-yak, a crossbreed of cattle and yak, has evident heterosis but F1 male cattle-yak is unable to generate sperm and is sterile, which limits the fixation of heterosis. This study analyzed the differences in testicular tissue development between four-year-old yak and cattle-yak from the perspective of histomorphological changes and sequenced the testicular tissue of the two using RNA-seq technology, examining the differential gene expression related to spermatogenesis and apoptosis. These findings offer a theoretical explanation for the sterility in F1 male cattle-yak that can help yak hybridization.

Abstract

Male-derived sterility in cattle-yaks, a hybrid deriving from yak and cattle, is a challenging problem. This study compared and analyzed the histomorphological differences in testis between sexually mature yak and cattle-yak, and examined the transcriptome differences employing RNA-seq. The study found that yak seminiferous tubules contained spermatogenic cells at all levels, while cattle-yak seminiferous tubules had reduced spermatogonia (SPG) and primary spermatocyte (Pri-SPC), fewer secondary spermatocytes (Sec-SPC), an absence of round spermatids (R-ST) and sperms (S), and possessed large vacuoles. All of these conditions could have significantly reduced the volume and weight of cattle-yak testis compared to that of yak. RNA-seq analysis identified 8473 differentially expressed genes (DEGs; 3580 upregulated and 4893 downregulated). GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment evaluations for DEGs found their relation mostly to spermatogenesis and apoptosis. Among the DEGs, spermatogonia stem cell (SSCs) marker genes (Gfra1, CD9, SOHLH1, SALL4, ID4, and FOXO1) and genes involved in apoptosis (Fas, caspase3, caspase6, caspase7, caspase8, CTSK, CTSB and CTSC) were significantly upregulated, while differentiation spermatogenic cell marker genes (Ccna1, PIWIL1, TNP1, and TXNDC2) and meiosis-related genes (TEX14, TEX15, MEIOB, STAG3 and M1AP) were significantly downregulated in cattle-yak. Furthermore, the alternative splicing events in cattle-yak were substantially decreased than in yak, suggesting that the lack of protein subtypes could be another reason for spermatogenic arrest in cattle-yak testis.

Regulatory role of phosphoproteins in the development of bovine small intestine during early life

The small intestine is the primary site of nutrient digestion and absorption, which plays a key role in the survival of neonatal calves. A comprehensive assessment of the phosphoproteomic changes in the small intestine of neonatal calves is unavailable; therefore, we used phosphopeptide enrichment coupled with liquid chromatography-tandem mass spectrometry to investigate the changes in the phosphoproteome profile in the bovine small intestine during the first 36 h of life. Twelve neonatal male calves were assigned to one of the following groups: (1) calves not fed colostrum and slaughtered approximately 2 h postpartum (n = 3), (2) calves fed colostrum at 1 to 2 h and slaughtered 8 h postpartum (n = 3), (3) calves fed 2 colostrum meals (at 1-2 and 10-12 h) and slaughtered 24 h postpartum (n = 3), (4) calves fed 3 colostrum meals (at 1-2, 10-12, and 22-24 h) and slaughtered 36 h postpartum (n = 3). Mid-duodenal, jejunal, and ileal samples of the calves were collected after slaughter. We identified 1,678 phosphoproteins with approximately 3,080 phosphosites, which were mainly Ser (89.9%), Thr (9.8%), and Tyr (0.3%) residues; they belonged to the prodirected (52.9%), basic (20.4%), acidic (16.6%), and Tyr-directed (1.7%) motif categories. The regional differentially expressed phosphoproteins included zonula occludens 2, sorting nexin 12, and protein kinase C, which are mainly associated with developmental processes, intracellular transport, vesicle-mediated transport, and immune system process. They are enriched in the endocytosis, tight junction, insulin signaling, and focal adhesion pathways. The temporal differentially expressed phosphoproteins included occludin, epsin 1, and bridging integrator 1, which were mainly associated with macromolecule metabolic process, cell adhesion, and growth. They were enriched in the spliceosomes, adherens junctions, and tight junctions. The observed changes in the phosphoproteins in the tissues of small intestine suggest the protein phosphorylation plays an important role in nutrient transport and immune response of calves during early life, which needs to be confirmed in a larger study.

Inflammatory and Proliferative Pathway Activation in Human Esophageal Myofibroblasts Treated with Acidic Bile Salts

Subepithelial human esophageal myofibroblasts (HEMFs) in gastroesophageal reflux disease (GERD) are exposed to luminal contents via impaired squamous epithelium barrier integrity. The supernatant of HEMFs treated with acidic bile salts reflective of in vivo reflux increases squamous epithelial thickness. We aimed to identify the involved mechanisms using an unbiased approach. Acidic-bile-salt-treated primary HEMF cultures (n = 4) were submitted for RNA-Seq and analyzed with Partek Flow followed by Ingenuity Pathway Analysis (IPA). A total of 1165 molecules (579 downregulated, 586 upregulated) were differentially expressed, with most top regulated molecules either extracellular or in the plasma membrane. Increases in HEMF CXCL-8, IL-6, AREG, and EREG mRNA, and protein secretion were confirmed. Top identified canonical pathways were agranulocyte and granulocyte adhesion and diapedesis, PI3K/AKT signaling, CCR5 signaling in macrophages, and the STAT3 pathway. Top diseases and biological functions were cellular growth and development, hematopoiesis, immune cell trafficking, and cell-mediated response. The targets of the top upstream regulator ErbB2 included CXCL-8, IL-6, and AREG and the inhibition of CXCL-8 in the HEMF supernatant decreased squamous epithelial proliferation. Our work shows an inflammatory/immune cell and proliferative pathways activation in HEMFs in the GERD environment and identifies CXCL-8 as a HEMF-derived chemokine with paracrine proliferative effects on squamous epithelium.

A widespread length-dependent splicing dysregulation in cancer

Dysregulation of alternative splicing is a key molecular hallmark of cancer. However, the common features and underlying mechanisms remain unclear. Here, we report an intriguing length-dependent splicing regulation in cancers. By systematically analyzing the transcriptome of thousands of cancer patients, we found that short exons are more likely to be mis-spliced and preferentially excluded in cancers. Compared to other exons, cancer-associated short exons (CASEs) are more conserved and likely to encode in-frame low-complexity peptides, with functional enrichment in GTPase regulators and cell adhesion. We developed a CASE-based panel as reliable cancer stratification markers and strong predictors for survival, which is clinically useful because the detection of short exon splicing is practical. Mechanistically, mis-splicing of CASEs is regulated by elevated transcription and alteration of certain RNA binding proteins in cancers. Our findings uncover a common feature of cancer-specific splicing dysregulation with important clinical implications in cancer diagnosis and therapies.

Combined treatment with epigenetic agents enhances anti-tumor activity of MAGE-D4 peptide-specific T cells by upregulating the MAGE-D4 expression in glioma

Objective

The study evaluated the efficacy of combined epigenetic drugs of decitabine (DAC), valproic acid (VPA), and trichostatin A (TSA) on immunotherapy against glioma.

Methods

The expression and prognosis of MAGE-D4 in glioma were analyzed online, and the expression of MAGE-D4 and HLA-A2 in glioma induced by epigenetic drugs was detected by qRT-PCR, Western blot, and flow cytometry. The methylation status of the MAGE-D4 promoter was determined by pyrosequencing. An HLA-A2 restricted MAGE-D4 peptide was predicted and synthesized. An affinity assay and a peptide/HLA complex stability assay were performed to determine the affinity between peptide and HLA. CCK8 assay, CFSE assay, ELISA and ELISPOT were performed to detect the function of MAGE-D4 peptide-specific T cells. Flow cytometry, ELISA, and cytotoxicity assays were used to detect the cytotoxicity effect of MAGE-D4 peptide-specific T cells combined with epigenetic drugs against glioma in vitro. Finally, the glioma-loaded mouse model was applied to test the inhibitory effect of specific T cells on gliomas in vivo.

Results

MAGE-D4 was highly expressed in glioma and correlated with poor prognosis. Glioma cells could be induced to express MAGE-D4 and HLA-A2 by epigenetic drugs. MAGE-D4-associated peptides were found that induce DCs to stimulate the highest T-cell activities of proliferation, IL-2 excretion, and IFN-γ secretion. MAGE-D4 peptide-specific T cells treated with TSA only or combining TSA and DAC had the most cytotoxicity effect, and its cytotoxicity effect on glioma cells decreased significantly after HLA blocking. In vivo experiments also confirmed that MAGE-D4-specific T cells inhibit TSA-treated glioma.

Conclusion

MAGE-D4 is highly expressed in glioma and correlated with the prognosis of glioma. The novel MAGE-D4 peptide identified was capable of inducing MAGE-D4-specific T cells that can effectively inhibit glioma growth, and the epigenetic drug application can enhance this inhibition.

Myocardin regulates exon usage in smooth muscle cells through induction of splicing regulatory factors

Differentiation of smooth muscle cells (SMCs) depends on serum response factor (SRF) and its co-activator myocardin (MYOCD). The role of MYOCD for the SMC program of gene transcription is well established. In contrast, the role of MYOCD in control of SMC-specific alternative exon usage, including exon splicing, has not been explored. In the current work we identified four splicing factors (MBNL1, RBPMS, RBPMS2, and RBFOX2) that correlate with MYOCD across human SMC tissues. Forced expression of MYOCD family members in human coronary artery SMCs in vitro upregulated expression of these splicing factors. For global profiling of transcript diversity, we performed RNA-sequencing after MYOCD transduction. We analyzed alternative transcripts with three different methods. Exon-based analysis identified 1637 features with differential exon usage. For example, usage of 3´ exons in MYLK that encode telokin increased relative to 5´ exons, as did the 17 kDa telokin to 130 kDa MYLK protein ratio. Dedicated event-based analysis identified 239 MYOCD-driven splicing events. Events involving MBNL1, MCAM, and ACTN1 were among the most prominent, and this was confirmed using variant-specific PCR analyses. In support of a role for RBPMS and RBFOX2 in MYOCD-driven splicing we found enrichment of their binding motifs around differentially spliced exons. Moreover, knockdown of either RBPMS or RBFOX2 antagonized splicing events stimulated by MYOCD, including those involving ACTN1, VCL, and MBNL1. Supporting an in vivo role of MYOCD-SRF-driven splicing, we demonstrate altered Rbpms expression and splicing in inducible and SMC-specific Srf knockout mice. We conclude that MYOCD-SRF, in part via RBPMS and RBFOX2, induce a program of differential exon usage and alternative splicing as part of the broader program of SMC differentiation.

EventPointer 3.0: flexible and accurate splicing analysis that includes studying the differential usage of protein-domains

Alternative splicing (AS) plays a key role in cancer: all its hallmarks have been associated with different mechanisms of abnormal AS. The improvement of the human transcriptome annotation and the availability of fast and accurate software to estimate isoform concentrations has boosted the analysis of transcriptome profiling from RNA-seq. The statistical analysis of AS is a challenging problem not yet fully solved. We have included in EventPointer (EP), a Bioconductor package, a novel statistical method that can use the bootstrap of the pseudoaligners. We compared it with other state-of-the-art algorithms to analyze AS. Its performance is outstanding for shallow sequencing conditions. The statistical framework is very flexible since it is based on design and contrast matrices. EP now includes a convenient tool to find the primers to validate the discoveries using PCR. We also added a statistical module to study alteration in protein domain related to AS. Applying it to 9514 patients from TCGA and TARGET in 19 different tumor types resulted in two conclusions: i) aberrant alternative splicing alters the relative presence of Protein domains and, ii) the number of enriched domains is strongly correlated with the age of the patients.

Osmoregulation and the Hypothalamic Supraoptic Nucleus: From Genes to Functions

Due to the relatively high permeability to water of the plasma membrane, water tends to equilibrate its chemical potential gradient between the intra and extracellular compartments. Because of this, changes in osmolality of the extracellular fluid are accompanied by changes in the cell volume. Therefore, osmoregulatory mechanisms have evolved to keep the tonicity of the extracellular compartment within strict limits. This review focuses on the following aspects of osmoregulation: 1) the general problems in adjusting the "milieu interieur" to challenges imposed by water imbalance, with emphasis on conceptual aspects of osmosis and cell volume regulation; 2) osmosensation and the hypothalamic supraoptic nucleus (SON), starting with analysis of the electrophysiological responses of the magnocellular neurosecretory cells (MNCs) involved in the osmoreception phenomenon; 3) transcriptomic plasticity of SON during sustained hyperosmolality, to pinpoint the genes coding membrane channels and transporters already shown to participate in the osmosensation and new candidates that may have their role further investigated in this process, with emphasis on those expressed in the MNCs, discussing the relationships of hydration state, gene expression, and MNCs electrical activity; and 4) somatodendritic release of neuropeptides in relation to osmoregulation. Finally, we expect that by stressing the relationship between gene expression and the electrical activity of MNCs, studies about the newly discovered plastic-regulated genes that code channels and transporters in the SON may emerge.

New Insights on the Regulatory Gene Network Disturbed in Central Areolar Choroidal Dystrophy-Beyond Classical Gene Candidates

Central areolar choroidal dystrophy (CACD) is a rare hereditary disease that mainly affects the macula, resulting in progressive and usually profound visual loss. Being part of congenital retinal dystrophies, it may have an autosomal dominant or recessive inheritance and, until now, has no effective treatment. Given the shortage of genotypic information about the disease, this work systematically reviews the literature for CACD-causing genes. Three independent researchers selected 33 articles after carefully searching and filtering the Scielo, Pubmed, Lilacs, Web of Science, Scopus, and Embase databases. Mutations of six genes (PRPH2, GUCA1A, GUCY2D, CDHR1, ABCA4, and TTLL5) are implicated in the monogenic dominant inheritance of CACD. They are functionally related to photoreceptors (either in the phototransduction process, as in the case of GUCY2D, or the recovery of retinal photodegradation in photoreceptors for GUCA1A, or the formation and maintenance of specific structures within photoreceptors for PRPH2). The identified genetic variants do not explain all observed clinical features, calling for further whole-genome and functional studies for this disease. A network analysis with the CACD-related genes identified in the systematic review resulted in the identification of another 20 genes that may influence CACD onset and symptoms. Furthermore, an enrichment analysis allowed the identification of 13 transcription factors and 4 long noncoding RNAs interacting with the products of the previously mentioned genes. If mutated or dysregulated, they may be directly involved in CACD development and related disorders. More than half of the genes identified by bioinformatic tools do not appear in commercial gene panels, calling for more studies about their role in the maintenance of the retina and phototransduction process, as well as for a timely update of these gene panels.

Identifying genes with conserved splicing structure and orthologous isoforms in human, mouse and dog

BACKGROUND

In eukaryote transcriptomes, a significant amount of transcript diversity comes from genes' capacity to generate different transcripts through alternative splicing. Identifying orthologous alternative transcripts across multiple species is of particular interest for genome annotators. However, there is no formal definition of transcript orthology based on the splicing structure conservation. Likewise there is no public dataset benchmark providing groups of orthologous transcripts sharing a conserved splicing structure.

RESULTS

We introduced a formal definition of splicing structure orthology and we predicted transcript orthologs in human, mouse and dog. Applying a selective strategy, we analyzed 2,167 genes and their 18,109 known transcripts and identified a set of 253 gene orthologs that shared a conserved splicing structure in all three species. We predicted 6,861 transcript CDSs (coding sequence), mainly for dog, an emergent model species. Each predicted transcript was an ortholog of a known transcript: both share the same CDS splicing structure. Evidence for the existence of the predicted CDSs was found in external data.

CONCLUSIONS

We generated a dataset of 253 gene triplets, structurally conserved and sharing all their CDSs in human, mouse and dog, which correspond to 879 triplets of spliced CDS orthologs. We have released the dataset both as an SQL database and as tabulated files. The data consists of the 879 CDS orthology groups with their detailed splicing structures, and the predicted CDSs, associated with their experimental evidence. The 6,861 predicted CDSs are provided in GTF files. Our data may contribute to compare highly conserved genes across three species, for comparative transcriptomics at the isoform level, or for benchmarking splice aligners and methods focusing on the identification of splicing orthologs. The data is available at https://data-access.cesgo.org/index.php/s/V97GXxOS66NqTkZ .

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

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

Striking a balance: PIP2 and PIP3 signaling in neuronal health and disease

Phosphoinositides are membrane phospholipids involved in a variety of cellular processes like growth, development, metabolism, and transport. This review focuses on the maintenance of cellular homeostasis of phosphatidylinositol 4,5-bisphosphate (PIP₂), and phosphatidylinositol 3,4,5-trisphosphate (PIP₃). The critical balance of these PIPs is crucial for regulation of neuronal form and function. The activity of PIP₂ and PIP₃ can be regulated through kinases, phosphatases, phospholipases and cholesterol microdomains. PIP₂ and PIP₃ carry out their functions either indirectly through their effectors activating integral signaling pathways, or through direct regulation of membrane channels, transporters, and cytoskeletal proteins. Any perturbations to the balance between PIP₂ and PIP₃ signaling result in neurodevelopmental and neurodegenerative disorders. This review will discuss the upstream modulators and downstream effectors of the PIP₂ and PIP₃ signaling, in the context of neuronal health and disease.

In silico assessment of DNA damage response gene variants associated with head and neck cancer

Head and neck cancer (HNC), the sixth most common cancer globally, stands first in India, especially Northeast India, where tobacco usage is predominant, which introduces various carcinogens leading to malignancies by accumulating DNA damages. Consequently, the present work aimed to predict the impact of significant germline variants in DNA repair and Tumour Suppressor genes on HNC development. WES in Ion Proton^TM platform on 'discovery set' (n = 15), followed by recurrence assessment of the observed variants on 'confirmation set' (n = 40) using Sanger Sequencing was performed on the HNC-prevalent NE Indian populations. Initially, 53 variants were identified, of which seven HNC-linked DNA damage response gene variants were frequent in the studied populations. Different tools ascertained the biological consequences of these variants, of which the non-coding variants viz. EXO1_rs4150018, RAD52_rs6413436, CHD5_rs2746066, HACE1_rs6918700 showed risk, while FLT3_rs2491227 and BMPR1A_rs7074064 conferred protection against HNC by affecting transcriptional regulation and splicing mechanism. Molecular Dynamics Simulation of the full-length p53 model predicted that the observed coding TP53_rs1042522 variant conferred HNC-risk by altering the structural dynamics of the protein, which displayed difficulty in the transition between active and inactive conformations due to high-energy barrier. Subsequent pathway and gene ontology analysis revealed that EXO1, RAD52 and TP53 variants affected the Double-Strand Break Repair pathway, whereas CHD5 and HACE1 variants inactivated DNA repair cascade, facilitating uncontrolled cell proliferation, impaired apoptosis and malignant transformation. Conversely, FLT3 and BMPR1A variants protected against HNC by controlling tumorigenesis, which requires experimental validation. These findings may serve as prognostic markers for developing preventive measures against HNC.

Pax-5 Protein Expression Is Regulated by Transcriptional 3'UTR Editing

The Pax-5 gene encodes a transcription factor that is essential for B-cell commitment and maturation. However, Pax-5 deregulation is associated with various cancer lesions, notably hematopoietic cancers. Mechanistically, studies have characterized genetic alterations within the Pax-5 locus that result in either dominant oncogenic function or haploinsufficiency-inducing mutations leading to oncogenesis. Apart from these mutations, some examples of aberrant Pax-5 expression cannot be associated with genetic alterations. In the present study, we set out to elucidate potential alterations in post-transcriptional regulation of Pax-5 expression and establish that Pax-5 transcript editing represents an important means to aberrant expression. Upon the profiling of Pax-5 mRNA in leukemic cells, we found that the 3′end of the Pax-5 transcript is submitted to alternative polyadenylation (APA) and alternative splicing events. Using rapid amplification of cDNA ends (3′RACE) from polysomal fractions, we found that Pax-5 3′ untranslated region (UTR) shortening correlates with increased ribosomal occupancy for translation. These observations were also validated using reporter gene assays with truncated 3′UTR regions cloned downstream of a luciferase gene. We also showed that Pax-5 3′UTR editing has direct repercussions on regulatory elements such as miRNAs, which in turn impact Pax-5 protein expression. More importantly, we found that advanced staging of various hematopoietic cancer lesions relates to shorter Pax-5 3′UTRs. Altogether, our findings identify novel molecular mechanisms that account for aberrant expression and function of the Pax-5 oncogene in cancer cells. These findings also present new avenues for strategic intervention in Pax-5-mediated cancers.

Genome-wide Analysis of Alternative Gene Splicing Associated with Virulence in the Brown Planthopper Nilaparvata lugens (Hemiptera: Delphacidae)

Alternative splicing of protein coding genes plays a profound role in phenotypic variation for many eukaryotic organisms. The development of high-throughput sequencing and bioinformatics algorithms provides the possibility of genome-wide identification of alternative splicing events in eukaryotes. However, for the brown planthopper Nilaparvata lugens, a destructive pest of rice crops, whole-genome distribution of alternative splicing events and the role of alternative splicing in the phenotypic plasticity of virulence have not previously been estimated. Here, we developed an analysis pipeline to identify alternative splicing events in the genome of N. lugens. Differential expression analysis and functional annotation were performed on datasets related to different virulence phenotypes. In total, 27,880 alternative splicing events corresponding to 9,787 multi-exon genes were detected in N. lugens. Among them, specifically expressed alternative splicing transcripts in the virulent Mudgo population were enriched in metabolic process categories, while transcripts in the avirulent TN1 population were enriched in regulator activity categories. In addition, genes encoding odorant receptor, secreted saliva protein and xenobiotic metabolic P450 monooxygenase showed different splicing patterns between Mudgo population and TN1 population. Host change experiment also revealed that an isoform of a P450 gene could be specially induced by the stimulation of resistant rice variety Mudgo. This research pioneered a genome-wide study of alternative gene splicing in the rice brown planthopper. Differences in alternative splicing between virulent and avirulent populations indicated that alternative splicing might play an important role in the formation of virulence phenotypes in N. lugens.

Biological function and molecular mechanism of SRSF3 in cancer and beyond

Serine/arginine-rich splicing factor 3 (SRSF3; also known as SRp20), an important member of the family of SRSFs, is abnormally expressed in tumors, resulting in aberrant splicing of hub genes, such as CD44, HER2, MDM4, Rac family small GTPase 1 and tumor protein p53. Under normal conditions, the splicing and expression of SRSF3 are strictly regulated. However, the splicing, expression and phosphorylation of SRSF3 are abnormal in tumors. SRSF3 plays important roles in the occurrence and development of tumors, including the promotion of tumorigenesis, cellular proliferation, the cell cycle and metastasis, as well as inhibition of cell senescence, apoptosis and autophagy. SRSF3-knockdown significantly inhibits the proliferation and metastatic characteristics of tumor cells. Therefore, SRSF3 may be suggested as a novel anti-tumor target. The other biological functions of SRSF3 and its regulatory mechanisms are also summarized in the current review.

Evolution of the Early Spliceosomal Complex-From Constitutive to Regulated Splicing

Pre-mRNA splicing is a major process in the regulated expression of genes in eukaryotes, and alternative splicing is used to generate different proteins from the same coding gene. Splicing is a catalytic process that removes introns and ligates exons to create the RNA sequence that codifies the final protein. While this is achieved in an autocatalytic process in ancestral group II introns in prokaryotes, the spliceosome has evolved during eukaryogenesis to assist in this process and to finally provide the opportunity for intron-specific splicing. In the early stage of splicing, the RNA 5' and 3' splice sites must be brought within proximity to correctly assemble the active spliceosome and perform the excision and ligation reactions. The assembly of this first complex, termed E-complex, is currently the least understood process. We focused in this review on the formation of the E-complex and compared its composition and function in three different organisms. We highlight the common ancestral mechanisms in S. cerevisiae, S. pombe, and mammals and conclude with a unifying model for intron definition in constitutive and regulated co-transcriptional splicing.

Antisense Oligonucleotide-Mediated Splice Switching: Potential Therapeutic Approach for Cancer Mitigation

Splicing is an essential process wherein precursor messenger RNA (pre-mRNA) is reshaped into mature mRNA. In alternative splicing, exons of any pre-mRNA get rearranged to form mRNA variants and subsequently protein isoforms, which are distinct both by structure and function. On the other hand, aberrant splicing is the cause of many disorders, including cancer. In the past few decades, developments in the understanding of the underlying biological basis for cancer progression and therapeutic resistance have identified many oncogenes as well as carcinogenic splice variants of essential genes. These transcripts are involved in various cellular processes, such as apoptosis, cell signaling and proliferation. Strategies to inhibit these carcinogenic isoforms at the mRNA level are promising. Antisense oligonucleotides (AOs) have been developed to inhibit the production of alternatively spliced carcinogenic isoforms through splice modulation or mRNA degradation. AOs can also be used to induce splice switching, where the expression of an oncogenic protein can be inhibited by the induction of a premature stop codon. In general, AOs are modified chemically to increase their stability and binding affinity. One of the major concerns with AOs is efficient delivery. Strategies for the delivery of AOs are constantly being evolved to facilitate the entry of AOs into cells. In this review, the different chemical modifications employed and delivery strategies applied are discussed. In addition to that various AOs in clinical trials and their efficacy are discussed herein with a focus on six distinct studies that use AO-mediated exon skipping as a therapeutic strategy to combat cancer.

Aberrant Post-Transcriptional Regulation of Protein Expression in the Development of Chronic Obstructive Pulmonary Disease

Chronic obstructive pulmonary disease (COPD) is an incurable and prevalent respiratory disorder that is characterized by chronic inflammation and emphysema. COPD is primarily caused by cigarette smoke (CS). CS alters numerous cellular processes, including the post-transcriptional regulation of mRNAs. The identification of RNA-binding proteins (RBPs), microRNAs (miRNAs), and long non-coding RNAs (lncRNAs) as main factors engaged in the regulation of RNA biology opens the door to understanding their role in coordinating physiological cellular processes. Dysregulation of post-transcriptional regulation by foreign particles in CS may lead to the development of diseases such as COPD. Here we review current knowledge about post-transcriptional events that may be involved in the pathogenesis of COPD.

A Novel Humanized Model of NASH and Its Treatment With META4, A Potent Agonist of MET

BACKGROUND & AIMS

Nonalcoholic fatty liver disease is a frequent cause of hepatic dysfunction and is now a global epidemic. This ailment can progress to an advanced form called nonalcoholic steatohepatitis (NASH) and end-stage liver disease. Currently, the molecular basis of NASH pathogenesis is poorly understood, and no effective therapies exist to treat NASH. These shortcomings are due to the paucity of experimental NASH models directly relevant to humans.

METHODS

We used chimeric mice with humanized liver to investigate nonalcoholic fatty liver disease in a relevant model. We carried out histologic, biochemical, and molecular approaches including RNA-Seq. For comparison, we used side-by-side human NASH samples.

RESULTS

Herein, we describe a "humanized" model of NASH using transplantation of human hepatocytes into fumarylacetoacetate hydrolase-deficient mice. Once fed a high-fat diet, these mice develop NAFLD faithfully, recapitulating human NASH at the histologic, cellular, biochemical, and molecular levels. Our RNA-Seq analyses uncovered that a variety of important signaling pathways that govern liver homeostasis are profoundly deregulated in both humanized and human NASH livers. Notably, we made the novel discovery that hepatocyte growth factor (HGF) function is compromised in human and humanized NASH at several levels including a significant increase in the expression of the HGF antagonists known as NK1/NK2 and marked decrease in HGF activator. Based on these observations, we generated a potent, human-specific, and stable agonist of human MET that we have named META4 (Metaphor) and used it in the humanized NASH model to restore HGF function.

CONCLUSIONS

Our studies revealed that the humanized NASH model recapitulates human NASH and uncovered that HGF-MET function is impaired in this disease. We show that restoring HGF-MET function by META4 therapy ameliorates NASH and reinstates normal liver function in the humanized NASH model. Our results show that the HGF-MET signaling pathway is a dominant regulator of hepatic homeostasis.

Alternative RNA splicing in stem cells and cancer stem cells: Importance of transcript-based expression analysis

Alternative ribonucleic acid (RNA) splicing can lead to the assembly of different protein isoforms with distinctive functions. The outcome of alternative splicing (AS) can result in a complete loss of function or the acquisition of new functions. There is a gap in knowledge of abnormal RNA splice variants promoting cancer stem cells (CSCs), and their prospective contribution in cancer progression. AS directly regulates the self-renewal features of stem cells (SCs) and stem-like cancer cells. Notably, octamer-binding transcription factor 4A spliced variant of octamer-binding transcription factor 4 contributes to maintaining stemness properties in both SCs and CSCs. The epithelial to mesenchymal transition pathway regulates the AS events in CSCs to maintain stemness. The alternative spliced variants of CSCs markers, including cluster of differentiation 44, aldehyde dehydrogenase, and doublecortin-like kinase, α6β1 integrin, have pivotal roles in increasing self-renewal properties and maintaining the pluripotency of CSCs. Various splicing analysis tools are considered in this study. LeafCutter software can be considered as the best tool for differential splicing analysis and identification of the type of splicing events. Additionally, LeafCutter can be used for efficient mapping splicing quantitative trait loci. Altogether, the accumulating evidence re-enforces the fact that gene and protein expression need to be investigated in parallel with alternative splice variants.

Expression of Myosin 5a splice variants in murine stomach

BACKGROUND

The motor protein, Myosin 5a (Myo5a) is known to play a role in inhibitory neurotransmission in gastric fundus. However, there is no information regarding the relative expression of total Myo5a, or of its alternative exon splice variants, across the stomach. This study investigated the differential distribution of Myo5a variants expressed within distinct anatomical regions of murine stomach.

METHODS

The distribution of Myo5a protein and mRNA in the stomach was assessed by immunofluorescence microscopy and fluorescent in situ hybridization. Quantitative PCR, restriction enzyme analysis, and electrophoresis were used to identify Myo5a splice variants and quantify their expression levels in the fundus, body, antrum, and pylorus.

KEY RESULTS

Myo5a protein colocalized with βIII-Tubulin in the myenteric plexus, and with synaptophysin in nerve fibers. Total Myo5a mRNA expression was lower in pylorus than in antrum, body, or fundus (p < 0.001), which expressed equivalent amounts of Myo5a. However, Myo5a splice variants were differentially expressed across the stomach. While the ABCE splice variant predominated in the antrum and body regions, the ACEF/ACDEF variants were enriched in fundus and pylorus.

CONCLUSIONS AND INFERENCES

Myo5a splice variants varied in their relative expression across anatomically distinguishable stomach regions and might mediate distinct physiological functions in gastric neurotransmission.

Identification of novel OCRL isoforms associated with phenotypic differences between Dent disease-2 and Lowe syndrome

BACKGROUND

Although Lowe syndrome and Dent disease-2 are both caused by OCRL mutations, their clinical severities differ substantially, and their molecular mechanisms remain unclear. Truncating mutations in OCRL exons 1 through 7 lead to Dent disease-2, whereas those in exons 8 through 24 lead to Lowe syndrome. Herein, we identified the mechanism underlying the action of novel OCRL protein isoforms.

METHODS

mRNA samples extracted from cultured urine-derived cells from a healthy control and the Dent disease-2 patient were examined to detect the 5' end of the OCRL isoform. For protein expression and functional analysis, vectors containing (1) the full-length OCRL transcripts, (2) the isoform transcripts, and (3) transcripts with truncating mutations detected in Lowe syndrome and Dent disease-2 patients were transfected into HeLa cells.

RESULTS

We successfully cloned the novel isoform transcripts from OCRL exons 6-24, including the translation-initiation codons present in exon 8. In vitro protein-expression analysis detected proteins of two different sizes (105 and 80 kDa) translated from full-length OCRL, whereas only one protein (80 kDa) was found from the isoform and Dent disease-2 variants. No protein expression was observed for the Lowe syndrome variants. The isoform enzyme activity was equivalent to that of full-length OCRL; the Dent disease-2 variants retained > 50% enzyme activity, whereas the Lowe syndrome variants retained < 20% activity.

CONCLUSIONS

We elucidated the molecular mechanism underlying the two different phenotypes in OCRL-related diseases; the functional OCRL isoform translated starting at exon 8 was associated with this mechanism.

A Novel Splice Site Variant in the LDLRAP1 Gene Causes Familial Hypercholesterolemia

Background:

FH, a hereditary disorder, is caused by pathogenic variants in the LDLR, APOB, and PCSK9 genes. This study has assessed genetic variants in a family, clinically diagnosed with FH.

Methods:

A family was recruited from MASHAD study in Iran with possible FH based on the Simon Broom criteria. The DNA sample of an affected individual (proband) was analyzed using WES, followed by bioinformatics and segregation analyses.

Results:

A novel splice site variant (c.345-2A>G) was detected in the LDLRAP1 gene, which was segregated in all affected family members. Moreover, HMGCR rs3846662 g.23092A>G was found to be homozygous (G/G) in the proband, probably leading to reduced response to simvastatin and pravastatin.

Conclusion:

LDLRAP1 c.345-2A>G could alter the PTB, which acts as an important part of biological pathways related to lipid metabolism.

Detection of alternative splicing in western corn rootworm (Diabrotica virgifera virgifera LeConte) in association with eCry3.1Ab resistance using RNA-seq and PacBio Iso-Seq

Alternative splicing is a common feature in eukaryotes that not only increases the transcript diversity, but also has functional consequences. In insects, alternative splicing has been found associated with resistance to pesticides and Bt toxins. Up to date, the alternative splicing in western corn rootworm (Diabrotica virgifera virgifera LeConte) has not been studied. To investigate its alternative splicing pattern and relation to Bt resistance, we carried out single-molecule real-time (SMRT) transcript sequencing and Iso-seq analysis on resistant, eCry3.1Ab-selected and susceptible, unselected, western corn rootworm neonate midguts which fed on seedling maize with and without eCry3.1Ab for 12 and 24 h. We present transcriptome-wide alternative splicing patterns of western corn rootworm midgut in response to feeding on eCry3.1Ab-expressing corn using a comprehensive approach that combines both RNA-seq and SMRT transcript sequencing techniques. The results showed genes in western corn rootworm are highly alternatively spliced, which happens on 67.73% of multi-exon genes. One of the alternative splicing events we identified was a novel peritrophic matrix protein with two alternative splicing isoforms. Analysis of differential exon usage between resistant and susceptible colonies showed that in eCry3.1Ab-resistant western corn rootworm, expression of one isoform was significantly higher than in the susceptible colony, while no significant differences between colonies were observed with the other isoform. Our results provide the first survey of alternative splicing in western corn rootworm and suggest that the observed alternatively spliced isoforms of peritrophic matrix protein may be associated with eCry3.1Ab resistance in western corn rootworm.

The clinical importance of tandem exon duplication-derived substitutions

Most coding genes in the human genome are annotated with multiple alternative transcripts. However, clear evidence for the functional relevance of the protein isoforms produced by these alternative transcripts is often hard to find. Alternative isoforms generated from tandem exon duplication-derived substitutions are an exception. These splice events are rare, but have important functional consequences. Here, we have catalogued the 236 tandem exon duplication-derived substitutions annotated in the GENCODE human reference set. We find that more than 90% of the events have a last common ancestor in teleost fish, so are at least 425 million years old, and twenty-one can be traced back to the Bilateria clade. Alternative isoforms generated from tandem exon duplication-derived substitutions also have significantly more clinical impact than other alternative isoforms. Tandem exon duplication-derived substitutions have >25 times as many pathogenic and likely pathogenic mutations as other alternative events. Tandem exon duplication-derived substitutions appear to have vital functional roles in the cell and may have played a prominent part in metazoan evolution.

Alternative splicing of BCL-X and implications for treating hematological malignancies

BCL-X is a member of the BCL-2 family. It regulates apoptosis and plays a critical role in hematological malignancies. It is well-known that >90% of human genes undergo alternative splicing. A total of 10 distinct splicing transcripts of the BCL-X gene have been identified, including transcript variants 1–9 and ABALON. Different transcripts from the same gene have different functions. The present review discusses the progress in understanding the different alternative splicing transcripts of BCL-X, including their characteristics, functions and expression patterns. The potential use of BCL-X in targeted therapies for hematological malignancies is also discussed.

Single-molecule, full-length transcript isoform sequencing reveals disease-associated RNA isoforms in cardiomyocytes

Alternative splicing generates differing RNA isoforms that govern phenotypic complexity of eukaryotes. Its malfunction underlies many diseases, including cancer and cardiovascular diseases. Comparative analysis of RNA isoforms at the genome-wide scale has been difficult. Here, we establish an experimental and computational pipeline that performs de novo transcript annotation and accurately quantifies transcript isoforms from cDNA sequences with a full-length isoform detection accuracy of 97.6%. We generate a searchable, quantitative human transcriptome annotation with 31,025 known and 5,740 novel transcript isoforms ( http://steinmetzlab.embl.de/iBrowser/ ). By analyzing the isoforms in the presence of RNA Binding Motif Protein 20 (RBM20) mutations associated with aggressive dilated cardiomyopathy (DCM), we identify 121 differentially expressed transcript isoforms in 107 cardiac genes. Our approach enables quantitative dissection of complex transcript architecture instead of mere identification of inclusion or exclusion of individual exons, as exemplified by the discovery of IMMT isoforms mis-spliced by RBM20 mutations. Thereby we achieve a path to direct differential expression testing independent of an existing annotation of transcript isoforms, providing more immediate biological interpretation and higher resolution transcriptome comparisons.

Double the Fun, Double the Trouble: Paralogs and Homologs Functioning in the Endoplasmic Reticulum

The evolution of eukaryotic genomes has been propelled by a series of gene duplication events, leading to an expansion in new functions and pathways. While duplicate genes may retain some functional redundancy, it is clear that to survive selection they cannot simply serve as a backup but rather must acquire distinct functions required for cellular processes to work accurately and efficiently. Understanding these differences and characterizing gene-specific functions is complex. Here we explore different gene pairs and families within the context of the endoplasmic reticulum (ER), the main cellular hub of lipid biosynthesis and the entry site for the secretory pathway. Focusing on each of the ER functions, we highlight specificities of related proteins and the capabilities conferred to cells through their conservation. More generally, these examples suggest why related genes have been maintained by evolutionary forces and provide a conceptual framework to experimentally determine why they have survived selection.

Consequences for Pancreatic β-Cell Identity and Function of Unregulated Transcript Processing

Mounting evidence suggests a role for alternative splicing (AS) of transcripts in the normal physiology and pathophysiology of the pancreatic β-cell. In the apparent absence of RNA repair systems, RNA decay pathways are likely to play an important role in controlling the stability, distribution and diversity of transcript isoforms in these cells. Around 35% of alternatively spliced transcripts in human cells contain premature termination codons (PTCs) and are targeted for degradation via nonsense-mediated decay (NMD), a vital quality control process. Inflammatory cytokines, whose levels are increased in both type 1 (T1D) and type 2 (T2D) diabetes, stimulate alternative splicing events and the expression of NMD components, and may or may not be associated with the activation of the NMD pathway. It is, however, now possible to infer that NMD plays a crucial role in regulating transcript processing in normal and stress conditions in pancreatic β-cells. In this review, we describe the possible role of Regulated Unproductive Splicing and Translation (RUST), a molecular mechanism embracing NMD activity in relationship to AS and translation of damaged transcript isoforms in these cells. This process substantially reduces the abundance of non-functional transcript isoforms, and its dysregulation may be involved in pancreatic β-cell failure in diabetes.

A comparison of alternative mRNA splicing in the CD4 and CD8 T cell lineages

T cells can be subdivided into a number of different subsets that are defined by their distinct functions. While the specialization of different T cell subsets is partly achieved by the expression of specific genes, the overall transcriptional profiles of all T cells appear very similar. Alternative mRNA splicing is a mechanism that facilitates greater transcript/protein diversity from a limited number of genes, which may contribute to the functional specialization of distinct T cell subsets. In this study we employ a combination of short-read and long-read sequencing technologies to compare alternative mRNA splicing between the CD4 and CD8 T cell lineages. While long-read technology was effective at assembling full-length alternatively spliced transcripts, the low sequencing depth did not facilitate accurate quantitation. On the other hand, short-read technology was ineffective at assembling full-length transcripts but was highly accurate for quantifying expression. We show that integrating long-read and short-read data together achieves a more complete view of transcriptomic diversity. We found that while the overall usage of transcript isoforms was very similar between the CD4 and CD8 lineages, there were numerous alternative spliced mRNA isoforms that were preferentially used by one lineage over the other. These alternative spliced isoforms included ones with different exon usage, exon exclusion or intron inclusion, all of which are expected to significantly alter the protein sequence.

A mechanistic view on the neurotoxic effects of air pollution on central nervous system: risk for autism and neurodegenerative diseases

Many reports have shown a strong association between exposure to neurotoxic air pollutants like heavy metal and particulate matter (PM) as an active participant and neurological disorders. While the effects of these toxic pollutants on cardiopulmonary morbidity have principally been studied, growing evidence has shown that exposure to polluted air is associated with memory impairment, communication deficits, and anxiety/depression among all ages. So, these toxic pollutants in the environment increase the risk of neurodegenerative disease, ischemia, and autism spectrum disorders (ASD). The precise mechanisms in which air pollutants lead to communicative inability, social inability, and declined cognition have remained unknown. Various animal model studies show that amyloid precursor protein (APP), processing, oxidant/antioxidant balance, and inflammation pathways change following the exposure to constituents of polluted air. In the present review study, we collect the probable molecular mechanisms of deleterious CNS effects in response to various air pollutants.

Personalised medicine in hypercholesterolaemia: the role of pharmacogenetics in statin therapy

Statins are the first-line choice in Lipid-lowering therapy to reduce cardiovascular risk. In a continuous attempt to optimise treatment success, there is a need for additional research on genes and related molecular pathways that can determine the efficacy and toxicity of lipid-lowering drugs. Several variations within genes associated with lipid metabolism, including those involved in uptake, distribution and metabolism of statins have been reported. The purpose of this study was to evaluate the effect of genetic variations in the key genes responsible for statins' metabolism and their role in personalised medicine and pharmacogenetic testing (PGx) in patients treated with such drugs. Genetic assessment for specific known SNPs within the most known genes such as ABCG2, SLCO1B1, CYP3A4, and HMGCR, appears likely to predict the efficacy of statin therapy and prevent their side effects but does not necessarily reduce the risk of cardiovascular events. Key Messages Hypercholesterolaemia patients show different response to statin therapy. Several variations within genes associated with statin metabolism have been investigated. Genetic assessment for specific known SNPs within the most known genes may improve the efficacy of statins treatment and prevent their side effects.

Identification and Characterization of Hypothalamic Alternative Splicing Events and Variants in Ovine Fecundity-Related Genes

Simple Summary

Previous studies revealed that alternative splicing (AS) events and gene variants played key roles in reproduction. However, their location and distribution in hypothalamic fecundity-related genes in sheep without the FecB mutation remain largely unknown. In this study, we performed a correlation analysis of transcriptomics and proteomics, and the results suggested several differentially expressed genes (DEGs)/differentially expressed proteins (DEPs), including galectin 3 (LGALS3), aspartoacylase (ASPA) and transthyretin (TTR), could be candidate genes influencing ovine litter size. Further analysis suggested that AS events, single nucleotide polymorphisms (SNPs) and microRNA (miRNA)-binding sites existed in key DEGs/DEPs, such as ASPA and TTR. This study provides a new insight into ovine and even other mammalian reproduction.

Abstract

Previous studies revealed that alternative splicing (AS) events and gene variants played key roles in reproduction; however, their location and distribution in hypothalamic fecundity-related genes in sheep without the FecB mutation remain largely unknown. Therefore, in this study, we described the hypothalamic AS events and variants in differentially expressed genes (DEGs) in Small Tail Han sheep without the FecB mutation at polytocous sheep in the follicular phase vs. monotocous sheep in the follicular phase (PF vs. MF) and polytocous sheep in the luteal phase vs. monotocous sheep in the luteal phase (PL vs. ML) via an RNA-seq study for the first time. We found 39 DEGs with AS events (AS DEGs) in PF vs. MF, while 42 AS DEGs were identified in PL vs. ML. No DEGs with single nucleotide polymorphisms (SNPs) were observed in PF vs. MF, but five were identified in PL vs. ML. We also performed a correlation analysis of transcriptomics and proteomics, and the results suggested several key DEGs/differentially expressed proteins (DEPs), such as galectin 3 (LGALS3) in PF vs. MF and aspartoacylase (ASPA) and transthyretin (TTR) in PL vs. ML, could be candidate genes influencing ovine litter size. In addition, further analyses suggested that AS events, SNPs and miRNA-binding sites existed in key DEGs/DEPs, such as ASPA and TTR. All in all, this study provides a new insight into ovine and even other mammalian reproduction.

Next-generation sequencing reveals alternative L-DOPA decarboxylase (DDC) splice variants bearing novel exons, in human hepatocellular and lung cancer cells

The human L-DOPA decarboxylase (DDC) is an enzyme that displays a pivotal role in metabolic processes. It is implicated in various human disorders, including hepatocellular and lung cancer. Several splice variants of DDC have previously been described, most of which encode for protein isoforms of this enzyme. In the present study, we used next-generation sequencing (NGS) technology along with nested touchdown PCR and Sanger sequencing to identify new splice variants bearing novel exons of the DDC gene, in hepatocellular and lung cancer cell lines. Using an in-house-developed algorithm, we discovered seven novel DDC exons. Next, we determined the structure of ten novel DDC transcripts, three of which contain an open reading frame (ORF) and probably encode for three previously unknown protein isoforms of this enzyme. Future studies should focus on the elucidation of their role in cellular physiology and cancer pathobiology.