The Landscape of Isoform Switches in Human Cancers

Direct and indirect effects of spliceosome disruption compromise gene regulation by Nonsense-Mediated mRNA Decay

Pre-mRNA splicing, carried out in the nucleus by a large ribonucleoprotein machine known as the spliceosome, is functionally and physically coupled to the mRNA surveillance pathway in the cytoplasm called nonsense mediated mRNA decay (NMD). The NMD pathway monitors for premature translation termination signals, which can result from alternative splicing, by relying on the exon junction complex (EJC) deposited on exon-exon junctions by the spliceosome. Recently, multiple genetic screens in human cell lines have identified numerous spliceosome components as putative NMD factors. Using publicly available RNA-seq datasets from K562 and HepG2 cells depleted of 18 different spliceosome components, we find that natural NMD targeted mRNA isoforms are upregulated when members of the catalytic spliceosome are reduced. While some of this increase could be due to widespread pleiotropic effects of spliceosome dysfunction (e.g., reduced expression of NMD factors due to mis-splicing of their mRNAs), we identify that AQR, SF3B1, SF3B4 and CDC40 may have a more direct role in NMD. We also test the hypothesis that increased production of novel NMD substrates may overwhelm the pathway to find a direct correlation between the amount of novel NMD substrates detected and the degree of NMD inhibition observed. Finally, similar transcriptome alterations and NMD substrate upregulation are also observed in cells treated with spliceosome inhibitors and in cells derived from retinitis pigmentosa patients with mutations in PRPF8 and PRPF31. Overall, our results show that regardless of the cause, spliceosome disruption upregulates a broad set of NMD targets, which could contribute to cellular dysfunction in spliceosomopathies.

Tumor-suppressive activities for pogo transposable element derived with KRAB domain via ribosome biogenesis restriction

Transposable elements (TEs) are indispensable for human development, with critical functions in pluripotency and embryogenesis. TE sequences also contribute to human pathologies, especially cancer, with documented activities as cis/trans transcriptional regulators, as sources of non-coding RNAs, and as mutagens that disrupt tumor suppressors. Despite this knowledge, little is known regarding the involvement of TE-derived genes (TEGs) in tumor pathogenesis. Here, systematic analyses of TEG expression across human cancer reveal a prominent role for pogo TE derived with KRAB domain (POGK). We show that POGK acts as a tumor suppressor in triple-negative breast cancer (TNBC) cells and that it couples with the co-repressor TRIM28 to directly block the transcription of ribosomal genes RPS16 and RPS29, in turn causing widespread inhibition of ribosomal biogenesis. We report that POGK undergoes deactivation by isoform switching in clinical TNBC, altogether revealing its exapted activities in tumor growth control.

Alternative splicing modulates chromatin interactome and phase separation of the RIF1 C-terminal domain

RIF1 (RAP1 interacting factor) fulfills diverse roles in DNA double-strand break repair, DNA replication, and nuclear organization. RIF1 is expressed as two splice variants, RIF1-Long (RIF1-L) and RIF1-Short (RIF1-S), from the alternative splicing (AS) of Exon 32 (Ex32) which encodes a 26 aa Ser/Lys-rich cassette peptide in the RIF1 C-terminal domain (CTD). Here we demonstrate that Ex32 inclusion was repressed by DNA damage and oncogenesis but peaked at G2/M phase of the cell cycle. Ex32 splice-in was catalyzed by positive regulators including SRSF1, which bound to Ex32 directly, and negative regulators such as PTBP1 and SRSF3. Isoform proteomics revealed enhanced association of RIF1-L with MDC1, whose recruitment to IR-induced foci was strengthened by RIF1-L. RIF1-L and RIF1-S also exhibited unique phase separation and chromatin-binding characteristics that were regulated by CDK1-dependent CTD phosphorylation. These combined findings suggest that regulated AS affects multiple aspects of RIF1 function in genome protection and organization.

Sex-specific gene expression differences in the prefrontal cortex of major depressive disorder individuals

Major depressive disorder (MDD) is a leading global cause of disability, being more prevalent in females, possibly due to molecular and neuronal pathway differences between females and males. However, the connection between transcriptional changes and MDD remains unclear. We identified transcriptionally altered genes (TAGs) in MDD through gene and transcript expression analyses, focusing on sex-specific differences. Analyzing 263 brain samples from both sexes, we conducted differential gene expression, differential transcript expression, and differential transcript usage analyses, revealing 1169 unique TAGs, primarily in the prefrontal areas, with nearly half exhibiting transcript-level alterations. Females showed notable RNA splicing and export process disruptions in the orbitofrontal cortex, alongside altered DDX39B gene expression in five of the six brain regions in both sexes. Our findings suggest that disruptions in RNA processing pathways may play a vital role in MDD.

Epstein-Barr virus induces host shutoff extensively via BGLF5-independent mechanisms

Epstein-Barr virus (EBV) is a ubiquitous oncogenic virus associated with multiple cancers and autoimmune diseases. Unlike most herpesviruses, EBV reactivation from latency occurs asymptomatically, allowing it to spread efficiently to other hosts. However, available models are limited by the inefficient and asynchronous reactivation from latency into lytic replication. To address this problem, we develop a dual-fluorescent lytic reporter (DFLR) EBV that specifically labels cells in the early and late stages of replication. Using lymphoblastoid cell lines transformed by DFLR EBV as a model for EBV reactivation in B cells, we observe extensive reprogramming of the host cell transcriptome during lytic cycle progression. This includes widespread shutoff of host gene expression and disruption of mRNA processing. Unexpectedly, host shutoff remains extensive even in cells infected with DFLR EBV deleted for the BGLF5 nuclease. These findings implicate BGLF5-independent mechanisms as the primary drivers of host transcriptome remodeling during EBV lytic replication.

Isoform-Level Transcriptome Analysis of Peripheral Blood Mononuclear Cells from Breast Cancer Patients Identifies a Disease-Associated RASGEF1A Isoform

Background: Breast cancer (BC) comprises multiple subtypes with distinct molecular features, which differ in their interplay with host immunity, prognosis, and treatment. Non-invasive blood analyses can provide valuable insights into systemic immunity during cancer. The aim of this study was to analyze the expression of transcriptional isoforms in peripheral blood mononuclear cells (PBMCs) from BC patients and healthy women to identify potential BC immune biomarkers. Methods: RNA sequencing and isoform-level bioinformatics were performed on PBMCs from 12 triple-negative and 13 luminal A patients. Isoform expression validation by qRT-PCR and clinicopathological correlations were performed in a larger cohort (156 BC patients and 32 healthy women). Results: Transcriptional analyses showed a significant (p < 0.001) decrease in the ENST00000374459 RASGEF1A isoform in PBMCs of BC compared to healthy subjects, indicating disease-related expression changes. The decrease was associated with higher ctDNA and Ki-67 values. Conclusions: The levels of the RASGEF1A transcriptional isoform ENST00000374459 may have the potential to distinguish between BC and healthy subjects. The downregulation of ENST00000374459 in breast cancer is associated with higher proliferation and ctDNA shedding. Specialized bioinformatics analyses such as isoform analyses hold significant promise in the detection of biomarkers, since standard RNA sequencing analyses may overlook specific transcriptional changes that may be disease-associated and biologically important.

In silico and in cellulo approaches for functional annotation of human protein splice variants

The elegance of pre-mRNA splicing mechanisms continues to interest scientists even after over a half century, since the discovery of the fact that coding regions in genes are interrupted by non-coding sequences. The vast majority of human genes have several mRNA variants, coding structurally and functionally different protein isoforms in a tissue-specific manner and with a linkage to specific developmental stages of the organism. Alteration of splicing patterns shifts the balance of functionally distinct proteins in living systems, distorts normal molecular pathways, and may trigger the onset and progression of various pathologies. Over the past two decades, numerous studies have been conducted in various life sciences disciplines to deepen our understanding of splicing mechanisms and the extent of their impact on the functioning of living systems. This review aims to summarize experimental and computational approaches used to elucidate the functions of splice variants of a single gene based on our experience accumulated in the laboratory of interactomics of proteoforms at the Institute of Biomedical Chemistry (IBMC) and best global practices.

Comparison of Peptidomes Extracted from Healthy Tissue and Tumor Tissue of the Parotid Glands and Saliva Samples

