Divergent molecular networks program functionally distinct CD8+ skin-resident memory T cells.

Skin-resident CD8+ T cells include distinct interferon-γ-producing [tissue-resident memory T type 1 (TRM1)] and interleukin-17 (IL-17)-producing (TRM17) subsets that differentially contribute to immune responses. However, whether these populations use common mechanisms to establish tissue residence is unknown. In this work, we show that TRM1 and TRM17 cells navigate divergent trajectories to acquire tissue residency in the skin. TRM1 cells depend on a T-bet-Hobit-IL-15 axis, whereas TRM17 cells develop independently of these factors. Instead, c-Maf commands a tissue-resident program in TRM17 cells parallel to that induced by Hobit in TRM1 cells, with an ICOS-c-Maf-IL-7 axis pivotal to TRM17 cell commitment. Accordingly, by targeting this pathway, skin TRM17 cells can be ablated without compromising their TRM1 counterparts. Thus, skin-resident T cells rely on distinct molecular circuitries, which can be exploited to strategically modulate local immunity.

Establishing the Link between X-Chromosome Aberrations and TP53 Status, with Breast Cancer Patient Outcomes.

Ubiquitous to normal female human somatic cells, X-chromosome inactivation (XCI) tightly regulates the transcriptional silencing of a single X chromosome from each pair. Some genes escape XCI, including crucial tumour suppressors. Cancer susceptibility can be influenced by the variability in the genes that escape XCI. The mechanisms of XCI dysregulation remain poorly understood in complex diseases, including cancer. Using publicly available breast cancer next-generation sequencing data, we show that the status of the major tumour suppressor TP53 from Chromosome 17 is highly associated with the genomic integrity of the inactive X (Xi) and the active X (Xa) chromosomes. Our quantification of XCI and XCI escape demonstrates that aberrant XCI is linked to poor survival. We derived prognostic gene expression signatures associated with either large deletions of Xi; large amplifications of Xa; or abnormal X-methylation. Our findings expose a novel insight into female cancer risks, beyond those associated with the standard molecular subtypes.

sccomp: Robust differential composition and variability analysis for single-cell data.

Cellular omics such as single-cell genomics, proteomics, and microbiomics allow the characterization of tissue and microbial community composition, which can be compared between conditions to identify biological drivers. This strategy has been critical to revealing markers of disease progression, such as cancer and pathogen infection. A dedicated statistical method for differential variability analysis is lacking for cellular omics data, and existing methods for differential composition analysis do not model some compositional data properties, suggesting there is room to improve model performance. Here, we introduce sccomp, a method for differential composition and variability analyses that jointly models data count distribution, compositionality, group-specific variability, and proportion mean-variability association, being aware of outliers. sccomp provides a comprehensive analysis framework that offers realistic data simulation and cross-study knowledge transfer. Here, we demonstrate that mean-variability association is ubiquitous across technologies, highlighting the inadequacy of the very popular Dirichlet-multinomial distribution. We show that sccomp accurately fits experimental data, significantly improving performance over state-of-the-art algorithms. Using sccomp, we identified differential constraints and composition in the microenvironment of primary breast cancer.

Spatial analysis with SPIAT and spaSim to characterize and simulate tissue microenvironments.

Spatial proteomics technologies have revealed an underappreciated link between the location of cells in tissue microenvironments and the underlying biology and clinical features, but there is significant lag in the development of downstream analysis methods and benchmarking tools. Here we present SPIAT (spatial image analysis of tissues), a spatial-platform agnostic toolkit with a suite of spatial analysis algorithms, and spaSim (spatial simulator), a simulator of tissue spatial data. SPIAT includes multiple colocalization, neighborhood and spatial heterogeneity metrics to characterize the spatial patterns of cells. Ten spatial metrics of SPIAT are benchmarked using simulated data generated with spaSim. We show how SPIAT can uncover cancer immune subtypes correlated with prognosis in cancer and characterize cell dysfunction in diabetes. Our results suggest SPIAT and spaSim as useful tools for quantifying spatial patterns, identifying and validating correlates of clinical outcomes and supporting method development.

Early immune pressure initiated by tissue-resident memory T cells sculpts tumor evolution in non-small cell lung cancer.

