The Landscape of Isoform Switches in Human Cancers

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.

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.

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.

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.

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.

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.

Characterizing isoform switching events in esophageal adenocarcinoma

Isoform switching events with predicted functional consequences are common in many cancers, but characterization of switching events in esophageal adenocarcinoma (EAC) is lacking. Next-generation sequencing was used to detect levels of RNA transcripts and identify specific isoforms in treatment-naïve esophageal tissues ranging from premalignant Barrett’s esophagus (BE), BE with low- or high-grade dysplasia (BE.LGD, BE.HGD), and EAC. Samples were stratified by histopathology and TP53 mutation status, identifying significant isoform switching events with predicted functional consequences. Comparing BE.LGD with BE.HGD, a histopathology linked to cancer progression, isoform switching events were identified in 75 genes including KRAS, RNF128, and WRAP53. Stratification based on TP53 status increased the number of significant isoform switches to 135, suggesting switching events affect cellular functions based on TP53 mutation and tissue histopathology. Analysis of isoforms agnostic, exclusive, and shared with mutant TP53 revealed unique signatures including demethylation, lipid and retinoic acid metabolism, and glucuronidation, respectively. Nearly half of isoform switching events were identified without significant gene-level expression changes. Importantly, two TP53-interacting isoforms, RNF128 and WRAP53, were significantly linked to patient survival. Thus, analysis of isoform switching events may provide new insight for the identification of prognostic markers and inform new potential therapeutic targets for EAC.

The landscape of isoform switches in sepsis: a multicenter cohort study

Sepsis is caused by an uncontrolled inflammatory response, whose underlying mechanisms are not fully understood. It is well known that the majority of human genes can be expressed as alternative isoforms. While isoform switching is implicated in many diseases and is particularly prominent in cancer, it has never been reported in the context of sepsis. Patients presented to the emergency department of three tertiary care hospitals from January 2020 to December 2020 were enrolled. Clinical variables and genome-wide transcriptome of peripheral blood mononuclear cells (PBMC) were obtained. Isoform switching analysis were performed to identify significant isoform switches and relevant biological consequences. A total of 48 subjects with sepsis, involving 42 survivors and 6 non-survivors, admitted to the emergency department of three tertiary care hospitals were enrolled in this study. PBMCs were extracted for RNA sequencing (RNA-seq). Patients (n = 4) with mild stroke or acute coronary syndrome without infection were enrolled in this study as controls. The most frequent functional changes resulting from isoform switching were changes affecting the open reading frame, protein domains and intron retention. Many genes without differences in gene expression showed significant isoform switching. Many genes with significant isoform switches ([Formula: see text]> 0.1) were associated with higher mortality risk, including PIGS, CASP3, LITAF, HBB and RUVBL2. The study for the first time described the landscape of isoform switching in sepsis, including differentially expressed isoform fractions between patients with and without sepsis and survivors and nonsurvivors. The biological consequences of isoform switching, including protein domain loss, signal peptide gain, and intron retention, were identified.

In-Depth Temporal Transcriptome Profiling of an Alphaherpesvirus Using Nanopore Sequencing

In this work, a long-read sequencing (LRS) technique based on the Oxford Nanopore Technology MinION platform was used for quantifying and kinetic characterization of the poly(A) fraction of bovine alphaherpesvirus type 1 (BoHV-1) lytic transcriptome across a 12-h infection period. Amplification-based LRS techniques frequently generate artefactual transcription reads and are biased towards the production of shorter amplicons. To avoid these undesired effects, we applied direct cDNA sequencing, an amplification-free technique. Here, we show that a single promoter can produce multiple transcription start sites whose distribution patterns differ among the viral genes but are similar in the same gene at different timepoints. Our investigations revealed that the circ gene is expressed with immediate–early (IE) kinetics by utilizing a special mechanism based on the use of the promoter of another IE gene (bicp4) for the transcriptional control. Furthermore, we detected an overlap between the initiation of DNA replication and the transcription from the bicp22 gene, which suggests an interaction between the two molecular machineries. This study developed a generally applicable LRS-based method for the time-course characterization of transcriptomes of any organism.

