KLF5 governs sphingolipid metabolism and barrier function of the skin.

Stem cells are fundamental units of tissue remodeling whose functions are dictated by lineage-specific transcription factors. Home to epidermal stem cells and their upward-stratifying progenies, skin relies on its secretory functions to form the outermost protective barrier, of which a transcriptional orchestrator has been elusive. KLF5 is a Krüppel-like transcription factor broadly involved in development and regeneration whose lineage specificity, if any, remains unclear. Here we report KLF5 specifically marks the epidermis, and its deletion leads to skin barrier dysfunction in vivo. Lipid envelopes and secretory lamellar bodies are defective in KLF5-deficient skin, accompanied by preferential loss of complex sphingolipids. KLF5 binds to and transcriptionally regulates genes encoding rate-limiting sphingolipid metabolism enzymes. Remarkably, skin barrier defects elicited by KLF5 ablation can be rescued by dietary interventions. Finally, we found that KLF5 is widely suppressed in human diseases with disrupted epidermal secretion, and its regulation of sphingolipid metabolism is conserved in human skin. Altogether, we established KLF5 as a disease-relevant transcription factor governing sphingolipid metabolism and barrier function in the skin, likely representing a long-sought secretory lineage-defining factor across tissue types.

Oncogenic Kras drives invasion and maintains metastases in colorectal cancer.

Human colorectal cancer (CRC) is a major cause of cancer mortality and frequently harbors activating mutations in the KRAS gene. To understand the role of oncogenic KRAS in CRC, we engineered a mouse model of metastatic CRC that harbors an inducible oncogenic Kras allele (Krasmut ) and conditional null alleles of Apc and Trp53 (iKAP). The iKAP model recapitulates tumor progression from adenoma through metastases. Whole-exome sequencing revealed that the Krasmut allele was heterogenous in primary tumors yet homogenous in metastases, a pattern consistent with activated Krasmut signaling being a driver of progression to metastasis. System-level and functional analyses revealed the TGF-ß pathway as a key mediator of Krasmut -driven invasiveness. Genetic extinction of Krasmut resulted in specific elimination of the Krasmut subpopulation in primary and metastatic tumors, leading to apoptotic elimination of advanced invasive and metastatic disease. This faithful CRC model provides genetic evidence that Krasmut drives CRC invasion and maintenance of metastases.

Experimental models of undifferentiated pleomorphic sarcoma and malignant peripheral nerve sheath tumor.

Undifferentiated pleomorphic sarcoma (UPS) and malignant peripheral nerve sheath tumor (MPNST) are aggressive soft tissue sarcomas that do not respond well to current treatment modalities. The limited availability of UPS and MPNST cell lines makes it challenging to identify potential therapeutic targets in a laboratory setting. Understanding the urgent need for improved treatments for these tumors and the limited cellular models available, we generated additional cell lines to study these rare cancers. Patient-derived tumors were used to establish 4 new UPS models, including one radiation-associated UPS-UPS271.1, UPS511, UPS0103, and RIS620, one unclassified spindle cell sarcoma-USC060.1, and 3 new models of MPNST-MPNST007, MPNST3813E, and MPNST4970. This study examined the utility of the new cell lines as sarcoma models by assessing their tumorigenic potential and mutation status for known sarcoma-related genes. All the cell lines formed colonies and migrated in vitro. The in vivo tumorigenic potential of the cell lines and corresponding xenografts was determined by subcutaneous injection or xenograft re-passaging into immunocompromised mice. USC060.1 and UPS511 cells formed tumors in mice upon subcutaneous injection. UPS0103 and RIS620 tumor implants formed tumors in vivo, as did MPNST007 and MPNST3813E tumor implants. Targeted sequencing analysis of a panel of genes frequently mutated in sarcomas identified TP53, RB1, and ATRX mutations in a subset of the cell lines. These new cellular models provide the scientific community with powerful tools for detailed studies of tumorigenesis and for investigating novel therapies for UPS and MPNST.

Metagenomes of rectal swabs in larger, advanced stage cervical cancers have enhanced mucus degrading functionalities and distinct taxonomic structure.

