Single cell view of tumor microenvironment gradients in pleural mesothelioma.

Immunotherapies have shown great promise in pleural mesothelioma (PM), yet most patients still do not achieve significant clinical response, highlighting the importance of improving understanding of the tumor microenvironment (TME). Here, we utilized high-throughput, single-cell RNA-sequencing to de novo identify 54 expression programs and construct a comprehensive cellular catalogue of the PM TME. We found four cancer-intrinsic programs associated with poor disease outcome and a novel fetal-like, endothelial cell population that likely responds to VEGF signaling and promotes angiogenesis. Throughout cellular compartments, we observe substantial difference in the TME associated with a cancer-intrinsic sarcomatoid signature, including enrichment in fetal-like endothelial cells, CXCL9+ macrophages, cytotoxic, exhausted, and regulatory T cells, which we validated using imaging and bulk deconvolution analyses on independent cohorts. Finally, we show, both computationally and experimentally, that NKG2A-HLA-E interaction between NK and tumor cells represents an important new therapeutic axis in PM, especially for epithelioid cases.

A PRMT5-ZNF326 axis mediates innate immune activation upon replication stress.

DNA replication stress (RS) is a widespread phenomenon in carcinogenesis, causing genomic instability and extensive chromatin alterations. DNA damage leads to activation of innate immune signaling, but little is known about transcriptional regulators mediating such signaling upon RS. Using a chemical screen, we identified protein arginine methyltransferase 5 (PRMT5) as a key mediator of RS-dependent induction of interferon-stimulated genes (ISGs). This response is also associated with reactivation of endogenous retroviruses (ERVs). Using quantitative mass spectrometry, we identify proteins with PRMT5-dependent symmetric dimethylarginine (SDMA) modification induced upon RS. Among these, we show that PRMT5 targets and modulates the activity of ZNF326, a zinc finger protein essential for ISG response. Our data demonstrate a role for PRMT5-mediated SDMA in the context of RS-induced transcriptional induction, affecting physiological homeostasis and cancer therapy.

Bat pluripotent stem cells reveal unusual entanglement between host and viruses.

Bats are distinctive among mammals due to their ability to fly, use laryngeal echolocation, and tolerate viruses. However, there are currently no reliable cellular models for studying bat biology or their response to viral infections. Here, we created induced pluripotent stem cells (iPSCs) from two species of bats: the wild greater horseshoe bat (Rhinolophus ferrumequinum) and the greater mouse-eared bat (Myotis myotis). The iPSCs from both bat species showed similar characteristics and had a gene expression profile resembling that of cells attacked by viruses. They also had a high number of endogenous viral sequences, particularly retroviruses. These results suggest that bats have evolved mechanisms to tolerate a large load of viral sequences and may have a more intertwined relationship with viruses than previously thought. Further study of bat iPSCs and their differentiated progeny will provide insights into bat biology, virus host relationships, and the molecular basis of bats' special traits.

Hedgehog-interacting protein acts in the habenula to regulate nicotine intake.

Hedgehog-interacting protein (HHIP) sequesters Hedgehog ligands to repress Smoothened (SMO)-mediated recruitment of the GLI family of transcription factors. Allelic variation in HHIP confers risk of chronic obstructive pulmonary disease and other smoking-related lung diseases, but underlying mechanisms are unclear. Using single-cell and cell-type-specific translational profiling, we show that HHIP expression is highly enriched in medial habenula (MHb) neurons, particularly MHb cholinergic neurons that regulate aversive behavioral responses to nicotine. HHIP deficiency dysregulated the expression of genes involved in cholinergic signaling in the MHb and disrupted the function of nicotinic acetylcholine receptors (nAChRs) through a PTCH-1/cholesterol-dependent mechanism. Further, CRISPR/Cas9-mediated genomic cleavage of the Hhip gene in MHb neurons enhanced the motivational properties of nicotine in mice. These findings suggest that HHIP influences vulnerability to smoking-related lung diseases in part by regulating the actions of nicotine on habenular aversion circuits.

NKG2A and HLA-E define an alternative immune checkpoint axis in bladder cancer.

Programmed cell death protein 1 (PD-1)/programmed death-ligand 1 (PD-L1)-blockade immunotherapies have limited efficacy in the treatment of bladder cancer. Here, we show that NKG2A associates with improved survival and responsiveness to PD-L1 blockade immunotherapy in bladder tumors that have high abundance of CD8+ T cells. In bladder tumors, NKG2A is acquired on CD8+ T cells later than PD-1 as well as other well-established immune checkpoints. NKG2A+ PD-1+ CD8+ T cells diverge from classically defined exhausted T cells through their ability to react to human leukocyte antigen (HLA) class I-deficient tumors using T cell receptor (TCR)-independent innate-like mechanisms. HLA-ABC expression by bladder tumors is progressively diminished as disease progresses, framing the importance of targeting TCR-independent anti-tumor functions. Notably, NKG2A+ CD8+ T cells are inhibited when HLA-E is expressed by tumors and partly restored upon NKG2A blockade in an HLA-E-dependent manner. Overall, our study provides a framework for subsequent clinical trials combining NKG2A blockade with other T cell-targeted immunotherapies, where tumors express higher levels of HLA-E.

