Long RNA Sequencing and Ribosome Profiling of Inflamed β-Cells Reveal an Extensive Translatome Landscape

Comprehensive assessment of mRNA isoform detection methods for long-read sequencing data

The advancement of Long-Read Sequencing (LRS) techniques has significantly increased the length of sequencing to several kilobases, thereby facilitating the identification of alternative splicing events and isoform expressions. Recently, numerous computational tools for isoform detection using long-read sequencing data have been developed. Nevertheless, there remains a deficiency in comparative studies that systemically evaluate the performance of these tools, which are implemented with different algorithms, under various simulations that encompass potential influencing factors. In this study, we conducted a benchmark analysis of thirteen methods implemented in nine tools capable of identifying isoform structures from long-read RNA-seq data. We evaluated their performances using simulated data, which represented diverse sequencing platforms generated by an in-house simulator, RNA sequins (sequencing spike-ins) data, as well as experimental data. Our findings demonstrate IsoQuant as a highly effective tool for isoform detection with LRS, with Bambu and StringTie2 also exhibiting strong performance. These results offer valuable guidance for future research on alternative splicing analysis and the ongoing improvement of tools for isoform detection using LRS data.

Pharmacological inhibition of tyrosine protein-kinase 2 reduces islet inflammation and delays type 1 diabetes onset in mice

Tyrosine protein-kinase 2 (TYK2), a member of the Janus kinase family, mediates inflammatory signaling through multiple cytokines, including interferon-α (IFNα), interleukin (IL)-12, and IL-23. Missense mutations in TYK2 are associated with protection against type 1 diabetes (T1D), and inhibition of TYK2 shows promise in the management of other autoimmune conditions. Here, we evaluated the effects of specific TYK2 inhibitors (TYK2is) in pre-clinical models of T1D. First, human β cells, cadaveric donor islets, and iPSC-derived islets were treated in vitro with IFNα in combination with a small molecule TYK2i (BMS-986165 or a related molecule BMS-986202). TYK2 inhibition prevented IFNα-induced β cell HLA class I up-regulation, endoplasmic reticulum stress, and chemokine production. In co-culture studies, pre-treatment of β cells with a TYK2i prevented IFNα-induced activation of T cells targeting an epitope of insulin. In vivo administration of BMS-986202 in two mouse models of T1D (RIP-LCMV-GP mice and NOD mice) reduced systemic and tissue-localized inflammation, prevented β cell death, and delayed T1D onset. Transcriptional phenotyping of pancreatic islets, pancreatic lymph nodes (PLN), and spleen during early disease pathogenesis highlighted a role for TYK2 inhibition in modulating signaling pathways associated with inflammation, translational control, stress signaling, secretory function, immunity, and diabetes. Additionally, TYK2i treatment changed the composition of innate and adaptive immune cell populations in the blood and disease target tissues, resulting in an immune phenotype with a diminished capacity for β cell destruction. Overall, these findings indicate that TYK2i has beneficial effects in both the immune and endocrine compartments in models of T1D, thus supporting a path forward for testing TYK2 inhibitors in human T1D.

RedRibbon: A new rank-rank hypergeometric overlap for gene and transcript expression signatures

High-throughput omics technologies have generated a wealth of large protein, gene, and transcript datasets that have exacerbated the need for new methods to analyse and compare big datasets. Rank-rank hypergeometric overlap is an important threshold-free method to combine and visualize two ranked lists of P-values or fold-changes, usually from differential gene expression analyses. Here, we introduce a new rank-rank hypergeometric overlap-based method aimed at gene level and alternative splicing analyses at transcript or exon level, hitherto unreachable as transcript numbers are an order of magnitude larger than gene numbers. We tested the tool on synthetic and real datasets at gene and transcript levels to detect correlation and anticorrelation patterns and found it to be fast and accurate, even on very large datasets thanks to an evolutionary algorithm-based minimal P-value search. The tool comes with a ready-to-use permutation scheme allowing the computation of adjusted P-values at low time cost. The package compatibility mode is a drop-in replacement to previous packages. RedRibbon holds the promise to accurately extricate detailed information from large comparative analyses.

A computational method for small molecule-RNA binding sites identification by utilizing position specificity and complex network information

