Genome-wide transcriptome analysis identifies alternative splicing regulatory network and key splicing factors in mouse and human psoriasis

RNA-binding proteins potentially regulate alternative splicing of immune/inflammatory-associated genes during the progression of generalized pustular psoriasis

Generalized pustular psoriasis (GPP) is a rare but severe form of psoriasis. However, the pathogenesis of GPP has not been fully elucidated. Although RNA-binding proteins (RBPs) and the alternative splicing (AS) process are essential for regulating post-transcriptional gene expression, their roles in GPP are still unclear. We aimed to elucidate the regulatory mechanisms to identify potential new therapeutic targets. Here, We analyzed an RNA sequencing (RNA-seq) dataset (GSE200977) of peripheral blood mononuclear cells (PBMCs) of 24 patients with GPP, psoriasis vulgaris (PV), and healthy controls (HCs) from the Gene Expression Omnibus (GEO) database. We found that the abnormal alternative splicing (AS) events associated with GPP were mainly "alt3p/alt5p", and 15 AS genes were differentially expressed. Notably, the proportions of different immune cell types were correlated with the expression levels of regulatory alternatively spliced genes (RASGs): significant differences were observed in expression levels of DTD2, NDUFAF3, NBPF15, and FBLN7 in B cells and ARFIP1, IPO11, and RP11-326L24.9 in neutrophils in the GPP samples. Furthermore, We identified 32 differentially expressed RNA-binding proteins (RBPs) (18 up-regulated and 14 down-regulated). Co-expression networks between 14 pairs of differentially expressed RBPs and RASGs were subsequently constructed, demonstrating that these differentially expressed RBPs may affect the progression of GPP by regulating the AS of downstream immune/inflammatory-related genes such as LINC00989, ENC1 and MMP25-AS1. Our results were innovative in revealing the involvement of inflammation-related RBPs and RASGs in the development of GPP from the perspective of RBP-regulated AS.

Secukinumab and Dead Sea Climatotherapy Impact Resolved Psoriasis Skin Differently Potentially Affecting Disease Memory

Secukinumab and Dead Sea treatment result in clear skin for many psoriasis patients, through distinct mechanisms. However, recurrence in the same areas after treatments suggests the existence of a molecular scar. We aimed to compare the molecular and genetic differences in psoriasis patients who achieved complete response from secukinumab and Dead Sea climatotherapy treatments. We performed quantitative immunohistochemical and transcriptomic analysis, in addition to digital spatial profiling of skin punch biopsies. Histologically, both treatments resulted in a normalization of the lesional skin to a level resembling nonlesional skin. Interestingly, the transcriptome was not normalized by either treatments. We revealed 479 differentially expressed genes between secukinumab and Dead Sea climatotherapy at the end of treatment, with a psoriasis panel identifying SERPINB4, SERPINB13, IL36G, IL36RN, and AKR1B10 as upregulated in Dead Sea climatotherapy compared with secukinumab. Using digital spatial profiling, pan-RAS was observed to be differentially expressed in the microenvironment surrounding CD103+ cells, and IDO1 was differentially expressed in the dermis when comparing the two treatments. The differences observed between secukinumab and Dead Sea climatotherapy suggest the presence of a molecular scar, which may stem from mechanistically different pathways and potentially contribute to disease recurrence. This may be important for determining treatment response duration and disease memory.

Single-cell sequencing analysis and multiple machine-learning models revealed the cellular crosstalk of dendritic cells and identified FABP5 and KLRB1 as novel biomarkers for psoriasis

Background

Psoriasis is an immune-mediated disorder influenced by environmental factors on a genetic basis. Despite advancements, challenges persist, including the diminishing efficacy of biologics and small-molecule targeted agents, alongside managing recurrence and psoriasis-related comorbidities. Unraveling the underlying pathogenesis and identifying valuable biomarkers remain pivotal for diagnosing and treating psoriasis.

Methods

We employed a series of bioinformatics (including single-cell sequencing data analysis and machine learning techniques) and statistical methods to integrate and analyze multi-level data. We observed the cellular changes in psoriatic skin tissues, screened the key genes Fatty acid binding protein 5 (FABP5) and The killer cell lectin-like receptor B1 (KLRB1), evaluated the efficacy of six widely prescribed drugs on psoriasis treatment in modulating the dendritic cell-associated pathway, and assessed their overall efficacy. Finally, RT-qPCR, immunohistochemistry, and immunofluorescence assays were used to validate.

