Defective structural RNA processing in relapsing-remitting multiple sclerosis

The Osgin Gene Family: Underexplored Yet Essential Mediators of Oxidative Stress

The Osgin gene family consists of two members, Osgin1 and Osgin2, involved in the cellular oxidative stress response. While many members of this essential cellular pathway have been extensively characterized, the Osgin gene family, despite its broad phylogenetic distribution, has received far less attention. Here, we review published articles and open-source databases to synthesize the current research on the evolutionary history, structure, biochemical and physiological functions, expression patterns, and role in disease of the Osgin gene family. Although Osgin displays broad spatiotemporal expression during development and adulthood, there is ambiguity regarding the cellular functions of the OSGIN proteins. A recent study identified OSGIN-1 as a flavin-dependent monooxygenase, but the biochemical role of OSGIN-2 has not yet been defined. Moreover, while the Osgin genes are implicated as mediators of cell proliferation, apoptosis, and autophagy, these functions have not been connected to the enzymatic classification of OSGIN. Misregulation of Osgin expression has long been associated with various disease states, yet recent analyses highlight the mechanistic role of OSGIN in pathogenesis and disease progression, underscoring the therapeutic potential of targeting OSGIN. In light of these findings, we suggest further avenues of research to advance our understanding of this essential, yet underexplored, gene family.

Defective Integrator activity shapes the transcriptome of patients with multiple sclerosis

HP1α/CBX5 is an epigenetic regulator with a suspected role in multiple sclerosis (MS). Here, using high-depth RNA sequencing on monocytes, we identified a subset of MS patients with reduced CBX5 expression, correlating with progressive stages of the disease and extensive transcriptomic alterations. Examination of rare non-coding RNA species in these patients revealed impaired maturation/degradation of U snRNAs and enhancer RNAs, indicative of reduced activity of the Integrator, a complex with suspected links to increased MS risk. At protein-coding genes, compromised Integrator activity manifested in reduced pre-mRNA splicing efficiency and altered expression of genes regulated by RNA polymerase II pause-release. Inactivation of Cbx5 in the mouse mirrored most of these transcriptional defects and resulted in hypersensitivity to experimental autoimmune encephalomyelitis. Collectively, our observations suggested a major contribution of the Integrator complex in safeguarding against transcriptional anomalies characteristic of MS, with HP1α/CBX5 emerging as an unexpected regulator of this complex's activity. These findings bring novel insights into the transcriptional aspects of MS and provide potential new criteria for patient stratification.

Ro60-Roles in RNA Processing, Inflammation, and Rheumatic Autoimmune Diseases

The Ro60/SSA2 autoantigen is an RNA-binding protein and a core component of nucleocytoplasmic ribonucleoprotein (RNP) complexes. Ro60 is essential in RNA metabolism, cell stress response pathways, and cellular homeostasis. It stabilises and mediates the quality control and cellular distribution of small RNAs, including YRNAs (for the 'y' in 'cytoplasmic'), retroelement transcripts, and misfolded RNAs. Ro60 transcriptional dysregulation or loss of function can result in the generation and release of RNA fragments from YRNAs and other small RNAs. Small RNA fragments can instigate an inflammatory cascade through endosomal toll-like receptors (TLRs) and cytoplasmic RNA sensors, which typically sense pathogen-associated molecular patterns, and mount the first line of defence against invading pathogens. However, the recognition of host-originating RNA moieties from Ro60 RNP complexes can activate inflammatory response pathways and compromise self-tolerance. Autoreactive B cells may produce antibodies targeting extracellular Ro60 RNP complexes. Ro60 autoantibodies serve as diagnostic markers for various autoimmune diseases, including Sjögren's disease (SjD) and systemic lupus erythematosus (SLE), and they may also act as predictive markers for anti-drug antibody responses among rheumatic patients. Understanding Ro60's structure, function, and role in self-tolerance can enhance our understanding of the underlying molecular mechanisms of autoimmune conditions.

Transcriptome driven discovery of novel candidate genes for human neurological disorders in the telomer-to-telomer genome assembly era

