An atlas of epithelial cell states and plasticity in lung adenocarcinoma.

Understanding the cellular processes that underlie early lung adenocarcinoma (LUAD) development is needed to devise intervention strategies1. Here we studied 246,102 single epithelial cells from 16 early-stage LUADs and 47 matched normal lung samples. Epithelial cells comprised diverse normal and cancer cell states, and diversity among cancer cells was strongly linked to LUAD-specific oncogenic drivers. KRAS mutant cancer cells showed distinct transcriptional features, reduced differentiation and low levels of aneuploidy. Non-malignant areas surrounding human LUAD samples were enriched with alveolar intermediate cells that displayed elevated KRT8 expression (termed KRT8+ alveolar intermediate cells (KACs) here), reduced differentiation, increased plasticity and driver KRAS mutations. Expression profiles of KACs were enriched in lung precancer cells and in LUAD cells and signified poor survival. In mice exposed to tobacco carcinogen, KACs emerged before lung tumours and persisted for months after cessation of carcinogen exposure. Moreover, they acquired Kras mutations and conveyed sensitivity to targeted KRAS inhibition in KAC-enriched organoids derived from alveolar type 2 (AT2) cells. Last, lineage-labelling of AT2 cells or KRT8+ cells following carcinogen exposure showed that KACs are possible intermediates in AT2-to-tumour cell transformation. This study provides new insights into epithelial cell states at the root of LUAD development, and such states could harbour potential targets for prevention or intervention.

Androgen receptor blockade promotes response to BRAF/MEK-targeted therapy.

Treatment with therapy targeting BRAF and MEK (BRAF/MEK) has revolutionized care in melanoma and other cancers; however, therapeutic resistance is common and innovative treatment strategies are needed1,2. Here we studied a group of patients with melanoma who were treated with neoadjuvant BRAF/MEK-targeted therapy ( NCT02231775 , n = 51) and observed significantly higher rates of major pathological response (MPR; ≤10% viable tumour at resection) and improved recurrence-free survival (RFS) in female versus male patients (MPR, 66% versus 14%, P = 0.001; RFS, 64% versus 32% at 2 years, P = 0.021). The findings were validated in several additional cohorts2-4 of patients with unresectable metastatic melanoma who were treated with BRAF- and/or MEK-targeted therapy (n = 664 patients in total), demonstrating improved progression-free survival and overall survival in female versus male patients in several of these studies. Studies in preclinical models demonstrated significantly impaired anti-tumour activity in male versus female mice after BRAF/MEK-targeted therapy (P = 0.006), with significantly higher expression of the androgen receptor in tumours of male and female BRAF/MEK-treated mice versus the control (P = 0.0006 and P = 0.0025). Pharmacological inhibition of androgen receptor signalling improved responses to BRAF/MEK-targeted therapy in male and female mice (P = 0.018 and P = 0.003), whereas induction of androgen receptor signalling (through testosterone administration) was associated with a significantly impaired response to BRAF/MEK-targeted therapy in male and female patients (P = 0.021 and P < 0.0001). Together, these results have important implications for therapy.

B cells and tertiary lymphoid structures promote immunotherapy response.

Treatment with immune checkpoint blockade (ICB) has revolutionized cancer therapy. Until now, predictive biomarkers1-10 and strategies to augment clinical response have largely focused on the T cell compartment. However, other immune subsets may also contribute to anti-tumour immunity11-15, although these have been less well-studied in ICB treatment16. A previously conducted neoadjuvant ICB trial in patients with melanoma showed via targeted expression profiling17 that B cell signatures were enriched in the tumours of patients who respond to treatment versus non-responding patients. To build on this, here we performed bulk RNA sequencing and found that B cell markers were the most differentially expressed genes in the tumours of responders versus non-responders. Our findings were corroborated using a computational method (MCP-counter18) to estimate the immune and stromal composition in this and two other ICB-treated cohorts (patients with melanoma and renal cell carcinoma). Histological evaluation highlighted the localization of B cells within tertiary lymphoid structures. We assessed the potential functional contributions of B cells via bulk and single-cell RNA sequencing, which demonstrate clonal expansion and unique functional states of B cells in responders. Mass cytometry showed that switched memory B cells were enriched in the tumours of responders. Together, these data provide insights into the potential role of B cells and tertiary lymphoid structures in the response to ICB treatment, with implications for the development of biomarkers and therapeutic targets.