Results

The regulatory influence of dendritic cells (DCs) on T cells through the CD70/CD27 signaling pathway may emerge as a significant facet of the inflammatory response in psoriasis. Notably, FABP5 and KLRB1 exhibited up-regulation and co-localization in psoriatic skin tissues and M5-induced HaCaT cells, serving as potential biomarkers influencing psoriasis development.

Conclusion

Our study analyzed the impact of DC-T cell crosstalk in psoriasis, elucidated the characterization of two biomarkers, FABP5 and KLRB1, in psoriasis, and highlighted the promise and value of tofacitinib in psoriasis therapy targeting DCs.

Covalent Targeting of Splicing in T Cells

Despite significant interest in therapeutic targeting of splicing, few chemical probes are available for the proteins involved in splicing. Here, we show that elaborated stereoisomeric acrylamide chemical probe EV96 and its analogues lead to a selective T cell state-dependent loss of interleukin 2-inducible T cell kinase (ITK) by targeting one of the core splicing factors SF3B1. Mechanistic investigations suggest that the state-dependency stems from a combination of differential protein turnover rates and availability of functional mRNA pools that can be depleted due to extensive alternative splicing. We further introduce a comprehensive list of proteins involved in splicing and leverage both cysteine- and protein-directed activity-based protein profiling (ABPP) data with electrophilic scout fragments to demonstrate covalent ligandability for many classes of splicing factors and splicing regulators in primary human T cells. Taken together, our findings show how chemical perturbation of splicing can lead to immune state-dependent changes in protein expression and provide evidence for the broad potential to target splicing factors with covalent chemistry.

IL-17D-induced inhibition of DDX5 expression in keratinocytes amplifies IL-36R-mediated skin inflammation

Aberrant RNA splicing in keratinocytes drives inflammatory skin disorders. In the present study, we found that the RNA helicase DDX5 was downregulated in keratinocytes from the inflammatory skin lesions in patients with atopic dermatitis and psoriasis, and that mice with keratinocyte-specific deletion of Ddx5 (Ddx5^∆KC) were more susceptible to cutaneous inflammation. Inhibition of DDX5 expression in keratinocytes was induced by the cytokine interleukin (IL)-17D through activation of the CD93-p38 MAPK-AKT-SMAD2/3 signaling pathway and led to pre-messenger RNA splicing events that favored the production of membrane-bound, intact IL-36 receptor (IL-36R) at the expense of soluble IL-36R (sIL-36R) and to the selective amplification of IL-36R-mediated inflammatory responses and cutaneous inflammation. Restoration of sIL-36R in Ddx5^∆KC mice with experimental atopic dermatitis or psoriasis suppressed skin inflammation and alleviated the disease phenotypes. These findings indicate that IL-17D modulation of DDX5 expression controls inflammation in keratinocytes during inflammatory skin diseases.

DAZAP1 overexpression promotes growth of HCC cell lines: a primary study using CEUS

PURPOSE

Hepatocellular carcinoma (HCC) is one of the most common types of hepatic carcinoma. The overall prognosis is poor. DAZAP1, a regulator of alternative splicing (AS) events, may participate in tumor growth.

METHODS

We collected 105 HCC patients and tissue samples from the Department of Hepatological Surgery in the Second Affiliated Hospital of Qiqihar Medical University. TCGA datasets were downloaded and operated using the R project. DAZAP1 expressions were examined by quantitative RT-PCR and western blotting. CCK8 assay was used to investigate the cell proliferation, and transwell assay was employed to examine the ability of migration and invasion in vitro. Contrast-enhanced ultrasound (CEUS) was used to evaluate images and parameters of the tumor.

RESULTS

DAZAP1 is highly expressed in the tissue samples of HCC. The peak intensity (PI) and area under the curve (AUC) of the tumor is higher than that of liver parenchyma, and correlated with high DAZAP1 expression. Parameters of CEUS in the tumor are correlated with TNM stage, tumor size, and vascularity. High DAZAP1 expression correlates with a shorter survival time and advanced histologic grade (G3-G4). Bioinformatical analysis revealed that downregulation of DAZAP1 identified differentiated expressed genes (DEGs) involved in the tumor growth process.

