Macrophage development and activation involve coordinated intron retention in key inflammatory regulators

Atherosis-associated lnc_000048 activates PKR to enhance STAT1-mediated polarization of THP-1 macrophages to M1 phenotype

JOURNAL/nrgr/04.03/01300535-202419110-00029/figure1/v/2024-03-08T184507Z/r/image-tiff Our previous study has demonstrated that lnc_000048 is upregulated in large-artery atherosclerotic stroke and promotes atherosclerosis in ApoE-/- mice. However, little is known about the role of lnc_000048 in classically activated macrophage (M1) polarization. In this study, we established THP-1-derived testing state macrophages (M0), M1 macrophages, and alternately activated macrophages (M2). Real-time fluorescence quantitative PCR was used to verify the expression of marker genes and the expression of lnc_000048 in macrophages. Flow cytometry was used to detect phenotypic proteins (CD11b, CD38, CD80). We generated cell lines with lentivirus-mediated upregulation or downregulation of lnc_000048. Flow cytometry, western blot, and real-time fluorescence quantitative PCR results showed that down-regulation of lnc_000048 reduced M1 macrophage polarization and the inflammation response, while over-expression of lnc_000048 led to the opposite effect. Western blot results indicated that lnc_000048 enhanced the activation of the STAT1 pathway and mediated the M1 macrophage polarization. Moreover, catRAPID prediction, RNA-pull down, and mass spectrometry were used to identify and screen the protein kinase RNA-activated (PKR), then catRAPID and RPIseq were used to predict the binding ability of lnc_000048 to PKR. Immunofluorescence (IF)-RNA fluorescence in situ hybridization (FISH) double labeling was performed to verify the subcellular colocalization of lnc_000048 and PKR in the cytoplasm of M1 macrophage. We speculate that lnc_000048 may form stem-loop structure-specific binding and activate PKR by inducing its phosphorylation, leading to activation of STAT1 phosphorylation and thereby enhancing STAT1 pathway-mediated polarization of THP-1 macrophages to M1 and inflammatory factor expression. Taken together, these results reveal that the lnc_000048/PKR/STAT1 axis plays a crucial role in the polarization of M1 macrophages and may be a novel therapeutic target for atherosclerosis alleviation in stroke.

Reduced protein-coding transcript diversity in severe dengue emphasises the role of alternative splicing

Dengue fever, a neglected tropical arboviral disease, has emerged as a global health concern in the past decade. Necessitating a nuanced comprehension of the intricate dynamics of host-virus interactions influencing disease severity, we analysed transcriptomic patterns using bulk RNA-seq from 112 age- and gender-matched NS1 antigen-confirmed hospital-admitted dengue patients with varying severity. Severe cases exhibited reduced platelet count, increased lymphocytosis, and neutropenia, indicating a dysregulated immune response. Using bulk RNA-seq, our analysis revealed a minimal overlap between the differentially expressed gene and transcript isoform, with a distinct expression pattern across the disease severity. Severe patients showed enrichment in retained intron and nonsense-mediated decay transcript biotypes, suggesting altered splicing efficiency. Furthermore, an up-regulated programmed cell death, a haemolytic response, and an impaired interferon and antiviral response at the transcript level were observed. We also identified the potential involvement of the RBM39 gene among others in the innate immune response during dengue viral pathogenesis, warranting further investigation. These findings provide valuable insights into potential therapeutic targets, underscoring the importance of exploring transcriptomic landscapes between different disease sub-phenotypes in infectious diseases.

PRPF40A induces inclusion of exons in GC-rich regions important for human myeloid cell differentiation

We characterized the regulatory mechanisms and role in human myeloid cell survival and differentiation of PRPF40A, a splicing factor lacking a canonical RNA Binding Domain. Upon PRPF40A knockdown, HL-60 cells displayed increased cell death, decreased proliferation and slight differentiation phenotype with upregulation of immune activation genes. Suggestive of both redundant and specific functions, cell death but not proliferation was rescued by overexpression of its paralog PRPF40B. Transcriptomic analysis revealed the predominant role of PRPF40A as an activator of cassette exon inclusion of functionally relevant splicing events. Mechanistically, the exons exclusively upregulated by PRPF40A are flanked by short and GC-rich introns which tend to localize to nuclear speckles in the nucleus center. These PRPF40A regulatory features are shared with other splicing regulators such as SRRM2, SON, PCBP1/2, and to a lesser extent TRA2B and SRSF2, as a part of a functional network that regulates splicing partly via co-localization in the nucleus.

Characterised intron retention profiles in muscle tissue of idiopathic inflammatory myopathy subtypes

OBJECTIVES

Idiopathic inflammatory myopathies (IIMs) are a group of heterogeneous autoimmune diseases. Intron retention (IR) serves as an important post-transcriptional and translational regulatory mechanism. This study aims to identify changes in IR profiles in IIM subtypes, investigating their influence on proteins and their correlations with clinical features.

METHODS

RNA sequencing and liquid chromatography-tandem mass spectrometry were performed on muscle tissues obtained from 174 patients with IIM and 19 controls, following QC procedures. GTFtools and iREAD software were used for IR identification. An analysis of differentially expressed IRs (DEIs), exons and proteins was carried out using edgeR or DEP. Functional analysis was performed with clusterProfiler, and SPIRON was used to assess splicing factors.

RESULTS

A total of 6783 IRs located in 3111 unique genes were identified in all IIM subtypes compared with controls. IIM subtype-specific DEIs were associated with the pathogenesis of respective IIM subtypes. Splicing factors YBX1 and HSPA2 exhibited the most changes in dermatomyositis and immune-mediated necrotising myopathy. Increased IR was associated with reduced protein expression. Some of the IIM-specific DEIs were correlated with clinical parameters (skin rash, MMT-8 scores and muscle enzymes) and muscle histopathological features (myofiber necrosis, regeneration and inflammation). IRs in IFIH1 and TRIM21 were strongly correlated with anti-MDA5+ antibody, while IRs in SRP14 were associated with anti-SRP+ antibody.

CONCLUSION

This study revealed distinct IRs and specific splicing factors associated with IIM subtypes, which might be contributing to the pathogenesis of IIM. We also emphasised the potential impact of IR on protein expression in IIM muscles.

Dynamic changes in RNA m6A and 5 hmC influence gene expression programs during macrophage differentiation and polarisation