Some computational methods have been given for small molecule-RNA binding site identification due to that it plays a significant role in revealing biology function researches. However, it is still challenging to design an accurate model, especially for MCC. We designed a feature extraction technology from two aspects (position specificity and complex network information). Specifically, complex network was employed to express the space topological structure and sequence position information for improving prediction effect. Then, the features fused position specificity and complex network information were input into random forest classifier for model construction. The AUC of 88.22%, 77.92% and 81.46% were obtained on three independent datasets (RB19, CS71, RB78). Compared with the existing method, the best MCC were obtained on three datasets, which were 8.19%, 0.59% and 4.35% higher than the state-of-the-art prediction methods, respectively. The outstanding performances show that our method is a powerful tool to identify RNA binding sites, helping to the design RNA-targeting small molecule drugs. The data and resource codes are available at https://github.com/Kangxiaoneuq/PCN_RNAsite.

Non-mutational neoantigens in disease

The ability of mammals to mount adaptive immune responses culminating with the establishment of immunological memory is predicated on the ability of the mature T cell repertoire to recognize antigenic peptides presented by syngeneic MHC class I and II molecules. Although it is widely believed that mature T cells are highly skewed towards the recognition of antigenic peptides originating from genetically diverse (for example, foreign or mutated) protein-coding regions, preclinical and clinical data rather demonstrate that novel antigenic determinants efficiently recognized by mature T cells can emerge from a variety of non-mutational mechanisms. In this Review, we describe various mechanisms that underlie the formation of bona fide non-mutational neoantigens, such as epitope mimicry, upregulation of cryptic epitopes, usage of non-canonical initiation codons, alternative RNA splicing, and defective ribosomal RNA processing, as well as both enzymatic and non-enzymatic post-translational protein modifications. Moreover, we discuss the implications of the immune recognition of non-mutational neoantigens for human disease.

IFNɣ but not IFNα increases recognition of insulin defective ribosomal product-derived antigen to amplify islet autoimmunity

AIMS/HYPOTHESIS

The inflammatory milieu characteristic of insulitis affects translation fidelity and generates defective ribosomal products (DRiPs) that participate in autoimmune beta cell destruction in type 1 diabetes. Here, we studied the role of early innate cytokines (IFNα) and late immune adaptive events (IFNɣ) in insulin DRiP-derived peptide presentation to diabetogenic CD8+ T cells.

METHODS

Single-cell transcriptomics of human pancreatic islets was used to study the composition of the (immuno)proteasome. Specific inhibition of the immunoproteasome catalytic subunits was achieved using siRNA, and antigenic peptide presentation at the cell surface of the human beta cell line EndoC-βH1 was monitored using peptide-specific CD8 T cells.

RESULTS

We found that IFNγ induces the expression of the PSMB10 transcript encoding the β2i catalytic subunit of the immunoproteasome in endocrine beta cells, revealing a critical role in insulin DRiP-derived peptide presentation to T cells. Moreover, we showed that PSMB10 is upregulated in a beta cell subset that is preferentially destroyed in the pancreases of individuals with type 1 diabetes.

CONCLUSIONS/INTERPRETATION

Our data highlight the role of the degradation machinery in beta cell immunogenicity and emphasise the need for evaluation of targeted immunoproteasome inhibitors to limit beta cell destruction in type 1 diabetes.

DATA AVAILABILITY

The single-cell RNA-seq dataset is available from the Gene Expression Omnibus (GEO) using the accession number GSE218316 ( https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE218316 ).

Heimdall, an alternative protein issued from a ncRNA related to kappa light chain variable region of immunoglobulins from astrocytes: a new player in neural proteome

The dogma "One gene, one protein" is clearly obsolete since cells use alternative splicing and generate multiple transcripts which are translated into protein isoforms, but also use alternative translation initiation sites (TISs) and termination sites on a given transcript. Alternative open reading frames for individual transcripts give proteins originate from the 5'- and 3'-UTR mRNA regions, frameshifts of mRNA ORFs or from non-coding RNAs. Longtime considered as non-coding, recent in-silico translation prediction methods enriched the protein databases allowing the identification of new target structures that have not been identified previously. To gain insight into the role of these newly identified alternative proteins in the regulation of cellular functions, it is crucial to assess their dynamic modulation within a framework of altered physiological modifications such as experimental spinal cord injury (SCI). Here, we carried out a longitudinal proteomic study on rat SCI from 12 h to 10 days. Based on the alternative protein predictions, it was possible to identify a plethora of newly predicted protein hits. Among these proteins, some presented a special interest due to high homology with variable chain regions of immunoglobulins. We focus our interest on the one related to Kappa variable light chains which is similarly highly produced by B cells in the Bence jones disease, but here expressed in astrocytes. This protein, name Heimdall is an Intrinsically disordered protein which is secreted under inflammatory conditions. Immunoprecipitation experiments showed that the Heimdall interactome contained proteins related to astrocyte fate keepers such as "NOTCH1, EPHA3, IPO13" as well as membrane receptor protein including "CHRNA9; TGFBR, EPHB6, and TRAM". However, when Heimdall protein was neutralized utilizing a specific antibody or its gene knocked out by CRISPR-Cas9, sprouting elongations were observed in the corresponding astrocytes. Interestingly, depolarization assays and intracellular calcium measurements in Heimdall KO, established a depolarization effect on astrocyte membranes KO cells were more likely that the one found in neuroprogenitors. Proteomic analyses performed under injury conditions or under lipopolysaccharides (LPS) stimulation, revealed the expression of neuronal factors, stem cell proteins, proliferation, and neurogenesis of astrocyte convertor factors such as EPHA4, NOTCH2, SLIT3, SEMA3F, suggesting a role of Heimdall could regulate astrocytic fate. Taken together, Heimdall could be a novel member of the gatekeeping astrocyte-to-neuroprogenitor conversion factors.