CONCLUSIONS

DAZAP1 is highly expressed in hepatic carcinoma and related to the blood flow, and high DAZAP1 expression predicts poor prognosis. DAZAP1 may promote liver carcinoma cell proliferation, migration, and invasion of HEPG2 cells. CEUS parameters are related to the high DAZAP1 expression, and will help to differentiate the HCC tumor.

Transcriptional Basis of Psoriasis from Large Scale Gene Expression Studies: The Importance of Moving towards a Precision Medicine Approach

Transcriptome profiling techniques, such as microarrays and RNA sequencing (RNA-seq), are valuable tools for deciphering the regulatory network underlying psoriasis and have revealed large number of differentially expressed genes in lesional and non-lesional skin. Such approaches provide a more precise measurement of transcript levels and their isoforms than any other methods. Large cohort transcriptomic analyses have greatly improved our understanding of the physiological and molecular mechanisms underlying disease pathogenesis and progression. Here, we mostly review the findings of some important large scale psoriatic transcriptomic studies, and the benefits of such studies in elucidating potential therapeutic targets and biomarkers for psoriasis treatment. We also emphasised the importance of looking into the alternatively spliced RNA isoforms/transcripts in psoriasis, rather than focussing only on the gene-level annotation. The neutrophil and blood transcriptome signature in psoriasis is also briefly reviewed, as it provides the immune status information of patients and is a less invasive platform. The application of precision medicine in current management of psoriasis, by combining transcriptomic data, improves the clinical response outcome in individual patients. Drugs tailored to individual patient's genetic profile will greatly improve patient outcome and cost savings for the healthcare system.

Single-cell characterization of CRISPR-modified transcript isoforms with nanopore sequencing

We developed a single-cell approach to detect CRISPR-modified mRNA transcript structures. This method assesses how genetic variants at splicing sites and splicing factors contribute to alternative mRNA isoforms. We determine how alternative splicing is regulated by editing target exon-intron segments or splicing factors by CRISPR-Cas9 and their consequences on transcriptome profile. Our method combines long-read sequencing to characterize the transcript structure and short-read sequencing to match the single-cell gene expression profiles and gRNA sequence and therefore provides targeted genomic edits and transcript isoform structure detection at single-cell resolution.

MiR-221-3p as a Potential Biomarker for Patients with Psoriasis and Its Role in Inflammatory Responses in Keratinocytes

INTRODUCTION

This study investigated serum miR-221-3p levels in psoriatic patients and the characterization of serum miR-221-3p in keratinocyte inflammatory responses was further assessed.

METHODS

qRT-PCR was used to detect the expression level of miR-221-3p in the serum of 46 patients with psoriasis and 42 healthy controls. The receiver operating characteristic curve evaluated the diagnostic ability of miR-221-3p in psoriasis. The effect of miR-221-3p on HaCaT cell proliferation was detected by using a cell counting Kit-8 and Transwell. ELISA was used to detect serum and keratinocyte pro-inflammatory factors.

RESULTS

miR-221-3p was significantly increased in the serum of patients with psoriasis. The area under the curve was 0.861, the sensitivity was 80.4%, and the specificity was 85.7%. Serum miR-221-3p was positively correlated with the expression levels of tumor necrosis factor-α, interleukin (IL)-17A, and IL-22. Cell experiments showed that reducing the expression of miR-221-3p could significantly inhibit cell proliferation. Additionally, miR-221-3p downregulation also inhibited the release of some inflammatory factors in the HaCaT cells.

DISCUSSION/CONCLUSION

MiR-221-3p is a latent biomarker of psoriasis patients. Lower expression of miR-221-3p inhibits the cell proliferation and inflammatory responses of HaCaT cells, which offers a possible target for the therapeutic interventions of psoriasis.

Identification of protein/mRNA network involving the PSORS1 locus gene CCHCR1 and the PSORS4 locus gene HAX1