Tissue-resident memory T (TRM) cells provide immune defense against local infection and can inhibit cancer progression. However, it is unclear to what extent chronic inflammation impacts TRM activation and whether TRM cells existing in tissues before tumor onset influence cancer evolution in humans. We performed deep profiling of healthy lungs and lung cancers in never-smokers (NSs) and ever-smokers (ESs), finding evidence of enhanced immunosurveillance by cells with a TRM-like phenotype in ES lungs. In preclinical models, tumor-specific or bystander TRM-like cells present prior to tumor onset boosted immune cell recruitment, causing tumor immune evasion through loss of MHC class I protein expression and resistance to immune checkpoint inhibitors. In humans, only tumors arising in ES patients underwent clonal immune evasion, unrelated to tobacco-associated mutagenic signatures or oncogenic drivers. These data demonstrate that enhanced TRM-like activity prior to tumor development shapes the evolution of tumor immunogenicity and can impact immunotherapy outcomes.

Intratumoral CD8+ T cells with a tissue-resident memory phenotype mediate local immunity and immune checkpoint responses in breast cancer.

CD8+ tumor-infiltrating lymphocytes with a tissue-resident memory T (TRM) cell phenotype are associated with favorable prognosis in patients with triple-negative breast cancer (TNBC). However, the relative contribution of CD8+ TRM cells to anti-tumor immunity and immune checkpoint blockade efficacy in breast cancer remains unknown. Here, we show that intratumoral CD8+ T cells in murine mammary tumors transcriptionally resemble those from TNBC patients. Phenotypic and transcriptional studies established two intratumoral sub-populations: one more enriched in markers of terminal exhaustion (TEX-like) and the other with a bona fide resident phenotype (TRM-like). Treatment with anti-PD-1 and anti-CTLA-4 therapy resulted in expansion of these intratumoral populations, with the TRM-like subset displaying significantly enhanced cytotoxic capacity. TRM-like CD8+ T cells could also provide local immune protection against tumor rechallenge and a TRM gene signature extracted from tumor-free tissue was significantly associated with improved clinical outcomes in TNBC patients treated with checkpoint inhibitors.

Removing unwanted variation from large-scale RNA sequencing data with PRPS.

Accurate identification and effective removal of unwanted variation is essential to derive meaningful biological results from RNA sequencing (RNA-seq) data, especially when the data come from large and complex studies. Using RNA-seq data from The Cancer Genome Atlas (TCGA), we examined several sources of unwanted variation and demonstrate here how these can significantly compromise various downstream analyses, including cancer subtype identification, association between gene expression and survival outcomes and gene co-expression analysis. We propose a strategy, called pseudo-replicates of pseudo-samples (PRPS), for deploying our recently developed normalization method, called removing unwanted variation III (RUV-III), to remove the variation caused by library size, tumor purity and batch effects in TCGA RNA-seq data. We illustrate the value of our approach by comparing it to the standard TCGA normalizations on several TCGA RNA-seq datasets. RUV-III with PRPS can be used to integrate and normalize other large transcriptomic datasets coming from multiple laboratories or platforms.

Estimation of tumor cell total mRNA expression in 15 cancer types predicts disease progression.

Single-cell RNA sequencing studies have suggested that total mRNA content correlates with tumor phenotypes. Technical and analytical challenges, however, have so far impeded at-scale pan-cancer examination of total mRNA content. Here we present a method to quantify tumor-specific total mRNA expression (TmS) from bulk sequencing data, taking into account tumor transcript proportion, purity and ploidy, which are estimated through transcriptomic/genomic deconvolution. We estimate and validate TmS in 6,590 patient tumors across 15 cancer types, identifying significant inter-tumor variability. Across cancers, high TmS is associated with increased risk of disease progression and death. TmS is influenced by cancer-specific patterns of gene alteration and intra-tumor genetic heterogeneity as well as by pan-cancer trends in metabolic dysregulation. Taken together, our results indicate that measuring cell-type-specific total mRNA expression in tumor cells predicts tumor phenotypes and clinical outcomes.

A mouse-specific retrotransposon drives a conserved Cdk2ap1 isoform essential for development.