BACKGROUND

With the first complete draft of a human genome, the Telomere-to-Telomere Consortium unlocked previously concealed genomic regions for genetic analyses. These regions harbour nearly 2000 potential novel genes with unknown function. In order to uncover candidate genes associated with human neurological pathologies, a comparative transcriptome study using the T2T-CHM13 and the GRCh38 genome assemblies was conducted on previously published datasets for eight distinct human neurological disorders.

RESULTS

The analysis of differential expression in RNA sequencing data led to the identification of 336 novel candidate genes linked to human neurological disorders. Additionally, it was revealed that, on average, 3.6% of the differentially expressed genes detected with the GRCh38 assembly may represent potential false positives. Among the noteworthy findings, two novel genes were discovered, one encoding a pore-structured protein and the other a highly ordered β-strand-rich protein. These genes exhibited upregulation in multiple epilepsy datasets and hold promise as candidate genes potentially modulating the progression of the disease. Furthermore, an analysis of RNA derived from white matter lesions in multiple sclerosis patients indicated significant upregulation of 26 rRNA encoding genes. Additionally, putative pathology related genes were identified for Alzheimer's disease, amyotrophic lateral sclerosis, glioblastoma, glioma, and conditions resulting from the m.3242 A > G mtDNA mutation.

CONCLUSION

The results presented here underline the potential of the T2T-CHM13 assembly in facilitating the discovery of candidate genes from transcriptome data in the context of human disorders. Moreover, the results demonstrate the value of remapping sequencing data to a superior genome assembly. Numerous potential pathology related genes, either as causative factors or related elements, have been unveiled, warranting further experimental validation.

The circular RNA landscape in multiple sclerosis: Disease-specific associated variants and exon methylation shape circular RNA expression profile

BACKGROUND

Circular RNAs (circRNAs) are a class of non-coding RNAs increasingly emerging as crucial actors in the pathogenesis of human diseases, including autoimmune and neurological disorders as multiple sclerosis (MS). Despite several efforts, the mechanisms regulating circRNAs expression are still largely unknown and the circRNA profile and regulation in MS-relevant cell models has not been completely investigated. In this work, we aimed at exploring the global landscape of circRNA expression in MS patients, also evaluating a possible correlation with their genetic and epigenetic background.

METHODS

We performed RNA-seq experiments on circRNA-enriched samples, derived from peripheral blood mononuclear cells (PBMCs) of 10 MS patients and 10 matched controls and performed differential circRNA expression. The genetic background was evaluated using array genotyping, and an expression quantitative trait loci (eQTL) analysis was carried out.

RESULTS

Expression analysis revealed 166 differentially expressed circRNAs in MS patients, 125 of which are downregulated. One of the top dysregulated circRNAs, hsa_circ_0007990, derives from the PGAP3 gene, encoding a protein relevant for the control of autoimmune responses. The downregulation of this circRNA was confirmed in two independent replication cohorts, suggesting its implementation as a possible RNA-based biomarker. The eQTL analysis evidenced a significant association between 89 MS-associated loci and the expression of at least one circRNA, suggesting that MS-associated variants could impact on disease pathogenesis by altering circRNA profiles. Finally, we found a significant correlation between exon methylation and circRNA expression levels, supporting the hypothesis that epigenetic features may play an important role in the definition of the cell circRNA pool.

CONCLUSION

We described the circRNA expression profile of PBMCs in MS patients, suggesting that MS-associated variants may tune the expression levels of circRNAs acting as "circ-QTLs", and proposing a role for exon-based DNA methylation in regulating circRNA expression.

A transcriptome atlas and interactive analysis platform for autoimmune disease