TIGIT is the central player in T-cell suppression associated with CAR T-cell relapse in mantle cell lymphoma.

BACKGROUND: Chimeric antigen receptor (CAR) T-cell therapy using brexucabtagene autoleucel (BA) induces remission in many patients with mantle cell lymphoma (MCL), and BA is the only CAR T-cell therapy approved by the FDA for MCL. However, development of relapses to BA is recognized with poor patient outcomes. Multiple CAR T-cell therapies have been approved for other lymphomas and the resistance mechanisms have been investigated. However, the mechanisms underlying BA relapse in MCL have not been investigated and whether any previously reported resistance mechanisms apply to BA-relapsed patients with MCL is unknown. METHODS: To interrogate BA resistance mechanisms in MCL, we performed single-cell RNA sequencing on 39 longitudinally collected samples from 15 BA-treated patients, and multiplex cytokine profiling on 80 serial samples from 20 patients. RESULTS: We demonstrate that after BA relapse, the proportion of T cells, especially cytotoxic T cells (CTLs), decreased among non-tumor cells, while the proportion of myeloid cells correspondingly increased. TIGIT, LAG3, and CD96 were the predominant checkpoint molecules expressed on exhausted T cells and CTLs; only TIGIT was significantly increased after relapse. CTLs expanded during remission, and then contracted during relapse with upregulated TIGIT expression. Tumor cells also acquired TIGIT expression after relapse, leading to the enhanced interaction of tumor cell TIGIT with monocyte CD155/PVR. In myeloid cells, post-relapse HLA-II expression was reduced relative to pretreatment and during remission. Myeloid-derived suppressor cells (MDSCs) were enriched after relapse with elevated expression of activation markers, including CLU (clusterin) and VCAN (versican). Extracellular chemokines (CCL4, CXCL9, CXCL13), soluble checkpoint inhibitors (sPD-L1, sTIM3, s4-1BB), and soluble receptors (sIL-2R, sTNFRII) were decreased during remission but elevated after relapse. CONCLUSIONS: Our data demonstrate that multiple tumor-intrinsic and -extrinsic factors are associated with T-cell suppression and BA relapse. Among these, TIGIT appears to be the central player given its elevated expression after BA relapse in not only CTLs but also MCL cells. The acquisition of TIGIT expression on tumor cells is MCL-specific and has not been reported in other CAR T-treated diseases. Together, our data suggest that co-targeting TIGIT may prevent CAR T relapses and thus promote long-term progression-free survival in MCL patients.

Integrated genomic profiling and modelling for risk stratification in patients with advanced oesophagogastric adenocarcinoma.

OBJECTIVE: Prognosis of patients with advanced oesophagogastric adenocarcinoma (mEGAC) is poor and molecular determinants of shorter or longer overall survivors are lacking. Our objective was to identify molecular features and develop a prognostic model by profiling the genomic features of patients with mEGAC with widely varying outcomes. DESIGN: We profiled 40 untreated mEGACs (20 shorter survivors <13 months and 20 longer survivors >36 months) with whole-exome sequencing (WES) and RNA sequencing and performed an integrated analysis of exome, transcriptome, immune profile and pathological phenotypes to identify the molecular determinants, developing an integrated model for prognosis and comparison with The Cancer Genome Atlas (TCGA) cohorts. RESULTS: KMT2C alterations were exclusively observed in shorter survivors together with high level of intratumour heterogeneity and complex clonal architectures, whereas the APOBEC mutational signatures were significantly enriched in longer survivors. Notably, the loss of heterozygosity in chromosome 4 (Chr4) was associated with shorter survival and 'cold' immune phenotype characterised by decreased B, CD8, natural killer cells and interferon-gamma responses. Unsupervised transcriptomic clustering revealed a shorter survivor subtype with distinct expression features (eg, upregulated druggable targets JAK2, MAP3K13 and MECOM). An integrated model was then built based on clinical variables and the identified molecular determinants, which significantly segregated shorter and longer survivors. All the above features and the integrated model have been validated independently in multiple TCGA cohorts. CONCLUSION: This study discovered novel molecular features prognosticating overall survival in patients with mEGAC and identified potential novel targets in shorter survivors.