Retrotransposons mediate gene regulation in important developmental and pathological processes. Here, we characterized the transient retrotransposon induction during preimplantation development of eight mammals. Induced retrotransposons exhibit similar preimplantation profiles across species, conferring gene regulatory activities, particularly through long terminal repeat (LTR) retrotransposon promoters. A mouse-specific MT2B2 retrotransposon promoter generates an N-terminally truncated Cdk2ap1ΔN that peaks in preimplantation embryos and promotes proliferation. In contrast, the canonical Cdk2ap1 peaks in mid-gestation and represses cell proliferation. This MT2B2 promoter, whose deletion abolishes Cdk2ap1ΔN production, reduces cell proliferation and impairs embryo implantation, is developmentally essential. Intriguingly, Cdk2ap1ΔN is evolutionarily conserved in sequence and function yet is driven by different promoters across mammals. The distinct preimplantation Cdk2ap1ΔN expression in each mammalian species correlates with the duration of its preimplantation development. Hence, species-specific transposon promoters can yield evolutionarily conserved, alternative protein isoforms, bestowing them with new functions and species-specific expression to govern essential biological divergence.

The Ratio of Exhausted to Resident Infiltrating Lymphocytes Is Prognostic for Colorectal Cancer Patient Outcome.

Immunotherapy success in colorectal cancer is mainly limited to patients whose tumors exhibit high microsatellite instability (MSI). However, there is variability in treatment outcomes within this group, which is in part driven by the frequency and characteristics of tumor-infiltrating immune cells. Indeed, the presence of specific infiltrating immune-cell subsets has been shown to correlate with immunotherapy response and is in many cases prognostic of treatment outcome. Tumor-infiltrating lymphocytes (TIL) can undergo distinct differentiation programs, acquiring features of tissue-residency or exhaustion, a process during which T cells upregulate inhibitory receptors, such as PD-1, and lose functionality. Although residency and exhaustion programs of CD8+ T cells are relatively well studied, these programs have only recently been appreciated in CD4+ T cells and remain largely unknown in tumor-infiltrating natural killer (NK) cells. In this study, we used single-cell RNA sequencing (RNA-seq) data to identify signatures of residency and exhaustion in colorectal cancer-infiltrating lymphocytes, including CD8+, CD4+, and NK cells. We then tested these signatures in independent single-cell data from tumor and normal tissue-infiltrating immune cells. Furthermore, we used versions of these signatures designed for bulk RNA-seq data to explore tumor-intrinsic mutations associated with residency and exhaustion from TCGA data. Finally, using two independent transcriptomic datasets from patients with colon adenocarcinoma, we showed that combinations of these signatures, in particular combinations of NK-cell activity signatures, together with tumor-associated signatures, such as TGFß signaling, were associated with distinct survival outcomes in patients with colon adenocarcinoma.

Finding a suitable library size to call variants in RNA-Seq.

BACKGROUND: RNA sequencing allows the study of both gene expression changes and transcribed mutations, providing a highly effective way to gain insight into cancer biology. When planning the sequencing of a large cohort of samples, library size is a fundamental factor affecting both the overall cost and the quality of the results. Here we specifically address how overall library size influences the detection of somatic mutations in RNA-seq data in two acute myeloid leukaemia datasets. RESULTS : We simulated shallower sequencing depths by downsampling 45 acute myeloid leukaemia samples (100 bp PE) that are part of the Leucegene project, which were originally sequenced at high depth. We compared the sensitivity of six methods of recovering validated mutations on the same samples. The methods compared are a combination of three popular callers (MuTect, VarScan, and VarDict) and two filtering strategies. We observed an incremental loss in sensitivity when simulating libraries of 80M, 50M, 40M, 30M and 20M fragments, with the largest loss detected with less than 30M fragments (below 90%, average loss of 7%). The sensitivity in recovering insertions and deletions varied markedly between callers, with VarDict showing the highest sensitivity (60%). Single nucleotide variant sensitivity is relatively consistent across methods, apart from MuTect, whose default filters need adjustment when using RNA-Seq. We also analysed 136 RNA-Seq samples from the TCGA-LAML cohort (50 bp PE) and assessed the change in sensitivity between the initial libraries (average 59M fragments) and after downsampling to 40M fragments. When considering single nucleotide variants in recurrently mutated myeloid genes we found a comparable performance, with a 6% average loss in sensitivity using 40M fragments. CONCLUSIONS: Between 30M and 40M 100 bp PE reads are needed to recover 90-95% of the initial variants on recurrently mutated myeloid genes. To extend this result to another cancer type, an exploration of the characteristics of its mutations and gene expression patterns is suggested.

Correction to: Serum microRNA is a biomarker for post-operative monitoring in glioma.

For the reference citation '[57]' in the second paragraph of the Results section of the original article there was no corresponding entry in the References section. It should have referred to the below mentioned article by Ebrahimkhani et al. (2018).