With the rapid development of next-generation sequencing technology, many laboratories have produced a large amount of single-cell transcriptome data of blood and tissue samples from patients with autoimmune diseases, which enables in-depth studies of the relationship between gene transcription and autoimmune diseases. However, there is still a lack of a database that integrates the large amount of autoimmune disease transcriptome sequencing data and conducts effective analysis. In this study, we developed a user-friendly web database tool, Interactive Analysis and Atlas for Autoimmune disease (IAAA), which integrates bulk RNA-seq data of 929 samples of 10 autoimmune diseases and single-cell RNA-seq data of 783 203 cells in 96 samples of 6 autoimmune diseases. IAAA also provides customizable analysis modules, including gene expression, difference, correlation, similar gene detection and cell–cell interaction, and can display results in three formats (plot, table and pdf) through custom parameters. IAAA provides valuable data resources for researchers studying autoimmune diseases and helps users deeply explore the potential value of the current transcriptome data. IAAA is available.

Database URL: http://galaxy.ustc.edu.cn/IAAA

Application of Genetic Algorithm-Based Support Vector Machine in Identification of Gene Expression Signatures for Psoriasis Classification: A Hybrid Model

BACKGROUND

Psoriasis is a chronic autoimmune disease impairing significantly the quality of life of the patient. The diagnosis of the disease is done via a visual inspection of the lesional skin by dermatologists. Classification of psoriasis using gene expression is an important issue for the early and effective treatment of the disease. Therefore, gene expression data and selection of suitable gene signatures are effective sources of information.

METHODS

We aimed to develop a hybrid classifier for the diagnosis of psoriasis based on two machine learning models of the genetic algorithm and support vector machine (SVM). The method also conducts gene signature selection. A publically available gene expression dataset was used to test the model.

RESULTS

A number of 181 probe sets were selected among the original 54,675 probes using the hybrid model with a prediction accuracy of 100% over the test set. A number of 10 hub genes were identified using the protein-protein interaction network. Nine out of 10 identified genes were found in significant modules.

CONCLUSIONS

The results showed that the genetic algorithm improved the SVM classifier performance significantly implying the ability of the proposed model in terms of detecting relevant gene expression signatures as the best features.

Small noncoding RNA profiling across cellular and biofluid compartments and their implications for multiple sclerosis immunopathology

Dysregulation of microRNAs (miRNAs), a type of small noncoding RNAs (sncRNAs), has frequently been associated with multiple sclerosis (MS). However, most studies have focused on peripheral blood, and few investigated other classes of sncRNAs. To address this, we analyzed all classes of sncRNAs in matching peripheral blood mononuclear cells, plasma, cerebrospinal fluid (CSF) cells, and cell-free CSF from MS patients and controls. We demonstrate widespread alterations of small nuclear (snRNA)–derived RNAs, small nucleolar-derived RNAs (sdRNAs), transfer RNA–derived fragments, and miRNAs, particularly in CSF cells. The striking contrast between the periphery and central nervous system and between relapse and remission phases of disease highlights the importance of sncRNA-mediated mechanisms in MS, in particular alternative splicing and mRNA translation.

Multiple sclerosis (MS) is a chronic inflammatory demyelinating disease affecting the central nervous system (CNS). Small non-coding RNAs (sncRNAs) and, in particular, microRNAs (miRNAs) have frequently been associated with MS. Here, we performed a comprehensive analysis of all classes of sncRNAs in matching samples of peripheral blood mononuclear cells (PBMCs), plasma, cerebrospinal fluid (CSF) cells, and cell-free CSF from relapsing-remitting (RRMS, n = 12 in relapse and n = 11 in remission) patients, secondary progressive (SPMS, n = 6) MS patients, and noninflammatory and inflammatory neurological disease controls (NINDC, n = 11; INDC, n = 5). We show widespread changes in miRNAs and sncRNA-derived fragments of small nuclear, nucleolar, and transfer RNAs. In CSF cells, 133 out of 133 and 115 out of 117 differentially expressed sncRNAs were increased in RRMS relapse compared to remission and RRMS compared to NINDC, respectively. In contrast, 65 out of 67 differentially expressed PBMC sncRNAs were decreased in RRMS compared to NINDC. The striking contrast between the periphery and CNS suggests that sncRNA-mediated mechanisms, including alternative splicing, RNA degradation, and mRNA translation, regulate the transcriptome of pathogenic cells primarily in the CNS target organ.