Presence of immunogenic alternatively spliced insulin gene product in human pancreatic delta cells

AIMS/HYPOTHESIS

Transcriptome analyses revealed insulin-gene-derived transcripts in non-beta endocrine islet cells. We studied alternative splicing of human INS mRNA in pancreatic islets.

METHODS

Alternative splicing of insulin pre-mRNA was determined by PCR analysis performed on human islet RNA and single-cell RNA-seq analysis. Antisera were generated to detect insulin variants in human pancreatic tissue using immunohistochemistry, electron microscopy and single-cell western blot to confirm the expression of insulin variants. Cytotoxic T lymphocyte (CTL) activation was determined by MIP-1β release.

RESULTS

We identified an alternatively spliced INS product. This variant encodes the complete insulin signal peptide and B chain and an alternative C-terminus that largely overlaps with a previously identified defective ribosomal product of INS. Immunohistochemical analysis revealed that the translation product of this INS-derived splice transcript was detectable in somatostatin-producing delta cells but not in beta cells; this was confirmed by light and electron microscopy. Expression of this alternatively spliced INS product activated preproinsulin-specific CTLs in vitro. The exclusive presence of this alternatively spliced INS product in delta cells may be explained by its clearance from beta cells by insulin-degrading enzyme capturing its insulin B chain fragment and a lack of insulin-degrading enzyme expression in delta cells.

CONCLUSIONS/INTERPRETATION

Our data demonstrate that delta cells can express an INS product derived from alternative splicing, containing both the diabetogenic insulin signal peptide and B chain, in their secretory granules. We propose that this alternative INS product may play a role in islet autoimmunity and pathology, as well as endocrine or paracrine function or islet development and endocrine destiny, and transdifferentiation between endocrine cells. INS promoter activity is not confined to beta cells and should be used with care when assigning beta cell identity and selectivity.

DATA AVAILABILITY

The full EM dataset is available via www.nanotomy.org (for review: http://www.nanotomy.org/OA/Tienhoven2021SUB/6126-368/ ). Single-cell RNA-seq data was made available by Segerstolpe et al [13] and can be found at https://sandberglab.se/pancreas . The RNA and protein sequence of INS-splice was uploaded to GenBank (BankIt2546444 INS-splice OM489474).

The Hidden Enemy Within: Non-canonical Peptides in Virus-Induced Autoimmunity

Viruses play a key role in explaining the pathogenesis of various autoimmune disorders, whose underlying principle is defined by the activation of autoreactive T-cells. In many cases, T-cells escape self-tolerance due to the failure in encountering certain MHC-I self-peptide complexes at substantial levels, whose peptides remain invisible from the immune system. Over the years, contribution of unstable defective ribosomal products (DRiPs) in immunosurveillance has gained prominence. A class of unstable products emerge from non-canonical translation and processing of unannotated mammalian and viral ORFs and their peptides are cryptic in nature. Indeed, high throughput sequencing and proteomics have revealed that a substantial portion of our genomes comprise of non-canonical ORFs, whose generation is significantly modulated during disease. Many of these ORFs comprise short ORFs (sORFs) and upstream ORFs (uORFs) that resemble DRiPs and may hence be preferentially presented. Here, we discuss how such products, normally “hidden” from the immune system, become abundant in viral infections activating autoimmune T-cells, by discussing their emerging role in infection and disease. Finally, we provide a perspective on how these mechanisms can explain several autoimmune disorders in the wake of the COVID-19 pandemic.