Salivary gland tumors are highly variable in clinical presentation and histology. The World Health Organization (WHO) classifies 22 types of malignant and 11 types of benign tumors of the salivary glands. Diagnosis of salivary gland tumors is based on imaging (ultrasound, magnetic resonance imaging) and fine-needle aspiration biopsy, but the final diagnosis is based on histopathological examination of the removed tumor tissue. In this pilot study, we are testing a new approach to identifying peptide biomarkers in saliva that can be used to diagnose salivary gland tumors. The research material for the peptidomic studies was extracts from washings of neoplastic tissues and healthy tissues (control samples). At the same time, saliva samples from patients and healthy individuals were analyzed. The comparison of the peptidome composition of tissue extracts and saliva samples may allow the identification of potential peptide markers of salivary gland tumors in patients' saliva. The peptidome compositions extracted from 18 tumor and 18 healthy tissue samples, patients' saliva samples (11 samples), and healthy saliva samples (8 samples) were analyzed by LC-MS tandem mass spectrometry. A group of 109 peptides was identified that were present only in the tumor tissue extracts and in the patients' saliva samples. Some of the identified peptides were derived from proteins previously suggested as potential biomarkers of salivary gland tumors (ANXA1, BPIFA2, FGB, GAPDH, HSPB1, IGHG1, VIM) or tumors of other tissues or organs (SERPINA1, APOA2, CSTB, GSTP1, S100A8, S100A9, TPI1). Unfortunately, none of the identified peptides were present in all samples analyzed. This may be due to the high heterogeneity of this type of cancer. The surprising result was that extracts from tumor tissue did not contain peptides derived from salivary gland-specific proteins (STATH, SMR3B, HTN1, HTN3). These results could suggest that the developing tumor suppresses the production of proteins that are essential components of saliva.

Identifying transcriptomic profiles in ovine spleen after repetitive vaccination

Aluminum hydroxide has long been employed as a vaccine adjuvant for its safety profile, although its precise mechanism of action remains elusive. In this study, we investigated the transcriptomic responses in sheep spleen following repetitive vaccination with aluminum adjuvanted vaccines and aluminum hydroxide alone. Notably, this work represents the first exploration of the sheep spleen transcriptome in such conditions. Animals were splitted in 3 treatment groups: vaccine group, adjuvant alone group and control group. A total of 18 high-depth RNA-seq libraries were sequenced, resulting in a rich dataset which also allowed isoform-level analysis. The comparisons between vaccine-treated and control groups (V vs C) as well as between vaccine-treated and adjuvant-alone groups (V vs A) revealed significant alterations in gene expression profiles, including protein coding genes and long non-coding RNAs. Among the differentially expressed genes, many were associated with processes such as endoplasmic reticulum (ER) stress, immune response and cell cycle. The analysis of co-expression modules further indicated a correlation between vaccine treatment and genes related to ER stress and unfolded protein response. Surprisingly, adjuvant-alone treatment had little impact on the spleen transcriptome. Additionally, the role of alternative splicing in the immune response was explored. We identified isoform switches in genes associated with immune regulation and inflammation, potentially influencing protein function. In conclusion, this study provides valuable insights into the transcriptomic changes in sheep spleen following vaccination with aluminum adjuvanted vaccines and aluminum hydroxide alone. These findings shed light on the molecular mechanisms underlying vaccine-induced immune responses and emphasize the significance of antigenic components in aluminum adjuvant mechanism of action. Furthermore, the analysis of alternative splicing revealed an additional layer of complexity in the immune response to vaccination in a livestock species.

Modulation of metabolic and immunoregulatory pathways in the gut transcriptome of Atlantic salmon (Salmo salar L.) after early nutritional programming during first feeding with plant-based diet

Introduction

Plant-based nutritional programming is the concept of exposing fish at very early life stages to a plant-based diet for a short duration to improve physiological responses when exposed to a similar plant-rich diet at a later developmental stage. The mechanisms of action underlying nutritional programming have not been fully deciphered, and the responses may be controlled at multiple levels.

Methods

This 22-week study examines gut transcriptional changes after nutritional programming. Triplicate groups of Atlantic salmon were fed with a plant (V) vs. a marine-rich (M, control) diet for 2 weeks (stimulus phase) at the first exogenous feeding. Both stimulus fish groups (M and V fish) were then fed the M diet for 12 weeks (intermediate phase) and lastly fed the V diet (challenge phase) for 6 weeks, generating two dietary regimes (MMV and VMV) across phases. This study used a whole-transcriptome approach to analyse the effects of the V diet at the end of stimulus (short-term effects) and 22 weeks post-first feeding (long-term effects). After the stimulus, due to its developmental stage, the whole intestine was used, whereas, after the challenge, pyloric caeca and middle and distal intestines were examined.

Results and discussion

At the stimulus end, genes with increased expression in V fish enriched pathways including regulatory epigenetic responses and lipid metabolism, and genes involved in innate immune response were downregulated. In the middle intestine at the end of the challenge, expression levels of genes of lipid, carbohydrate, and energy metabolism were increased in V fish, while M fish revealed increased expression of genes associated with autoimmune and acute adaptive immune response. The distal intestine of V fish showed increased expression of genes associated with immune response and potential immune tolerance. Conversely, the distal intestine of M fish at challenge revealed upregulation of lipid and carbohydrate metabolic pathways, tissue degeneration, and apoptotic responses. The present study demonstrated nutritional programming-associated changes in the intestinal transcriptome, with altered expression of genes involved in both immune responses and different metabolic processes. While there were limited changes in growth between the groups, the results show that there were transcriptional differences, suggesting a programming response, although the mechanism of this response still requires to be fully elucidated.

Role of the histone variant H2A.Z.1 in memory, transcription, and alternative splicing is mediated by lysine modification

Creating long-lasting memories requires learning-induced changes in gene expression, which are impacted by epigenetic modifications of DNA and associated histone proteins. Post-translational modifications (PTMs) of histones are key regulators of transcription, with different PTMs producing unique effects on gene activity and behavior. Although recent studies implicate histone variants as novel regulators of memory, effects of PTMs on the function of histone variants are rarely considered. We previously showed that the histone variant H2A.Z suppresses memory, but it is unclear if this role is impacted by H2A.Z acetylation, a PTM that is typically associated with positive effects on transcription and memory. To answer this question, we used a mutation approach to manipulate acetylation on H2A.Z without impacting acetylation of other histone types. Specifically, we used adeno-associated virus (AAV) constructs to overexpress mutated H2A.Z.1 isoforms that either mimic acetylation (acetyl-mimic) by replacing lysines 4, 7 and 11 with glutamine (KQ), or H2A.Z.1 with impaired acetylation (acetyl-defective) by replacing the same lysines with alanine (KA). Expressing the H2A.Z.1 acetyl-mimic (H2A.Z.1KQ) improved memory under weak learning conditions, whereas expressing the acetyl-defective H2A.Z.1KA generally impaired memory, indicating that the effect of H2A.Z.1 on memory depends on its acetylation status. RNA sequencing showed that H2A.Z.1KQ and H2A.Z.1KA uniquely impact the expression of different classes of genes in both females and males. Specifically, H2A.Z.1KA preferentially impacts genes involved in synaptic function, suggesting that acetyl-defective H2A.Z.1 impairs memory by altering synaptic regulation. Finally, we describe, for the first time, that H2A.Z is also involved in alternative splicing of neuronal genes, whereby H2A.Z depletion, as well as expression of H2A.Z.1 lysine mutants influence transcription and splicing of different gene targets, suggesting that H2A.Z.1 can impact behavior through effects on both splicing and gene expression. This is the first study to demonstrate that direct manipulation of H2A.Z post-translational modifications regulates memory, whereby acetylation adds another regulatory layer by which histone variants can fine tune higher brain functions through effects on gene expression and splicing.

Modifications of lipid pathways restrict SARS-CoV-2 propagation in human induced pluripotent stem cell-derived 3D airway organoids

BACKGROUND

Modifications of lipid metabolism were closely associated with the manifestations and prognosis of coronavirus disease of 2019 (COVID-19). Pre-existing metabolic conditions exacerbated the severity of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection while modulations of aberrant lipid metabolisms alleviated the manifestations. To elucidate the underlying mechanisms, an experimental platform that reproduces human respiratory physiology is required.

METHODS