Exon junction complex-associated multi-adapter RNPS1 nucleates splicing regulatory complexes to maintain transcriptome surveillance

The exon junction complex (EJC) is an RNA-binding multi-protein complex with critical functions in post-transcriptional gene regulation. It is deposited on the mRNA during splicing and regulates diverse processes including pre-mRNA splicing and nonsense-mediated mRNA decay (NMD) via various interacting proteins. The peripheral EJC-binding protein RNPS1 was reported to serve two insufficiently characterized functions: suppressing mis-splicing of cryptic splice sites and activating NMD in the cytoplasm. The analysis of transcriptome-wide effects of EJC and RNPS1 knockdowns in different human cell lines supports the conclusion that RNPS1 can moderately influence NMD activity, but is not a globally essential NMD factor. However, numerous aberrant splicing events strongly suggest that the main function of RNPS1 is splicing regulation. Rescue analyses revealed that the RRM and C-terminal domain of RNPS1 both contribute partially to regulate RNPS1-dependent splicing events. We defined the RNPS1 core interactome using complementary immunoprecipitations and proximity labeling, which identified interactions with splicing-regulatory factors that are dependent on the C-terminus or the RRM domain of RNPS1. Thus, RNPS1 emerges as a multifunctional splicing regulator that promotes correct and efficient splicing of different vulnerable splicing events via the formation of diverse splicing-promoting complexes.

The HLA-G immune checkpoint: a new immuno-stimulatory role for the α1-domain-deleted isoform

The heterogeneity of cancer cells, in part maintained via the expression of multiple isoforms, introduces significant challenges in designing effective therapeutic approaches. In this regard, isoforms of the immune checkpoint HLA-G have been found in most of the tumors analyzed, such as ccRCC, the most common human renal malignancy. In particular, HLA-G∆α1, which is the only HLA-G isoform described that lacks the α1 extracellular domain, has been newly identified in ccRCC and now here in trophoblasts. Using a cellular model expressing HLA-G∆α1, we have uncovered its specific and overlapping functional roles, relative to the main HLA-G isoform, i.e., the full-length HLA-G1. We found that HLA-G∆α1 has several particular features: (i) although possessing the α3 domain, it does not associate with β2-microglobulin; (ii) it may not present peptides to T cells due to absence of the peptide-binding groove; and (iii) it exerts immune-stimulatory activity towards peripheral blood NK and T cells, while all known isoforms of HLA-G are immune-inhibitory checkpoint molecules. Such immune-stimulatory properties of HLA-G∆α1 on the cytotoxic function of peripheral blood NK cells are individual dependent and are not exerted through the interaction with the known HLA-G receptor, ILT2. Importantly, we are faced here with a potential antitumor effect of an HLA-G isoform, opposed to the pro-tumor properties described for all other HLA-G isoforms, which should be taken into account in future therapeutic designs aimed at blocking this immune checkpoint.

ERBB3 methylation and immune infiltration in tumor microenvironment of cervical cancer

ERBB3, a member of the ERBB family of receptor tyrosine kinases, plays an important role in cancer, despite its lack of intrinsic carcinogenic mechanism of cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC). Research on bioinformatics methods through multi-omics, this work proves that ERBB3 gene mutation, methylation modification have extensive regulatory mechanisms on the CESC microenvironment. We found that ERBB3 is involved in carcinogenesis of cervical cancer and is not associated with its prognosis. The carcinogenic mechanism is mainly related to the suppression of the immune system between tumor infiltrating lymphocytes (TILs) and the methylation of the RNA level. Our study indicated ERBB3 is more likely to be a carcinogenic factor than a key prognostic factor for cervical cancer. Methylation of ERBB3 may work as a checkpoint immunotherapy target in CESC, DNA methylation modification of the 4480 base pair downstream of ERBB3 transcription initiation site was the highest.