RNA modifications are essential for the establishment of cellular identity. Although increasing evidence indicates that RNA modifications regulate the innate immune response, their role in monocyte-to-macrophage differentiation and polarisation is unclear. While m6A has been widely studied, other RNA modifications, including 5 hmC, remain poorly characterised. We profiled m6A and 5 hmC epitranscriptomes, transcriptomes, translatomes and proteomes of monocytes and macrophages at rest and pro- and anti-inflammatory states. Transcriptome-wide mapping of m6A and 5 hmC reveals enrichment of m6A and/or 5 hmC on specific categories of transcripts essential for macrophage differentiation. Our analyses indicate that m6A and 5 hmC modifications are present in transcripts with critical functions in pro- and anti-inflammatory macrophages. Notably, we also discover the co-occurrence of m6A and 5 hmC on alternatively-spliced isoforms and/or opposing ends of the untranslated regions (UTR) of mRNAs with key roles in macrophage biology. In specific examples, RNA 5 hmC controls the decay of transcripts independently of m6A. This study provides (i) a comprehensive dataset to interrogate the role of RNA modifications in a plastic system (ii) a resource for exploring different layers of gene expression regulation in the context of human monocyte-to-macrophage differentiation and polarisation, (iii) new insights into RNA modifications as central regulators of effector cells in innate immunity.

Ribonuclease activity undermines immune sensing of naked extracellular RNA

The plasma membrane and the membrane of endosomal vesicles are considered physical barriers preventing extracellular RNA uptake. While naked RNA can be spontaneously internalized by certain cells types, functional delivery of naked RNA into the cytosol has been rarely observed. Here we show that extracellular ribonucleases, mainly derived from cell culture supplements, have so far hindered the study of extracellular RNA functionality. In the presence of active ribonuclease inhibitors (RI), naked bacterial RNA is pro-inflammatory when spiked in the media of dendritic cells and macrophages. In murine cells, this response mainly depends on the action of endosomal Toll-like receptors. However, we also show that naked RNA can perform endosomal escape and engage with cytosolic RNA sensors and ribosomes. For example, naked mRNAs encoding reporter proteins can be spontaneously internalized and translated by a variety of cell types, in an RI-dependent manner. In vivo, RI co-injection enhances the activation induced by naked extracellular RNA on splenic lymphocytes and myeloid-derived leukocytes. Furthermore, naked extracellular RNA is inherently pro-inflammatory in ribonuclease-poor compartments such as the peritoneal cavity. Overall, these results demonstrate that naked RNA is bioactive and does not need encapsulation inside synthetic or biological lipid vesicles for functional uptake, making a case for nonvesicular extracellular RNA-mediated intercellular communication.

Serine/arginine-rich splicing factor 7 promotes the type I interferon response by activating Irf7 transcription

Tight regulation of macrophage immune gene expression is required to fight infection without risking harmful inflammation. The contribution of RNA-binding proteins (RBPs) to shaping the macrophage response to pathogens remains poorly understood. Transcriptomic analysis reveals that a member of the serine/arginine-rich (SR) family of mRNA processing factors, SRSF7, is required for optimal expression of a cohort of interferon-stimulated genes in macrophages. Using genetic and biochemical assays, we discover that in addition to its canonical role in regulating alternative splicing, SRSF7 drives transcription of interferon regulatory transcription factor 7 (IRF7) to promote antiviral immunity. At the Irf7 promoter, SRSF7 maximizes STAT1 transcription factor binding and RNA polymerase II elongation via cooperation with the H4K20me1 histone methyltransferase KMT5a (SET8). These studies define a role for an SR protein in activating transcription and reveal an RBP-chromatin network that orchestrates macrophage antiviral gene expression.

A multiomics dataset for the study of RNA modifications in human macrophage differentiation and polarisation

RNA modifications have emerged as central regulators of gene expression programs. Amongst RNA modifications are N6-methyladenosine (m6A) and RNA 5-hydroxymethylcytosine (5hmC). While m6A is established as a versatile regulator of RNA metabolism, the functions of RNA 5hmC are unclear. Despite some evidence linking RNA modifications to immunity, their implications in gene expression control in macrophage development and functions remain unclear. Here we present a multi-omics dataset capturing different layers of the gene expression programs driving macrophage differentiation and polarisation. We obtained mRNA-Seq, m6A-IP-Seq, 5hmC-IP-Seq, Polyribo-Seq and LC-MS/MS data from monocytes and resting-, pro- and anti-inflammatory-like macrophages. We present technical validation showing high quality and correlation between samples for all datasets, and evidence of biological consistency of modelled macrophages at the transcriptomic, epitranscriptomic, translational and proteomic levels. This multi-omics dataset provides a resource for the study of RNA m6A and 5hmC in the context of macrophage biology and spans the gene expression process from transcripts to proteins.

Replicative aging in yeast involves dynamic intron retention patterns associated with mRNA processing/export and protein ubiquitination