Multiplex profiling of peritoneal metastases from gastric adenocarcinoma identified novel targets and molecular subtypes that predict treatment response.

OBJECTIVE: Peritoneal carcinomatosis (PC) occurs frequently in patients with gastric adenocarcinoma (GAC) and confers a poor prognosis. Multiplex profiling of primary GACs has been insightful but the underpinnings of PC's development/progression remain largely unknown. We characterised exome/transcriptome/immune landscapes of PC cells from patients with GAC aiming to identify novel therapeutic targets. DESIGN: We performed whole-exome sequencing (WES) and whole transcriptome sequencing (RNA-seq) on 44 PC specimens (43 patients with PC) including an integrative analysis of WES, RNA-seq, immune profile, clinical and pathological phenotypes to dissect the molecular pathogenesis, identifying actionable targets and/or biomarkers and comparison with TCGA primary GACs. RESULTS: We identified distinct alterations in PC versus primary GACs, such as more frequent CDH1 and TAF1 mutations, 6q loss and chr19 gain. Alterations associated with aggressive PC phenotypes emerged with increased mutations in TP53, CDH1, TAF1 and KMT2C, higher level of 'clock-like' mutational signature, increase in whole-genome doublings, chromosomal instability (particularly, copy number losses), reprogrammed microenvironment, enriched cell cycle pathways, MYC activation and impaired immune response. Integrated analysis identified two main molecular subtypes: 'mesenchymal-like' and 'epithelial-like' with discriminating response to chemotherapy (31% vs 71%). Patients with the less responsive 'mesenchymal-like' subtype had high expression of immune checkpoint T-Cell Immunoglobulin And Mucin Domain-Containing Protein 3 (TIM-3), its ligand galectin-9, V-domain Ig suppressor of T cell activation (VISTA) and transforming growth factor-ß as potential therapeutic immune targets. CONCLUSIONS: We have uncovered the unique mutational landscape, copy number alteration and gene expression profile of PC cells and defined PC molecular subtypes, which correlated with PC therapy resistance/response. Novel targets and immune checkpoint proteins have been identified with a potential to be translated into clinics.

Efficacy of venetoclax in high risk relapsed mantle cell lymphoma (MCL) - outcomes and mutation profile from venetoclax resistant MCL patients.

Venetoclax is effective in relapsed patients with mantle cell lymphoma (MCL). Mechanisms of resistance to venetoclax in MCL are poorly understood. We describe the clinical outcomes and genomic characteristics of 24 multiply relapsed patients (median of five prior lines of therapy) who received venetoclax-based therapies; 67% had progressed on BTK inhibitors (BTKi) and 54% had blastoid or pleomorphic histology. Median follow up after venetoclax treatment was 17 months. The overall response rate was 50% and complete response (CR) rate was 21%, 16 patients had progressed and 15 died. The median progression free, overall and post venetoclax survival were 8, 13.5 and 7.3 months respectively. Whole-exome sequencing (WES) was performed on samples collected from seven patients (including five pairs; before starting venetoclax and after progression on venetoclax). The SMARCA4 and BCL2 alterations were noted only after progression, while TP53, CDKN2A, KMT2D, CELSR3, CCND1, NOTCH2 and ATM were altered 2-4-fold more frequently after progression. In two patients with serial samples, we demonstrated clonal evolution of novel SMARCA4 and KMT2C/D mutations at progression. Mutation dynamics in venetoclax resistant MCL is demonstrated. Our data indicates that venetoclax resistance in MCL is predominantly associated with non-BCL2 gene mutations. Further studies are ongoing in MCL patients to evaluate the efficacy of venetoclax in combination with other agents and understand the biology of venetoclax resistance in MCL.