RETRACTED: Parallel bimodal single-cell sequencing of transcriptome and methylome provides molecular and translational insights on oocyte maturation and maternal aging

This article has been retracted: please see Elsevier Policy on Article Withdrawal (http://www.elsevier.com/locate/withdrawalpolicy). This article has been retracted at the request of the Editor-in-Chief. It has been brought to our attention that the authors of the article "Parallel bimodal single-cell sequencing of transcriptome and methylome provides molecular and translational insights on oocyte maturation and maternal aging" cannot agree on who should be listed as an author of the article. Further inquiry by the journal revealed that the authorship was also changed at the revision stages of the article without notifying the handling Editor, which is contrary to the journal policy on changes to authorship. The journal considers this unacceptable practice, and the Editor-in-Chief decided to retract the article.

Human UPF3A and UPF3B enable fault-tolerant activation of nonsense-mediated mRNA decay

The paralogous human proteins UPF3A and UPF3B are involved in recognizing mRNAs targeted by nonsense‐mediated mRNA decay (NMD). UPF3B has been demonstrated to support NMD, presumably by bridging an exon junction complex (EJC) to the NMD factor UPF2. The role of UPF3A has been described either as a weak NMD activator or an NMD inhibitor. Here, we present a comprehensive functional analysis of UPF3A and UPF3B in human cells using combinatory experimental approaches. Overexpression or knockout of UPF3A as well as knockout of UPF3B did not substantially change global NMD activity. In contrast, the co‐depletion of UPF3A and UPF3B resulted in a marked NMD inhibition and a transcriptome‐wide upregulation of NMD substrates, demonstrating a functional redundancy between both NMD factors. In rescue experiments, UPF2 or EJC binding‐deficient UPF3B largely retained NMD activity. However, combinations of different mutants, including deletion of the middle domain, showed additive or synergistic effects and therefore failed to maintain NMD. Collectively, UPF3A and UPF3B emerge as fault‐tolerant, functionally redundant NMD activators in human cells.

Target isoforms are an overlooked challenge and opportunity in chimeric antigen receptor cell therapy

The development of novel chimeric antigen receptor (CAR) cell therapies is rapidly growing, with 299 new agents being reported and 109 new clinical trials initiated so far this year. One critical lesson from approved CD19-specific CAR therapies is that target isoform switching has been shown to cause tumour relapse, but little is known about the isoforms of CAR targets in solid cancers. Here we assess the protein isoform landscape and identify both the challenges and opportunities protein isoform switching present as CAR therapy is applied to solid cancers.

Placental Transcription Profiling in 6-23 Weeks' Gestation Reveals Differential Transcript Usage in Early Development

The human placenta is a rapidly developing transient organ that is key to pregnancy success. Early development of the conceptus occurs in a low oxygen environment before oxygenated maternal blood begins to flow into the placenta at ~10-12 weeks' gestation. This process is likely to substantially affect overall placental gene expression. Transcript variability underlying gene expression has yet to be profiled. In this study, accurate transcript expression profiles were identified for 84 human placental chorionic villus tissue samples collected across 6-23 weeks' gestation. Differential gene expression (DGE), differential transcript expression (DTE) and differential transcript usage (DTU) between 6-10 weeks' and 11-23 weeks' gestation groups were assessed. In total, 229 genes had significant DTE yet no significant DGE. Integration of DGE and DTE analyses found that differential expression patterns of individual transcripts were commonly masked upon aggregation to the gene-level. Of the 611 genes that exhibited DTU, 534 had no significant DGE or DTE. The four most significant DTU genes ADAM10, VMP1, GPR126, and ASAH1, were associated with hypoxia-responsive pathways. Transcript usage is a likely regulatory mechanism in early placentation. Identification of functional roles will facilitate new insight in understanding the origins of pregnancy complications.