Saccharomyces cerevisiae (baker's yeast) has yielded relevant insights into some of the basic mechanisms of organismal aging. Among these are genomic instability, oxidative stress, caloric restriction and mitochondrial dysfunction. Several genes are known to have an impact on the aging process, with corresponding mutants exhibiting short- or long-lived phenotypes. Research dedicated to unraveling the underlying cellular mechanisms can support the identification of conserved mechanisms of aging in other species. One of the hitherto less studied fields in yeast aging is how the organism regulates its gene expression at the transcriptional level. To our knowledge, we present the first investigation into alternative splicing, particularly intron retention, during replicative aging of S. cerevisiae. This was achieved by utilizing the IRFinder algorithm on a previously published RNA-seq data set by Janssens et al. (2015). In the present work, 44 differentially retained introns in 43 genes were identified during replicative aging. We found that genes with altered intron retention do not display significant changes in overall transcript levels. It was possible to functionally assign distinct groups of these genes to the cellular processes of mRNA processing and export (e.g., YRA1) in early and middle-aged yeast, and protein ubiquitination (e.g., UBC5) in older cells. In summary, our work uncovers a previously unexplored layer of the transcriptional program of yeast aging and, more generally, expands the knowledge on the occurrence of alternative splicing in baker's yeast.

A comparative study of platelet storage lesion in platelet-rich plasma under cryopreservation

Platelet-rich plasma (PRP) has significant potential for various applications and holds clinical value in regenerative medicine. Cryopreservation is used to extend the preservation period of PRP, facilitating its clinical application. However, the potential negative effects of long-term cryopreservation on platelet storage lesion are still uncertain. In this study, PRP was stored at - 30 °C or - 80 °C. Platelet count, apoptosis, reactive oxygen species (ROS) content, and CD62P expression were assessed on the 14th and 28th days. The study also evaluated platelet mitochondria morphology and function, serotonin (5-HT) secretion by platelets, and the inflammatory activating effect of cryopreserved platelets in PRP. The results showed that there were no significant differences in platelet count, the content of 5-HT, and inflammatory effects between fresh PRP and PRP cryopreserved at both - 30 °C and - 80 °C. However, there was an increase in ROS level, apoptosis, and CD62P level after cryopreservation at both temperatures. Additionally, the levels of ROS, apoptosis, and CD62P in platelets were similar after storage at - 30 °C and - 80 °C. The main difference observed was that the morphology and function of mitochondria were severely damaged after storage at - 30 °C, while they were less affected at - 80 °C. Based on these findings, it can be concluded that storing PRP at - 80 °C is more suitable for achieving a better therapeutic effect in clinical applications, but cryopreservation could not replace the current standard.

SpliceWiz: interactive analysis and visualization of alternative splicing in R

Alternative splicing (AS) is a crucial mechanism for regulating gene expression and isoform diversity in eukaryotes. However, the analysis and visualization of AS events from RNA sequencing data remains challenging. Most tools require a certain level of computer literacy and the available means of visualizing AS events, such as coverage and sashimi plots, have limitations and can be misleading. To address these issues, we present SpliceWiz, an R package with an interactive Shiny interface that allows easy and efficient AS analysis and visualization at scale. A novel normalization algorithm is implemented to aggregate splicing levels within sample groups, thereby allowing group differences in splicing levels to be accurately visualized. The tool also offers downstream gene ontology enrichment analysis, highlighting ASEs belonging to functional pathways of interest. SpliceWiz is optimized for speed and efficiency and introduces a new file format for coverage data storage that is more efficient than BigWig. Alignment files are processed orders of magnitude faster than other R-based AS analysis tools and on par with command-line tools. Overall, SpliceWiz streamlines AS analysis, enabling reliable identification of functionally relevant AS events for further characterization. SpliceWiz is a Bioconductor package and is also available on GitHub (https://github.com/alexchwong/SpliceWiz).

Simultaneous detection of DNA variation and methylation at HLA class II locus and immune gene promoters using targeted SureSelect Methyl-Sequencing

The Human Leukocyte Antigen (HLA) locus associates with a variety of complex diseases, particularly autoimmune and inflammatory conditions. The HLA-DR15 haplotype, for example, confers the major risk for developing Multiple Sclerosis in Caucasians, pinpointing an important role in the etiology of this chronic inflammatory disease of the central nervous system. In addition to the protein-coding variants that shape the functional HLA-antigen-T cell interaction, recent studies suggest that the levels of HLA molecule expression, that are epigenetically controlled, also play a role in disease development. However, deciphering the exact molecular mechanisms of the HLA association has been hampered by the tremendous genetic complexity of the locus and a lack of robust approaches to investigate it. Here, we developed a method to specifically enrich the genomic DNA from the HLA class II locus (chr6:32,426,802-34,167,129) and proximal promoters of 2,157 immune-relevant genes, utilizing the Agilent RNA-based SureSelect Methyl-Seq Capture related method, followed by sequencing to detect genetic and epigenetic variation. We demonstrated successful simultaneous detection of the genetic variation and quantification of DNA methylation levels in HLA locus. Moreover, by the detection of differentially methylated positions in promoters of immune-related genes, we identified relevant pathways following stimulation of cells. Taken together, we present a method that can be utilized to study the interplay between genetic variance and epigenetic regulation in the HLA class II region, potentially, in a wide disease context.

Serine arginine-rich splicing factor (SRSF7) cooperates with the histone methyltransferase KMT5a to promote the type I interferon response via transcriptional activation of IRF7

Tight regulation of macrophage immune gene expression is required to fight infection without risking harmful inflammation. The contribution of RNA binding proteins (RBPs) to shaping the macrophage response to pathogens remains poorly understood. Transcriptomic analysis revealed that a member of the serine/arginine-rich (SR) family of mRNA processing factors, SRSF7, is required for optimal expression of a cohort of interferon stimulated genes (ISGs) in macrophages. Using genetic and biochemical assays, we discovered that in addition to its canonical role in regulating alternative splicing, SRSF7 drives transcription of interferon regulatory transcription factor 7 (IRF7) to promote antiviral immunity. At the Irf7 promoter, SRSF7 maximizes STAT1 transcription factor binding and RNA polymerase II elongation via cooperation with the H4K20me1 histone methyltransferase KMT5a (SET8). These studies define an unorthodox role for an SR protein in activating transcription and reveal an unappreciated RNA binding protein-chromatin network that orchestrates macrophage antiviral gene expression.

The central inflammatory regulator IκBζ: induction, regulation and physiological functions

IκBζ (encoded by NFKBIZ) is the most recently identified IkappaB family protein. As an atypical member of the IkappaB protein family, NFKBIZ has been the focus of recent studies because of its role in inflammation. Specifically, it is a key gene in the regulation of a variety of inflammatory factors in the NF-KB pathway, thereby affecting the progression of related diseases. In recent years, investigations into NFKBIZ have led to greater understanding of this gene. In this review, we summarize the induction of NFKBIZ and then elucidate its transcription, translation, molecular mechanism and physiological function. Finally, the roles played by NFKBIZ in psoriasis, cancer, kidney injury, autoimmune diseases and other diseases are described. NFKBIZ functions are universal and bidirectional, and therefore, this gene may exert a great influence on the regulation of inflammation and inflammation-related diseases.

The NR_109/FUBP1/c-Myc axis regulates TAM polarization and remodels the tumor microenvironment to promote cancer development

BACKGROUND

Tumor-associated macrophages (TAMs) are a major component of the tumor microenvironment (TME) and exert an important role in tumor progression. Due to the heterogeneity and plasticity of TAMs, modulating the polarization states of TAMs is considered as a potential therapeutic strategy for tumors. Long noncoding RNAs (lncRNAs) have been implicated in various physiological and pathological processes, yet the underlying mechanism on how lncRNAs manipulate the polarization states of TAMs is still unclear and remains to be further investigated.

METHODS

Microarray analyses were employed to characterize the lncRNA profile involved in THP-1-induced M0, M1 and M2-like macrophage. Among those differentially expressed lncRNAs, NR_109 was further studied, for its function in M2-like macrophage polarization and the effects of the condition medium or macrophages mediated by NR_109 on tumor proliferation, metastasis and TME remodeling both in vitro and in vivo. Moreover, we revealed how NR_109 interacted with far upstream element-binding protein 1 (FUBP1) to regulate the protein stability through hindering ubiquitination modification by competitively binding with JVT-1. Finally, we examined sections of tumor patients to probe the correlation among the expression of NR_109 and related proteins, showing the clinical significance of NR_109.

RESULTS

We found that lncRNA NR_109 was highly expressed in M2-like macrophages. Knockdown NR_109 impeded IL-4 induced M2-like macrophage polarization and significantly reduced the activity of M2-like macrophages to support the proliferation and metastasis of tumor cells in vitro and in vivo. Mechanistically, NR_109 competed with JVT-1 to bind FUBP1 at its C-terminus domain, impeded the ubiquitin-mediated degradation of FUBP1, activated c-Myc transcription and thus promoted M2-like macrophages polarization. Meanwhile, as a transcription factor, c-Myc could bind to the promoter of NR_109 and enhance the transcription of NR_109. Clinically, high NR_109 expression was found in CD163⁺ TAMs from tumor tissues and was positively correlated with poor clinical stages of patients with gastric cancer and breast cancer.

CONCLUSIONS

Our work revealed for the first time that NR_109 exerted a crucial role in regulating the phenotype-remodeling and function of M2-like macrophages via a NR_109/FUBP1/c-Myc positive feedback loop. Thus, NR_109 has great translational potentials in the diagnosis, prognosis and immunotherapy of cancer.

A Potential Therapy Using Antisense Oligonucleotides to Treat Autosomal Recessive Polycystic Kidney Disease

(1) Background: Autosomal recessive polycystic kidney disease (ARPKD) is a rare ciliopathy characterized by progressively enlarged kidneys with fusiform dilatation of the collecting ducts. Loss-of-function mutations in the PKHD1 gene, which encodes fibrocystin/polyductin, cause ARPKD; however, an efficient treatment method and drug for ARPKD have yet to be found. Antisense oligonucleotides (ASOs) are short special oligonucleotides which function to regulate gene expression and alter mRNA splicing. Several ASOs have been approved by the FDA for the treatment of genetic disorders, and many are progressing at present. We designed ASOs to verify whether ASOs mediate the correction of splicing further to treat ARPKD arising from splicing defects and explored them as a potential treatment option. (2) Methods: We screened 38 children with polycystic kidney disease for gene detection using whole-exome sequencing (WES) and targeted next-generation sequencing. Their clinical information was investigated and followed up. The PKHD1 variants were summarized and analyzed, and association analysis was carried out to analyze the relationship between genotype and phenotype. Various bioinformatics tools were used to predict pathogenicity. Hybrid minigene analysis was performed as part of the functional splicing analysis. Moreover, the de novo protein synthesis inhibitor cycloheximide was selected to verify the degraded pathway of abnormal pre-mRNAs. ASOs were designed to rescue aberrant splicing, and this was verified. (3) Results: Of the 11 patients with PKHD1 variants, all of them exhibited variable levels of complications of the liver and kidneys. We found that patients with truncating variants and variants in certain regions had a more severe phenotype. Two splicing variants of the PKHD1 genotypes were studied via the hybrid minigene assay: variants c.2141-3T>C and c.11174+5G>A. These cause aberrant splicing, and their strong pathogenicity was confirmed. We demonstrated that the abnormal pre-mRNAs produced from the variants escaped from the NMD pathway with the use of the de novo protein synthesis inhibitor cycloheximide. Moreover, we found that the splicing defects were rescued by using ASOs, which efficiently induced the exclusion of pseudoexons. (4) Conclusion: Patients with truncating variants and variants in certain regions had a more severe phenotype. ASOs are a potential drug for treating ARPKD patients harboring splicing mutations of the PKHD1 gene by correcting the splicing defects and increasing the expression of the normal PKHD1 gene.

A multiplex inhalation platform to model in situ like aerosol delivery in a breathing lung-on-chip

Prolonged exposure to environmental respirable toxicants can lead to the development and worsening of severe respiratory diseases such as asthma, chronic obstructive pulmonary disease (COPD) and fibrosis. The limited number of FDA-approved inhaled drugs for these serious lung conditions has led to a shift from in vivo towards the use of alternative in vitro human-relevant models to better predict the toxicity of inhaled particles in preclinical research. While there are several inhalation exposure models for the upper airways, the fragile and dynamic nature of the alveolar microenvironment has limited the development of reproducible exposure models for the distal lung. Here, we present a mechanistic approach using a new generation of exposure systems, the Cloud α AX12. This novel in vitro inhalation tool consists of a cloud-based exposure chamber (VITROCELL) that integrates the breathing ^AXLung-on-chip system (AlveoliX). The ultrathin and porous membrane of the AX12 plate was used to create a complex multicellular model that enables key physiological culture conditions: the air-liquid interface (ALI) and the three-dimensional cyclic stretch (CS). Human-relevant cellular models were established for a) the distal alveolar-capillary interface using primary cell-derived immortalized alveolar epithelial cells (^AXiAECs), macrophages (THP-1) and endothelial (HLMVEC) cells, and b) the upper-airways using Calu3 cells. Primary human alveolar epithelial cells (^AXhAEpCs) were used to validate the toxicity results obtained from the immortalized cell lines. To mimic in vivo relevant aerosol exposures with the Cloud α AX12, three different models were established using: a) titanium dioxide (TiO2) and zinc oxide nanoparticles b) polyhexamethylene guanidine a toxic chemical and c) an anti-inflammatory inhaled corticosteroid, fluticasone propionate (FL). Our results suggest an important synergistic effect on the air-blood barrier sensitivity, cytotoxicity and inflammation, when air-liquid interface and cyclic stretch culture conditions are combined. To the best of our knowledge, this is the first time that an in vitro inhalation exposure system for the distal lung has been described with a breathing lung-on-chip technology. The Cloud α AX12 model thus represents a state-of-the-art pre-clinical tool to study inhalation toxicity risks, drug safety and efficacy.

