(A) Overall study design. The numbers indicate sample counts in each group. (B) Frequency of CH variant detection in each gene in NSCLC samples before (light blue) and after (blue) ICI treatment. (C)…
We collected numbers indicate count of variants after each variant filtering scheme (all samples/excluding controls).
(A) Variant depth distributions from raw calls and putative clonal hematopoiesis of indeterminate potential (CHIP) variants. (B) Distribution of variant annotation from each variant filtering step.
(A) Gene distribution in before and after immune checkpoint inhibitor (ICI) treatment. Color indicates variant annotation. (B) Waterfall plot showing CHIP mutations from baseline, cohort samples.
(A–D) Discovery cohort (n = 100), (E–G) replication cohort (n = 180; see Methods for description). (A, E) CH prevalence by pathology (lung squamous cell carcinoma [LUSC] and lung adenocarcinoma …
CHIP status was displayed using (A) binned variant allele frequency (VAF) and (B) variant count per patient. Colors indicate binned VAF and variants count.
CHIP variant counts and major clinical parameters were presented. Color indicates the binned number of CHIP variants.
(A) Variant depth distribution from raw and CHIP variant. Colored bar indicates variants cutoff, with depth = 40. (B) Distribution of variant annotation from each variant filtering step. (C) All …
(A) Uniform manifold approximation and projection (UMAP) plot of scRNA-seq data from the discovery cohort (n = 63). (B) Effect of pathology (left) and variant allele frequency (VAF) (right) on cell …
(A) UMAP plot of scRNA-seq data, with CH variant allele frequency (VAF) bins. (B) Cell-type-specific marker gene expression from each cluster. (C) Cell composition plot from each sample.
Color represents normalized effect score (NES), and dot size represents adjusted p-values.
(A) Results of GSEA analysis in response to ICI and CHIP VAF bins, from before treatment samples. Color represents normalized effect score (NES), and dot size represents adjusted p-values. (B) …
(A) WGCNA dendrogram, derived from myeloid gene expressions. (B) Cluster-wise ‘black’ module score distribution. (C) Gene ontology (GO) enrichment analysis of the black module genes. The colored bar …
(A) Clustered, representative GRNs from scRNA-seq dataset. 1500 cells from each variant allele frequency (VAF) bin were presented for this visualization. Color represents GRN AUC scores. (B) …
Heatmap illustrating the numbers of cell–cell interactions using (A) high variant allele frequency (VAF) clonal hematopoiesis of indeterminate potential (CHIP) and (B) CHIP negative samples. (C) …
Survival curve using progression-free survival in clonal hematopoiesis of indeterminate potential (CHIP) negative (N = 71) and high-burden CHIP (N = 11), using survplot in R.
Hypoxia and TNF-ɑ signaling pathways across clusters. The bars represent the count of clusters in which each gene is included in the leading edge.
Supplementary tables provided in an excel file.
(a) Clinical information of non-small cell lung cancer (NSCLC) cohort with immune checkpoint inhibitor (ICI) treatment. (b) Information of control samples. (c) List of clonal hematopoiesis of indeterminate potential (CHIP) variants detected from targeted sequencing. (d) List of CHIP variants detected from whole exome sequencing (WES). (e) Identified DEGs from annotated clusters of single-cell RNA sequencing (scRNA-seq). (f) Gene set enrichment analysis (GSEA) results related to myeloid lineage, related to Figure 3C. (g) Identified DEGs from each cluster by variant allele frequency (VAF) bin. (h) Gene-module annotation from weighted correlation network analysis (WGCNA).