BACKGROUND: Gut microbiome community composition differs between cervical cancer (CC) patients and healthy controls, and increased gut diversity is associated with improved outcomes after treatment. We proposed that functions of specific microbial species adjoining the mucus layer may directly impact the biology of CC. METHOD: Metagenomes of rectal swabs in 41 CC patients were examined by whole-genome shotgun sequencing to link taxonomic structures, molecular functions, and metabolic pathway to patient's clinical characteristics. RESULTS: Significant association of molecular functions encoded by the metagenomes was found with initial tumor size and stage. Profiling of the molecular function abundances and their distributions identified 2 microbial communities co-existing in each metagenome but having distinct metabolism and taxonomic structures. Community A (Clostridia and Proteobacteria predominant) was characterized by high activity of pathways involved in stress response, mucus glycan degradation and utilization of degradation byproducts. This community was prevalent in patients with larger, advanced stage tumors. Conversely, community B (Bacteroidia predominant) was characterized by fast growth, active oxidative phosphorylation, and production of vitamins. This community was prevalent in patients with smaller, early-stage tumors. CONCLUSIONS: In this study, enrichment of mucus degrading microbial communities in rectal metagenomes of CC patients was associated with larger, more advanced stage tumors.

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.

Multi-omic molecular profiling reveals potentially targetable abnormalities shared across multiple histologies of brain metastasis.

The deadly complication of brain metastasis (BM) is largely confined to a relatively narrow cross-section of systemic malignancies, suggesting a fundamental role for biological mechanisms shared across commonly brain metastatic tumor types. To identify and characterize such mechanisms, we performed genomic, transcriptional, and proteomic profiling using whole-exome sequencing, mRNA-seq, and reverse-phase protein array analysis in a cohort of the lung, breast, and renal cell carcinomas consisting of BM and patient-matched primary or extracranial metastatic tissues. While no specific genomic alterations were associated with BM, correlations with impaired cellular immunity, upregulated oxidative phosphorylation (OXPHOS), and canonical oncogenic signaling pathways including phosphoinositide 3-kinase (PI3K) signaling, were apparent across multiple tumor histologies. Multiplexed immunofluorescence analysis confirmed significant T cell depletion in BM, indicative of a fundamentally altered immune microenvironment. Moreover, functional studies using in vitro and in vivo modeling demonstrated heightened oxidative metabolism in BM along with sensitivity to OXPHOS inhibition in murine BM models and brain metastatic derivatives relative to isogenic parentals. These findings demonstrate that pathophysiological rewiring of oncogenic signaling, cellular metabolism, and immune microenvironment broadly characterizes BM. Further clarification of this biology will likely reveal promising targets for therapeutic development against BM arising from a broad variety of systemic cancers.

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.

Clinical implications of cancer gene mutations in patients with chronic lymphocytic leukemia treated with lenalidomide.

Lenalidomide is clinically active in chronic lymphocytic leukemia (CLL), but its effectiveness in the context of the CLL mutational landscape is unknown. We performed targeted capture sequencing of 295 cancer genes in specimens from 102 CLL patients with treatment-naïve disease (TN patients) and 186 CLL patients with relapsed/refractory disease (R/R patients) who received lenalidomide-based therapy at our institution. The most frequently mutated gene was SF3B1 (15%), followed by NOTCH1 (14%) and TP53 (14%), with R/R patients having significantly more TP53 mutations than did TN patients. Among all lenalidomide-treated patients, del(17p) (P ≤ .001), del(11q) (P = .032), and complex karyotype (P = .022), along with mutations in TP53 (P ≤ .001), KRAS (P = .034), and DDX3X (P ≤ .001), were associated with worse overall response (OR). R/R patients with SF3B1 and MGA mutations had significantly worse OR (P = .025 and .035, respectively). TN and R/R patients with del(17p) and TP53 mutations had worse overall survival (OS) and progression-free survival (PFS). In R/R patients, complex karyotype and SF3B1 mutations were associated with worse OS and PFS; DDX3X mutations were associated with worse PFS only. Weibull regression multivariate analysis revealed that TP53 aberrations (del(17p), TP53 mutation, or both), along with complex karyotype and SF3B1 mutations, were associated with worse OS in the R/R cohort. Taken together, cancer gene mutations in CLL contribute to the already comprehensive risk stratification and add to prognosis and response to treatment. The related trials were registered at www.clinicaltrials.gov as #NCT00267059, #NCT00535873, #NCT00759603, #NCT01446133, and #NCT01002755.

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.