Single-cell and long-read sequencing to enhance modelling of splicing and cell-fate determination

Single-cell sequencing technologies have revolutionised the life sciences and biomedical research. Single-cell sequencing provides high-resolution data on cell heterogeneity, allowing high-fidelity cell type identification, and lineage tracking. Computational algorithms and mathematical models have been developed to make sense of the data, compensate for errors and simulate the biological processes, which has led to breakthroughs in our understanding of cell differentiation, cell-fate determination and tissue cell composition. The development of long-read (a.k.a. third-generation) sequencing technologies has produced powerful tools for investigating alternative splicing, isoform expression (at the RNA level), genome assembly and the detection of complex structural variants (at the DNA level). In this review, we provide an overview of the recent advancements in single-cell and long-read sequencing technologies, with a particular focus on the computational algorithms that help in correcting, analysing, and interpreting the resulting data. Additionally, we review some mathematical models that use single-cell and long-read sequencing data to study cell-fate determination and alternative splicing, respectively. Moreover, we highlight the emerging opportunities in modelling cell-fate determination that result from the combination of single-cell and long-read sequencing technologies.

Expression and prognostic value of MIP-1α in neonatal acute respiratory distress syndrome

OBJECTIVE

To investigate the expression and prognostic value of macrophage inflammatory protein 1α (MIP-1α) in neonatal acute respiratory distress syndrome (NRDS).

