Cancer cells are embedded in the tumor microenvironment (TME), a complex ecosystem of stromal cells. Here, we present a 52,698-cell catalog of the TME transcriptome in human lung tumors at single-cell resolution, validated in independent samples where 40,250 additional cells were sequenced. By comparing with matching non-malignant lung samples, we reveal a highly complex TME that profoundly molds stromal cells. We identify 52 stromal cell subtypes, including novel subpopulations in cell types hitherto considered to be homogeneous, as well as transcription factors underlying their heterogeneity. For instance, we discover fibroblasts expressing different collagen sets, endothelial cells downregulating immune cell homing and genes coregulated with established immune checkpoint transcripts and correlating with T-cell activity. By assessing marker genes for these cell subtypes in bulk RNA-sequencing data from 1,572 patients, we illustrate how these correlate with survival, while immunohistochemistry for selected markers validates them as separate cellular entities in an independent series of lung tumors. Hence, in providing a comprehensive catalog of stromal cells types and by characterizing their phenotype and co-optive behavior, this resource provides deeper insights into lung cancer biology that will be helpful in advancing lung cancer diagnosis and therapy.

TERT-locus single nucleotide polymorphisms (SNPs) and leucocyte telomere measures are reportedly associated with risks of multiple cancers. Using the iCOGs chip, we analysed ~480 TERT-locus SNPs in breast (n=103,991), ovarian (n=39,774) and BRCA1 mutation carrier (11,705) cancer cases and controls. 53,724 participants have leucocyte telomere measures. Most associations cluster into three independent peaks. Peak 1 SNP rs2736108 minor allele associates with longer telomeres (P=5.8×10 −7 ), reduced estrogen receptor negative (ER-negative) (P=1.0×10 −8 ) and BRCA1 mutation carrier (P=1.1×10 −5 ) breast cancer risks, and altered promoter-assay signal. Peak 2 SNP rs7705526 minor allele associates with longer telomeres (P=2.3×10 −14 ), increased low malignant potential ovarian cancer risk (P=1.3×10 −15 ) and increased promoter activity. Peak 3 SNPs rs10069690 and rs2242652 minor alleles increase ER-negative (P=1.2×10 −12 ) and BRCA1 mutation carrier (P=1.6×10 −14 ) breast and invasive ovarian (P=1.3×10 −11 ) cancer risks, but not via altered telomere length. The cancer-risk alleles of rs2242652 and rs10069690 respectively increase silencing and generate a truncated TERT splicevariant.

The stromal compartment of the tumor microenvironment consists of a heterogeneous set of tissue-resident and tumor-infiltrating cells, which are profoundly moulded by cancer cells. An outstanding question is to what extent this heterogeneity is similar between cancers affecting different organs. Here, we profile 233,591 single cells from patients with lung, colorectal, ovary and breast cancer ( n = 36) and construct a pan-cancer blueprint of stromal cell heterogeneity using different single-cell RNA and protein-based technologies. We identify 68 stromal cell populations, of which 46 are shared between cancer types and 22 are unique. We also characterise each population phenotypically by highlighting its marker genes, transcription factors, metabolic activities and tissue-specific expression differences. Resident cell types are characterised by substantial tissue specificity, while tumor-infiltrating cell types are largely shared across cancer types. Finally, by applying the blueprint to melanoma tumors treated with checkpoint immunotherapy and identifying a naïve CD4 + T-cell phenotype predictive of response to checkpoint immunotherapy, we illustrate how it can serve as a guide to interpret scRNA-seq data. In conclusion, by providing a comprehensive blueprint through an interactive web server, we generate the first panoramic view on the shared complexity of stromal cells in different cancers.

Mutations or amplification of the MET proto-oncogene are involved in the pathogenesis of several tumours1-4, which rely on the constitutive engagement of this pathway for their growth and survival1,5. However, MET is expressed not only by cancer cells but also by tumour-associated stromal cells although its precise role in this compartment is not well characterized6-11. Here, we show that MET is required for neutrophil chemoattraction and cytotoxicity in response to its ligand HGF. Met deletion in neutrophils enhances tumour growth and metastasis. This phenotype correlates with reduced neutrophil infiltration to both primary tumour and metastatic site. Similarly, Met is necessary for neutrophil transudation during colitis, skin rash or peritonitis. Mechanistically, Met is induced by tumour-derived TNF-α or other inflammatory stimuli in both mouse and human neutrophils. This induction is instrumental for neutrophil transmigration across an activated endothelium and iNOS production upon HGF stimulation. Consequently, HGF/MET-dependent nitric oxide release by neutrophils promotes cancer cell killing, which abates tumour growth and metastasis. Following systemic administration of a MET kinase inhibitor, we prove that the therapeutic benefit of MET targeting in cancer cells is partly countered by the pro-tumoural effect rising from MET blockade in neutrophils. Our work identifies an unprecedented role of MET in neutrophils, suggests a potential “Achilles’ heel” of MET-targeted therapies in cancer, and supports the rationale for evaluating anti-MET drugs in certain inflammatory diseases.