TissGDB: tissue-specific gene database in cancer.

Tissue-specific gene expression is critical in understanding biological processes, physiological conditions, and disease. The identification and appropriate use of tissue-specific genes (TissGenes) will provide important insights into disease mechanisms and organ-specific therapeutic targets. To better understand the tissue-specific features for each cancer type and to advance the discovery of clinically relevant genes or mutations, we built TissGDB (Tissue specific Gene DataBase in cancer) available at http://zhaobioinfo.org/TissGDB. We collected and curated 2461 tissue specific genes (TissGenes) across 22 tissue types that matched the 28 cancer types of The Cancer Genome Atlas (TCGA) from three representative tissue-specific gene expression resources: The Human Protein Atlas (HPA), Tissue-specific Gene Expression and Regulation (TiGER), and Genotype-Tissue Expression (GTEx). For these 2461 TissGenes, we performed gene expression, somatic mutation, and prognostic marker-based analyses across 28 cancer types using TCGA data. Our analyses identified hundreds of TissGenes, including genes that universally kept or lost tissue-specific gene expression, with other features: cancer type-specific isoform expression, fusion with oncogenes or tumor suppressor genes, and markers for protective or risk prognosis. TissGDB provides seven categories of annotations: TissGeneSummary, TissGeneExp, TissGene-miRNA, TissGeneMut, TissGeneNet, TissGeneProg, TissGeneClin.

LncRNA PVT1 up-regulation is a poor prognosticator and serves as a therapeutic target in esophageal adenocarcinoma.

BACKGROUND: PVT1 has emerged as an oncogene in many tumor types. However, its role in Barrett's esophagus (BE) and esophageal adenocarcinoma (EAC) is unknown. The aim of this study was to assess the role of PVT1 in BE/EAC progression and uncover its therapeutic value against EAC. METHODS: PVT1 expression was assessed by qPCR in normal, BE, and EAC tissues and statistical analysis was performed to determine the association of PVT1 expression and EAC (stage, metastases, and survival). PVT1 antisense oligonucleotides (ASOs) were tested for their antitumor activity in vitro and in vivo. RESULTS: PVT1 expression was up-regulated in EACs compared with paired BEs, and normal esophageal tissues. High expression of PVT1 was associated with poor differentiation, lymph node metastases, and shorter survival. Effective knockdown of PVT1 in EAC cells using PVT1 ASOs resulted in decreased cell proliferation, invasion, colony formation, tumor sphere formation, and reduced proportion of ALDH1A1+ cells. Mechanistically, we discovered mutual regulation of PVT1 and YAP1 in EAC cells. Inhibition of PVT1 by PVT1 ASOs suppressed YAP1 expression through increased phosphor-LATS1and phosphor-YAP1 while knockout of YAP1 in EAC cells significantly suppressed PVT1 levels indicating a positive regulation of PVT1 by YAP1. Most importantly, we found that targeting both PVT1 and YAP1 using their specific ASOs led to better antitumor activity in vitro and in vivo. CONCLUSIONS: Our results provide strong evidence that PVT1 confers an aggressive phenotype to EAC and is a poor prognosticator. Combined targeting of PVT1 and YAP1 provided the highest therapeutic index and represents a novel therapeutic strategy.

SZGR 2.0: a one-stop shop of schizophrenia candidate genes.