METHODS

In this retrospective analysis, 96 newborns with NRDS in Affiliated Lianyungang Hospital of Xuzhou Medical University from January 2018 to June 2021 were included in the experimental group (EG), while the other 60 normal neonates were included as the control group (CG). The concentration of MIP-1α in umbilical cord blood was tested by Elisa method. The clinical value of MIP-1α in diagnosing NRDS was assessed via receiver operating characteristic (ROC) curve. According to the 28-day survival data, children were divided into a survival group and a death group. The prognostic factors were assessed by Cox regression analysis. The correlation between MIP-1α and IL-1β, IL-6, TNF-α, SNAPPE-II scores were evaluated by Pearson test. The relationship between the MIP-1α level and the severity of the disease was assessed.

RESULTS

The MIF-1α level in cord blood of children in the EG was dramatically higher than that in the CG (P<0.05). Besides, ROC curve further found that the area of MIF-1α under the curve of diagnosing NRSD was 0.949. MIF-1α was positively correlated with the levels of interleukin-1β (IL-1β), interleukin-6 (IL-6), tumor necrosis factor α (TNF-α) and SNAPPE-II score (P<0.001). With the increase of NRDS, the serum MIF-1α level increased, showing a positive association (P<0.05). Cox regression analysis revealed that the severity and MIF-1α level were independent prognostic factors of survival (P<0.001). The survival rate of children with MIF-1α <281.58 pg/mL as well as children with I-II grade was higher than those with MIF-1α >281.58 pg/mL as well as children with III-IV grade (P<0.05).

CONCLUSION

The increase of serum MIP-1α level is relevant to the condition and prognosis of NRDS children. The level of cord blood MIP-1α level is expected to become a potential outcome measure.

Increased chromatin accessibility facilitates intron retention in specific cell differentiation states

Dynamic intron retention (IR) in vertebrate cells is of widespread biological importance. Aberrant IR is associated with numerous human diseases including several cancers. Despite consistent reports demonstrating that intrinsic sequence features can help introns evade splicing, conflicting findings about cell type- or condition-specific IR regulation by trans-regulatory and epigenetic mechanisms demand an unbiased and systematic analysis of IR in a controlled experimental setting. We integrated matched mRNA sequencing (mRNA-Seq), whole-genome bisulfite sequencing (WGBS), nucleosome occupancy methylome sequencing (NOMe-Seq) and chromatin immunoprecipitation sequencing (ChIP-Seq) data from primary human myeloid and lymphoid cells. Using these multi-omics data and machine learning, we trained two complementary models to determine the role of epigenetic factors in the regulation of IR in cells of the innate immune system. We show that increased chromatin accessibility, as revealed by nucleosome-free regions, contributes substantially to the retention of introns in a cell-specific manner. We also confirm that intrinsic characteristics of introns are key for them to evade splicing. This study suggests an important role for chromatin architecture in IR regulation. With an increasing appreciation that pathogenic alterations are linked to RNA processing, our findings may provide useful insights for the development of novel therapeutic approaches that target aberrant splicing.

Stimulus-specific remodeling of the neuronal transcriptome through nuclear intron-retaining transcripts

The nuclear envelope has long been considered primarily a physical barrier separating nuclear and cytosolic contents. More recently, nuclear compartmentalization has been shown to have additional regulatory functions in controlling gene expression. A sizeable proportion of protein-coding mRNAs is more prevalent in the nucleus than in the cytosol, suggesting regulated mRNA trafficking to the cytosol, but the mechanisms underlying controlled nuclear mRNA retention remain unclear. Here, we provide a comprehensive map of the subcellular localization of mRNAs in mature mouse cortical neurons, and reveal that transcripts retained in the nucleus comprise the majority of stable intron-retaining mRNAs. Systematically probing the fate of nuclear transcripts upon neuronal stimulation, we found opposite effects on sub-populations of transcripts: while some are targeted for degradation, others complete splicing to generate fully mature mRNAs that are exported to the cytosol and mediate rapid increases in protein levels. Finally, different forms of stimulation mobilize distinct groups of intron-retaining transcripts, with this selectivity arising from the activation of specific signaling pathways. Overall, our findings uncover a cue-specific control of intron retention as a major regulator of acute remodeling of the neuronal transcriptome.

Transcriptional regulatory networks of circulating immune cells in type 1 diabetes: A community knowledgebase

Investigator-generated transcriptomic datasets interrogating circulating immune cell (CIC) gene expression in clinical type 1 diabetes (T1D) have underappreciated re-use value. Here, we repurposed these datasets to create an open science environment for the generation of hypotheses around CIC signaling pathways whose gain or loss of function contributes to T1D pathogenesis. We firstly computed sets of genes that were preferentially induced or repressed in T1D CICs and validated these against community benchmarks. We then inferred and validated signaling node networks regulating expression of these gene sets, as well as differentially expressed genes in the original underlying T1D case:control datasets. In a set of three use cases, we demonstrated how informed integration of these networks with complementary digital resources supports substantive, actionable hypotheses around signaling pathway dysfunction in T1D CICs. Finally, we developed a federated, cloud-based web resource that exposes the entire data matrix for unrestricted access and re-use by the research community.

•

Re-use of transcriptomic type 1 diabetes (T1D) circulating immune cells (CICs) datasets

•

We generated transcriptional regulatory networks for T1D CICs

•

Use cases generate substantive hypotheses around signaling pathway dysfunction in T1D CICs

•

Networks are freely accessible on the web for re-use by the research community

3’UTR-Seq analysis of chicken abdominal adipose tissue reveals widespread intron retention in 3’UTR and provides insight into molecular basis of feed efficiency