Here we generated induced pluripotent stem cell-derived airway organoids (iPSC-AOs) that resemble the human native airway. Single-cell sequencing (ScRNAseq) and microscopic examination verified the cellular heterogeneity and microstructures of iPSC-AOs, respectively. We subjected iPSC-AOs to SARS-CoV-2 infection and investigated the treatment effect of lipid modifiers statin drugs on viral pathogenesis, gene expression, and the intracellular trafficking of the SARS-CoV-2 entry receptor angiotensin-converting enzyme-2 (ACE-2).

RESULTS

In SARS-CoV-2-infected iPSC-AOs, immunofluorescence staining detected the SARS-CoV-2 spike (S) and nucleocapsid (N) proteins and bioinformatics analysis further showed the aberrant enrichment of lipid-associated pathways. In addition, SARS-CoV-2 hijacked the host RNA replication machinery and generated the new isoforms of a high-density lipoprotein constituent apolipoprotein A1 (APOA1) and the virus-scavenging protein deleted in malignant brain tumors 1 (DMBT1). Manipulating lipid homeostasis using cholesterol-lowering drugs (e.g. Statins) relocated the viral entry receptor angiotensin-converting enzyme-2 (ACE-2) and decreased N protein expression, leading to the reduction of SARS-CoV-2 entry and replication. The same lipid modifications suppressed the entry of luciferase-expressing SARS-CoV-2 pseudoviruses containing the S proteins derived from different SARS-CoV-2 variants, i.e. wild-type, alpha, delta, and omicron.

CONCLUSIONS

Together, our data demonstrated that modifications of lipid pathways restrict SARS-CoV-2 propagation in the iPSC-AOs, which the inhibition is speculated through the translocation of ACE2 from the cell membrane to the cytosol. Considering the highly frequent mutation and generation of SARS-CoV-2 variants, targeting host metabolisms of cholesterol or other lipids may represent an alternative approach against SARS-CoV-2 infection.

Upstream open reading frames: new players in the landscape of cancer gene regulation

The translation of RNA by ribosomes represents a central biological process and one of the most dysregulated processes in cancer. While translation is traditionally thought to occur exclusively in the protein-coding regions of messenger RNAs (mRNAs), recent transcriptome-wide approaches have shown abundant ribosome activity across diverse stretches of RNA transcripts. The most common type of this kind of ribosome activity occurs in gene leader sequences, also known as 5' untranslated regions (UTRs) of the mRNA, that precede the main coding sequence. Translation of these upstream open reading frames (uORFs) is now known to occur in upwards of 25% of all protein-coding genes. With diverse functions from RNA regulation to microprotein generation, uORFs are rapidly igniting a new arena of cancer biology, where they are linked to cancer genetics, cancer signaling, and tumor-immune interactions. This review focuses on the contributions of uORFs and their associated 5'UTR sequences to cancer biology.

Mapping medically relevant RNA isoform diversity in the aged human frontal cortex with deep long-read RNA-seq

Determining whether the RNA isoforms from medically relevant genes have distinct functions could facilitate direct targeting of RNA isoforms for disease treatment. Here, as a step toward this goal for neurological diseases, we sequenced 12 postmortem, aged human frontal cortices (6 Alzheimer disease cases and 6 controls; 50% female) using one Oxford Nanopore PromethION flow cell per sample. We identified 1,917 medically relevant genes expressing multiple isoforms in the frontal cortex where 1,018 had multiple isoforms with different protein-coding sequences. Of these 1,018 genes, 57 are implicated in brain-related diseases including major depression, schizophrenia, Parkinson's disease and Alzheimer disease. Our study also uncovered 53 new RNA isoforms in medically relevant genes, including several where the new isoform was one of the most highly expressed for that gene. We also reported on five mitochondrially encoded, spliced RNA isoforms. We found 99 differentially expressed RNA isoforms between cases with Alzheimer disease and controls.

Detecting differential transcript usage in complex diseases with SPIT

Differential transcript usage (DTU) plays a crucial role in determining how gene expression differs among cells, tissues, and developmental stages, contributing to the complexity and diversity of biological systems. In abnormal cells, it can also lead to deficiencies in protein function and underpin disease pathogenesis. Analyzing DTU via RNA sequencing (RNA-seq) data is vital, but the genetic heterogeneity in populations with complex diseases presents an intricate challenge due to diverse causal events and undetermined subtypes. Although the majority of common diseases in humans are categorized as complex, state-of-the-art DTU analysis methods often overlook this heterogeneity in their models. We therefore developed SPIT, a statistical tool that identifies predominant subgroups in transcript usage within a population along with their distinctive sets of DTU events. This study provides comprehensive assessments of SPIT's methodology and applies it to analyze brain samples from individuals with schizophrenia, revealing previously unreported DTU events in six candidate genes.

Quantifying transcriptome diversity: a review

Following the central dogma of molecular biology, gene expression heterogeneity can aid in predicting and explaining the wide variety of protein products, functions and, ultimately, heterogeneity in phenotypes. There is currently overlapping terminology used to describe the types of diversity in gene expression profiles, and overlooking these nuances can misrepresent important biological information. Here, we describe transcriptome diversity as a measure of the heterogeneity in (1) the expression of all genes within a sample or a single gene across samples in a population (gene-level diversity) or (2) the isoform-specific expression of a given gene (isoform-level diversity). We first overview modulators and quantification of transcriptome diversity at the gene level. Then, we discuss the role alternative splicing plays in driving transcript isoform-level diversity and how it can be quantified. Additionally, we overview computational resources for calculating gene-level and isoform-level diversity for high-throughput sequencing data. Finally, we discuss future applications of transcriptome diversity. This review provides a comprehensive overview of how gene expression diversity arises, and how measuring it determines a more complete picture of heterogeneity across proteins, cells, tissues, organisms and species.

Functional requirement of alternative splicing in epithelial-mesenchymal transition of pancreatic circulating tumor

Circulating tumor cells (CTCs) that undergo epithelial-to-mesenchymal transition (EMT) can provide valuable information regarding metastasis and potential therapies. However, current studies on the EMT overlook alternative splicing. Here, we used single-cell full-length transcriptome data and mRNA sequencing of CTCs to identify stage-specific alternative splicing of partial EMT and mesenchymal states during pancreatic cancer metastasis. We classified definitive tumor and normal epithelial cells via genetic aberrations and demonstrated dynamic changes in the epithelial-mesenchymal continuum in both epithelial cancer cells and CTCs. We provide the landscape of alternative splicing in CTCs at different stages of EMT, uncovering cell-type-specific splicing patterns and splicing events in cell surface proteins suitable for therapies. We show that MBNL1 governs cell fate through alternative splicing independently of changes in gene expression and affects the splicing pattern during EMT. We found a high frequency of events that contained multiple premature termination codons and were enriched with C and G nucleotides in close proximity, which influence the likelihood of stop codon readthrough and expand the range of potential therapeutic targets. Our study provides insights into the EMT transcriptome's dynamic changes and identifies potential diagnostic and therapeutic targets in pancreatic cancer.

Isoform alterations in the ubiquitination machinery impacting gastrointestinal malignancies

The advancement of RNAseq and isoform-specific expression platforms has led to the understanding that isoform changes can alter molecular signaling to promote tumorigenesis. An active area in cancer research is uncovering the roles of ubiquitination on spliceosome assembly contributing to transcript diversity and expression of alternative isoforms. However, the effects of isoform changes on functionality of ubiquitination machineries (E1, E2, E3, E4, and deubiquitinating (DUB) enzymes) influencing onco- and tumor suppressor protein stabilities is currently understudied. Characterizing these changes could be instrumental in improving cancer outcomes via the identification of novel biomarkers and targetable signaling pathways. In this review, we focus on highlighting reported examples of direct, protein-coded isoform variation of ubiquitination enzymes influencing cancer development and progression in gastrointestinal (GI) malignancies. We have used a semi-automated system for identifying relevant literature and applied established systems for isoform categorization and functional classification to help structure literature findings. The results are a comprehensive snapshot of known isoform changes that are significant to GI cancers, and a framework for readers to use to address isoform variation in their own research. One of the key findings is the potential influence that isoforms of the ubiquitination machinery have on oncoprotein stability.

Downregulation of Ribosomal Protein Genes Is Revealed in a Model of Rat Hippocampal Neuronal Culture Activation with GABA(A)R/GlyRa2 Antagonist Picrotoxin