CCHCR1 (Coiled-Coil alpha-Helical Rod 1), maps to chromosomal region 6p21.3, within the major psoriasis susceptibility locus PSORS1. CCHCR1 itself is a plausible psoriasis candidate gene, however its role in psoriasis pathogenesis remains unclear. We previously demonstrated that CCHCR1 protein acts as a cytoplasmic docking site for RNA polymerase II core subunit 3 (RPB3) in cycling cells, suggesting a role for CCHCR1 in vesicular trafficking between cellular compartments. Here, we report a novel interaction between CCHCR1 and the RNA binding protein HAX1. HAX1 maps to chromosomal region 1q21.3 within the PSORS4 locus and is over-expressed in psoriasis. Both CCHCR1 and HAX1 share subcellular co-localization with mitochondria, nuclei and cytoplasmic vesicles as P-bodies. By a series of ribonucleoprotein immunoprecipitation (RIP) assays, we isolated a pool of mRNAs complexed with HAX1 and/or CCHCR1 proteins. Among the mRNAs complexed with both CCHCR1 and HAX1 proteins, there are Vimentin mRNA, previously described to be bound by HAX1, and CAMP/LL37 mRNA, whose gene product is over-expressed in psoriasis.

Integrated analysis of a compendium of RNA-Seq datasets for splicing factors

A vast amount of public RNA-sequencing datasets have been generated and used widely to study transcriptome mechanisms. These data offer precious opportunity for advancing biological research in transcriptome studies such as alternative splicing. We report the first large-scale integrated analysis of RNA-Seq data of splicing factors for systematically identifying key factors in diseases and biological processes. We analyzed 1,321 RNA-Seq libraries of various mouse tissues and cell lines, comprising more than 6.6 TB sequences from 75 independent studies that experimentally manipulated 56 splicing factors. Using these data, RNA splicing signatures and gene expression signatures were computed, and signature comparison analysis identified a list of key splicing factors in Rett syndrome and cold-induced thermogenesis. We show that cold-induced RNA-binding proteins rescue the neurite outgrowth defects in Rett syndrome using neuronal morphology analysis, and we also reveal that SRSF1 and PTBP1 are required for energy expenditure in adipocytes using metabolic flux analysis. Our study provides an integrated analysis for identifying key factors in diseases and biological processes and highlights the importance of public data resources for identifying hypotheses for experimental testing.

Review-Current Concepts in Inflammatory Skin Diseases Evolved by Transcriptome Analysis: In-Depth Analysis of Atopic Dermatitis and Psoriasis

During the last decades, high-throughput assessment of gene expression in patient tissues using microarray technology or RNA-Seq took center stage in clinical research. Insights into the diversity and frequency of transcripts in healthy and diseased conditions provide valuable information on the cellular status in the respective tissues. Growing with the technique, the bioinformatic analysis toolkit reveals biologically relevant pathways which assist in understanding basic pathophysiological mechanisms. Conventional classification systems of inflammatory skin diseases rely on descriptive assessments by pathologists. In contrast to this, molecular profiling may uncover previously unknown disease classifying features. Thereby, treatments and prognostics of patients may be improved. Furthermore, disease models in basic research in comparison to the human disease can be directly validated. The aim of this article is not only to provide the reader with information on the opportunities of these techniques, but to outline potential pitfalls and technical limitations as well. Major published findings are briefly discussed to provide a broad overview on the current findings in transcriptomics in inflammatory skin diseases.

A developmentally regulated spliced variant of PTBP1 is upregulated in type 1 diabetic hearts

Alternative splicing (AS) is dysregulated in Type 1 diabetic (T1D) hearts but mechanisms responsible are unclear. Here, we provide evidence that the RNA binding protein (RBP) PTBP1 is modulated in adult T1D hearts contributing to AS changes. We show that a spliced variant of PTBP1 that is highly expressed in normal newborn mouse hearts is aberrantly expressed in adult T1D mouse hearts. Comparing known PTBP1-target datasets to our T1D mouse transcriptome datasets, we discovered a group of genes with PTBP1 binding sites in their pre-mRNAs that are differentially spliced in T1D mouse hearts. We demonstrated that inducible expression of diabetes-induced PTBP1 spliced variant has less repressive splicing function. Notably, PTBP1 regulates AS of some of its targets antagonistically to RBFOX2. In sum, our results indicate that diabetic conditions disrupt developmental regulation of PTBP1 leading to differential AS of PTBP1 target genes. Identification of PTBP1 and PTBP1-regulated RNA networks can provide RNA-based therapies for the treatment of diabetes cardiac complications.