Feed efficiency (FE) is an important trait in the broiler industry due to its direct correlation to efficient muscle growth instead of fat deposition. The present study characterized and compared gene expression profiles in abdominal fat from broiler chickens of different FE levels to enhance the understanding of FE biology. Specifically, traditional whole-transcript RNA-sequencing (RNA-seq) and 3’ UTR-sequencing (3’ UTR-seq) were applied to 22 and 61 samples, respectively. Overall, these two sequencing techniques shared a high correlation (0.76) between normalized counts, although 3’ UTR-seq showed a higher variance in sequencing and mapping performance statistics across samples and a lower rate of uniquely mapped reads. A higher percentage of 3’ UTR-seq reads mapped to introns suggested the frequent presence of cleavage sites in introns, thus warranting future research to study its regulatory function. Differential expression analysis identified 1198 differentially expressed genes (DEGs) between high FE (HFE) and intermediate FE (IFE) chickens with False Discovery Rate < 0.05 and fold change > 1.2. The processes that were significantly enriched by the DEGs included extracellular matrix remodeling and mechanisms impacting gene expression at the transcriptional and translational levels. Gene ontology enrichment analysis suggested that the divergence in fat deposition and FE in broiler chickens could be associated with peroxisome and lipid metabolism possibly regulated by G0/G1 switch gene 2 (G0S2).

The role of glycolysis and lactate in the induction of tumor-associated macrophages immunosuppressive phenotype

Growing evidence highlights that glycolysis and tumor-derived lactate could skew tumor-associated macrophages (TAMs) toward an immunosuppressive phenotype. However, the updated research has not been systematically summarized yet. TAMs are educated by the tumor microenvironment (TME) and exert immunosuppressive functions and tumorigenic effects via multiple biological processes. It is well known that lactate generated by aerobic glycolysis is significantly accumulated in TME and promotes tumor progression in solid tumors. Moreover, some recent research demonstrated that glycolysis is activated in TAMs to support M2-like polarization, which is absolutely in contrast with the metabolic profile of M2 macrophages in inflammation. Notably, lactate produced by high levels of glycolysis is not only a metabolic by-product but also an oncometabolite. TAMs could access the biological information delivered by lactate and further enhance protumor functions such as immunosuppression and angiogenesis. Here, we outline the connection between glycolysis and TAM phenotype to elucidate the metabolic characteristics of TAMs. Further, insights into the specific molecular mechanisms of lactate-induced TAM polarization and potential therapeutic targets are summarized. We sought to discuss the reciprocal interaction between tumor cells and TAMs mediated by lactate, which will lay a foundation for the research aiming to elucidate the complex functions of TAMs.

Comment on: The m6A Reader IGF2BP2 Regulates Macrophage Phenotypic Activation and Inflammatory Diseases by Stabilizing TSC1 and PPARγ

Recently, first insights into the regulation and the role of the RNA-binding protein IMP2 in macrophage activation have been published by Wang et al. This study addresses differences in the regulation of IMP2 between the human and murine system. While the expression of IMP2 in anti-inflammatory macrophages is synchronous in mice and men, IMP2 expression is regulated differently in inflammatory macrophages.

Towards resolution of the intron retention paradox in breast cancer

BACKGROUND

After many years of neglect in the field of alternative splicing, the importance of intron retention (IR) in cancer has come into focus following landmark discoveries of aberrant IR patterns in cancer. Many solid and liquid tumours are associated with drastic increases in IR, and such patterns have been pursued as both biomarkers and therapeutic targets. Paradoxically, breast cancer (BrCa) is the only tumour type in which IR is reduced compared to adjacent normal breast tissue.

METHODS

In this study, we have conducted a pan-cancer analysis of IR with emphasis on BrCa and its subtypes. We explored mechanisms that could cause aberrant and pathological IR and clarified why normal breast tissue has unusually high IR.

RESULTS

Strikingly, we found that aberrantly decreasing IR in BrCa can be largely attributed to normal breast tissue having the highest occurrence of IR events compared to other healthy tissues. Our analyses suggest that low numbers of IR events in breast tumours are associated with poor prognosis, particularly in the luminal B subtype. Interestingly, we found that IR frequencies negatively correlate with cell proliferation in BrCa cells, i.e. rapidly dividing tumour cells have the lowest number of IR events. Aberrant RNA-binding protein expression and changes in tissue composition are among the causes of aberrantly decreasing IR in BrCa.

CONCLUSIONS

Our results suggest that IR should be considered for therapeutic manipulation in BrCa patients with aberrantly low IR levels and that further work is needed to understand the cause and impact of high IR in other tumour types.

Selection and validation of reference genes for qPCR analysis of differentiation and maturation of THP-1 cells into M1 macrophage-like cells

For cell-based assays studying monocytes and macrophages, the immortalized monocyte cell line THP-1 is widely used and stimulated with phorbol 12-myristate 13-acetate, lipopolysaccharide (LPS) and/or interferon-γ (IFN-γ), after which it differentiates and polarizes into proinflammatory M1-like macrophages. For the quantification of this and the effect of different factors affecting these processes, the expression levels of various maturation markers are determined using reverse transcription-quantitative PCR. For this purpose, stably expressed reference genes are crucial. However, no studies evaluating the stability of reference genes in THP-1 cells stimulated with LPS and IFN-γ have been performed. Therefore, this paper describes the selection of the most used reference genes [RPL37A (ribosomal protein L37a), GAPDH (glyceraldehyde-3-phosphate dehydrogenase), UBC (ubiquitin C), B2M (0β2-microbulin), ACTB (β-actin) and PPIA (cyclophilin A)], the in silico primer design, the analysis and the validation of these in accordance with the MIQE (Minimum Information for Publication of Quantitative Real-Time PCR Experiments) guidelines and more recent recommendations for the validation of the stability of reference genes. Using the RefFinder platform, including the four most popular algorithms for reference gene validation, the Delta CT, BestKeeper, NormFinder and geNorm, we find the reference genes GAPDH and UBC to be the most stable. Furthermore, we demonstrate that the normalization of gene expression data using the least stable reference genes, ACTB and B2M, dramatically affects the interpretation of experimental data. Taken together, it is vital to validate the stability of reference genes under the specific experimental conditions used when utilizing the THP-1 monocyte model system.

Dynamic intron retention modulates gene expression in the monocytic differentiation pathway

Monocytes play a crucial role in maintaining homeostasis and mediating a successful innate immune response. They also act as central players in diverse pathological conditions, thus making them an attractive therapeutic target. Within the bone marrow, monocytes arise from a committed precursor termed cMoP (Common Monocyte Progenitor). However, molecular mechanisms that regulate the differentiation of cMoP to various monocytic subsets remain unclear. Herein, we purified murine myeloid precursors for deep poly-A enriched RNA sequencing to understand the role of alternative splicing in the development and differentiation of monocytes under homeostasis. Our analyses revealed intron retention to be the major alternative splicing mechanism involved in the monocyte differentiation cascade, especially in the differentiation of Ly6C^hi monocytes to Ly6C^lo monocytes. Furthermore, we found that the key genes regulated by intron retention in the differentiation of murine Ly6C^hi to Ly6C^lo monocytes were also conserved in humans. Our data highlight the unique role of intron retention in the regulation of the monocytic differentiation pathway.