Serum microRNA is a biomarker for post-operative monitoring in glioma.

PURPOSE: A circulating biomarker has potential to provide more accurate information for glioma progression post treatment, however no such biomarker is currently available. We aimed to discover a microRNA serum biomarker for longitudinal monitoring of glioma patients. METHODS: A prospectively collected cohort of 91 glioma patients and 17 healthy controls underwent pre and post-operative serum miRNA profiling using Nanostring®. Differentially expressed miRNAs were discovered using a machine learning random forest analysis. Candidate miRNAs were then assessed by droplet digital PCR in 11 patients with multiple follow up samples and compared to tumor volume based on magnetic resonance imaging. RESULTS: A 9-gene miRNA signature was identified that could distinguish between glioma and healthy controls with 99.8% accuracy. Two miRNAs miR-223 and miR-320e, best demonstrated dynamic changes that correlated closely with tumor volume in LGG and GBM respectively. Importantly, miRNA levels did not increase in two cases of pseudo-progression, indicating the potential utility of this test in guiding treatment decisions. CONCLUSIONS: We identified a highly accurate 9-miRNA signature associated with glioma serum. Additionally, we observed dynamic changes in specific miRNAs correlating with tumor volume over long-term follow up. These results support a large prospective validation study of serum miRNA biomarkers in glioma.

Removing unwanted variation with CytofRUV to integrate multiple CyTOF datasets.

Mass cytometry (CyTOF) is a technology that has revolutionised single-cell biology. By detecting over 40 proteins on millions of single cells, CyTOF allows the characterisation of cell subpopulations in unprecedented detail. However, most CyTOF studies require the integration of data from multiple CyTOF batches usually acquired on different days and possibly at different sites. To date, the integration of CyTOF datasets remains a challenge due to technical differences arising in multiple batches. To overcome this limitation, we developed an approach called CytofRUV for analysing multiple CyTOF batches, which includes an R-Shiny application with diagnostic plots. CytofRUV can correct for batch effects and integrate data from large numbers of patients and conditions across batches, to confidently compare cellular changes and correlate these with clinically relevant outcomes.

Strategies to enable large-scale proteomics for reproducible research.

Reproducible research is the bedrock of experimental science. To enable the deployment of large-scale proteomics, we assess the reproducibility of mass spectrometry (MS) over time and across instruments and develop computational methods for improving quantitative accuracy. We perform 1560 data independent acquisition (DIA)-MS runs of eight samples containing known proportions of ovarian and prostate cancer tissue and yeast, or control HEK293T cells. Replicates are run on six mass spectrometers operating continuously with varying maintenance schedules over four months, interspersed with ~5000 other runs. We utilise negative controls and replicates to remove unwanted variation and enhance biological signal, outperforming existing methods. We also design a method for reducing missing values. Integrating these computational modules into a pipeline (ProNorM), we mitigate variation among instruments over time and accurately predict tissue proportions. We demonstrate how to improve the quantitative analysis of large-scale DIA-MS data, providing a pathway toward clinical proteomics.

Changes in ferrous iron and glutathione promote ferroptosis and frailty in aging Caenorhabditis elegans.

All eukaryotes require iron. Replication, detoxification, and a cancer-protective form of regulated cell death termed ferroptosis, all depend on iron metabolism. Ferrous iron accumulates over adult lifetime in Caenorhabditis elegans. Here, we show that glutathione depletion is coupled to ferrous iron elevation in these animals, and that both occur in late life to prime cells for ferroptosis. We demonstrate that blocking ferroptosis, either by inhibition of lipid peroxidation or by limiting iron retention, mitigates age-related cell death and markedly increases lifespan and healthspan. Temporal scaling of lifespan is not evident when ferroptosis is inhibited, consistent with this cell death process acting at specific life phases to induce organismal frailty, rather than contributing to a constant aging rate. Because excess age-related iron elevation in somatic tissue, particularly in brain, is thought to contribute to degenerative disease, post-developmental interventions to limit ferroptosis may promote healthy aging.

Identification of cancer sex-disparity in the functional integrity of p53 and its X chromosome network.