Incident Cancer Risk and Signatures Among Older MUTYH Carriers: Analysis of Population-Based and Genomic Cohorts.

MUTYH carriers have an increased colorectal cancer risk in case-control studies, with loss of heterozygosity (LOH) as the presumed mechanism. We evaluated cancer risk among carriers in a prospective, population-based cohort of older adults. In addition, we assessed if cancers from carriers demonstrated mutational signatures (G:C>T:A transversions) associated with early LOH. We calculated incident risk of cancer and colorectal cancer among 13,131 sequenced study participants of the ASPirin in Reducing Events in the Elderly cohort, stratified by sex and adjusting for age, smoking, alcohol use, BMI, polyp history, history of cancer, and aspirin use. MUTYH carriers were identified among 13,033 participants in The Cancer Genome Atlas and International Cancer Genome Consortium, and somatic signatures of cancers were analyzed. Male MUTYH carriers demonstrated an increased risk for overall cancer incidence [multivariable HR, 1.66; 95% confidence interval (CI), 1.03-2.68; P = 0.038] driven by increased colorectal cancer incidence (multivariable HR, 3.55; 95% CI, 1.42-8.78; P = 0.007), as opposed to extracolonic cancer incidence (multivariable HR, 1.40; 95% CI, 0.81-2.44; P = 0.229). Female carriers did not demonstrate increased risk of cancer, colorectal cancer, or extracolonic cancers. Analysis of mutation signatures from cancers of MUTYH carriers revealed no significant contribution toward early mutagenesis from widespread G:C>T:A transversions among gastrointestinal epithelial cancers. Among cancers from carriers, somatic transversions associated with base-excision repair deficiency are uncommon, suggestive of diverse mechanisms of carcinogenesis in carriers compared with those who inherit biallelic MUTYH mutations. PREVENTION RELEVANCE: Despite absence of loss of heterozygosity in colorectal cancers, elderly male MUTYH carriers appeared to be at increased of colorectal cancer.

Dietary fructose improves intestinal cell survival and nutrient absorption.

Fructose consumption is linked to the rising incidence of obesity and cancer, which are two of the leading causes of morbidity and mortality globally1,2. Dietary fructose metabolism begins at the epithelium of the small intestine, where fructose is transported by glucose transporter type 5 (GLUT5; encoded by SLC2A5) and phosphorylated by ketohexokinase to form fructose 1-phosphate, which accumulates to high levels in the cell3,4. Although this pathway has been implicated in obesity and tumour promotion, the exact mechanism that drives these pathologies in the intestine remains unclear. Here we show that dietary fructose improves the survival of intestinal cells and increases intestinal villus length in several mouse models. The increase in villus length expands the surface area of the gut and increases nutrient absorption and adiposity in mice that are fed a high-fat diet. In hypoxic intestinal cells, fructose 1-phosphate inhibits the M2 isoform of pyruvate kinase to promote cell survival5-7. Genetic ablation of ketohexokinase or stimulation of pyruvate kinase prevents villus elongation and abolishes the nutrient absorption and tumour growth that are induced by feeding mice with high-fructose corn syrup. The ability of fructose to promote cell survival through an allosteric metabolite thus provides additional insights into the excess adiposity generated by a Western diet, and a compelling explanation for the promotion of tumour growth by high-fructose corn syrup.

Myeloid Cell-associated Resistance to PD-1/PD-L1 Blockade in Urothelial Cancer Revealed Through Bulk and Single-cell RNA Sequencing.