Long-read transcriptome sequencing provides us with a convenient tool for the thorough study of biological processes such as neuronal plasticity. Here, we aimed to perform transcriptional profiling of rat hippocampal primary neuron cultures after stimulation with picrotoxin (PTX) to further understand molecular mechanisms of neuronal activation. To overcome the limitations of short-read RNA-Seq approaches, we performed an Oxford Nanopore Technologies MinION-based long-read sequencing and transcriptome assembly of rat primary hippocampal culture mRNA at three time points after the PTX activation. We used a specific approach to exclude uncapped mRNAs during sample preparation. Overall, we found 23,652 novel transcripts in comparison to reference annotations, out of which ~6000 were entirely novel and mostly transposon-derived loci. Analysis of differentially expressed genes (DEG) showed that 3046 genes were differentially expressed, of which 2037 were upregulated and 1009 were downregulated at 30 min after the PTX application, with only 446 and 13 genes differentially expressed at 1 h and 5 h time points, respectively. Most notably, multiple genes encoding ribosomal proteins, with a high basal expression level, were downregulated after 30 min incubation with PTX; we suggest that this indicates redistribution of transcriptional resources towards activity-induced genes. Novel loci and isoforms observed in this study may help us further understand the functional mRNA repertoire in neuronal plasticity processes. Together with other NGS techniques, differential gene expression analysis of sequencing data obtained using MinION platform might provide a simple method to optimize further study of neuronal plasticity.

Identification of unique genomic signatures in patients with fibromyalgia and chronic pain

Fibromyalgia (FM) is a chronic pain syndrome characterized by widespread pain. The pathophysiology of fibromyalgia is not clearly understood and there are no specific biomarkers available for accurate diagnosis. Here we define genomic signatures using high throughput RNA sequencing on 96 fibromyalgia and 93 control cases. Our findings revealed three major fibromyalgia-associated expression signatures. The first group included 43 patients with a signature enriched for gene expression associated with extracellular matrix and downregulation of RhoGDI signaling pathway. The second group included 30 patients and showed a profound reduction in the expression of inflammatory mediators with an increased expression of genes involved in the CLEAR signaling pathway. These results suggest defective tissue homeostasis associated with the extra-cellular matrix and cellular program that regulates lysosomal biogenesis and participates in macromolecule clearance in fibromyalgia. The third group of 17 FM patients showed overexpression of pathways that control acute inflammation and dysfunction of the global transcriptional process. The result of this study indicates that FM is a heterogeneous and complex disease. Further elucidation of these pathways will lead to the development of accurate diagnostic markers, and effective therapeutic options for fibromyalgia.

Rbm3 deficiency leads to transcriptome-wide splicing alterations

Rbm3 (RNA-binding motif protein 3) is a stress responsive gene, which maintains cellular homeostasis and promotes survival upon various harmful cellular stimuli. Rbm3 protein shows conserved structural and molecular similarities to heterogeneous nuclear ribonucleoproteins (hnRNPs), which regulate all steps of the mRNA metabolism. Growing evidence is pointing towards a broader role of Rbm3 in various steps of gene expression. Here, we demonstrate that Rbm3 deficiency is linked to transcriptome-wide pre-mRNA splicing alterations, which can be reversed through Rbm3 co-expression from a cDNA. Using an MS2 tethering assay, we show that Rbm3 regulates splice site selection similar to other hnRNP proteins when recruited between two competing 5  ' splice sites. Furthermore, we show that the N-terminal part of Rbm3 encompassing the RNA recognition motif (RRM), is sufficient to elicit changes in splice site selection. On the basis of these findings, we propose a novel, undescribed function of Rbm3 in RNA splicing that contributes to the preservation of transcriptome integrity.

Long read sequencing to reveal the full complexity of a plant transcriptome by targeting both standard and long workflows

BACKGROUND

Long read sequencing allows the analysis of full-length transcripts in plants without the challenges of reliable transcriptome assembly. Long read sequencing of transcripts from plant genomes has often utilized sized transcript libraries. However, the value of including libraries of differing sizes has not been established.

METHODS

A comprehensive transcriptome of the leaves of Jojoba (Simmondsia chinensis) was generated from two different PacBio library preparations: standard workflow (SW) and long workflow (LW).

RESULTS

The importance of using both transcript groups in the analysis was demonstrated by the high proportion of unique sequences (74.6%) that were not shared between the groups. A total of 37.8% longer transcripts were only detected in the long dataset. The completeness of the combined transcriptome was indicated by the presence of 98.7% of genes predicted in the jojoba male reference genome. The high coverage of the transcriptome was further confirmed by BUSCO analysis showing the presence of 96.9% of the genes from the core viridiplantae_odb10 lineage. The high-quality isoforms post Cd-Hit merged dataset of the two workflows had a total of 167,866 isoforms. Most of the transcript isoforms were protein-coding sequences (71.7%) containing open reading frames (ORFs) ≥ 100 amino acids (aa). Alternative splicing and intron retention were the basis of most transcript diversity when analysed at the whole genome level and by specific analysis of the apetala2 gene families.

CONCLUSION

This suggests the need to specifically target the capture of longer transcripts to provide more comprehensive genome coverage in plant transcriptome analysis and reveal the high level of alternative splicing.

Expression and splicing mediate distinct biological signals

BACKGROUND

Through alternative splicing, most human genes produce multiple isoforms in a cell-, tissue-, and disease-specific manner. Numerous studies show that alternative splicing is essential for development, diseases, and their treatments. Despite these important examples, the extent and biological relevance of splicing are currently unknown.

RESULTS

To solve this problem, we developed pairedGSEA and used it to profile transcriptional changes in 100 representative RNA-seq datasets. Our systematic analysis demonstrates that changes in splicing, on average, contribute to 48.1% of the biological signal in expression analyses. Gene-set enrichment analysis furthermore indicates that expression and splicing both convey shared and distinct biological signals.

CONCLUSIONS

These findings establish alternative splicing as a major regulator of the human condition and suggest that most contemporary RNA-seq studies likely miss out on critical biological insights. We anticipate our results will contribute to the transition from a gene-centric to an isoform-centric research paradigm.

Isoform switching leads to downregulation of cytokine producing genes in estrogen receptor positive breast cancer

Objective: Estrogen receptor breast cancer (BC) is characterized by the expression of estrogen receptors. It is the most common cancer among women, with an incidence rate of 2.26 million cases worldwide. The aim of this study was to identify differentially expressed genes and isoform switching between estrogen receptor positive and triple negative BC samples. Methods: The data were collected from ArrayExpress, followed by preprocessing and subsequent mapping from HISAT2. Read quantification was performed by StringTie, and then R package ballgown was used to perform differential expression analysis. Functional enrichment analysis was conducted using Enrichr, and then immune genes were shortlisted based on the ScType marker database. Isoform switch analysis was also performed using the IsoformSwitchAnalyzeR package. Results: A total of 9,771 differentially expressed genes were identified, of which 86 were upregulated and 117 were downregulated. Six genes were identified as mainly associated with estrogen receptor positive BC, while a novel set of ten genes were found which have not previously been reported in estrogen receptor positive BC. Furthermore, alternative splicing and subsequent isoform usage in the immune system related genes were determined. Conclusion: This study identified the differential usage of isoforms in the immune system related genes in cancer cells that suggest immunosuppression due to the dysregulation of CXCR chemokine receptor binding, iron ion binding, and cytokine activity.

The Role of NCS1 in Immunotherapy and Prognosis of Human Cancer

The Neural Calcium Sensor1 (NCS1) is a crucial protein that binds to Ca2+ and is believed to play a role in regulating tumor invasion and cell proliferation. However, the role of NCS1 in immune infiltration and cancer prognosis is still unknown. Our study aimed to explore the expression profile, immune infiltration pattern, prognostic value, biological function, and potential compounds targeting NCS1 using public databases. High expression of NCS1 was detected by immune histochemical staining in LIHC (Liver hepatocellular carcinoma), BRCA (Breast invasive carcinoma), KIRC (Kidney renal clear cell carcinoma), and SKCM (Skin Cutaneous Melanoma). The expression of NCS1 in cancer was determined by TCGA (The Cancer Genome Atlas Program), GTEx (The Genotype-Tissue Expression), the Kaplan-Meier plotter, GEO (Gene Expression Omnibus), GEPIA2.0 (Gene Expression Profiling Interactive Analysis 2.0), HPA (The Human Protein Atlas), UALCAN, TIMER2.0, TISIDB, Metascape, Drugbank, chEMBL, and ICSDB databases. NCS1 has genomic mutations as well as aberrant DNA methylation in multiple cancers compared to normal tissues. Also, NCS1 was significantly different in the immune microenvironment, tumor mutational burden (TMB), microsatellite instability (MSI), and immune infiltrate-associated cells in different cancers, which could be used for the typing of immune and molecular subtypes of cancer and the presence of immune checkpoint resistance in several cancers. Univariate regression analysis, multivariate regression analysis, and gene enrichment analysis to construct prognostic models revealed that NCS1 is involved in immune regulation and can be used as a prognostic biomarker for SKCM, LIHC, BRCA, COAD, and KIRC. These results provide clues from a bioinformatic perspective and highlight the importance of NCS1 in a variety of cancers.