Characterization of HPV and host genome interactions in primary head and neck cancers.

Previous studies have established that a subset of head and neck tumors contains human papillomavirus (HPV) sequences and that HPV-driven head and neck cancers display distinct biological and clinical features. HPV is known to drive cancer by the actions of the E6 and E7 oncoproteins, but the molecular architecture of HPV infection and its interaction with the host genome in head and neck cancers have not been comprehensively described. We profiled a cohort of 279 head and neck cancers with next generation RNA and DNA sequencing and show that 35 (12.5%) tumors displayed evidence of high-risk HPV types 16, 33, or 35. Twenty-five cases had integration of the viral genome into one or more locations in the human genome with statistical enrichment for genic regions. Integrations had a marked impact on the human genome and were associated with alterations in DNA copy number, mRNA transcript abundance and splicing, and both inter- and intrachromosomal rearrangements. Many of these events involved genes with documented roles in cancer. Cancers with integrated vs. nonintegrated HPV displayed different patterns of DNA methylation and both human and viral gene expressions. Together, these data provide insight into the mechanisms by which HPV interacts with the human genome beyond expression of viral oncoproteins and suggest that specific integration events are an integral component of viral oncogenesis.

The complete genome of an individual by massively parallel DNA sequencing.

The association of genetic variation with disease and drug response, and improvements in nucleic acid technologies, have given great optimism for the impact of 'genomic medicine'. However, the formidable size of the diploid human genome, approximately 6 gigabases, has prevented the routine application of sequencing methods to deciphering complete individual human genomes. To realize the full potential of genomics for human health, this limitation must be overcome. Here we report the DNA sequence of a diploid genome of a single individual, James D. Watson, sequenced to 7.4-fold redundancy in two months using massively parallel sequencing in picolitre-size reaction vessels. This sequence was completed in two months at approximately one-hundredth of the cost of traditional capillary electrophoresis methods. Comparison of the sequence to the reference genome led to the identification of 3.3 million single nucleotide polymorphisms, of which 10,654 cause amino-acid substitution within the coding sequence. In addition, we accurately identified small-scale (2-40,000 base pair (bp)) insertion and deletion polymorphism as well as copy number variation resulting in the large-scale gain and loss of chromosomal segments ranging from 26,000 to 1.5 million base pairs. Overall, these results agree well with recent results of sequencing of a single individual by traditional methods. However, in addition to being faster and significantly less expensive, this sequencing technology avoids the arbitrary loss of genomic sequences inherent in random shotgun sequencing by bacterial cloning because it amplifies DNA in a cell-free system. As a result, we further demonstrate the acquisition of novel human sequence, including novel genes not previously identified by traditional genomic sequencing. This is the first genome sequenced by next-generation technologies. Therefore it is a pilot for the future challenges of 'personalized genome sequencing'.

Targeting MCL1-driven anti-apoptotic pathways overcomes blast progression after hypomethylating agent failure in chronic myelomonocytic leukemia.

RAS pathway mutations, which are present in 30% of patients with chronic myelomonocytic leukemia (CMML) at diagnosis, confer a high risk of resistance to and progression after hypomethylating agent (HMA) therapy, the current standard of care for the disease. Here, using single-cell, multi-omics technologies, we seek to dissect the biological mechanisms underlying the initiation and progression of RAS pathway-mutated CMML. We identify that RAS pathway mutations induce transcriptional reprogramming of hematopoietic stem and progenitor cells (HSPCs) and downstream monocytic populations in response to cell-intrinsic and -extrinsic inflammatory signaling that also impair the functions of immune cells. HSPCs expand at disease progression after therapy with HMA or the BCL2 inhibitor venetoclax and rely on the NF-κB pathway effector MCL1 to maintain survival. Our study has implications for the development of therapies to improve the survival of patients with RAS pathway-mutated CMML.

Intratumoral microbiome of adenoid cystic carcinomas and comparison with other head and neck cancers.