Disease severity-specific neutrophil signatures in blood transcriptomes stratify COVID-19 patients

BACKGROUND

The SARS-CoV-2 pandemic is currently leading to increasing numbers of COVID-19 patients all over the world. Clinical presentations range from asymptomatic, mild respiratory tract infection, to severe cases with acute respiratory distress syndrome, respiratory failure, and death. Reports on a dysregulated immune system in the severe cases call for a better characterization and understanding of the changes in the immune system.

METHODS

In order to dissect COVID-19-driven immune host responses, we performed RNA-seq of whole blood cell transcriptomes and granulocyte preparations from mild and severe COVID-19 patients and analyzed the data using a combination of conventional and data-driven co-expression analysis. Additionally, publicly available data was used to show the distinction from COVID-19 to other diseases. Reverse drug target prediction was used to identify known or novel drug candidates based on finding from data-driven findings.

RESULTS

Here, we profiled whole blood transcriptomes of 39 COVID-19 patients and 10 control donors enabling a data-driven stratification based on molecular phenotype. Neutrophil activation-associated signatures were prominently enriched in severe patient groups, which was corroborated in whole blood transcriptomes from an independent second cohort of 30 as well as in granulocyte samples from a third cohort of 16 COVID-19 patients (44 samples). Comparison of COVID-19 blood transcriptomes with those of a collection of over 3100 samples derived from 12 different viral infections, inflammatory diseases, and independent control samples revealed highly specific transcriptome signatures for COVID-19. Further, stratified transcriptomes predicted patient subgroup-specific drug candidates targeting the dysregulated systemic immune response of the host.

CONCLUSIONS

Our study provides novel insights in the distinct molecular subgroups or phenotypes that are not simply explained by clinical parameters. We show that whole blood transcriptomes are extremely informative for COVID-19 since they capture granulocytes which are major drivers of disease severity.

Disease severity-specific neutrophil signatures in blood transcriptomes stratify COVID-19 patients

BACKGROUND

The SARS-CoV-2 pandemic is currently leading to increasing numbers of COVID-19 patients all over the world. Clinical presentations range from asymptomatic, mild respiratory tract infection, to severe cases with acute respiratory distress syndrome, respiratory failure, and death. Reports on a dysregulated immune system in the severe cases call for a better characterization and understanding of the changes in the immune system.

METHODS

In order to dissect COVID-19-driven immune host responses, we performed RNA-seq of whole blood cell transcriptomes and granulocyte preparations from mild and severe COVID-19 patients and analyzed the data using a combination of conventional and data-driven co-expression analysis. Additionally, publicly available data was used to show the distinction from COVID-19 to other diseases. Reverse drug target prediction was used to identify known or novel drug candidates based on finding from data-driven findings.

RESULTS

Here, we profiled whole blood transcriptomes of 39 COVID-19 patients and 10 control donors enabling a data-driven stratification based on molecular phenotype. Neutrophil activation-associated signatures were prominently enriched in severe patient groups, which was corroborated in whole blood transcriptomes from an independent second cohort of 30 as well as in granulocyte samples from a third cohort of 16 COVID-19 patients (44 samples). Comparison of COVID-19 blood transcriptomes with those of a collection of over 3100 samples derived from 12 different viral infections, inflammatory diseases, and independent control samples revealed highly specific transcriptome signatures for COVID-19. Further, stratified transcriptomes predicted patient subgroup-specific drug candidates targeting the dysregulated systemic immune response of the host.

CONCLUSIONS

Our study provides novel insights in the distinct molecular subgroups or phenotypes that are not simply explained by clinical parameters. We show that whole blood transcriptomes are extremely informative for COVID-19 since they capture granulocytes which are major drivers of disease severity.

Adenosine-to-Inosine RNA Editing of Alu Double-Stranded (ds)RNAs Is Markedly Decreased in Multiple Sclerosis and Unedited Alu dsRNAs Are Potent Activators of Proinflammatory Transcriptional Responses