SZGR 2.0 is a comprehensive resource of candidate variants and genes for schizophrenia, covering genetic, epigenetic, transcriptomic, translational and many other types of evidence. By systematic review and curation of multiple lines of evidence, we included almost all variants and genes that have ever been reported to be associated with schizophrenia. In particular, we collected ∼4200 common variants reported in genome-wide association studies, ∼1000 de novo mutations discovered by large-scale sequencing of family samples, 215 genes spanning rare and replication copy number variations, 99 genes overlapping with linkage regions, 240 differentially expressed genes, 4651 differentially methylated genes and 49 genes as antipsychotic drug targets. To facilitate interpretation, we included various functional annotation data, especially brain eQTL, methylation QTL, brain expression featured in deep categorization of brain areas and developmental stages and brain-specific promoter and enhancer annotations. Furthermore, we conducted cross-study, cross-data type and integrative analyses of the multidimensional data deposited in SZGR 2.0, and made the data and results available through a user-friendly interface. In summary, SZGR 2.0 provides a one-stop shop of schizophrenia variants and genes and their function and regulation, providing an important resource in the schizophrenia and other mental disease community. SZGR 2.0 is available at https://bioinfo.uth.edu/SZGR/.

Unique protein expression signatures of survival time in kidney renal clear cell carcinoma through a pan-cancer screening.

BACKGROUND: In 2016, it is estimated that there will be 62,700 new cases of kidney cancer in the United States, and 14,240 patients will die from the disease. Because the incidence of kidney renal clear cell carcinoma (KIRC), the most common type of kidney cancer, is expected to continue to increase in the US, there is an urgent need to find effective diagnostic biomarkers for KIRC that could help earlier detection of and customized treatment strategies for the disease. Accordingly, in this study we systematically investigated KIRC's prognostic biomarkers for survival using the reverse phase protein array (RPPA) data and the high throughput sequencing data from The Cancer Genome Atlas (TCGA). RESULTS: With comprehensive data available in TCGA, we systematically screened protein expression based survival biomarkers in 10 major cancer types, among which KIRC presented many protein prognostic biomarkers of survival time. This is in agreement with a previous report that expression level changes (mRNAs, microRNA and protein) may have a better performance for prognosis of KIRC. In this study, we also identified 52 prognostic genes for KIRC, many of which are involved in cell-cycle and cancer signaling, as well as 15 tumor-stage-specific prognostic biomarkers. Notably, we found fewer prognostic biomarkers for early-stage than for late-stage KIRC. Four biomarkers (the RPPA protein IDs: FASN, ACC1, Cyclin_B1 and Rad51) were found to be prognostic for survival based on both protein and mRNA expression data. CONCLUSIONS: Through pan-cancer screening, we found that many protein biomarkers were prognostic for patients' survival in KIRC. Stage-specific survival biomarkers in KIRC were also identified. Our study indicated that these protein biomarkers might have potential clinical value in terms of predicting survival in KIRC patients and developing individualized treatment strategies. Importantly, we found many biomarkers in KIRC at both the mRNA expression level and the protein expression level. These biomarkers shared a significant overlap, indicating that they were technically replicable.

Distinctive RNA expression profiles in blood associated with Alzheimer disease after accounting for white matter hyperintensities.

BACKGROUND: Defining the RNA transcriptome in Alzheimer Disease (AD) will help understand the disease mechanisms and provide biomarkers. Though the AD blood transcriptome has been studied, effects of white matter hyperintensities (WMH) were not considered. This study investigated the AD blood transcriptome and accounted for WMH. METHODS: RNA from whole blood was processed on whole-genome microarrays. RESULTS: A total of 293 probe sets were differentially expressed in AD versus controls, 5 of which were significant for WMH status. The 288 AD-specific probe sets classified subjects with 87.5% sensitivity and 90.5% specificity. They represented 188 genes of which 29 have been reported in prior AD blood and 89 in AD brain studies. Regulated blood genes included MMP9, MME (Neprilysin), TGFß1, CA4, OCLN, ATM, TGM3, IGFR2, NOV, RNF213, BMX, LRRN1, CAMK2G, INSR, CTSD, SORCS1, SORL1, and TANC2. CONCLUSIONS: RNA expression is altered in AD blood irrespective of WMH status. Some genes are shared with AD brain.

Proteogenomic landscape of gastric adenocarcinoma peritoneal metastases.