The disproportionately high prevalence of male cancer is poorly understood. We tested for sex-disparity in the functional integrity of the major tumor suppressor p53 in sporadic cancers. Our bioinformatics analyses expose three novel levels of p53 impact on sex-disparity in 12 non-reproductive cancer types. First, TP53 mutation is more frequent in these cancers among US males than females, with poorest survival correlating with its mutation. Second, numerous X-linked genes are associated with p53, including vital genomic regulators. Males are at unique risk from alterations of their single copies of these genes. High expression of X-linked negative regulators of p53 in wild-type TP53 cancers corresponds with reduced survival. Third, females exhibit an exceptional incidence of non-expressed mutations among p53-associated X-linked genes. Our data indicate that poor survival in males is contributed by high frequencies of TP53 mutations and an inability to shield against deregulated X-linked genes that engage in p53 networks.

A new normalization for Nanostring nCounter gene expression data.

The Nanostring nCounter gene expression assay uses molecular barcodes and single molecule imaging to detect and count hundreds of unique transcripts in a single reaction. These counts need to be normalized to adjust for the amount of sample, variations in assay efficiency and other factors. Most users adopt the normalization approach described in the nSolver analysis software, which involves background correction based on the observed values of negative control probes, a within-sample normalization using the observed values of positive control probes and normalization across samples using reference (housekeeping) genes. Here we present a new normalization method, Removing Unwanted Variation-III (RUV-III), which makes vital use of technical replicates and suitable control genes. We also propose an approach using pseudo-replicates when technical replicates are not available. The effectiveness of RUV-III is illustrated on four different datasets. We also offer suggestions on the design and analysis of studies involving this technology.

Dysregulated IL-1ß-GM-CSF Axis in Acute Rheumatic Fever That Is Limited by Hydroxychloroquine.

BACKGROUND: Acute rheumatic fever (ARF) and rheumatic heart disease are autoimmune consequences of group A streptococcus infection and remain major causes of cardiovascular morbidity and mortality around the world. Improved treatment has been stymied by gaps in understanding key steps in the immunopathogenesis of ARF and rheumatic heart disease. This study aimed to identify (1) effector T cell cytokine(s) that might be dysregulated in the autoimmune response of patients with ARF by group A streptococcus, and (2) an immunomodulatory agent that suppresses this response and could be clinically translatable to high-risk patients with ARF. METHODS: The immune response to group A streptococcus was analyzed in peripheral blood mononuclear cells from an Australian Aboriginal ARF cohort by a combination of multiplex cytokine array, flow cytometric analysis, and global gene expression analysis by RNA sequencing. The immunomodulatory drug hydroxychloroquine was tested for effects on this response. RESULTS: We found a dysregulated interleukin-1ß-granulocyte-macrophage colony-stimulating factor (GM-CSF) cytokine axis in ARF peripheral blood mononuclear cells exposed to group A streptococcus in vitro, whereby persistent interleukin-1ß production is coupled to overproduction of GM-CSF and selective expansion of CXCR3+CCR4-CCR6- CD4 T cells. CXCR3+CCR4-CCR6- CD4 T cells are the major source of GM-CSF in human CD4 T cells and CXCL10, a CXCR3 ligand and potent T helper 1 chemoattractant, was elevated in sera from patients with ARF. GM-CSF has recently emerged as a key T cell-derived effector cytokine in numerous autoimmune diseases, including myocarditis, and the production of CXCL10 may explain selective trafficking of these cells to the heart. We provide evidence that interleukin-1ß amplifies the expansion of GM-CSF-expressing CD4 T cells, which is effectively suppressed by hydroxychloroquine. RNA sequencing showed shifts in gene expression profiles and differentially expressed genes in peripheral blood mononuclear cells derived from patients at different clinical stages of ARF. CONCLUSIONS: Given the safety profile of hydroxychloroquine and its clinical pedigree in treating autoimmune diseases such as rheumatoid arthritis, where GM-CSF plays a pivotal role, we propose that hydroxychloroquine could be repurposed to reduce the risk of rheumatic heart disease after ARF.

Publisher Correction: Single-cell profiling of breast cancer T cells reveals a tissue-resident memory subset associated with improved prognosis.

In the version of this article originally published, the institution in affiliation 10 was missing. Affiliation 10 was originally listed as Department of Surgery, Royal Melbourne Hospital and Royal Womens' Hospital, Melbourne, Victoria, Australia. It should have been Department of Surgery, Royal Melbourne Hospital and Royal Womens' Hospital, University of Melbourne, Melbourne, Victoria, Australia. The error has been corrected in the HTML and PDF versions of this article.