Sensors that detect dsRNA stimulate IFN responses as a defense against viral infection. IFN responses are also well documented in a variety of human autoimmune diseases, including relapsing-remitting multiple sclerosis (MS), in which increased IFN responses result from increased levels of double-stranded endogenous Alu RNAs. Mechanisms underlying increases in double-stranded Alu RNAs in MS are obscure. We find widespread loss of adenosine-to-inosine editing of Alu RNAs in MS. Unedited Alu RNAs are potent activators of both IFN and NF-κB responses via the dsRNA sensors, RIG-I, and TLR3. Minor editing of highly active Alu elements abrogates the ability to activate both transcriptional responses. Thus, adenosine-to-inosine editing may also represent an important defense against autoimmune diseases such as MS.

Not only cancer: the long non-coding RNA MALAT1 affects the repertoire of alternatively spliced transcripts and circular RNAs in multiple sclerosis

Long non-coding RNAs (lncRNAs) are post-transcriptional and epigenetic regulators, whose implication in neurodegenerative and autoimmune diseases remains poorly understood. We analyzed publicly available microarray data sets to identify dysregulated lncRNAs in multiple sclerosis (MS), a neuroinflammatory autoimmune disease. We found a consistent upregulation in MS of the lncRNA MALAT1 (2.7-fold increase; meta-analysis, P = 1.3 × 10-8; 190 cases, 182 controls), known to regulate alternative splicing (AS). We confirmed MALAT1 upregulation in two independent MS cohorts (1.5-fold increase; P < 0.01; 59 cases, 50 controls). We hence performed MALAT1 overexpression/knockdown in cell lines, demonstrating that its modulation impacts on endogenous expression of splicing factors (HNRNPF and HNRNPH1) and on AS of MS-associated genes (IL7R and SP140). Minigene-based splicing assays upon MALAT1 modulation recapitulated IL7R and SP140 isoform unbalances observed in patients. RNA-sequencing of MALAT1-knockdown Jurkat cells further highlighted MALAT1 role in splicing (approximately 1100 significantly-modulated AS events) and revealed its contribution to backsplicing (approximately 50 differentially expressed circular RNAs). Our study proposes a possible novel role for MALAT1 dysregulation and the consequent AS alteration in MS pathogenesis, based on anomalous splicing/backsplicing profiles of MS-relevant genes.

Aberrant expression of alternative splicing variants in multiple sclerosis - A systematic review

OBJECTIVE

Alternative splicing is an important form of RNA processing that affects nearly all human genes. The differential expression of specific transcript and protein isoforms holds the potential of novel biomarkers for complex diseases. In this systematic review, we compiled the existing literature on aberrant alternative splicing events in multiple sclerosis (MS).

METHODS

A systematic literature search in the PubMed database was carried out and supplemented by screening the reference lists of the identified articles. We selected only MS-related original research studies which compared the levels of different isoforms of human protein-coding genes. A narrative synthesis of the research findings was conducted. Additionally, we performed a case-control analysis using high-density transcriptome microarray data to reevaluate the genes that were examined in the reviewed studies.

RESULTS

A total of 160 records were screened. Of those, 36 studies from the last two decades were included. Most commonly, peripheral blood samples were analyzed (32 studies), and PCR-based techniques were usually employed (27 studies) for measuring the expression of selected genes. Two studies used an exploratory genome-wide approach. Overall, 27 alternatively spliced genes were investigated. Nine of these genes appeared in at least two studies (CD40, CFLAR, FOXP3, IFNAR2, IL7R, MOG, PTPRC, SP140 and TNFRSF1A). The microarray data analysis confirmed differential alternative pre-mRNA splicing for 19 genes.

CONCLUSIONS

An altered RNA processing of genes mediating immune signaling pathways has been repeatedly implicated in MS. The analysis of individual exon-level expression patterns is stimulated by the advancement of transcriptome profiling technologies. In particular, the examination of genes encoded in MS-associated genetic regions may provide important insights into the pathogenesis of the disease and help to identify new biomarkers.

Endogenous double-stranded Alu RNA elements stimulate IFN-responses in relapsing remitting multiple sclerosis