PURPOSE: To define dominant molecular and cellular features associated with PD-1/PD-L1 blockade resistance in metastatic urothelial cancer. EXPERIMENTAL DESIGN: We pursued an unbiased approach using bulk RNA sequencing data from two clinical trials to discover (IMvigor 210) and validate (CheckMate 275) pretreatment molecular features associated with resistance to PD-1/PD-L1 blockade in metastatic urothelial cancer. We then generated single-cell RNA sequencing (scRNA-seq) data from muscle-invasive bladder cancer specimens to dissect the cellular composition underlying the identified gene signatures. RESULTS: We identified an adaptive immune response gene signature associated with response and a protumorigenic inflammation gene signature associated with resistance to PD-1/PD-L1 blockade. The adaptive immune response:protumorigenic inflammation signature expression ratio, coined the 2IR score, best correlated with clinical outcomes, and was externally validated. Mapping these bulk gene signatures onto scRNA-seq data uncovered their underlying cellular diversity, with prominent expression of the protumorigenic inflammation signature by myeloid phagocytic cells. However, heterogeneity in expression of adaptive immune and protumorigenic inflammation genes was observed among single myeloid phagocytic cells, quantified as the myeloid single cell immune:protumorigenic inflammation ratio (Msc2IR) score. Single myeloid phagocytic cells with low Msc2IR scores demonstrated upregulation of proinflammatory cytokines/chemokines and downregulation of antigen presentation genes, were unrelated to M1 versus M2 polarization, and were enriched in pretreatment blood samples from patients with PD-L1 blockade-resistant metastatic urothelial cancer. CONCLUSIONS: The balance of adaptive immunity and protumorigenic inflammation in individual tumor microenvironments is associated with PD-1/PD-L1 resistance in urothelial cancer with the latter linked to a proinflammatory cellular state of myeloid phagocytic cells detectable in tumor and blood.See related commentary by Drake, p. 4139.

A reference profile-free deconvolution method to infer cancer cell-intrinsic subtypes and tumor-type-specific stromal profiles.

BACKGROUND: Patient stratification based on molecular subtypes is an important strategy for cancer precision medicine. Deriving clinically informative cancer molecular subtypes from transcriptomic data generated on whole tumor tissue samples is a non-trivial task, especially given the various non-cancer cellular elements intertwined with cancer cells in the tumor microenvironment. METHODS: We developed a computational deconvolution method, DeClust, that stratifies patients into subtypes based on cancer cell-intrinsic signals identified by distinguishing cancer-type-specific signals from non-cancer signals in bulk tumor transcriptomic data. DeClust differs from most existing methods by directly incorporating molecular subtyping of solid tumors into the deconvolution process and outputting molecular subtype-specific tumor reference profiles for the cohort rather than individual tumor profiles. In addition, DeClust does not require reference expression profiles or signature matrices as inputs and estimates cancer-type-specific microenvironment signals from bulk tumor transcriptomic data. RESULTS: DeClust was evaluated on both simulated data and 13 solid tumor datasets from The Cancer Genome Atlas (TCGA). DeClust performed among the best, relative to existing methods, for estimation of cellular composition. Compared to molecular subtypes reported by TCGA or other similar approaches, the subtypes generated by DeClust had higher correlations with cancer-intrinsic genomic alterations (e.g., somatic mutations and copy number variations) and lower correlations with tumor purity. While DeClust-identified subtypes were not more significantly associated with survival in general, DeClust identified a poor prognosis subtype of clear cell renal cancer, papillary renal cancer, and lung adenocarcinoma, all of which were characterized by CDKN2A deletions. As a reference profile-free deconvolution method, the tumor-type-specific stromal profiles and cancer cell-intrinsic subtypes generated by DeClust were supported by single-cell RNA sequencing data. CONCLUSIONS: DeClust is a useful tool for cancer cell-intrinsic molecular subtyping of solid tumors. DeClust subtypes, together with the tumor-type-specific stromal profiles generated by this pan-cancer study, may lead to mechanistic and clinical insights across multiple tumor types.

Microevolution in response to transient heme-iron restriction enhances intracellular bacterial community development and persistence.

Bacterial pathogens must sense, respond and adapt to a myriad of dynamic microenvironmental stressors to survive. Adaptation is key for colonization and long-term ability to endure fluctuations in nutrient availability and inflammatory processes. We hypothesize that strains adapted to survive nutrient deprivation are more adept for colonization and establishment of chronic infection. In this study, we detected microevolution in response to transient nutrient limitation through mutation of icc. The mutation results in decreased 3',5'-cyclic adenosine monophosphate phosphodiesterase activity in nontypeable Haemophilus influenzae (NTHI). In a preclinical model of NTHI-induced otitis media (OM), we observed a significant decrease in the recovery of effusion from ears infected with the icc mutant strain. Clinically, resolution of OM coincides with the clearance of middle ear fluid. In contrast to this clinical paradigm, we observed that the icc mutant strain formed significantly more intracellular bacterial communities (IBCs) than the parental strain early during experimental OM. Although the number of IBCs formed by the parental strain was low at early stages of OM, we observed a significant increase at later stages that coincided with absence of recoverable effusion, suggesting the presence of a mucosal reservoir following resolution of clinical disease. These data provide the first insight into NTHI microevolution during nutritional limitation and provide the first demonstration of IBCs in a preclinical model of chronic OM.