IRFinder-S: a comprehensive suite to discover and explore intron retention

Accurate quantification and detection of intron retention levels require specialized software. Building on our previous software, we create a suite of tools called IRFinder-S, to analyze and explore intron retention events in multiple samples. Specifically, IRFinder-S allows a better identification of true intron retention events using a convolutional neural network, allows the sharing of intron retention results between labs, integrates a dynamic database to explore and contrast available samples, and provides a tested method to detect differential levels of intron retention.

Sequence determinants as key regulators in gene expression of T cells

T cell homeostasis, T cell differentiation, and T cell effector function rely on the constant fine-tuning of gene expression. To alter the T cell state, substantial remodeling of the proteome is required. This remodeling depends on the intricate interplay of regulatory mechanisms, including post-transcriptional gene regulation. In this review, we discuss how the sequence of a transcript influences these post-transcriptional events. In particular, we review how sequence determinants such as sequence conservation, GC content, and chemical modifications define the levels of the mRNA and the protein in a T cell. We describe the effect of different forms of alternative splicing on mRNA expression and protein production, and their effect on subcellular localization. In addition, we discuss the role of sequences and structures as binding hubs for miRNAs and RNA-binding proteins in T cells. The review thus highlights how the intimate interplay of post-transcriptional mechanisms dictate cellular fate decisions in T cells.

The Expanding Role of Alternative Splicing in Vascular Smooth Muscle Cell Plasticity

Vascular smooth muscle cells (VSMCs) display extraordinary phenotypic plasticity. This allows them to differentiate or dedifferentiate, depending on environmental cues. The ability to ‘switch’ between a quiescent contractile phenotype to a highly proliferative synthetic state renders VSMCs as primary mediators of vascular repair and remodelling. When their plasticity is pathological, it can lead to cardiovascular diseases such as atherosclerosis and restenosis. Coinciding with significant technological and conceptual innovations in RNA biology, there has been a growing focus on the role of alternative splicing in VSMC gene expression regulation. Herein, we review how alternative splicing and its regulatory factors are involved in generating protein diversity and altering gene expression levels in VSMC plasticity. Moreover, we explore how recent advancements in the development of splicing-modulating therapies may be applied to VSMC-related pathologies.

Activation of cGAS/STING pathway upon paramyxovirus infection

During inflammatory diseases, cancer, and infection, the cGAS/STING pathway is known to recognize foreign or self-DNA in the cytosol and activate an innate immune response. Here, we report that negative-strand RNA paramyxoviruses, Nipah virus (NiV), and measles virus (MeV), can also trigger the cGAS/STING axis. Although mice deficient for MyD88, TRIF, and MAVS still moderately control NiV infection when compared with wild-type mice, additional STING deficiency resulted in 100% lethality, suggesting synergistic roles of these pathways in host protection. Moreover, deletion of cGAS or STING resulted in decreased type I interferon production with enhanced paramyxoviral infection in both human and murine cells. Finally, the phosphorylation and ubiquitination of STING, observed during viral infections, confirmed the activation of cGAS/STING pathway by NiV and MeV. Our data suggest that cGAS/STING activation is critical in controlling paramyxovirus infection and possibly represents attractive targets to develop countermeasures against severe disease induced by these pathogens.

Integrative functional genomic analysis of intron retention in human and mouse brain with Alzheimer's disease

Intron retention (IR) has been implicated in the pathogenesis of complex diseases such as cancers; its association with Alzheimer's disease (AD) remains unexplored. We performed genome-wide analysis of IR through integrating genetic, transcriptomic, and proteomic data of AD subjects and mouse models from the Accelerating Medicines Partnership-Alzheimer's Disease project. We identified 4535 and 4086 IR events in 2173 human and 1736 mouse genes, respectively. Quantitation of IR enabled the identification of differentially expressed genes that conventional exon-level approaches did not reveal. There were significant correlations of intron expression within innate immune genes, like HMBOX1, with AD in humans. Peptides with a high probability of translation from intron-retained mRNAs were identified using mass spectrometry. Further, we established AD-specific intron expression Quantitative Trait Loci, and identified splicing-related genes that may regulate IR. Our analysis provides a novel resource for the search for new AD biomarkers and pathological mechanisms.

Reactive astrocytes in ALS display diminished intron retention

Reactive astrocytes are implicated in amyotrophic lateral sclerosis (ALS), although the mechanisms controlling reactive transformation are unknown. We show that decreased intron retention (IR) is common to human-induced pluripotent stem cell (hiPSC)-derived astrocytes carrying ALS-causing mutations in VCP, SOD1 and C9orf72. Notably, transcripts with decreased IR and increased expression are overrepresented in reactivity processes including cell adhesion, stress response and immune activation. This was recapitulated in public-datasets for (i) hiPSC-derived astrocytes stimulated with cytokines to undergo reactive transformation and (ii) in vivo astrocytes following selective deletion of TDP-43. We also re-examined public translatome sequencing (TRAP-seq) of astrocytes from a SOD1 mouse model, which revealed that transcripts upregulated in translation significantly overlap with transcripts exhibiting decreased IR. Using nucleocytoplasmic fractionation of VCP mutant astrocytes coupled with mRNA sequencing and proteomics, we identify that decreased IR in nuclear transcripts is associated with enhanced nonsense mediated decay and increased cytoplasmic expression of transcripts and proteins regulating reactive transformation. These findings are consistent with a molecular model for reactive transformation in astrocytes whereby poised nuclear reactivity-related IR transcripts are spliced, undergo nuclear-to-cytoplasmic translocation and translation. Our study therefore provides new insights into the molecular regulation of reactive transformation in astrocytes.

Aberrant functional connectivity between the suprachiasmatic nucleus and the superior temporal gyrus: Bridging RORA gene polymorphism with diurnal mood variation in major depressive disorder

Diurnal mood variation (DMV), a common symptom of major depressive disorder (MDD), is associated with circadian related genes and dysregulation of the suprachiasmatic nucleus (SCN). Previous research confirmed that the RORA gene is involved in the regulation of circadian rhythms. In this study, we hypothesized that polymorphisms of RORA may affect DMV symptoms of MDD through functional changes in the SCN. A total of 208 patients diagnosed with depression and 120 control subjects were enrolled and underwent a resting-state functional magnetic resonance imaging (rs-fMRI). Blood samples were collected and genotyping of 9 RORA gene SNPs were performed using next-generation sequencing technology. Patients were categorized as an AA genotype or C allele carriers based on RORA rs72752802 polymorphism. SCN-seed functional connectivity (FC) was compared between the two groups and correlation with severity of DMV was analyzed. Finally, a mediation analysis was performed to further determine FC intermediary effects. We observed that rs72752802 was significantly associated with patients' DMV symptoms. C allele carriers of rs72752802 showed significantly decreased FC between the right SCN and right superior temporal gyrus (rSTG). This was also correlated with DMV symptoms. In addition, the rs72752802 SNP influenced DMV symptoms through intermediary effects of SCN-rSTG connectivity. The study presented here provides a neurological and genetic basis for understanding depressed patients experiencing DMV.