Adenoid cystic carcinoma (ACC) is a rare, usually slow-growing yet aggressive head and neck malignancy. Despite its clinical significance, our understanding of the cellular evolution and microenvironment in ACC remains limited. We investigated the intratumoral microbiomes of 50 ACC tumor tissues and 33 adjacent normal tissues using 16S rRNA gene sequencing. This allowed us to characterize the bacterial communities within the ACC and explore potential associations between the bacterial community structure, patient clinical characteristics, and tumor molecular features obtained through RNA sequencing. The bacterial composition in the ACC was significantly different from that in adjacent normal salivary tissue, and the ACC exhibited diverse levels of species richness. We identified two main microbial subtypes within the ACC: oral-like and gut-like. Oral-like microbiomes, characterized by increased diversity and abundance of Neisseria, Leptotrichia, Actinomyces, Streptococcus, Rothia, and Veillonella (commonly found in healthy oral cavities), were associated with a less aggressive ACC-II molecular subtype and improved patient outcomes. Notably, we identified the same oral genera in oral cancer and head and neck squamous cell carcinomas. In both cancers, they were part of shared oral communities associated with a more diverse microbiome, less aggressive tumor phenotype, and better survival that reveal the genera as potential pancancer biomarkers for favorable microbiomes in ACC and other head and neck cancers. Conversely, gut-like intratumoral microbiomes, which feature low diversity and colonization by gut mucus layer-degrading species, such as Bacteroides, Akkermansia, Blautia, Bifidobacterium, and Enterococcus, were associated with poorer outcomes. Elevated levels of Bacteroides thetaiotaomicron were independently associated with significantly worse survival and positively correlated with tumor cell biosynthesis of glycan-based cell membrane components.

Tsyn-Seq: a T-cell Synapse-Based Antigen Identification Platform.

Tools for genome-wide rapid identification of peptide-major histocompatibility complex targets of T-cell receptors (TCR) are not yet universally available. We present a new antigen screening method, the T-synapse (Tsyn) reporter system, which includes antigen-presenting cells (APC) with a Fas-inducible NF-κB reporter and T cells with a nuclear factor of activated T cells (NFAT) reporter. To functionally screen for target antigens from a cDNA library, productively interacting T cell-APC aggregates were detected by dual-reporter activity and enriched by flow sorting followed by antigen identification quantified by deep sequencing (Tsyn-seq). When applied to a previously characterized TCR specific for the E7 antigen derived from human papillomavirus type 16 (HPV16), Tsyn-seq successfully enriched the correct cognate antigen from a cDNA library derived from an HPV16-positive cervical cancer cell line. Tsyn-seq provides a method for rapidly identifying antigens recognized by TCRs of interest from a tumor cDNA library. See related Spotlight by Makani and Joglekar, p. 515.

A Modular Trial of Androgen Signaling Inhibitor Combinations Testing a Risk-Adapted Strategy in Patients with Metastatic Castration-Resistant Prostate Cancer.

PURPOSE: To determine the efficacy and safety of risk-adapted combinations of androgen signaling inhibitors and inform disease classifiers for metastatic castration-resistant prostate cancers. PATIENTS AND METHODS: In a modular, randomized phase II trial, 192 men were treated with 8 weeks of abiraterone acetate, prednisone, and apalutamide (AAPA; module 1) and then allocated to modules 2 or 3 based on satisfactory (≥50% PSA decline from baseline and <5 circulating tumor cell/7.5 mL) versus unsatisfactory status. Men in the former were randomly assigned to continue AAPA alone (module 2A) or with ipilimumab (module 2B). Men in the latter group had carboplatin + cabazitaxel added to AAPA (module 3). Optional baseline biopsies were subjected to correlative studies. RESULTS: Median overall survival (from allocation) was 46.4 [95% confidence interval (CI), 39.2-68.2], 41.4 (95% CI, 33.3-49.9), and 18.7 (95% CI, 14.3-26.3) months in modules 2A (n = 64), 2B (n = 64), and 3 (n = 59), respectively. Toxicities were within expectations. Of 192 eligible patients, 154 (80.2%) underwent pretreatment metastatic biopsies. The aggressive-variant prostate cancer molecular profile (defects in ≥2 of p53, RB1, and PTEN) was associated with unsatisfactory status. Exploratory analyses suggested that secreted phosphoprotein 1-positive and insulin-like growth factor-binding protein 2-positive macrophages, druggable myeloid cell markers, and germline pathogenic mutations were enriched in the unsatisfactory group. CONCLUSIONS: Adding ipilimumab to AAPA did not improve outcomes in men with androgen-responsive metastatic castration-resistant prostate cancer. Despite the addition of carboplatin + cabazitaxel, men in the unsatisfactory group had shortened survivals. Adaptive designs can enrich for biologically and clinically relevant disease subgroups to contribute to the development of marker-informed, risk-adapted therapy strategies in men with prostate cancer.