Various sensors that detect double-stranded RNA, presumably of viral origin, exist in eukaryotic cells and induce IFN-responses. Ongoing IFN-responses have also been documented in a variety of human autoimmune diseases including relapsing-remitting multiple sclerosis (RRMS) but their origins remain obscure. We find increased IFN-responses in leukocytes in relapsing-remitting multiple sclerosis at distinct stages of disease. Moreover, endogenous RNAs isolated from blood cells of these same patients recapitulate this IFN-response if transfected into naïve cells. These endogenous RNAs are double-stranded RNAs, contain Alu and Line elements and are transcribed from leukocyte transcriptional enhancers. Thus, transcribed endogenous retrotransposon elements can co-opt pattern recognition sensors to induce IFN-responses in RRMS.

Longitudinal changes in the expression of IL-33 and IL-33 regulated genes in relapsing remitting MS

Objective

We tested the hypothesis that the expression of IL-33 in MS is dynamic and is likely to reflect the clinical and radiological changes during the course of RRMS.

Methods

MS with either clinical or radiological relapses were recruited for the study and followed for one year. IL-33 and a panel of genes was measured by q PCR and flow cytometry at different time points.

Results

Among 22 RRMS patients, 4 patients showed highest levels of IL-33 at the time they were recruited to the study (Month 0); in 14 patients highest levels of IL-33 were seen between 6–11 months after relapse and in 4 patients maximal levels of IL-33 were seen 12 months after relapse. A similar pattern of IL-33 kinetics was seen when IL-33 was measured by flow cytometry in an additional cohort of 12 patients. The timing of the improvement clinically did not correlate with IL-33 expression with highest expression levels either preceding or following clinical recovery. From our whole genome RNA-sequencing data we found a strong correlation between expression levels of IL-33 and a ~2000 mRNA genes. However, none of these genes encoded proteins involved in either innate or adaptive immunity. Rather, many of the genes that correlated highly with IL-33 encoded to proteins involved in DNA repair or mitochondrial function and mRNA splicing pathways.

Interpretation

Given the neuro-reparative and remodeling functions attributed to IL-33, it is likely that some of the novel genes we have uncovered may be involved in repair and recovery of the CNS in MS.

Genetic Variant rs755622 Regulates Expression of the Multiple Sclerosis Severity Modifier D-Dopachrome Tautomerase in a Sex-Specific Way

Multiple sclerosis (MS) is a sex-specific autoimmune disease involving central nervous system. Previous studies determined that macrophage migration inhibitory factor (MIF) and its homologue D-dopachrome tautomerase (DDT) sex-specifically affect MS progression. Moreover, other studies reported that rs755622 polymorphism in promoter region of MIF gene is associated with risk of MS and affects the promoter activity to regulate MIF expression in a sex-specific way. Given that MIF and DDT share a part of promoter sequence, we surmise that rs755622 can also regulate DDT expression in a sex-specific way. However, this has not yet been studied. Here, we used five large-scale expression quantitative trait loci (eQTLs) and two RNA-seq datasets from brain and blood to assess the potential influence of rs755622 variant on expression of DDT in different genders by the linear regression and differential expression analysis. The results show that the minor allele frequency of rs755622 and expression of DDT are significantly increased in males for MS subjects and this minor allele variant can significantly upregulate DDT expression for males but not females, which suggests that the regulation of DDT expression level by rs755622 can affect MS progression in males. These findings further support and expand conclusions of previous studies and may help to better understand the mechanisms of MS.

Activation of Peptidylarginine Deiminase in the Salivary Glands of Balb/c Mice Drives the Citrullination of Ro and La Ribonucleoproteins