Widespread Aberrant Alternative Splicing despite Molecular Remission in Chronic Myeloid Leukaemia Patients

Simple Summary

This study provides new insights into the changing transcriptomic and epigenomic landscapes in chronic myeloid leukaemia (CML) patients who are receiving tyrosine kinase inhibitor (TKI) therapy (often life-long). Alternative splicing, vital for cellular homeostasis, is dysregulated in human cancers. Remarkably, we found abnormal splicing patterns despite molecular remission in peripheral blood cells of chronic-phase CML patients. This phenomenon is independent of the TKI drug used and in striking contrast to the normalisation of gene expression and DNA methylation patterns.

Abstract

Vast transcriptomics and epigenomics changes are characteristic of human cancers, including leukaemia. At remission, we assume that these changes normalise so that omics-profiles resemble those of healthy individuals. However, an in-depth transcriptomic and epigenomic analysis of cancer remission has not been undertaken. A striking exemplar of targeted remission induction occurs in chronic myeloid leukaemia (CML) following tyrosine kinase inhibitor (TKI) therapy. Using RNA sequencing and whole-genome bisulfite sequencing, we profiled samples from chronic-phase CML patients at diagnosis and remission and compared these to healthy donors. Remarkably, our analyses revealed that abnormal splicing distinguishes remission samples from normal controls. This phenomenon is independent of the TKI drug used and in striking contrast to the normalisation of gene expression and DNA methylation patterns. Most remarkable are the high intron retention (IR) levels that even exceed those observed in the diagnosis samples. Increased IR affects cell cycle regulators at diagnosis and splicing regulators at remission. We show that aberrant splicing in CML is associated with reduced expression of specific splicing factors, histone modifications and reduced DNA methylation. Our results provide novel insights into the changing transcriptomic and epigenomic landscapes of CML patients during remission. The conceptually unanticipated observation of widespread aberrant alternative splicing after remission induction warrants further exploration. These results have broad implications for studying CML relapse and treating minimal residual disease.

Endoglin: An 'Accessory' Receptor Regulating Blood Cell Development and Inflammation

Transforming growth factor-β1 (TGF-β1) is a pleiotropic factor sensed by most cells. It regulates a broad spectrum of cellular responses including hematopoiesis. In order to process TGF-β1-responses in time and space in an appropriate manner, there is a tight regulation of its signaling at diverse steps. The downstream signaling is mediated by type I and type II receptors and modulated by the ‘accessory’ receptor Endoglin also termed cluster of differentiation 105 (CD105). Endoglin was initially identified on pre-B leukemia cells but has received most attention due to its high expression on activated endothelial cells. In turn, Endoglin has been figured out as the causative factor for diseases associated with vascular dysfunction like hereditary hemorrhagic telangiectasia-1 (HHT-1), pre-eclampsia, and intrauterine growth restriction (IUPR). Because HHT patients often show signs of inflammation at vascular lesions, and loss of Endoglin in the myeloid lineage leads to spontaneous inflammation, it is speculated that Endoglin impacts inflammatory processes. In line, Endoglin is expressed on progenitor/precursor cells during hematopoiesis as well as on mature, differentiated cells of the innate and adaptive immune system. However, so far only pro-monocytes and macrophages have been in the focus of research, although Endoglin has been identified in many other immune system cell subsets. These findings imply a functional role of Endoglin in the maturation and function of immune cells. Aside the functional relevance of Endoglin in endothelial cells, CD105 is differentially expressed during hematopoiesis, arguing for a role of this receptor in the development of individual cell lineages. In addition, Endoglin expression is present on mature immune cells of the innate (i.e., macrophages and mast cells) and the adaptive (i.e., T-cells) immune system, further suggesting Endoglin as a factor that shapes immune responses. In this review, we summarize current knowledge on Endoglin expression and function in hematopoietic precursors and mature hematopoietic cells of different lineages.

Holding on to Junk Bonds: Intron Retention in Cancer and Therapy

Intron retention (IR) in cancer was for a long time overlooked by the scientific community, as it was previously considered to be an artifact of a dysfunctional spliceosome. Technological advancements made in the last decade offer unique opportunities to explore the role of IR as a widespread phenomenon that contributes to the transcriptional diversity of many cancers. Numerous studies in cancer have shed light on dysregulation of cellular mechanisms that lead to aberrant and pathologic IR. IR is not merely a mechanism of gene regulation, but rather it can mediate cancer pathogenesis and therapeutic resistance in various human diseases. The burden of IR in cancer is governed by perturbations to mechanisms known to regulate this phenomenon and include epigenetic variation, mutations within the gene body, and splicing factor dysregulation. This review summarizes possible causes for aberrant IR and discusses the role of IR in therapy or as a consequence of disease treatment. As neoepitopes originating from retained introns can be presented on the cancer cell surface, the development of personalized cancer vaccines based on IR-derived neoepitopes should be considered. Ultimately, a deeper comprehension about the origins and consequences of aberrant IR may aid in the development of such personalized cancer vaccines.

Ctcf haploinsufficiency mediates intron retention in a tissue-specific manner

CTCF is a master regulator of gene transcription and chromatin organisation with occupancy at thousands of DNA target sites genome-wide. While CTCF is essential for cell survival, CTCF haploinsufficiency is associated with tumour development and hypermethylation. Increasing evidence demonstrates CTCF as a key player in several mechanisms regulating alternative splicing (AS), however, the genome-wide impact of Ctcf dosage on AS has not been investigated.

We examined the effect of Ctcf haploinsufficiency on gene expression and AS in five tissues from Ctcf hemizygous (Ctcf^+/-) mice. Reduced Ctcf levels caused distinct tissue-specific differences in gene expression and AS in all tissues. An increase in intron retention (IR) was observed in Ctcf^+/- liver and kidney. In liver, this specifically impacted genes associated with cytoskeletal organisation, splicing and metabolism. Strikingly, most differentially retained introns were short, with a high GC content and enriched in Ctcf binding sites in their proximal upstream genomic region. This study provides new insights into the effects of CTCF haploinsufficiency on organ transcriptomes and the role of CTCF in AS regulation.