Most protein domains exist as variants with distinct functions across cells, tissues and diseases

Protein domains are the active subunits that provide proteins with specific functions through precise three-dimensional structures. Such domains facilitate most protein functions, including molecular interactions and signal transduction. Currently, these protein domains are described and analyzed as invariable molecular building blocks with fixed functions. Here, I show that most human protein domains exist as multiple distinct variants termed 'domain isotypes'. Domain isotypes are used in a cell, tissue and disease-specific manner and have surprisingly different 3D structures. Accordingly, domain isotypes, compared to each other, modulate or abolish the functionality of protein domains. These results challenge the current view of protein domains as invariable building blocks and have significant implications for both wet- and dry-lab workflows. The extensive use of protein domain isotypes within protein isoforms adds to the literature indicating we need to transition to an isoform-centric research paradigm.

mRNA 3'UTR lengthening by alternative polyadenylation attenuates inflammatory responses and correlates with virulence of Influenza A virus

Changes of mRNA 3'UTRs by alternative polyadenylation (APA) have been associated to numerous pathologies, but the mechanisms and consequences often remain enigmatic. By combining transcriptomics, proteomics and recombinant viruses we show that all tested strains of IAV, including A/PR/8/34(H1N1) (PR8) and A/Cal/07/2009 (H1N1) (Cal09), cause APA. We mapped the effect to the highly conserved glycine residue at position 184 (G184) of the viral non-structural protein 1 (NS1). Unbiased mass spectrometry-based analyses indicate that NS1 causes APA by perturbing the function of CPSF4 and that this function is unrelated to virus-induced transcriptional shutoff. Accordingly, IAV strain PR8, expressing an NS1 variant with weak CPSF binding, does not induce host shutoff but only APA. However, recombinant IAV (PR8) expressing NS1(G184R) lacks binding to CPSF4 and thereby also the ability to cause APA. Functionally, the impaired ability to induce APA leads to an increased inflammatory cytokine production and an attenuated phenotype in a mouse infection model. Investigating diverse viral infection models showed that APA induction is a frequent ability of many pathogens. Collectively, we propose that targeting of the CPSF complex, leading to widespread alternative polyadenylation of host transcripts, constitutes a general immunevasion mechanism employed by a variety of pathogenic viruses.

Investigating Open Reading Frames in Known and Novel Transcripts using ORFanage

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

LP30K protein manifested in hemocytes of Bombyx mori larva on Nosema bombycis infection and showed functional evolution based on glucose- binding domain

Infection by microsporidian Nosema bombycis induced appearance of exclusive protein conjugate of 190 kDa in hemocytes of silkworm Bombyx mori L (Lepidoptera: Bombycidae). Mass spectrometry of the band showed peptides of low molecular weight 30 kDa lipoprotein (LP30K). Six accessions of LP30K identified from the hemocytes comprised 30 K lipoprotein 1, 30 K protein 1, 2, 6, 7 and 11. Two uncharacterised proteins (UCP) identified from the hemocytes showed 100% similarity with LP30K sequence, altogether showed abundance after the infection. The LP30K accessions H9J4F6 (Q00802), E5EVW2 and the UCP accessions D4QGC0 and D4QGB9 showed presence of glucose binding protein I domain "ADSDVPNDILEEQLYNSIVVADYDSAVEK" that binds with fungal glucans to inhibit infection. However glucose binding protein II domain "TLAPRTDDVLAEQLYMSVVIGEYETAIAK" is absent in LP30K accessions from hemocytes showed loss of DNA sequences encoding the domain. The accessions H9J4F5, H9B440, A7LIK7 and H9B444 showed 92% identity with B. mori LP30K protein (NP_001095198.2) however the glucose binding domain I is absent in these accessions suggesting isoform- specific restricted fungal defense activity. Phylogeny tree of the LP30K homologues showed four groups including microvitellogenin and 30 kDa proteins showing functional diversity endorsed with evolutionary diversity. LP30K accessions with glucose binding domain diverged from that without glucose binding domain exemplify co-evolution for domain- dependent functional roles like storage and immune reactions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03685-x.

Transcriptomic analysis of paired healthy human skeletal muscles to identify modulators of disease severity in DMD

Muscle damage and fibro-fatty replacement of skeletal muscles is a main pathologic feature of Duchenne muscular dystrophy (DMD) with more proximal muscles affected earlier and more distal affected later in the disease course, suggesting that different skeletal muscle groups possess distinctive characteristics that influence their susceptibility to disease. To explore transcriptomic factors driving differential gene expression and modulating DMD skeletal muscle severity, we characterized the transcriptome of vastus lateralis (VL), a more proximal and susceptible muscle, relative to tibialis anterior (TA), a more distal and protected muscle, in 15 healthy individuals using bulk RNA sequencing to identify gene expression differences that may mediate their relative susceptibility to damage with loss of dystrophin. Matching single nuclei RNA sequencing data was generated for 3 of the healthy individuals, to infer cell composition in the bulk RNA sequencing dataset and to improve mapping of differentially expressed genes to their cell source of expression. A total of 3,410 differentially expressed genes were identified and mapped to cell type using single nuclei RNA sequencing of muscle, including long non-coding RNAs and protein coding genes. There was an enrichment of genes involved in calcium release from the sarcoplasmic reticulum, particularly in the myofibers and these myofiber genes were higher in the VL. There was an enrichment of genes in "Collagen-Containing Extracellular Matrix" expressed by fibroblasts, endothelial, smooth muscle and pericytes, with most genes higher in the TA, as well as genes in "Regulation Of Apoptotic Process" expressed across all cell types. Previously reported genetic modifiers were also enriched within the differentially expressed genes. We also identify 6 genes with differential isoform usage between the VL and TA. Lastly, we integrate our findings with DMD RNA sequencing data from the TA, and identify "Collagen-Containing Extracellular Matrix" and "Negative Regulation Of Apoptotic Process" as differentially expressed between DMD compared to healthy. Collectively, these findings propose novel candidate mechanisms that may mediate differential muscle susceptibility in muscular dystrophies and provide new insight into potential therapeutic targets.

Detecting differential transcript usage in complex diseases with SPIT

Differential transcript usage (DTU) plays a crucial role in determining how gene expression differs among cells, tissues, and different developmental stages, thereby contributing to the complexity and diversity of biological systems. In abnormal cells, it can also lead to deficiencies in protein function, potentially leading to pathogenesis of diseases. Detecting such events for single-gene genetic traits is relatively uncomplicated; however, the heterogeneity of populations with complex diseases presents an intricate challenge due to the presence of diverse causal events and undetermined subtypes. SPIT is the first statistical tool that quantifies the heterogeneity in transcript usage within a population and identifies predominant subgroups along with their distinctive sets of DTU events. We provide comprehensive assessments of SPIT's methodology in both single-gene and complex traits and report the results of applying SPIT to analyze brain samples from individuals with schizophrenia. Our analysis reveals previously unreported DTU events in six candidate genes.

Quantitative Analysis of Isoform Switching in Cancer

Over the past 8 years, multiple studies examined the phenomenon of isoform switching in human cancers and discovered that isoform switching is widespread, with hundreds to thousands of such events per cancer type. Although all of these studies used slightly different definitions of isoform switching, which in part led to a rather poor overlap of their results, they all leveraged transcript usage, a proportion of the transcript's expression in the total expression level of the parent gene, to detect isoform switching. However, how changes in transcript usage correlate with changes in transcript expression is not sufficiently explored. In this article, we adopt the most common definition of isoform switching and use a state-of-the-art tool for the analysis of differential transcript usage, SatuRn, to detect isoform switching events in 12 cancer types. We analyze the detected events in terms of changes in transcript usage and the relationship between transcript usage and transcript expression on a global scale. The results of our analysis suggest that the relationship between changes in transcript usage and changes in transcript expression is far from straightforward, and that such quantitative information can be effectively used for prioritizing isoform switching events for downstream analyses.