Advanced gastric adenocarcinoma (GAC) often leads to peritoneal carcinomatosis (PC) and is associated with very poor outcome. Here we report the comprehensive proteogenomic study of ascites derived cells from a prospective GAC cohort (n = 26 patients with peritoneal carcinomatosis, PC). A total of 16,449 proteins were detected from whole cell extracts (TCEs). Unsupervised hierarchical clustering resulted in three distinct groups that reflected extent of enrichment in tumor cells. Integrated analysis revealed enriched biological pathways and notably, some druggable targets (cancer-testis antigens, kinases, and receptors) that could be exploited to develop effective therapies and/or tumor stratifications. Systematic comparison of expression levels of proteins and mRNAs revealed special expression patterns of key therapeutics target notably high mRNA and low protein expression of HAVCR2 (TIM-3), and low mRNA but high protein expression of cancer-testis antigens CTAGE1 and CTNNA2. These results inform strategies to target GAC vulnerabilities.

Single cell clonotypic and transcriptional evolution of multiple myeloma precursor disease.

Multiple myeloma remains an incurable disease, and the cellular and molecular evolution from precursor conditions, including monoclonal gammopathy of undetermined significance and smoldering multiple myeloma, is incompletely understood. Here, we combine single-cell RNA and B cell receptor sequencing from fifty-two patients with myeloma precursors in comparison with myeloma and normal donors. Our comprehensive analysis reveals early genomic drivers of malignant transformation, distinct transcriptional features, and divergent clonal expansion in hyperdiploid versus non-hyperdiploid samples. Additionally, we observe intra-patient heterogeneity with potential therapeutic implications and identify distinct patterns of evolution from myeloma precursor disease to myeloma. We also demonstrate distinctive characteristics of the microenvironment associated with specific genomic changes in myeloma cells. These findings add to our knowledge about myeloma precursor disease progression, providing valuable insights into patient risk stratification, biomarker discovery, and possible clinical applications.

Evolution of immune and stromal cell states and ecotypes during gastric adenocarcinoma progression.

Understanding tumor microenvironment (TME) reprogramming in gastric adenocarcinoma (GAC) progression may uncover novel therapeutic targets. Here, we performed single-cell profiling of precancerous lesions, localized and metastatic GACs, identifying alterations in TME cell states and compositions as GAC progresses. Abundant IgA+ plasma cells exist in the premalignant microenvironment, whereas immunosuppressive myeloid and stromal subsets dominate late-stage GACs. We identified six TME ecotypes (EC1-6). EC1 is exclusive to blood, while EC4, EC5, and EC2 are highly enriched in uninvolved tissues, premalignant lesions, and metastases, respectively. EC3 and EC6, two distinct ecotypes in primary GACs, associate with histopathological and genomic characteristics, and survival outcomes. Extensive stromal remodeling occurs in GAC progression. High SDC2 expression in cancer-associated fibroblasts (CAFs) is linked to aggressive phenotypes and poor survival, and SDC2 overexpression in CAFs contributes to tumor growth. Our study provides a high-resolution GAC TME atlas and underscores potential targets for further investigation.

Pan-cancer T cell atlas links a cellular stress response state to immunotherapy resistance.

Tumor-infiltrating T cells offer a promising avenue for cancer treatment, yet their states remain to be fully characterized. Here we present a single-cell atlas of T cells from 308,048 transcriptomes across 16 cancer types, uncovering previously undescribed T cell states and heterogeneous subpopulations of follicular helper, regulatory and proliferative T cells. We identified a unique stress response state, TSTR, characterized by heat shock gene expression. TSTR cells are detectable in situ in the tumor microenvironment across various cancer types, mostly within lymphocyte aggregates or potential tertiary lymphoid structures in tumor beds or surrounding tumor edges. T cell states/compositions correlated with genomic, pathological and clinical features in 375 patients from 23 cohorts, including 171 patients who received immune checkpoint blockade therapy. We also found significantly upregulated heat shock gene expression in intratumoral CD4/CD8+ cells following immune checkpoint blockade treatment, particularly in nonresponsive tumors, suggesting a potential role of TSTR cells in immunotherapy resistance. Our well-annotated T cell reference maps, web portal and automatic alignment/annotation tool could provide valuable resources for T cell therapy optimization and biomarker discovery.