Advances in Analyzing Virus-Induced Alterations of Host Cell Splicing

Alteration of host cell splicing is a common feature of many viral infections which is underappreciated because of the complexity and technical difficulty of studying alternative splicing (AS) regulation. Recent advances in RNA sequencing technologies revealed that up to several hundreds of host genes can show altered mRNA splicing upon viral infection. The observed changes in AS events can be either a direct consequence of viral manipulation of the host splicing machinery or result indirectly from the virus-induced innate immune response or cellular damage. Analysis at a higher resolution with single-cell RNAseq, and at a higher scale with the integration of multiple omics data sets in a systems biology perspective, will be needed to further comprehend this complex facet of virus-host interactions.

RBPMetaDB: a comprehensive annotation of mouse RNA-Seq datasets with perturbations of RNA-binding proteins

RNA-binding proteins (RBPs) may play a critical role in gene regulation in various diseases or biological processes by controlling post-transcriptional events such as polyadenylation, splicing and mRNA stabilization via binding activities to RNA molecules. Owing to the importance of RBPs in gene regulation, a great number of studies have been conducted, resulting in a large amount of RNA-Seq datasets. However, these datasets usually do not have structured organization of metadata, which limits their potentially wide use. To bridge this gap, the metadata of a comprehensive set of publicly available mouse RNA-Seq datasets with perturbed RBPs were collected and integrated into a database called RBPMetaDB. This database contains 292 mouse RNA-Seq datasets for a comprehensive list of 187 RBPs. These RBPs account for only ∼10% of all known RBPs annotated in Gene Ontology, indicating that most are still unexplored using high-throughput sequencing. This negative information provides a great pool of candidate RBPs for biologists to conduct future experimental studies. In addition, we found that DNA-binding activities are significantly enriched among RBPs in RBPMetaDB, suggesting that prior studies of these DNA- and RNA-binding factors focus more on DNA-binding activities instead of RNA-binding activities. This result reveals the opportunity to efficiently reuse these data for investigation of the roles of their RNA-binding activities. A web application has also been implemented to enable easy access and wide use of RBPMetaDB. It is expected that RBPMetaDB will be a great resource for improving understanding of the biological roles of RBPs.

Database URL: http://rbpmetadb.yubiolab.org

CELF1 contributes to aberrant alternative splicing patterns in the type 1 diabetic heart

Dysregulated alternative splicing (AS) that contributes to diabetes pathogenesis has been identified, but little is known about the RNA binding proteins (RBPs) involved. We have previously found that the RBP CELF1 is upregulated in the diabetic heart; however, it is unclear if CELF1 contributes to diabetes-induced AS changes. Utilizing genome wide approaches, we identified extensive changes in AS patterns in Type 1 diabetic (T1D) mouse hearts. We discovered that many aberrantly spliced genes in T1D hearts have CELF1 binding sites. CELF1-regulated AS affects key genes within signaling pathways relevant to diabetes pathogenesis. Disruption of CELF1 binding sites impairs AS regulation by CELF1. In sum, our results indicate that CELF1 target RNAs are aberrantly spliced in the T1D heart leading to abnormal gene expression. These discoveries pave the way for targeting RBPs and their RNA networks as novel therapies for cardiac complications of diabetes.

Specification of Drosophila neuropeptidergic neurons by the splicing component brr2