Adipose Triglyceride Lipase Is a Therapeutic Target in Advanced Prostate Cancer That Promotes Metabolic Plasticity.

Lipid metabolism plays a central role in prostate cancer. To date, the major focus has centered on de novo lipogenesis and lipid uptake in prostate cancer, but inhibitors of these processes have not benefited patients. A better understanding of how cancer cells access lipids once they are created or taken up and stored could uncover more effective strategies to perturb lipid metabolism and treat patients. Here, we identified that expression of adipose triglyceride lipase (ATGL), an enzyme that controls lipid droplet homeostasis and a previously suspected tumor suppressor, correlates with worse overall survival in men with advanced, castration-resistant prostate cancer (CRPC). Molecular, genetic, or pharmacologic inhibition of ATGL impaired human and murine prostate cancer growth in vivo and in cell culture or organoids under conditions mimicking the tumor microenvironment. Mass spectrometry imaging demonstrated that ATGL profoundly regulates lipid metabolism in vivo, remodeling membrane composition. ATGL inhibition induced metabolic plasticity, causing a glycolytic shift that could be exploited therapeutically by cotargeting both metabolic pathways. Patient-derived phosphoproteomics identified ATGL serine 404 as a target of CAMKK2-AMPK signaling in CRPC cells. Mutation of serine 404 did not alter the lipolytic activity of ATGL but did decrease CRPC growth, migration, and invasion, indicating that noncanonical ATGL activity also contributes to disease progression. Unbiased immunoprecipitation/mass spectrometry suggested that mutation of serine 404 not only disrupts existing ATGL protein interactions but also leads to new protein-protein interactions. Together, these data nominate ATGL as a therapeutic target for CRPC and provide insights for future drug development and combination therapies. SIGNIFICANCE: ATGL promotes prostate cancer metabolic plasticity and progression through both lipase-dependent and lipase-independent activity, informing strategies to target ATGL and lipid metabolism for cancer treatment.

Clonal evolution of hematopoietic stem cells after cancer chemotherapy.

Normal hematopoietic stem and progenitor cells (HSPCs) inherently accumulate somatic mutations and lose clonal diversity with age, processes implicated in the development of myeloid malignancies 1 . The impact of exogenous stressors, such as cancer chemotherapies, on the genomic integrity and clonal dynamics of normal HSPCs is not well defined. We conducted whole-genome sequencing on 1,032 single-cell-derived HSPC colonies from 10 patients with multiple myeloma (MM), who had undergone various chemotherapy regimens. Our findings reveal that melphalan treatment distinctly increases mutational burden with a unique mutation signature, whereas other MM chemotherapies do not significantly affect the normal mutation rate of HSPCs. Among these therapy-induced mutations were several oncogenic drivers such as TET2 and PPM1D . Phylogenetic analysis showed a clonal architecture in post-treatment HSPCs characterized by extensive convergent evolution of mutations in genes such as TP53 and PPM1D . Consequently, the clonal diversity and structure of post-treatment HSPCs mirror those observed in normal elderly individuals, suggesting an accelerated clonal aging due to chemotherapy. Furthermore, analysis of matched therapy-related myeloid neoplasm (t-MN) samples, which occurred 1-8 years later, enabled us to trace the clonal origin of t-MNs to a single HSPC clone among a group of clones with competing malignant potential, indicating the critical role of secondary mutations in dictating clonal dominance and malignant transformation. Our findings suggest that cancer chemotherapy promotes an oligoclonal architecture with multiple HSPC clones possessing competing leukemic potentials, setting the stage for the selective emergence of a singular clone that evolves into t-MNs after acquiring secondary mutations. These results underscore the importance of further systematic research to elucidate the long-term hematological consequences of cancer chemotherapy.

Multi-omics Analysis Reveals Immune Features Associated with Immunotherapy Benefit in Patients with Squamous Cell Lung Cancer from Phase III Lung-MAP S1400I Trial.