Coupling of Single Molecule, Long Read Sequencing with IMGT/HighV-QUEST Analysis Expedites Identification of SIV gp140-Specific Antibodies from scFv Phage Display Libraries.

The simian immunodeficiency virus (SIV)/macaque model of human immunodeficiency virus (HIV)/acquired immunodeficiency syndrome pathogenesis is critical for furthering our understanding of the role of antibody responses in the prevention of HIV infection, and will only increase in importance as macaque immunoglobulin (IG) gene databases are expanded. We have previously reported the construction of a phage display library from a SIV-infected rhesus macaque (Macaca mulatta) using oligonucleotide primers based on human IG gene sequences. Our previous screening relied on Sanger sequencing, which was inefficient and generated only a few dozen sequences. Here, we re-analyzed this library using single molecule, real-time (SMRT) sequencing on the Pacific Biosciences (PacBio) platform to generate thousands of highly accurate circular consensus sequencing (CCS) reads corresponding to full length single chain fragment variable. CCS data were then analyzed through the international ImMunoGeneTics information system® (IMGT®)/HighV-QUEST (www.imgt.org) to identify variable genes and perform statistical analyses. Overall the library was very diverse, with 2,569 different IMGT clonotypes called for the 5,238 IGHV sequences assigned to an IMGT clonotype. Within the library, SIV-specific antibodies represented a relatively limited number of clones, with only 135 different IMGT clonotypes called from 4,594 IGHV-assigned sequences. Our data did confirm that the IGHV4 and IGHV3 gene usage was the most abundant within the rhesus antibodies screened, and that these genes were even more enriched among SIV gp140-specific antibodies. Although a broad range of VH CDR3 amino acid (AA) lengths was observed in the unpanned library, the vast majority of SIV gp140-specific antibodies demonstrated a more uniform VH CDR3 length (20 AA). This uniformity was far less apparent when VH CDR3 were classified according to their clonotype (range: 9-25 AA), which we believe is more relevant for specific antibody identification. Only 174 IGKV and 588 IGLV clonotypes were identified within the VL sequences associated with SIV gp140-specific VH. Together, these data strongly suggest that the combination of SMRT sequencing with the IMGT/HighV-QUEST querying tool will facilitate and expedite our understanding of polyclonal antibody responses during SIV infection and may serve to rapidly expand the known scope of macaque V genes utilized during these responses.

Characterization of drug-induced splicing complexity in prostate cancer cell line using long read technology.

We characterize the transcriptional splicing landscape of a prostate cancer cell line treated with a previously identified synergistic drug combination. We use a combination of third generation long-read RNA sequencing technology and short-read RNAseq to create a high-fidelity map of expressed isoforms and fusions to quantify splicing events triggered by treatment. We find strong evidence for drug-induced, coherent splicing changes which disrupt the function of oncogenic proteins, and detect novel transcripts arising from previously unreported fusion events.

Identification of Differentially Expressed Splice Variants by the Proteogenomic Pipeline Splicify.

Proteogenomics, i.e. comprehensive integration of genomics and proteomics data, is a powerful approach identifying novel protein biomarkers. This is especially the case for proteins that differ structurally between disease and control conditions. As tumor development is associated with aberrant splicing, we focus on this rich source of cancer specific biomarkers. To this end, we developed a proteogenomic pipeline, Splicify, which can detect differentially expressed protein isoforms. Splicify is based on integrating RNA massive parallel sequencing data and tandem mass spectrometry proteomics data to identify protein isoforms resulting from differential splicing between two conditions. Proof of concept was obtained by applying Splicify to RNA sequencing and mass spectrometry data obtained from colorectal cancer cell line SW480, before and after siRNA-mediated downmodulation of the splicing factors SF3B1 and SRSF1. These analyses revealed 2172 and 149 differentially expressed isoforms, respectively, with peptide confirmation upon knock-down of SF3B1 and SRSF1 compared with their controls. Splice variants identified included RAC1, OSBPL3, MKI67, and SYK. One additional sample was analyzed by PacBio Iso-Seq full-length transcript sequencing after SF3B1 downmodulation. This analysis verified the alternative splicing identified by Splicify and in addition identified novel splicing events that were not represented in the human reference genome annotation. Therefore, Splicify offers a validated proteogenomic data analysis pipeline for identification of disease specific protein biomarkers resulting from mRNA alternative splicing. Splicify is publicly available on GitHub (https://github.com/NKI-TGO/SPLICIFY) and suitable to address basic research questions using pre-clinical model systems as well as translational research questions using patient-derived samples, e.g. allowing to identify clinically relevant biomarkers.