High numbers of tissue-resident memory T (TRM) cells are associated with better clinical outcomes in cancer patients. However, the molecular characteristics that drive their efficient immune response to tumors are poorly understood. Here, single-cell and bulk transcriptomic analysis of TRM and non-TRM cells present in tumor and normal lung tissue from patients with lung cancer revealed that PD-1–expressing TRM cells in tumors were clonally expanded and enriched for transcripts linked to cell proliferation and cytotoxicity when compared with PD-1–expressing non-TRM cells. This feature was more prominent in the TRM cell subset coexpressing PD-1 and TIM-3, and it was validated by functional assays ex vivo and also reflected in their chromatin accessibility profile. This PD-1+TIM-3+ TRM cell subset was enriched in responders to PD-1 inhibitors and in tumors with a greater magnitude of CTL responses. These data highlight that not all CTLs expressing PD-1 are dysfunctional; on the contrary, TRM cells with PD-1 expression were enriched for features suggestive of superior functionality.

Motivation: MicroRNAs (miRNAs) are small non-coding RNAs that are involved in posttranscriptional regulation of gene expression. In this high-throughput sequencing era, a tremendous amount of RNA-seq data is accumulating, and full utilization of publicly available miRNA data is an important challenge. These data are useful to determine expression values for each miRNA, but quantification pipelines are in a primitive stage and still evolving; there are many factors that affect expression values significantly. Results: We used 304 high-quality microRNA sequencing (miRNA-seq) datasets from NCBI-SRA and calculated expression profiles for different tissues and cell-lines. In each miRNA-seq dataset, we found an average of more than 500 miRNAs with higher than 5x coverage, and we explored the top five highly expressed miRNAs in each tissue and cell-line. This user-friendly miRmine database has options to retrieve expression profiles of single or multiple miRNAs for a specific tissue or cellline, either normal or with disease information. Results can be displayed in multiple interactive, graphical and downloadable formats.

BackgroundEvidence is accumulating that non-coding transcripts, previously thought to be functionally inert, play important roles in various cellular activities. High throughput techniques like next generation sequencing have resulted in the generation of vast amounts of sequence data. It is therefore desirable, not only to discriminate coding and non-coding transcripts, but also to assign the noncoding RNA (ncRNA) transcripts into respective classes (families). Although there are several algorithms available for this task, their classification performance remains a major concern. Acknowledging the crucial role that non-coding transcripts play in cellular processes, it is required to develop algorithms that are able to precisely classify ncRNA transcripts.ResultsIn this study, we initially develop prediction tools to discriminate coding or non-coding transcripts and thereafter classify ncRNAs into respective classes. In comparison to the existing methods that employed multiple features, our SVM-based method by using a single feature (tri-nucleotide composition), achieved MCC of 0.98. Knowing that the structure of a ncRNA transcript could provide insights into its biological function, we use graph properties of predicted ncRNA structures to classify the transcripts into 18 different non-coding RNA classes. We developed classification models using a variety of algorithms (BayeNet, NaiveBayes, MultilayerPerceptron, IBk, libSVM, SMO and RandomForest) and observed that model based on RandomForest performed better than other models. As compared to the GraPPLE study, the sensitivity (of 13 classes) and specificity (of 14 classes) was higher. Moreover, the overall sensitivity of 0.43 outperforms the sensitivity of GraPPLE (0.33) whereas the overall MCC measure of 0.40 (in contrast to MCC of 0.29 of GraPPLE) was significantly higher for our method. This clearly demonstrates that our models are more accurate than existing models.ConclusionsThis work conclusively demonstrates that a simple feature, tri-nucleotide composition, is sufficient to discriminate between coding and non-coding RNA sequences. Similarly, graph properties based feature set along with RandomForest algorithm are most suitable to classify different ncRNA classes. We have also developed an online and standalone tool-- RNAcon ( http://crdd.osdd.net/raghava/rnacon).

Spinocerebellar ataxia type 3 (SCA3) is a neurodegenerative disorder caused by a polyglutamine-encoding CAG repeat expansion in the ATXN3 gene. This expansion leads to misfolding and aggregation of mutant ataxin-3 (ATXN3) and degeneration of select brain regions. A key unanswered question in SCA3 and other polyglutamine diseases is the extent to which neurodegeneration is mediated through gain-of-function versus loss-of-function. To address this question in SCA3, we performed transcriptional profiling on the brainstem, a highly vulnerable brain region in SCA3, in a series of mouse models with varying degrees of ATXN3 expression and aggregation. We include two SCA3 knock-in mouse models: our previously published model that erroneously harbors a tandem duplicate of the CAG repeat-containing exon, and a corrected model, introduced here. Both models exhibit dose-dependent neuronal accumulation and aggregation of mutant ATXN3, but do not exhibit a behavioral phenotype. We identified a molecular signature that correlates with ATXN3 neuronal aggregation yet is primarily linked to oligodendrocytes, highlighting early white matter dysfunction in SCA3. Two robustly elevated oligodendrocyte transcripts, Acy3 and Tnfrsf13c, were confirmed as elevated at the protein level in SCA3 human disease brainstem. To determine if mutant ATXN3 acts on oligodendrocytes cell-autonomously, we manipulated the repeat expansion in the variant SCA3 knock-in mouse by cell-type specific Cre/LoxP recombination. Changes in oligodendrocyte transcripts are driven cell-autonomously and occur independent of neuronal ATXN3 aggregation. Our findings support a primary toxic gain of function mechanism and highlight a previously unrecognized role for oligodendrocyte dysfunction in SCA3 disease pathogenesis.

Rheumatoid arthritis (RA) affects millions world-wide. While anti-TNF treatment is widely used to reduce disease progression, treatment fails in ∼one-third of patients. No biomarker currently exists that identifies non-responders before treatment. A rigorous community-based assessment of the utility of SNP data for predicting anti-TNF treatment efficacy in RA patients was performed in the context of a DREAM Challenge (http://www.synapse.org/RA_Challenge). An open challenge framework enabled the comparative evaluation of predictions developed by 73 research groups using the most comprehensive available data and covering a wide range of state-of-the-art modelling methodologies. Despite a significant genetic heritability estimate of treatment non-response trait (h2=0.18, P value=0.02), no significant genetic contribution to prediction accuracy is observed. Results formally confirm the expectations of the rheumatology community that SNP information does not significantly improve predictive performance relative to standard clinical traits, thereby justifying a refocusing of future efforts on collection of other data.