Alterations in oligodendrocyte transcriptional networks reveal region-specific vulnerabilities to neurological disease

Neurological disease is characterized the by dysfunction of specific neuroanatomical regions. To determine whether region-specific vulnerabilities have a transcriptional basis at cell-type-specific resolution, we analyzed gene expression in mouse oligodendrocytes across various brain regions. Oligodendrocyte transcriptomes cluster in an anatomical arrangement along the rostrocaudal axis. Moreover, regional oligodendrocyte populations preferentially regulate genes implicated in diseases that target their region of origin. Systems-level analyses identify five region-specific co-expression networks representing distinct molecular pathways in oligodendrocytes. The cortical network exhibits alterations in mouse models of intellectual disability and epilepsy, the cerebellar network in ataxia, and the spinal network in multiple sclerosis. Bioinformatic analyses reveal potential molecular regulators of these networks, which were confirmed to modulate network expression in vitro in human oligodendroglioma cells, including reversal of the disease-associated transcriptional effects of a pathogenic Spinocerebellar ataxia type 1 allele. These findings identify targetable region-specific vulnerabilities to neurological disease mediated by oligodendrocytes.

ACTOR: a latent Dirichlet model to compare expressed isoform proportions to a reference panel

The relative proportion of RNA isoforms expressed for a given gene has been associated with disease states in cancer, retinal diseases, and neurological disorders. Examination of relative isoform proportions can help determine biological mechanisms, but such analyses often require a per-gene investigation of splicing patterns. Leveraging large public data sets produced by genomic consortia as a reference, one can compare splicing patterns in a data set of interest with those of a reference panel in which samples are divided into distinct groups, such as tissue of origin, or disease status. We propose A latent Dirichlet model to Compare expressed isoform proportions TO a Reference panel (ACTOR), a latent Dirichlet model with Dirichlet Multinomial observations to compare expressed isoform proportions in a data set to an independent reference panel. We use a variational Bayes procedure to estimate posterior distributions for the group membership of one or more samples. Using the Genotype-Tissue Expression project as a reference data set, we evaluate ACTOR on simulated and real RNA-seq data sets to determine tissue-type classifications of genes. ACTOR is publicly available as an R package at https://github.com/mccabes292/actor.

Investigating Open Reading Frames in Known and Novel Transcripts using ORFanage

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

Transcriptome analyses in infertile men reveal germ cell-specific expression and splicing patterns

The process of spermatogenesis-when germ cells differentiate into sperm-is tightly regulated, and misregulation in gene expression is likely to be involved in the physiopathology of male infertility. The testis is one of the most transcriptionally rich tissues; nevertheless, the specific gene expression changes occurring during spermatogenesis are not fully understood. To better understand gene expression during spermatogenesis, we generated germ cell-specific whole transcriptome profiles by systematically comparing testicular transcriptomes from tissues in which spermatogenesis is arrested at successive steps of germ cell differentiation. In these comparisons, we found thousands of differentially expressed genes between successive germ cell types of infertility patients. We demonstrate our analyses' potential to identify novel highly germ cell-specific markers (TSPY4 and LUZP4 for spermatogonia; HMGB4 for round spermatids) and identified putatively misregulated genes in male infertility (RWDD2A, CCDC183, CNNM1, SERF1B). Apart from these, we found thousands of genes showing germ cell-specific isoforms (including SOX15, SPATA4, SYCP3, MKI67). Our approach and dataset can help elucidate genetic and transcriptional causes for male infertility.

The landscape of differential splicing and transcript alternations in severe COVID-19 infection

Viral infections can modulate the widespread alternations of cellular splicing, favouring viral replication within the host cells by overcoming host immune responses. However, how SARS-CoV-2 induces host cell differential splicing and affects the landscape of transcript alternation in severe COVID-19 infection remains elusive. Understanding the differential splicing and transcript alternations in severe COVID-19 infection may improve our molecular insights into the SARS-CoV-2 pathogenesis. In this study, we analysed the publicly available blood and lung transcriptome data of severe COVID-19 patients, blood transcriptome data of recovered COVID-19 patients at 12-, 16- and 24-week postinfection and healthy controls. We identified a significant transcript isoform switching in the individual blood and lung RNA-seq data of severe COVID-19-infected patients and 25 common genes that alter their transcript isoform in both blood and lung samples. Altered transcripts show significant loss of the open reading frame, functional domains and switch from coding to noncoding transcript, impacting normal cellular functions. Furthermore, we identified the expression of several novel recurrent chimeric transcripts in the blood samples from severe COVID-19 patients. Moreover, the analysis of the isoform switching into blood samples from recovered COVID-19 patients highlights that there is no significant isoform switching in 16- and 24-week postinfection, and the levels of expressed chimeric transcripts are reduced. This finding emphasizes that SARS-CoV-2 severe infection induces widespread splicing in the host cells, which could help the virus alter the host immune responses and facilitate the viral replication within the host and the efficient translation of viral proteins.

Identifying transcript-level differential expression in primary human immune cells

BACKGROUND

Multipotential hematopoietic stem cells differentiate into a wide variety of immune cells with a diversity of functions, including the ability to respond to a variety of stimuli. Importantly, numerous studies have demonstrated the importance of gene transcription in defining cell identity and functions. While these studies have primarily been performed at the level of the gene, it is known that key immune genes such as CD44 and CD45 generate multiple different transcripts that are differentially expressed across different immune cells, and that encode proteins with different sequences and functions. Prior genomic surveys have shown that the mechanisms for generating diversity in expressed transcripts (alternate splicing, alternate transcription start sites, etc.) are very active in immune cells, but have been lacking in terms of identifying genes with multiple transcripts, that are differentially expressed, and likely to affect cell functions.

METHODS

We first identified the set of genes that had at least two transcripts expressed in our RNA sequencing dataset generated from purified populations of neutrophils, monocytes and five lymphocyte populations (B, NK, γδ T, CD4 + T and CD8 + T) from twelve healthy donors. Next, we developed a heuristic approach to identify genes where two or more transcripts have distinct expression patterns across lymphoid and/or myeloid populations. We then focused our annotation and interpretation on differentially expressed transcripts that affect the coding sequence. This process was repeated to identify transcripts that were differentially expressed between monocytes and populations of macrophages and LPS-stimulated macrophages derived from these monocytes in vitro.

RESULTS

We found that over 55 % of genes had two or more expressed transcripts, with an average ∼3 transcripts per gene, and that 70 % of these had at least two of the transcripts that encoded proteins with different sequences. As expected, we identified a complex pattern of differential expression for multiple transcripts encoding the CD45 transmembrane protein, but we also found similar evidence for ten other genes (CD300A, FYB1, GPI, LITAF, PSMA1, PTMA, RPL32, SEPTIN9, SH3BP2, SH3KBP1) when comparing the expression patterns of transcripts within myeloid and lymphoid cells. We also identified five genes with differentially expressed transcripts associated with the transition from monocytes to macrophages (FNBP1, KLF6, and SEPTIN9) or between macrophages and LPS-stimulated macrophages (CD44, OAZ2, and SEPTIN9). For the most part, we found that the different transcripts of these genes are expected to impact specific biological functions, for example the different transcripts of SEPTIN9 likely regulate the cytoskeleton in immune cells via their interactions with actins filaments and microtubules.

CONCLUSIONS

This analytic approach successfully identified multi-transcript genes that are differentially expressed across immune cells and could be applied to other transcriptomic data.

DATA AVAILABILITY STATEMENT

Researchers can request access to the individual-level data from the current study by contacting the Montreal Heart Institute ethics committee at the following institutional email address: cer.icm@icm-mhi.org.

Systematic analysis of alternative splicing in time course data using Spycone

MOTIVATION

During disease progression or organism development, alternative splicing may lead to isoform switches that demonstrate similar temporal patterns and reflect the alternative splicing co-regulation of such genes. Tools for dynamic process analysis usually neglect alternative splicing.

RESULTS