How the innate and adaptive host immune system miscommunicate to worsen COVID-19 immunopathology has not been fully elucidated. Here, we perform single-cell deep-immune profiling of bronchoalveolar lavage (BAL) samples from 5 patients with mild and 26 with critical COVID-19 in comparison to BALs from non-COVID-19 pneumonia and normal lung. We use pseudotime inference to build T-cell and monocyte-to-macrophage trajectories and model gene expression changes along them. In mild COVID-19, CD8+ resident-memory (TRM) and CD4+ T-helper-17 (TH17) cells undergo active (presumably antigen-driven) expansion towards the end of the trajectory, and are characterized by good effector functions, while in critical COVID-19 they remain more naïve. Vice versa, CD4+ T-cells with T-helper-1 characteristics (TH1-like) and CD8+ T-cells expressing exhaustion markers (TEX-like) are enriched halfway their trajectories in mild COVID-19, where they also exhibit good effector functions, while in critical COVID-19 they show evidence of inflammation-associated stress at the end of their trajectories. Monocyte-to-macrophage trajectories show that chronic hyperinflammatory monocytes are enriched in critical COVID-19, while alveolar macrophages, otherwise characterized by anti-inflammatory and antigen-presenting characteristics, are depleted. In critical COVID-19, monocytes contribute to an ATP-purinergic signaling-inflammasome footprint that could enable COVID-19 associated fibrosis and worsen disease-severity. Finally, viral RNA-tracking reveals infected lung epithelial cells, and a significant proportion of neutrophils and macrophages that are involved in viral clearance.

Estrogen receptor (ER)-negative breast cancer shows a higher incidence in women of African ancestry compared to women of European ancestry. In search of common risk alleles for ER-negative breast cancer, we combined genome-wide association study (GWAS) data from women of African ancestry (1,004 ER-negative cases and 2,745 controls) and European ancestry (1,718 ER-negative cases and 3,670 controls), with replication testing conducted in an additional 2,292 ER-negative cases and 16,901 controls of European ancestry. We identified a common risk variant for ER-negative breast cancer at the TERT-CLPTM1L locus on chromosome 5p15 (rs10069690: per-allele odds ratio (OR) = 1.18 per allele, P = 1.0 × 10−10). The variant was also significantly associated with triple-negative (ER-negative, progesterone receptor (PR)-negative and human epidermal growth factor-2 (HER2)-negative) breast cancer (OR = 1.25, P = 1.1 × 10−9), particularly in younger women (<50 years of age) (OR = 1.48, P = 1.9 × 10−9). Our results identify a genetic locus associated with estrogen receptor negative breast cancer subtypes in multiple populations.

In the original article, the surname of co-author Wouter Everaerts was spelled incorrectly as "Everaert." It appears correctly in this Correction, and the error has been corrected online.

Metabolic reprogramming by oncogenic signals promotes cancer initiation and progression. The oncogene KRAS and tumor suppressor STK11, which encodes the kinase LKB1, regulate metabolism and are frequently mutated in non-small cell lung cancer (NSCLC). Concurrent KRAS mutation and LKB1 loss (KL) specifies aggressive oncological behavior1,2. We show that KL cells and tumors share metabolomic signatures of perturbed nitrogen handling. KL cells express the urea cycle enzyme carbamoyl phosphate synthetase-1 (CPS1), which produces carbamoyl phosphate (CP) in the mitochondria from ammonia and bicarbonate, initiating nitrogen disposal. CPS1 transcription is suppressed by LKB1 via AMPK, and CPS1 expression anticorrelates with LKB1 in human NSCLC. Silencing CPS1 in KL cells induces cell death and reduces tumor growth. Surprisingly, cell death results from pyrimidine depletion rather than ammonia toxicity, as CPS1 enables an unconventional pathway of nitrogen flow from ammonia into pyrimidines. CPS1 loss reduces the pyrimidine/purine ratio, compromises S-phase progression, and induces DNA polymerase stalling and DNA damage. Exogenous pyrimidines reverse DNA damage and rescue growth. The data indicate that the KL oncogenotype imposes a novel metabolic vulnerability related to exquisite dependence on a cross-compartmental pathway of pyrimidine metabolism in an aggressive subset of NSCLC.