The goal of the present study was to determine whether peptidylarginine deiminase PAD2 and PAD4 enzymes are present in Balb/c mouse salivary glands and whether they are able to citrullinate Ro and La ribonucleoproteins. Salivary glands from Balb/c mice were cultured in DMEM and supplemented with one of the following stimulants: ATP, LPS, TNF, IFNγ, or IL-6. A control group without stimulant was also evaluated. PAD2, PAD4, citrullinated peptides, Ro60, and La were detected by immunohistochemistry and double immunofluorescence. PAD2 and PAD4 mRNAs and protein expression were detected by qPCR and Western blot analysis. PAD activity was assessed using an antigen capture enzyme-linked immunosorbent assay. LPS, ATP, and TNF triggered PAD2 and PAD4 expression; in contrast, no expression was detected in the control group (p < 0.001). PAD transcription slightly increased in response to stimulation. Additionally, PAD2/4 activity modified the arginine residues of a reporter protein (fibrinogen) in vitro. PADs citrullinated Ro60 and La ribonucleoproteins in vivo. Molecular stimulants induced apoptosis in ductal cells and the externalization of Ro60 and La ribonucleoproteins onto apoptotic membranes. PAD enzymes citrullinate Ro and La ribonucleoproteins, and this experimental approach may facilitate our understanding of the role of posttranslational modifications in the pathophysiology of Sjögren's syndrome.

The Characterization of GSDMB Splicing and Backsplicing Profiles Identifies Novel Isoforms and a Circular RNA That Are Dysregulated in Multiple Sclerosis

Abnormalities in alternative splicing (AS) are emerging as recurrent features in autoimmune diseases (AIDs). In particular, a growing body of evidence suggests the existence of a pathogenic association between a generalized defect in splicing regulatory genes and multiple sclerosis (MS). Moreover, several studies have documented an unbalance in alternatively-spliced isoforms in MS patients possibly contributing to the disease etiology. In this work, using a combination of PCR-based techniques (reverse-transcription (RT)-PCR, fluorescent-competitive, real-time, and digital RT-PCR assays), we investigated the alternatively-spliced gene encoding Gasdermin B, GSDMB, which was repeatedly associated with susceptibility to asthma and AIDs. The in-depth characterization of GSDMB AS and backsplicing profiles led us to the identification of an exonic circular RNA (ecircRNA) as well as of novel GSDMB in-frame and out-of-frame isoforms. The non-productive splicing variants were shown to be downregulated by the nonsense-mediated mRNA decay (NMD) in human cell lines, suggesting that GSDMB levels are significantly modulated by NMD. Importantly, both AS isoforms and the identified ecircRNA were significantly dysregulated in peripheral blood mononuclear cells of relapsing-remitting MS patients compared to controls, further supporting the notion that aberrant RNA metabolism is a characteristic feature of the disease.

A fuzzy method for RNA-Seq differential expression analysis in presence of multireads

Background

When the reads obtained from high-throughput RNA sequencing are mapped against a reference database, a significant proportion of them - known as multireads - can map to more than one reference sequence. These multireads originate from gene duplications, repetitive regions or overlapping genes. Removing the multireads from the mapping results, in RNA-Seq analyses, causes an underestimation of the read counts, while estimating the real read count can lead to false positives during the detection of differentially expressed sequences.

Results

We present an innovative approach to deal with multireads and evaluate differential expression events, entirely based on fuzzy set theory. Since multireads cause uncertainty in the estimation of read counts during gene expression computation, they can also influence the reliability of differential expression analysis results, by producing false positives. Our method manages the uncertainty in gene expression estimation by defining the fuzzy read counts and evaluates the possibility of a gene to be differentially expressed with three fuzzy concepts: over-expression, same-expression and under-expression. The output of the method is a list of differentially expressed genes enriched with information about the uncertainty of the results due to the multiread presence.

We have tested the method on RNA-Seq data designed for case-control studies and we have compared the obtained results with other existing tools for read count estimation and differential expression analysis.

Conclusions

The management of multireads with the use of fuzzy sets allows to obtain a list of differential expression events which takes in account the uncertainty in the results caused by the presence of multireads. Such additional information can be used by the biologists when they have to select the most relevant differential expression events to validate with laboratory assays. Our method can be used to compute reliable differential expression events and to highlight possible false positives in the lists of differentially expressed genes computed with other tools.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1195-2) contains supplementary material, which is available to authorized users.