Here, we propose Spycone, a splicing-aware framework for time course data analysis. Spycone exploits a novel IS detection algorithm and offers downstream analysis such as network and gene set enrichment. We demonstrate the performance of Spycone using simulated and real-world data of SARS-CoV-2 infection.

AVAILABILITY AND IMPLEMENTATION

The Spycone package is available as a PyPI package. The source code of Spycone is available under the GPLv3 license at https://github.com/yollct/spycone and the documentation at https://spycone.readthedocs.io/en/latest/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Expression Profile of Sphingosine Kinase 1 Isoforms in Human Cancer Tissues and Cells: Importance and Clinical Relevance of the Neglected 1b-Isoform

Background

Overexpression of sphingosine kinase 1 (SphK1) is casually associated with many types of cancer, and inhibitors of SphK1 sensitize tumors to chemotherapy. SphK1 is expressed as two major isoforms, SphK1a and SphK1b. To date, no information has been reported on the SphK1 isoform expression profile and its clinical relevance.

Objective

The objective is to examine the expression profile of the SphK1a and SPhK1b isoforms in human cancer and noncancer tissues and cell lines and explore their clinical relevance.

Methods

We used PCR to qualitatively examine the expression profile of these two isoforms in breast, liver, and prostate cancer tissues plus paired adjacent tissues and in 11 cancer and normal cell lines (breast, cervical, bone, prostate, colon, brain, mesothelioma tumor and benign, and human kidney cells).

Results

We found that SphK1a was ubiquitously expressed in all cancer cells and tissues tested; in contrast, SphK1b was only expressed in selective cell types in breast, prostate, and lung cancer.

Conclusions

Our data suggest that SphK1a is important for generic SphK1/S1P functions, and SphK1b mediates specialized and/or unique pathways in a specific type of tissue and could be a biomarker for cancer. This discovery is important for future SphK1-related cancer research and may have clinical implications in drug development associated with SphK1-directed cancer treatment.

Possible Contribution of Alternative Transcript Isoforms in Mature Biofilm Growth Phase of Candida glabrata

Expression of genome-wide alternative transcript isoforms and differential transcript isoform usage in different biological conditions (isoform switching) are responsible for the varied proteomic functional diversity in higher eukaryotic organisms. However, these mechanisms have not been studied in Candida glabrata, which is a potent eukaryotic opportunistic pathogen. Biofilm formation is an important virulence factor of C. glabrata that withstands antifungal drug stress and overcomes the host-immune response. Here, we present the genome-wide differential transcript isoform expression (DTE) and differential transcript isoform usage (DTU) in a mature biofilm growth phase of C. glabrata (clinical isolate; NCCPF 100,037) using the RNA sequencing approach. The DTE analysis generated 7837 transcript isoforms from the C. glabrata genome (5293 genes in total), and revealed that transcript isoforms generated from 292 genes showed significant DTU in the mature biofilm cells. Gene ontology, pathway analysis and protein-protein interactions of significant transcript isoforms, further substantiated that their specific expression and differential usage is required for transitioning the planktonic cells to biofilm in C. glabrata. The present study reported the possible role of expression of alternative transcript isoforms and differential transcript isoform usage in the mature biofilms of C. glabrata. The observation derived from the study may prove to be beneficial for making future antifungal therapeutic strategies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12088-022-01036-7.

Investigating isoform switching in RHBDF2 and its role in neoplastic growth in breast cancer

Background

Breast cancer is the second leading cause of cancer-related deaths globally, and its prevalence rates are increasing daily. In the past, studies predicting therapeutic drug targets for cancer therapy focused on the assumption that one gene is responsible for producing one protein. Therefore, there is always an immense need to find promising and novel anti-cancer drug targets. Furthermore, proteases have an integral role in cell proliferation and growth because the proteolysis mechanism is an irreversible process that aids in regulating cellular growth during tumorigenesis. Therefore, an inactive rhomboid protease known as iRhom2 encoded by the gene RHBDF2 can be considered an important target for cancer treatment. Speculatively, previous studies on gene expression analysis of RHBDF2 showed heterogenous behaviour during tumorigenesis. Consistent with this, several studies have reported the antagonistic role of iRhom2 in tumorigenesis, i.e., either they are involved in negative regulation of EGFR ligands via the ERAD pathway or positively regulate EGFR ligands via the EGFR signalling pathway. Additionally, different opinions suggest iRhom2 mediated cleavage of EGFR ligands takes place TACE dependently or TACE independently. However, reconciling these seemingly opposing roles is still unclear and might be attributed to more than one transcript isoform of iRhom2.

Methods

To observe the differences at isoform resolution, the current strategy identified isoform switching in RHBDF2 via differential transcript usage using RNA-seq data during breast cancer initiation and progression. Furthermore, interacting partners were found via correlation and enriched to explain their antagonistic role.

Results

Isoform switching was observed at DCIS, grade 2 and grade 3, from canonical to the cub isoform. Neither EGFR nor ERAD was found enriched. However, pathways leading to TACE-dependent EGFR signalling pathways were more observant, specifically MAPK signalling pathways, GPCR signalling pathways, and toll-like receptor pathways. Nevertheless, it was noteworthy that during CTCs, the cub isoform switches back to the canonical isoform, and the proteasomal degradation pathway and cytoplasmic ribosomal protein pathways were significantly enriched. Therefore, it could be inferred that cub isoform functions during cancer initiation in EGFR signalling. In contrast, during metastasis, where invasion is the primary task, the isoform switches back to the canonical isoform.

Multi-omics approach to identifying isoform variants as therapeutic targets in cancer patients

Cancer-specific alternatively spliced events (ASE) play a role in cancer pathogenesis and can be targeted by immunotherapy, oligonucleotide therapy, and small molecule inhibition. However, identifying actionable ASE targets remains challenging due to the uncertainty of its protein product, structure impact, and proteoform (protein isoform) function. Here we argue that an integrated multi-omics profiling strategy can overcome these challenges, allowing us to mine this untapped source of targets for therapeutic development. In this review, we will provide an overview of current multi-omics strategies in characterizing ASEs by utilizing the transcriptome, proteome, and state-of-art algorithms for protein structure prediction. We will discuss limitations and knowledge gaps associated with each technology and informatics analytics. Finally, we will discuss future directions that will enable the full integration of multi-omics data for ASE target discovery.

Comprehensive Transcriptome Analysis Reveals Genome-Wide Changes Associated with Endoplasmic Reticulum (ER) Stress in Potato (Solanum tuberosum L.)

We treated potato (Solanum tuberosum L.) plantlets with TM and performed gene expression studies to identify genome-wide changes associated with endoplasmic reticulum (ER) stress and the unfolded protein response (UPR). An extensive network of responses was identified, including chromatin remodeling, transcriptional reprogramming, as well as changes in the structural components of the endomembrane network system. Limited genome-wide changes in alternative RNA splicing patterns of protein-coding transcripts were also discovered. Significant changes in RNA metabolism, components of the translation machinery, as well as factors involved in protein folding and maturation occurred, which included a broader set of genes than expected based on Arabidopsis research. Antioxidant defenses and oxygen metabolic enzymes are differentially regulated, which is expected of cells that may be experiencing oxidative stress or adapting to protect proteins from oxidation. Surges in protein kinase expression indicated early signal transduction events. This study shows early genomic responses including an array of differentially expressed genes that have not been reported in Arabidopsis. These data describe novel ER stress responses in a solanaceous host.

Analysis of matched primary and recurrent BRCA1/2 mutation-associated tumors identifies recurrence-specific drivers

Recurrence is a major cause of death among BRCA1/2 mutation carriers with breast (BrCa) and ovarian cancers (OvCa). Herein we perform multi-omic sequencing on 67 paired primary and recurrent BrCa and OvCa from 27 BRCA1/2 mutation carriers to identify potential recurrence-specific drivers. PARP1 amplifications are identified in recurrences (False Discovery Rate q = 0.05), and PARP1 is significantly overexpressed across primary BrCa and recurrent BrCa and OvCa, independent of amplification status. RNA sequencing analysis finds two BRCA2 isoforms, BRCA2-201/Long and BRCA2-001/Short, respectively predicted to be sensitive and insensitive to nonsense-mediated decay. BRCA2-001/Short is expressed more frequently in recurrences and associated with reduced overall survival in breast cancer (87 vs. 121 months; Hazard Ratio = 2.5 [1.18-5.5]). Loss of heterozygosity (LOH) status is discordant in 25% of patient's primary and recurrent tumors, with switching between both LOH and lack of LOH found. Our study reveals multiple potential drivers of recurrent disease in BRCA1/2 mutation-associated cancer, improving our understanding of tumor evolution and suggesting potential biomarkers.