PURPOSE: Identifying molecular and immune features to guide immune checkpoint inhibitor (ICI)-based regimens remains an unmet clinical need. EXPERIMENTAL DESIGN: Tissue and longitudinal blood specimens from phase III trial S1400I in patients with metastatic squamous non-small cell carcinoma (SqNSCLC) treated with nivolumab monotherapy (nivo) or nivolumab plus ipilimumab (nivo+ipi) were subjected to multi-omics analyses including multiplex immunofluorescence (mIF), nCounter PanCancer Immune Profiling Panel, whole-exome sequencing, and Olink. RESULTS: Higher immune scores from immune gene expression profiling or immune cell infiltration by mIF were associated with response to ICIs and improved survival, except regulatory T cells, which were associated with worse overall survival (OS) for patients receiving nivo+ipi. Immune cell density and closer proximity of CD8+GZB+ T cells to malignant cells were associated with superior progression-free survival and OS. The cold immune landscape of NSCLC was associated with a higher level of chromosomal copy-number variation (CNV) burden. Patients with LRP1B-mutant tumors had a shorter survival than patients with LRP1B-wild-type tumors. Olink assays revealed soluble proteins such as LAMP3 increased in responders while IL6 and CXCL13 increased in nonresponders. Upregulation of serum CXCL13, MMP12, CSF-1, and IL8 were associated with worse survival before radiologic progression. CONCLUSIONS: The frequency, distribution, and clustering of immune cells relative to malignant ones can impact ICI efficacy in patients with SqNSCLC. High CNV burden may contribute to the cold immune microenvironment. Soluble inflammation/immune-related proteins in the blood have the potential to monitor therapeutic benefit from ICI treatment in patients with SqNSCLC.

A novel Bayesian model for assessing intratumor heterogeneity of tumor infiltrating leukocytes with multi-region gene expression sequencing.

Intratumor heterogeneity (ITH) of tumor-infiltrated leukocytes (TILs) is an important phenomenon of cancer biology with potentially profound clinical impacts. Multi-region gene expression sequencing data provide a promising opportunity that allows for explorations of TILs and their intratumor heterogeneity for each subject. Although several existing methods are available to infer the proportions of TILs, considerable methodological gaps exist for evaluating intratumor heterogeneity of TILs with multi-region gene expression data. Here, we develop ICeITH, immune cell estimation reveals intratumor heterogeneity, a Bayesian hierarchical model that borrows cell type profiles as prior knowledge to decompose mixed bulk data while accounting for the within-subject correlations among tumor samples. ICeITH quantifies intratumor heterogeneity by the variability of targeted cellular compositions. Through extensive simulation studies, we demonstrate that ICeITH is more accurate in measuring relative cellular abundance and evaluating intratumor heterogeneity compared with existing methods. We also assess the ability of ICeITH to stratify patients by their intratumor heterogeneity score and associate the estimations with the survival outcomes. Finally, we apply ICeITH to two multi-region gene expression datasets from lung cancer studies to classify patients into different risk groups according to the ITH estimations of targeted TILs that shape either pro- or anti-tumor processes. In conclusion, ICeITH is a useful tool to evaluate intratumor heterogeneity of TILs from multi-region gene expression data.

Endothelial-to-osteoblast transition in normal mouse bone development.

Metastatic prostate cancer (PCa) in bone induces bone-forming lesions. We have previously shown that PCa-induced bone originates from endothelial cells (ECs) that have undergone EC-to-osteoblast (OSB) transition. Here, we investigated whether EC-to-OSB transition also occurs during normal bone formation. We developed an EC and OSB dual-color reporter mouse (DRM) model that marks EC-OSB hybrid cells with red and green fluorescent proteins. We observed EC-to-OSB transition (RFP and GFP co-expression) in both endochondral and intramembranous bone formation during embryonic development and in adults. Co-expression was confirmed in cells isolated from DRM. Bone marrow- and lung-derived ECs underwent transition to OSBs and mineralization in osteogenic medium. RNA-sequencing revealed GATA family transcription factors were upregulated in EC-OSB hybrid cells and knockdown of GATA3 inhibited BMP4-induced mineralization. Our findings support that EC-to-OSB transition occurs during normal bone development and suggest a new paradigm regarding the endothelial origin of OSBs.