During embryonic development, a number of genetic cues act to generate neuronal diversity. While intrinsic transcriptional cascades are well-known to control neuronal sub-type cell fate, the target cells can also provide critical input to specific neuronal cell fates. Such signals, denoted retrograde signals, are known to provide critical survival cues for neurons, but have also been found to trigger terminal differentiation of neurons. One salient example of such target-derived instructive signals pertains to the specification of the Drosophila FMRFamide neuropeptide neurons, the Tv4 neurons of the ventral nerve cord. Tv4 neurons receive a BMP signal from their target cells, which acts as the final trigger to activate the FMRFa gene. A recent FMRFa-eGFP genetic screen identified several genes involved in Tv4 specification, two of which encode components of the U5 subunit of the spliceosome: brr2 (l(3)72Ab) and Prp8. In this study, we focus on the role of RNA processing during target-derived signaling. We found that brr2 and Prp8 play crucial roles in controlling the expression of the FMRFa neuropeptide specifically in six neurons of the VNC (Tv4 neurons). Detailed analysis of brr2 revealed that this control is executed by two independent mechanisms, both of which are required for the activation of the BMP retrograde signaling pathway in Tv4 neurons: (1) Proper axonal pathfinding to the target tissue in order to receive the BMP ligand. (2) Proper RNA splicing of two genes in the BMP pathway: the thickveins (tkv) gene, encoding a BMP receptor subunit, and the Medea gene, encoding a co-Smad. These results reveal involvement of specific RNA processing in diversifying neuronal identity within the central nervous system.

The nervous system displays daunting cellular diversity, largely generated through complex regulatory input operating on stem cells and their neural lineages during development. Most of the reported mechanisms acting to generate neural diversity pertain to transcriptional regulation. In contrast, little is known regarding the post-transcriptional mechanisms involved. Here, we use a specific group of neurons, Apterous neurons, in the ventral nerve cord of Drosophila melanogaster as our model, to analyze the function of two essential components of the spliceosome; Brr2 and Prp8. Apterous neurons require a BMP retrograde signal for terminal differentiation, and we find that brr2 and Prp8 play crucial roles during this process. brr2 is critical for two independent events; axon pathfinding and BMP signaling, both of which are required for the activation of the retrograde signaling pathway necessary for Apterous neurons. These results identify a post-transcriptional mechanism as key for specifying neuronal identity, by ensuring the execution of a retrograde signal.

A data mining paradigm for identifying key factors in biological processes using gene expression data

A large volume of biological data is being generated for studying mechanisms of various biological processes. These precious data enable large-scale computational analyses to gain biological insights. However, it remains a challenge to mine the data efficiently for knowledge discovery. The heterogeneity of these data makes it difficult to consistently integrate them, slowing down the process of biological discovery. We introduce a data processing paradigm to identify key factors in biological processes via systematic collection of gene expression datasets, primary analysis of data, and evaluation of consistent signals. To demonstrate its effectiveness, our paradigm was applied to epidermal development and identified many genes that play a potential role in this process. Besides the known epidermal development genes, a substantial proportion of the identified genes are still not supported by gain- or loss-of-function studies, yielding many novel genes for future studies. Among them, we selected a top gene for loss-of-function experimental validation and confirmed its function in epidermal differentiation, proving the ability of this paradigm to identify new factors in biological processes. In addition, this paradigm revealed many key genes in cold-induced thermogenesis using data from cold-challenged tissues, demonstrating its generalizability. This paradigm can lead to fruitful results for studying molecular mechanisms in an era of explosive accumulation of publicly available biological data.

Activity-dependent aberrations in gene expression and alternative splicing in a mouse model of Rett syndrome

Rett syndrome (RTT) is a severe neurodevelopmental disorder that affects about 1 in 10,000 female live births. The underlying cause of RTT is mutations in the X-linked gene, methyl-CpG-binding protein 2 (MECP2); however, the molecular mechanism by which these mutations mediate the RTT neuropathology remains enigmatic. Specifically, although MeCP2 is known to act as a transcriptional repressor, analyses of the RTT brain at steady-state conditions detected numerous differentially expressed genes, while the changes in transcript levels were mostly subtle. Here we reveal an aberrant global pattern of gene expression, characterized predominantly by higher levels of expression of activity-dependent genes, and anomalous alternative splicing events, specifically in response to neuronal activity in a mouse model for RTT. Notably, the specific splicing modalities of intron retention and exon skipping displayed a significant bias toward increased retained introns and skipped exons, respectively, in the RTT brain compared with the WT brain. Furthermore, these aberrations occur in conjunction with higher seizure susceptibility in response to neuronal activity in RTT mice. Our findings advance the concept that normal MeCP2 functioning is required for fine-tuning the robust and immediate changes in gene transcription and for proper regulation of alternative splicing induced in response to neuronal stimulation.