Magnetique: an interactive web application to explore transcriptome signatures of heart failure

Background

Despite a recent increase in the number of RNA-seq datasets investigating heart failure (HF), accessibility and usability remain critical issues for medical researchers. We address the need for an intuitive and interactive web application to explore the transcriptional signatures of heart failure with this work.

Methods

We reanalysed the Myocardial Applied Genomics Network RNA-seq dataset, one of the largest publicly available datasets of left ventricular RNA-seq samples from patients with dilated (DCM) or hypertrophic (HCM) cardiomyopathy, as well as unmatched non-failing hearts (NFD) from organ donors and patient characteristics that allowed us to model confounding factors. We analyse differential gene expression, associated pathway signatures and reconstruct signaling networks based on inferred transcription factor activities through integer linear programming. We additionally focus, for the first time, on differential RNA transcript isoform usage (DTU) changes and predict RNA-binding protein (RBP) to target transcript interactions using a Global test approach. We report results for all pairwise comparisons (DCM, HCM, NFD).

Results

Focusing on the DCM versus HCM contrast (DCMvsHCM), we identified 201 differentially expressed genes, some of which can be clearly associated with changes in ERK1 and ERK2 signaling. Interestingly, the signs of the predicted activity for these two kinases have been inferred to be opposite to each other: In the DCMvsHCM contrast, we predict ERK1 to be consistently less activated in DCM while ERK2 was more activated in DCM. In the DCMvsHCM contrast, we identified 149 differently used transcripts. One of the top candidates is the O-linked N-acetylglucosamine (GlcNAc) transferase (OGT), which catalyzes a common post-translational modification known for its role in heart arrhythmias and heart hypertrophy. Moreover, we reconstruct RBP – target interaction networks and showcase the examples of CPEB1, which is differentially expressed in the DCMvsHCM contrast.

Conclusion

Magnetique (https://shiny.dieterichlab.org/app/magnetique) is the first online application to provide an interactive view of the HF transcriptome at the RNA isoform level and to include transcription factor signaling and RBP:RNA interaction networks. The source code for both the analyses (https://github.com/dieterich-lab/magnetiqueCode2022) and the web application (https://github.com/AnnekathrinSilvia/magnetique) is available to the public. We hope that our application will help users to uncover the molecular basis of heart failure.

Progression of the pluripotent epiblast depends upon the NMD factor UPF2

Nonsense-mediated RNA decay (NMD) is a highly conserved RNA turnover pathway that degrades RNAs harboring in-frame stop codons in specific contexts. Loss of NMD factors leads to embryonic lethality in organisms spanning the phylogenetic scale, but the mechanism remains unknown. Here, we report that the core NMD factor, UPF2, is required for expansion of epiblast cells within the inner cell mass of mice in vivo. We identify NMD target mRNAs in mouse blastocysts - both canonical and alternatively processed mRNAs - including those encoding cell cycle arrest and apoptosis factors, raising the possibility that NMD is essential for embryonic cell proliferation and survival. In support, the inner cell mass of Upf2-null blastocysts rapidly regresses with outgrowth and is incompetent for embryonic stem cell derivation in vitro. In addition, we uncovered concordant temporal- and lineage-specific regulation of NMD factors and mRNA targets, indicative of a shift in NMD magnitude during peri-implantation development. Together, our results reveal developmental and molecular functions of the NMD pathway in the early embryo.

Regulatory dynamics of gene expression in the developing male gametophyte of Arabidopsis

Sexual reproduction in angiosperms requires the production and delivery of two male gametes by a three-celled haploid male gametophyte. This demands synchronized gene expression in a short developmental window to ensure double fertilization and seed set. While transcriptomic changes in developing pollen are known for Arabidopsis, no studies have integrated RNA and proteomic data in this model. Further, the role of alternative splicing has not been fully addressed, yet post-transcriptional and post-translational regulation may have a key role in gene expression dynamics during microgametogenesis. We have refined and substantially updated global transcriptomic and proteomic changes in developing pollen for two Arabidopsis accessions. Despite the superiority of RNA-seq over microarray-based platforms, we demonstrate high reproducibility and comparability. We identify thousands of long non-coding RNAs as potential regulators of pollen development, hundreds of changes in alternative splicing and provide insight into mRNA translation rate and storage in developing pollen. Our analysis delivers an integrated perspective of gene expression dynamics in developing Arabidopsis pollen and a foundation for studying the role of alternative splicing in this model.

The Glycosyltransferase Pathway: An Integrated Analysis of the Cell Metabolome

Nucleotide sugar-dependent glycosyltransferases (UGTs) are critical to the homeostasis of endogenous metabolites and the detoxification of xenobiotics. Their impact on the cell metabolome remains unknown. Cellular metabolic changes resulting from human UGT expression were profiled by untargeted metabolomics. The abundant UGT1A1 and UGT2B7 were studied as UGT prototypes along with their alternative (alt.) splicing-derived isoforms displaying structural differences. Nineteen biochemical routes were modified, beyond known UGT substrates. Significant variations in glycolysis and pyrimidine pathways, and precursors of the co-substrate UDP-glucuronic acid were observed. Bioactive lipids such as arachidonic acid and endocannabinoids were highly enriched by up to 13.3-fold (p < 0.01) in cells expressing the canonical enzymes. Alt. UGT2B7 induced drastic and unique metabolic perturbations, including higher glucose (18-fold) levels and tricarboxylic acid cycle (TCA) cycle metabolites and abrogated the effects of the UGT2B7 canonical enzyme when co-expressed. UGT1A1 proteins promoted the accumulation of branched-chain amino acids (BCAA) and TCA metabolites upstream of the mitochondrial oxoglutarate dehydrogenase complex (OGDC). Alt. UGT1A1 exacerbated these changes, likely through its interaction with the OGDC component oxoglutarate dehydrogenase-like (OGDHL). This study expands the breadth of biochemical pathways associated with UGT expression and establishes extensive connectivity between UGT enzymes, alt. proteins and other metabolic processes.

Progeria and Aging-Omics Based Comparative Analysis

Since ancient times aging has also been regarded as a disease, and humankind has always strived to extend the natural lifespan. Analyzing the genes involved in aging and disease allows for finding important indicators and biological markers for pathologies and possible therapeutic targets. An example of the use of omics technologies is the research regarding aging and the rare and fatal premature aging syndrome progeria (Hutchinson-Gilford progeria syndrome, HGPS). In our study, we focused on the in silico analysis of differentially expressed genes (DEGs) in progeria and aging, using a publicly available RNA-Seq dataset (GEO dataset GSE113957) and a variety of bioinformatics tools. Despite the GSE113957 RNA-Seq dataset being well-known and frequently analyzed, the RNA-Seq data shared by Fleischer et al. is far from exhausted and reusing and repurposing the data still reveals new insights. By analyzing the literature citing the use of the dataset and subsequently conducting a comparative analysis comparing the RNA-Seq data analyses of different subsets of the dataset (healthy children, nonagenarians and progeria patients), we identified several genes involved in both natural aging and progeria (KRT8, KRT18, ACKR4, CCL2, UCP2, ADAMTS15, ACTN4P1, WNT16, IGFBP2). Further analyzing these genes and the pathways involved indicated their possible roles in aging, suggesting the need for further in vitro and in vivo research. In this paper, we (1) compare "normal aging" (nonagenarians vs. healthy children) and progeria (HGPS patients vs. healthy children), (2) enlist genes possibly involved in both the natural aging process and progeria, including the first mention of IGFBP2 in progeria, (3) predict miRNAs and interactomes for WNT16 (hsa-mir-181a-5p), UCP2 (hsa-mir-26a-5p and hsa-mir-124-3p), and IGFBP2 (hsa-mir-124-3p, hsa-mir-126-3p, and hsa-mir-27b-3p), (4) demonstrate the compatibility of well-established R packages for RNA-Seq analysis for researchers interested but not yet familiar with this kind of analysis, and (5) present comparative proteomics analyses to show an association between our RNA-Seq data analyses and corresponding changes in protein expression.