Clonal Profiles Research Articles

Clonal and resistance profiles of fluoroquinolone-resistant uropathogenic Escherichia coli in countries with different practices of antibiotic prescription.

Antibiotic prescription practices differ between countries, influencing regional antimicrobial resistance prevalence. However, comparisons of clonal diversity among resistant bacteria in countries with different prescribing practices are rare. The rise of fluoroquinolone-resistant Escherichia coli (FQREC), often multidrug-resistant, exacerbates global antibiotic resistance. Unlike in the USA, antibiotics are commonly dispensed in Iraq without prescriptions, leading to widespread overuse and misuse. This study aimed to assess the impact of varying antibiotic use practices on FQREC diversity. We compared FQREC prevalence, multidrug resistance, and clonality of FQREC among E. coli isolated from urine submitted between 2017 and 2018 to three US hospitals and two Iraqi hospitals. All FQREC isolates were analyzed for QRDR mutations and the presence of PMQR genes. A subset of FQREC strains from the ST131-H30R/Rx subgroups underwent whole-genome sequencing (WGS) and phylogenetic analysis. E. coli from Iraq showed significantly higher resistance to all tested antibiotics compared to those from the USA, with 76.2% being FQREC versus 31.2% in the USA (p < 0.01). Iraqi FQREC strains were more frequently multidrug resistant. The predominant subgroup in both countries was ST131-H30, with the notable absence of ST1193 among Iraqi FQREC. Iraqi-origin ST131-H30 strains exhibited higher minimum inhibitory concentrations (MICs) for ciprofloxacin and greater resistance to third-generation cephalosporins (3GC), trimethoprim/sulfamethoxazole (TMP/STX), and imipenem (IMI) than those from the USA. Increased 3GC resistance in Iraqi strains was linked to a higher proportion of bla CTX-M-15-carrying H30Rx subclade isolates. Additionally, Iraqi H30 strains exhibited higher MICs for fluoroquinolones due to more frequent carriage of PMQR determinants compared to US strains. Whole-genome sequencing was performed on 46 Iraqi and 63 US H30 isolates. Phylogenetic analysis revealed two clades-H30R and H30Rx-present in both countries, with isolates from both regions distributed throughout, without the emergence of distinct new major subclones. However, Iraqi isolates tended to cluster in separate subclades, indicating endemic circulation of the strain groups. In regions like Iraq, where antibiotics are overused and misused, resistance among uropathogenic E. coli to various antibiotics is significantly higher. Most Iraqi resistant strains belong to well-known international groups, and no new highly successful strains have emerged. The absence of ST1193 in Iraq may reflect regional, socioeconomic, demographic, or cultural factors that hinder the success of certain strain groups in the country.

Understanding clonal heterogeneity and mechanisms of methotrexate resistance in two osteosarcoma PDX models.

e23512 Background: Osteosarcoma is a primary malignant bone tumor that is known to have significant intratumoral heterogeneity with branched cancer evolution. The lung is the most frequent site of metastatic disease and relapse occurs in more than 30 % children and young adults. It has been revealed that the metastatic ability and resistance of cancers including osteosarcoma might be due to the presence of a subset of cells within the intra-tumorally heterogeneous population. We investigated the emergence of resistant clonal subpopulations following targeted methotrexate treatment in two patient derived xenograft (PDX) models in vitro and in vivo. Our aim was to not only track the subclonal population of cells that establish in the lungs but also to identify clonal drivers and the mutational signature of the genes and pathways involved in MTX resistance. Methods: We previously established 20 barcoded osteosarcoma PDX models and demonstrated that the application of methotrexate (MTX-IC95) treatment as a selective pressure in seven models produces clonal subpopulations that were more resistant to MTX. In the current work we used barcoded next generation and RNA sequencing techniques to map the clonal and transcriptomic architecture of MTX resistant single clones in M36 and OS252 PDX models. We next injected the barcoded PDX cells intratibially in twenty severe combined immunodeficiency disease (SCID) mice to track the clonal population/subpopulation of cells that metastasize to the lungs. Following primary tumor growth, limb amputation was carried out and mice were treated with methotrexate or vehicle for two weeks. Primary tumor, blood and lung samples were collected after death from metastasis or euthanasia. All samples collected were snap frozen in liquid nitrogen and stored at -800C. DNA and RNA was extracted and analyzed by PCR, NGS or WES to map clonality, mutational and evolutionary profiles. Results: MTX treated single clones showed diversity and barcode abundance when compared to vehicle; barcode abundance reduced in the resistant clonal subpopulations when compared to the vehicles. Ten unique subclones were shared in both vehicle and MTX resistant clones in M36 while OS252 had only one shared subclone respectively. More than 50% of the injected mice developed primary tumors and matched lung micrometastases. We are currently analyzing both targeted and RNA sequencing results to further understand clonality, the mutational signature of the genes (DHFR, MTR, FPGS, GGH) and pathways involved in MTX resistance. Conclusions: We investigated mechanism of methotrexate resistance in two barcoded PDX models in vitro and in vivo using cytotoxicity, next generation and RNA sequencing techniques. Targeted amplicon sequencing of key genes involved in antifolate pathway will yield new information on clonal evolution and mechanisms of resistance.

Abstract 882: Harnessing single cell multi-omics data to identify predictors of clinical outcome in CD19 CAR T cell therapy

Abstract Background: Chimeric antigen receptor (CAR) T cell therapy is effective in treating B-cell malignancies, but factors influencing the persistence of functional CAR+ T cells, such as product composition, patients’ lymphodepletion, and immune reconstitution, are not well understood. To shed light on this issue, we conducted single-cell multi-omics analysis of transcriptional, clonal, and phenotypic profiles from pre- to 1-month post-infusion of CAR+ and CAR− T cells from patients who received a donor-derived 4-1BB CAR product targeting CD19. Methods: Longitudinal blood samples and infusion products were obtained from 10 individuals with large B cell lymphoma or B-cell acute lymphoblastic leukemia (CARTELL cohort in Sydney Australia1). CD45+ cells and enriched CAR+ T cells were single cell sorted for simultaneous scRNA-seq, 41 barcoded proteins (N=41) and T cell receptors, n= 115,000 cells across 21 samples. Machine learning methods were used to identify subsets of immune cells associated with clinically relevant parameters. Results: Following infusion, CAR+ and CAR− T cells showed similar differentiation profiles with clonally expanded populations across heterogeneous phenotypes, demonstrating clonal lineages and phenotypic plasticity. These findings were validated by performing the same analysis on two publicly available scRNA-seq data sets from two cohorts 2, 3 of paediatric (N=18) and adult (N=32) individual with large B cell lymphoma who received CAR19 T cells. Analysis of mass-cytometry data from our cohort and from 28 patients with B cell lymphoma treated with commercial CAR19 T cells (axi-cel)4, identified an association of CAR− and also CAR+ T cells, both expressing an 2B4− CD11b+ CD8+ phenotype, with complete remission or progressive disease at 3 and 6 months post-infusion. Finally, we report additional analysis on the immunogenomics of circulating CAR T cells in two patients who developed CAR T cell induced lymphoma and compared it with ascites and CAR-Lymphoma cells. Discussion: These results suggest that CAR+ and CAR− CD8+ T cells share a differentiation trajectory terminating in an NK-like phenotype that is associated with clinical outcome, and that non-CAR-derived signals may influence patients’ immune recovery and outcomes.

Abstract 148: Multiomic sequencing reveals diversity in clonal landscape and genomic alterations in a lung metastatic PDX osteosarcoma model

Abstract Introduction: Osteosarcoma is a primary malignant bone tumor characterized by the production of spindle cells resulting in immature bone formation. The lung is the most frequent site of metastatic disease and relapse occurs in more than 30 % of patients. Recent results show that relapse is due to a subset of cells with different phenotypic and genetic signatures conferring advantage to drive progression and drug resistance within the intra-tumorally heterogeneous population. Thus, developing a defined in vivo model that can identify the rare genetic subpopulation and recapitulate clonal evolution is crucial. In this study, we used whole exome, barcode and RNA sequencing to characterize the landscape of genomic alterations and also track and identify the clonal subpopulations, clonal and genetic drivers of lungs metastasis. Methods: We injected barcoded PDX (OS17) cells intravenously in ten SCID mice to track the clonal subpopulation of cells that metastasize to the lungs. Metastasized lung nodules were collected following death or euthanasia. All samples collected were snap frozen in liquid nitrogen and stored at -80°C. DNA and RNA was extracted and analyzed by PCR, NGS, WES and RNA sequencing to map clonality, mutational and evolutionary profiles. Results: Four (40%) of the mice developed lung metastases with three mice having multiple metastatic lung nodules. Surgery to remove the metastatic nodules was performed 216 ± 55.5 days following tumor cell injection. Nine metastatic lung nodules from three mice were used for sequencing and analysis. A Shannon-Weaver and Jaccard similarity index show a diversity in clonal architecture between the lung nodules. A total of 15,394 somatic variants was identified including 12,547 single nucleotide variants (SNV), 4054 synonymous, 2662 missense, 28 nonsense mutations and 1639 somatic indels. The SNVs accounted for 81% of all variants followed by CNV 11 % and indels 8% respectively. Synonymous and missense mutations were the most frequently observed in all samples. Copy number loss was two-fold compared to copy number gain across all tumors. Furthermore, the deletions span longer genomic distances compared to amplifications. The mutational spectrum showed that the highest substitutions were C&amp;gt;T/G&amp;gt;A followed by T&amp;gt;C/A&amp;gt;G; the least substitutions were T&amp;gt;A/A&amp;gt;T. Also, we identified a large set of genes with genomic aberrations. These include amplifications in VEGF, RB1, RUNX, PARP4, ICAM3, EGFR, BRCA2, COLA6A1, COLA6A2, CCND3, CDKN2D etc. Conclusions: In the current model, we identified multiple tumorigenic seeding clones and the diverse set of genes mutated in osteosarcoma. Further analysis of the mutational signature resulted in the identification of potential genes that act as drivers in the lung metastatic process. Acknowledgments: Swim Across America, the Foster Foundation and the Barbara Epstein Foundation. Citation Format: Sylvester Jusu, Wendong Zhang, Zhongting Zhang, Xu Zhaohui, Sankaranarayanan Kannan, Yifei Wang, Zhou Xin, Yi Yanhua, Michael Roth, Jonathan Gill, Richard Gorlick. Multiomic sequencing reveals diversity in clonal landscape and genomic alterations in a lung metastatic PDX osteosarcoma model [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 148.

Abstract 3693: Molecular and functional characterization of circulating CK+/CD45+ cells in breast cancer

Abstract Background: Circulating tumor cells (CTCs) are considered a key element in metastatic progression and represent a validated prognostic biomarker in breast cancer (BC). CTCs in BC are normally defined as epithelial cells (positive for EPCAM or cytokeratin, CK) lacking the expression of the leukocyte marker CD45. Circulating cells expressing both EPCAM/CK and CD45 (dual positive, DPcells) have however been reported in BC but not investigated. It is believed that DPcells derive from the heterotypic cell fusion of tumor and immune cells, but this hypothesis has not been demonstrated yet. We recently reported, for the first time, the association of DPcells with worse survival in a cohort of 341 BC patients (pts) [Reduzzi et al. ASCO 2022] supporting their tumor origin. Here, we further investigated their malignancy and metastatic potential by molecular profiling and functional experiments. Methods: The study consisted of 2 components: A) Single CTCs and DPcells were enriched from the blood of BC patients with the CellSearch platform, and isolated with the DEPArray system. Single cells underwent whole genome amplification and lowpass whole genome sequencing for single cell copy number alteration (CNA) profiling; B) DPcells and CTCs were investigated in the blood of immunodeficient (NSG) and immunocompetent (FVB) mouse models using various cell lines (LM2, 4T1, MVT1). To test their metastatic potential, DPcells from the primary tumor of FVB donor mice were transplanted into recipient NSG mice, and the metastatic burden in the lungs and liver quantified, in comparison with control CD45- tumor cells. Results: We analyzed 73 DPcells and 18 CTCs collected from 16 BC pts: Aberrant genomes were observed in 28% and 93% of evaluable cells, respectively. When collected from the same pt, CTCs showed clonal CNA profiles, while DPcells were highly heterogeneous. In vivo, DPcells could only be detected in the blood of FVB mice but not in NSG mice. They were very rare, comprising of only about 3% of the total CTC count. When DPcells and control tumor cells were isolated from the primary tumor of 10 donor mice and transplanted into the tail vein of 10 recipient mice, we observed no significant differences in the metastatic colonization of the lungs and liver by the two population, indicating that DPcells do in fact possess metastatic ability, and to a similar extent as control cancer cells. Conclusions: DPcells can have an aberrant genome, supporting their malignancy. Contrary to CTCs, DPcells’ CNA profiles are not clonal and, in vivo, they can only be observed in immunocompetent models; both observations support the tumor/immune fusion hypothesis. Finally, DPcells show metastatic potential in mice, in line with the association with worse survival we previously observed in BC pts. Overall DPcells seem a new CTC subpopulation that should not be neglected anymore. More studies are needed to better understand their origin and phenotypic/molecular features. Citation Format: Carolina Reduzzi, Eleonora Nicolo, Lauren L. Ozimski, Marco Silvestri, Lorenzo Gerratana, Mara S. Serafini, Elisabetta Molteni, Nadia Bayou, Amanda K. Strickland, Huiping Liu, Vera Cappelletti, Nicola Aceto, Massimo Cristofanilli. Molecular and functional characterization of circulating CK+/CD45+ cells in breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 3693.

Abstract A010: Longitudinal histology and spatial transcriptomics profiling in targeted treatment of melanoma patient-derived models interrogates persister and resistant cell populations

Abstract Over half of BRAF-mutant melanoma patients with initial response to targeted therapy will recur with resistant disease. It is thought that recurrence arises from the ability of resistant clonal populations to enter and exit a slow-cycling persister state, evading treatment in quiescence before returning to a proliferative state. However, clonal transcriptional profiles and their spatial relationships in the tumor environment are not well characterized. Through longitudinal profiling of tumors from pre-treatment through the growth of resistant tumors, we track clonal lineages and transcriptional state changes to identify recurrent evolutionary patterns. BRAF mutant patient-derived melanoma xenografts were treated with BRAF/MEK inhibitors through maximum tumor size reduction and re-growth. Spatial transcriptomics (ST) and H&amp;E histological staining were performed at five timepoints across 94 and 133 days of PDX models WM4237 and WM4007 treatment, tracking clonal populations with intact tissue structure. Deep learning feature extraction on high-resolution H&amp;E-stained images allowed the detection of histological features that co-localize with expression levels. A novel computational workflow performed integration of samples across time points, pathway enrichment analysis, and pseudotemporal ordering, allowing spatial recreation of copy number variation based clonal phylogenies and cell fate trajectory inference during the treatment course. We characterize the clonal populations by the burden of accumulated copy number variations and establish clonal phylogeny. With the added temporal resolution, ST profiling during melanoma treatment enables observation of a global shift where all lineages begin the transition into the persister transcriptional state. The drug-sensitive cells are eliminated due to treatment, and distinct transcriptional states arise due to metabolic plasticity regardless of the lineage of origin. The clones of the transient persister state constitute the minimal residual disease, where melanoma cells have an increase in invasive capacity and a slowing of the cell cycle. The ability to survive through the persister state and re-emerge into a proliferating clone where the cell cycle is upregulated varies by lineage. We observe that the transcriptional changes during emergence from dormancy show decreased invasiveness and continued reliance on oxidative respiration versus glycolysis. We compare the re-emerged resistant cell populations and pre-treatment samples to pinpoint the changes in the transcriptional state of the MAPK signaling pathway and find that sets of upregulated genes differ by the model, WM4237 and WM4007. We connect deep learning imaging features extracted from longitudinal histology profiling to ST profiles to associate histology data with the cell cycle activation, invasiveness, hypoxia state, and metabolic mode across the timepoints. We observe that the imaging features alone can be potentially used to track the evolution of the tumor cell population phenotypes during drug treatment. Citation Format: Sergii Domanskyi, Todd B. Sheridan, Brian J. Sanderson, SungHee Park, Jessica Kaster, Haiyin Li, Olga Anczukow, Meenhard Herlyn, Jeffrey H. Chuang, Jill C. Rubinstein. Longitudinal histology and spatial transcriptomics profiling in targeted treatment of melanoma patient-derived models interrogates persister and resistant cell populations [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Translating Cancer Evolution and Data Science: The Next Frontier; 2023 Dec 3-6; Boston, Massachusetts. Philadelphia (PA): AACR; Cancer Res 2024;84(3 Suppl_2):Abstract nr A010.

A case series exploring the human milk polyclonal IgA1 response to repeated SARS-CoV-2 vaccinations by LC-MS based fab profiling.

Upon vaccination against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) humans will start to produce antibodies targeting virus specific antigens that will end up in circulation. In lactating women such antibodies will also end up in breastmilk, primarily in the form of secretory immunoglobulin A1 (SIgA1), the most abundant immunoglobulin (Ig) in human milk. Here we set out to investigate the SIgA1 clonal repertoire response to repeated SARS-CoV-2 vaccination, using a LC-MS fragment antigen-binding (Fab) clonal profiling approach. We analyzed the breastmilk of six donors from a larger cohort of 109 lactating mothers who received one of three commonly used SARS-CoV-2 vaccines. We quantitatively monitored the SIgA1 Fab clonal profile over 16 timepoints, from just prior to the first vaccination until 15 days after the second vaccination. In all donors, we detected a population of 89-191 vaccine induced clones. These populations were unique to each donor and heterogeneous with respect to individual clonal concentrations, total clonal titer, and population size. The vaccine induced clones were dominated by persistent clones (68%) which came up after the first vaccination and were retained or reoccurred after the second vaccination. However, we also observe transient SIgA1 clones (16%) which dissipated before the second vaccination, and vaccine induced clones which uniquely emerged only after the second vaccination (16%). These distinct populations were observed in all analyzed donors, regardless of the administered vaccine. Our findings suggest that while individual donors have highly unique human milk SIgA1 clonal profiles and a highly personalized SIgA1 response to SARS-CoV-2 vaccination, there are also commonalities in vaccine induced responses.

Into the Dark Serum Proteome: Personalized Features of IgG1 and IgA1 Repertoires in Severe COVID-19 Patients

Serum proteomics has matured and is now able to monitor hundreds of proteins quantitatively in large cohorts of patients. However, the fine characteristics of some of the most dominant proteins in serum, the immunoglobulins, are in these studies often ignored, due to their vast, and highly personalized, diversity in sequences. Here, we focus exclusively on these personalized features in the serum proteome and distinctively chose to study individual samples from a low diversity population: elderly donors infected by severe acute respiratory syndrome corona virus 2 (SARS-CoV-2). By using mass spectrometry-based methods, immunoglobulin IgG1 and IgA1 clonal repertoires were monitored quantitatively and longitudinally in more than 50 individual serum samples obtained from 17 Corona virus disease 2019 patients admitted to intensive care units. These clonal profiles were used to examine how each patient reacted to a severe SARS-CoV-2 infection. All 17 donors revealed unique polyclonal repertoires and substantial changes over time, with several new clones appearing following the infection, in a few cases leading to a few, very high, abundant clones dominating their repertoire. Several of these clones were de novo sequenced through combinations of top-down, middle-down, and bottom-up proteomics approaches. This revealed sequence features in line with sequences deposited in the SARS-CoV-specific antibody database. In other patients, the serological Ig profiles revealed the treatment with tocilizumab, that subsequently dominated their serological IgG1 repertoire. Tocilizumab clearance could be monitored, and a half-life of approximately 6 days was established. Overall, our longitudinal monitoring of IgG1 and IgA1 repertoires of individual donors reveals that antibody responses are highly personalized traits of each patient, affected by the disease and the chosen clinical treatment. The impact of these observations argues for a more personalized and longitudinal approach in patients’ diagnostics, both in serum proteomics as well as in monitoring immune responses.

CAR+ and CAR− T cells share a differentiation trajectory into an NK-like subset after CD19 CAR T cell infusion in patients with B cell malignancies

Chimeric antigen receptor (CAR) T cell therapy is effective in treating B cell malignancies, but factors influencing the persistence of functional CAR+ T cells, such as product composition, patients’ lymphodepletion, and immune reconstitution, are not well understood. To shed light on this issue, here we conduct a single-cell multi-omics analysis of transcriptional, clonal, and phenotypic profiles from pre- to 1-month post-infusion of CAR+ and CAR− T cells from patients from a CARTELL study (ACTRN12617001579381) who received a donor-derived 4-1BB CAR product targeting CD19. Following infusion, CAR+ T cells and CAR− T cells shows similar differentiation profiles with clonally expanded populations across heterogeneous phenotypes, demonstrating clonal lineages and phenotypic plasticity. We validate these findings in 31 patients with large B cell lymphoma treated with CD19 CAR T therapy. For these patients, we identify using longitudinal mass-cytometry data an association between NK-like subsets and clinical outcomes at 6 months with both CAR+ and CAR− T cells. These results suggest that non-CAR-derived signals can provide information about patients’ immune recovery and be used as correlate of clinically relevant parameters.

Determination of T-cell clonality and expression profiles of Toll-like receptors signaling pathway genes and related miRNAs in patients with mycosis fungoides

Cutaneous T-cell lymphomas (CTCL) encompass a group of diseases characterized by the presence of malignant clonal CD4+ T lymphocytes in the skin. Mycosis fungoides (MF) is the most prevalent form of CTCL, accounting for approximately 60 % of cutaneous T-cell lymphomas and 50 % of all primary cutaneous lymphomas. Despite ongoing research, the precise pathogenesis of MF remains incompletely understood. Toll-like receptors (TLRs) have the ability to specifically recognize ligands, subsequently induce the expression of diverse genes and activate innate immunity within the cell. Furthermore, miRNAs play a crucial role in regulating various aspects of immune cell function. The aim of our study was to explore the potential roles of TLRs and the genes implicated in their signal transduction, along with the expression status of miRNAs in the mechanisms underlying MF. Additionally, we assessed the clonal status and compared it with clinicopathological data using a T-cell clonality assay. To determine the expression status of TLR pathway genes and miRNAs, we conducted RT-PCR analysis on 52 MF samples and 50 control paraffin block materials. Pathway analysis were conducted using the KEGG database. T-cell receptor (TCR) gamma clonality changes were evaluated. Results from the study revealed increased expressions of TLR-1, -4, -8, IRF7, TRAF3, MEK1, MEK2, Elk1, NFkB, hsa-miR-21-5p, and hsa-miR-155-5p, as well as decreased expressions of hsa-miR-130a-3p, hsa-miR-210-3p, and hsa-let-7e-5p in the MF group. TCR gamma clonal change analysis demonstrated that 55.5 % of the analysed DNAs exhibited monoclonal and biallelic patterns, while 45.5 % displayed polyclonality. These findings collectively suggest the potential influence and therapeutic possibilities of the TLR signalling pathway in the molecular pathogenesis of MF.

Clonal dynamics and copy number variants by single-cell analysis in leukemic evolution of myeloproliferative neoplasms.

Transformation from chronic (CP) to blast phase (BP) in myeloproliferative neoplasm (MPN) remains poorly characterized, and no specific mutation pattern has been highlighted. BP-MPN represents an unmet need, due to its refractoriness to treatment and dismal outcome. Taking advantage of the granularity provided by single-cell sequencing (SCS), we analyzed paired samples of CP and BP in 10 patients to map clonal trajectories and interrogate target copy number variants (CNVs). Already at diagnosis, MPN present as oligoclonal diseases with varying ratio of mutated and wild-type cells, including cases where normal hematopoiesis was entirely surmised by mutated clones. BP originated from increasing clonal complexity, either on top or independent of a driver mutation, through acquisition of novel mutations as well as accumulation of clones harboring multiple mutations, that were detected at CP by SCS but were missed by bulk sequencing. There were progressive copy-number imbalances from CP to BP, that configured distinct clonal profiles and identified recurrences in genes including NF1, TET2, and BCOR, suggesting an additional level of complexity and contribution to leukemic transformation. EZH2 emerged as the gene most frequently affected by single nucleotide and CNVs, that might result in EZH2/PRC2-mediated transcriptional deregulation, as supported by combined scATAC-seq and snRNA-seq analysis of the leukemic clone in a representative case. Overall, findings provided insights into the pathogenesis of MPN-BP, identified CNVs as a hitherto poorly characterized mechanism and point to EZH2 dysregulation as target. Serial assessment of clonal dynamics might potentially allow early detection of impending disease transformation, with therapeutic implications.

CAR +and CAR −T cells differentiate into an NK-like subset that is associated with increased inflammatory cytokines following infusion

Abstract Chimeric antigen receptor (CAR) T cells have demonstrable efficacy in treating B-cell malignancies. Factors such as product composition, lymphodepletion and immune reconstitution are known to influence CAR +T cell survival and persistence, however the determinants of CAR +T cells’ fate and differentiation and how the patients’ reconstituting immune system influence therapy outcomes remain poorly understood. We applied single cell multi-omics analysis of RNA, protein and clonal profiles of CAR +T and PBMC in patients receiving donor-derived products from a piggybac CAR19 clinical trial (N=8 patients). We used these data to reconstruct a differentiation trajectory, which explained the observed phenotype and fate of CAR +and CAR −T cells within the first 30 days of post-infusion. Following initial lympho-depletion, endogenous CAR −CD8 +, γδ T cells and CAR+ T cells clonally expanded, and differentiated across heterogenous phenotypes. They transitioned from a dominant resting or proliferating state into precursor of exhausted T cells, and finally into a terminal NK-like phenotype. The subset of terminal NK-like CAR +T cells was associated with increased serum concentrations of proinflammatory cytokines including IL-12 and IL-18. Validation on a published cohort (61 patients) showed that this subset was associated with clinical outcome at 6 months. Finally, we analysed in detail one patient that developed CAR T cell induced lymphoma and discussed the impact of early phase therapy in determining this unexpected outcome. These results demonstrate that CAR- and CAR+ T cells share a differentiation profile which may be determined by non-CAR derived signals and which may influence patients’ immune-recovery and outcome. National Health And Medical Research Funding (Australia) (APP1121643, APP1128416) BD Bioscience

Antimicrobial resistance profiles and clonal diversity of Staphylococcus epidermidis isolates from pig farms, slaughterhouses, and retail pork

Livestock-associated coagulase-negative staphylococci (CoNS), such as Staphylococcus (S.) epidermidis, have emerged as a significant reservoir of antimicrobial resistance (AMR). In the current study, the AMR profiles and genetic diversity of S. epidermidis isolates obtained from pig farms, slaughterhouses, and retail pork were analyzed. A total of 89 S. epidermidis isolates, comprising 22 methicillin-resistant (MRSE) and 67 methicillin-susceptible S. epidermidis (MSSE) isolates, were assessed to determine (i) the clonal lineages of the isolates [multilocus sequence (MLST), agr, and staphylococcal cassette chromosome mec (SCCmec) types], (ii) the profiles of AMR phenotypes, and (iii) the carriage of genetic factors associated with major AMR phenotypes and zinc chloride resistance. Two dominant clonal lineages of S. epidermidis, ST100 and ST570, were observed on pig farms, especially in healthy pigs. In addition, potential transmission of pig-associated ST100 MRSE-SCCmec V and ST100 MSSE to farm workers was identified. The high prevalence of ST100 and ST570 isolates in pig farms was associated with high levels of AMR and zinc chloride resistance. In relation to resistance phenotypes, higher carriage rates of resistance genes, such as β–lactams (mecA), phenicols (fexA), and tetracyclines [tet(K), tet(L), tet(S), tet(M), and tet(O)], were identified in pig farm-associated isolates. Furthermore, cfr-mediated linezolid resistance was detected in 14 MSSE isolates from pig farms. This study is the first to provide important insights into the clonal structures and AMR profiles of S. epidermidis isolates collected from healthy pigs, carcass/pork samples, and human workers in Korea.

Abstract 2257: Activation and clonotypic expansion of the native T cell repertoire identifies durable response to CD19 CAR T cell therapy

Abstract Background: Axicabtagene ciloleucel (Axi-cel), a CD19 directed CAR T cell therapy, results in durable response in a subset of patients with relapsed/refractory large B cell lymphoma (LBCL) in the absence of persistent circulating CAR T cells. Aim: We postulated that long-term efficacy of CAR T therapy depends on the downstream triggering of native T cell immunity. We performed single-cell transcriptomics of longitudinal peripheral blood (PB) samples and RNAseq of tumor samples from patients of the ZUMA 1 Axi-cel study. Methods: Single cell immunoprofiling (5’ expression + V(D)J, 10x Genomics) was performed on PB mononuclear cell samples from ZUMA-1 patients (N=32), collected at leukapheresis, 4 weeks, 6 and 12 months post Axi-cel infusion. RNAseq was performed on FFPE lymphoma samples (N=17). Patients were divided into 3 groups: non-responders, relapsed within 1 year from CAR T infusion, and long-term responders. A total of 405,775 cells passed quality check, capturing 73 cellular populations. Results: Long-term responders presented a distinct T cell landscape with increased CD8 T cells and CD8/CD4 ratios prior to CAR T therapy, compared to the other groups, with similar trends observed across time points. They also presented an increased abundance of 3 distinct CD8 T cell populations: (a) cells expressing cytotoxic and NK cell markers, (b) CD8 T effector memory cells characterized by CXCR4, TGFB1, and BCL3, and (c) proinflammatory CD8 T cells. In contrast, patients with early relapse showed increased levels of regulatory T cells pre-/post-CAR T infusion and lower abundance of CD4 cytotoxic T cells.Comparisons of the TCR repertoire pre-/post-CAR T demonstrated a greater clonal expansion of cytotoxic CD4 and CD8 T cell populations in the long-term responders, with high similarity of expanded clones post CAR T. Shared PB T cell clones and tumor antigen sequences derived from the RNAseq samples suggest an antigen-specific response driven by common epitopes. In contrast to the lack of expansion of T reg clones in responders, relapsed patients demonstrated high T reg-clonal expansion 6 months post treatment.Monocyte and NK cells were less prevalent in responders before treatment, driven by differences in phagocytic monocyte and inhibitory NK cell abundance. Modeling the interaction between monocyte and effector cell populations based on ligand receptor expression was suggestive of negative immunoregulatory impact on the T cell populations. Conclusion: The application of single-cell immunoprofiling on longitudinal samples from ZUMA-1 patients demonstrated distinct cellular and clonal profiles among long-term responders. These findings confirm our hypothesis of an important role for the native immune cell repertoire in response to CAR T therapeutics and can be utilized to further our understanding but also to potentially inform patient stratification, management, and treatment. Citation Format: Dimitra Karagkouni, Giulia Cheloni, Yered Pita-Juarez, Daniela Torres, Eleni Kanata, Zachary Avigan, Jessica Liegel, Dina Stroopinsky, Brodie Miles, Gayatri Tiwari, Jenny Kim, Mike Mattie, Jacalyn Rosenblatt, David Avigan, Ioannis Vlachos. Activation and clonotypic expansion of the native T cell repertoire identifies durable response to CD19 CAR T cell therapy [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 2257.

Abstract 5932: Clonal mapping phylogenetic and transcriptomic analysis of an M36 PDX amputation lung metastatic osteosarcoma model

Abstract Introduction: Osteosarcoma is a primary malignant bone tumor characterized by the production of spindle cells resulting in immature bone formation. It has been revealed that metastatic cancers including osteosarcoma have a subset of cells with different phenotypic and genetic signatures conferring advantage to drive progression and drug resistance within the intra-tumorally heterogeneous population. However, the lack of a defined in vivo model that can identify the rare genetic subpopulation and recapitulate clonal evolution has been a challenge. In this study, we used an established model to track and identify the clonal subpopulations and assess the genetic signatures that confer specific advantages selecting for uncontrolled growth, adaptation and colonization in the lungs. Methods: We injected barcoded M36 PDX osteosarcoma cells intratibially in ten severe combined immunodeficiency disease (SCID) mice to track the clonal population/subpopulation of cells that metastasize the established primary tumors to the lungs. Limb amputation was performed on the tumor bearing leg upon detection of palpable tumor. Lung tissue was harvested post-amputation following death from metastasis or euthanasia. DNA and RNA from frozen tissues was extracted and analyzed by PCR, NGS, WES and RNA sequencing to map clonality, mutational and evolutionary profiles. Results: Limb amputation was performed in 7 mice with all mice surviving the surgery. The mean mouse survival following amputation was 46.9 days and the longest survival was 167 days. Four (40%) of the primary tumor bearing mice had lung metastasis. A Shannon-Wiener and Jaccard similarity index show a diversity in clonal architecture between the primary tumors of all mice injected on the same day and a diversity in clonal population in both the primary tumors and their matched lung metastatic nodules. Clonal mapping further revealed 21 unique clones shared in all primary tumors examined. However, the clonal frequencies and distribution varied within each primary tumor and across the tumor samples. Interestingly, clonal frequencies increased in the lung metastatic nodules when compared to their matched primary tumors. Four distinct clones that were present in all the primary tumors were also found in all lung metastatic samples. These unique clones, here in called “seeding clones” may play an important role in the lung metastatic process. Conclusions: We previously showed that the barcoded osteosarcoma PDX models are inter and intratumorally heterogenous both in vitro and in vivo. In the current model, we identify the multiple tumorigenic seeding clones that potentially act as drivers in the lung metastatic process. Furthermore, we show the dynamic evolutionary relationships between the clones by reconstructing a phylogenetic map of the M36 PDX tumor model. Funding: Swim Across America, Foster Foundation and the Barbara Epstein Foundation. Citation Format: Sylvester Jusu, Zhongting Zhang, Zhang Wendong, Zhaohui Xu, Michael Roth, Jonathan Gill, Richard Gorlick. Clonal mapping phylogenetic and transcriptomic analysis of an M36 PDX amputation lung metastatic osteosarcoma model. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 5932.

Abstract 182: Bladder cancer organoids as a functional system to model different disease stages and therapy response

Abstract Introduction: Interpatient heterogeneity is one of the causes of treatment failure in bladder cancer (BLCa) patients, who would therefore profit from tailor-made therapies. Toward the direction of precision medicine, patient-derived organoids (PDOs), which retain parental tumor (PT) features, have been successfully tested in different cancer types and shown to predict patients’ responses. Therefore, our study aimed to explore the ability of BLCa PDOs to retain PT features and to determine patient drug sensitivity. We correlated PDO drug response profiles to their genomic profiles and the patient's clinical history. Methods: PDOs were derived from the suspension culture of single cells. Whole exome sequencing was performed on genomic DNA and scRNA sequencing on 3 PT/PDO pairs. Histological evaluation of PDO and PT marker expression and morphology was performed. PDOs were treated for 48h and drug response was evaluated by viability assay. Results: We derived PDOs from non-muscle-invasive and muscle-invasive BLCa. Genomic analysis of 15 representative PT/PDOs pairs showed high similarity in terms of copy number variations and single nucleotide variants (SNVs) with preservation of peculiar BLCa alterations. PT genomic heterogeneity was preserved in PDOs as demonstrated by the correlation between SNVs clonality profiles of PT/PDOs pairs. Moreover, marker analyses at the scRNA and protein levels supported that PDOs retained the main tumor phenotype. Interestingly, PDO grew accordingly to three morphologies linked to marker expression and PT stage.PDOs were used in drug screen assays for testing standard-of-care (SOC) and FDA-approved drugs. PDO responses were highly heterogeneous highlighting interpatient heterogeneity. In 2 case, PDO sensitivity mimicked the patient’s response to SOC and in a few cases correlated with the sample-specific genomic background. Due to the high relapse rate of BLCa, we performed 2 longitudinal studies. In the first one, PDOs from a relapse collected after epirubicin treatment showed less sensitivity to the drug compared to the baseline sample. Moreover, the SNVs' mean clonality of relapse PDOs was significantly higher than the baseline PDOs which were characterized by mutations in DNA damage repair genes (e.g., ATM, FANC1, ERCC2, and BRCA2). We hypothesize a drug-induced selection of a pre-existing epirubicin-resistant population. In the second study, PDOs from the cystectomy and the following relapse showed similar genomic and drug sensitivity profiles suggesting that the first sample could be informative for adjuvant therapy. Conclusion: We established PDOs from different BLCa stages and grades which preserve PT features and can be implemented in drug screening assays for personalized approaches. The power of PDOs in precision medicine was highlighted by the longitudinal studies that showed that PDOs mirrored patient response and tumor evolution in vitro. Citation Format: Martina Minoli, Thomas Cantore, Daniel Hanhart, Mirjam Kiener, Tarcisio Fedrizzi, Federico La Manna, Sofia Karkampouna, Panagiotis Chouvardas, Vera Genitsch, José Antonio Rodriguez-Calero, Eva Compérat, Irena Klima, Paola Gasperini, Bernhard Kiss, Roland Seiler-Blarer, Francesca Demichelis, George N. Thalmann, Marianna Kruithof-de Julio. Bladder cancer organoids as a functional system to model different disease stages and therapy response [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 182.

Reconstructing clonal tree for phylo-phenotypic characterization of cancer using single-cell transcriptomics

Functional characterization of the cancer clones can shed light on the evolutionary mechanisms driving cancer’s proliferation and relapse mechanisms. Single-cell RNA sequencing data provide grounds for understanding the functional state of cancer as a whole; however, much research remains to identify and reconstruct clonal relationships toward characterizing the changes in functions of individual clones. We present PhylEx that integrates bulk genomics data with co-occurrences of mutations from single-cell RNA sequencing data to reconstruct high-fidelity clonal trees. We evaluate PhylEx on synthetic and well-characterized high-grade serous ovarian cancer cell line datasets. PhylEx outperforms the state-of-the-art methods both when comparing capacity for clonal tree reconstruction and for identifying clones. We analyze high-grade serous ovarian cancer and breast cancer data to show that PhylEx exploits clonal expression profiles beyond what is possible with expression-based clustering methods and clear the way for accurate inference of clonal trees and robust phylo-phenotypic analysis of cancer.

Integrated single-cell multiomics reveals novel immune candidate markers for post-traumatic coagulopathy.

Post-traumatic coagulopathy (PTC) is a critical pathology in traumatic brain injury (TBI), however, its potential mechanism is not clear. To explore this in peripheral samples, we integrated single cell RNA-sequencing and T cell repertoire (TCR)-sequencing across a cohort of patients with TBI. Clinical samples from patients with more brain severity demonstrated overexpression of T cell receptor-encoding genes and less TCR diversity. By mapping TCR clonality, we found patients with PTC have less TCR clones, and the TCR clones are mainly distributed in cytotoxic effector CD8+T cell. In addition, the counts of CD8+ T cell and natural killer (NK) cells are associated with the coagulation parameter by WGCNA, and the granzyme and lectin-like receptor profiles are also decreased in the peripheral blood from TBI patients, suggesting that reduced peripheral CD8+ clonality and cytotoxic profiles may be involved in PTC after TBI. Our work systematically revealed the critical immune status in PTC patients at the single-cell level.

Advances and challenges in zebrafish models of T-cell acute lymphoblastic leukemia

T-cell acute lymphoblastic leukemia (T-ALL) presents a significant challenge in oncology, necessitating innovative research approaches for improved comprehension and therapeutic advancements. In recent years, zebrafish (Danio rerio) have emerged as a vital model organism in T-ALL research, conferring distinctive advantages over conventional models. Zebrafish, sharing substantial genetic homology with humans, accurately replicate key phenotypic and molecular features of human T-ALL, including clonal expansion, aneuploidy, organ infiltration, and gene expression profiles. They also develop chromosomal translocations similar to those in human T-ALL. This comprehensive review explores the diverse zebrafish T-ALL models, encompassing genetic alterations and inducible systems that facilitate the investigation of disease initiation, progression, and response to therapy. These models have revealed critical insights into the roles of specific genes and the identification of novel oncogenic drivers. Moreover, zebrafish models enable high-throughput drug screening, providing an expedited path toward discovering potential therapeutics. However, while zebrafish models present exceptional advantages, they are not without challenges, including genetic differences, biological disparities, and technical complexities. To exploit the potential of zebrafish in T-ALL research, ongoing efforts are directed toward CRISPR/Cas9 gene editing, improved antibody development, and enhanced experimental techniques.

Expansions of tumor-reactive Vdelta1 gamma-delta T cells in newly diagnosed patients with chronic myeloid leukemia

Recent studies have underscored the importance of gamma-delta (γδ) T cells in mediating potent MHC-unrestricted cytotoxicity in numerous malignancies. Here, we analyzed Vδ1 and Vδ2 γδ T cell subsets in newly diagnosed chronic myeloid leukemia (CML) patients (n = 40) who had initiated tyrosine kinase inhibitor (TKI) therapy including imatinib (n = 22), nilotinib (n = 14) and dasatinib (n = 4). Patient peripheral blood samples were analyzed at diagnosis and monitored prospectively at 3, 6, 12 and 18 months post-TKI. γδ T cells isolated from healthy donors and CML patients were used against K562, LAMA-84 and KYO-1 cell lines and against primary CML cells in cytotoxicity assays. We found large expansions of Vδ1 and Vδ2 T cells in patients at diagnosis compared to age-matched healthy donors (n = 40) (p < 0.0001). The γδ T cell reconstitution in patients on imatinib and also on nilotinib showed significant reductions of Vδ1 T cell and Vδ2 T cell absolute counts at 3 months compared to diagnosis. Importantly, Vδ1 and Vδ2 T absolute cell counts remained at normal levels from 3 months throughout the follow-up. Next, we observed susceptibility to specific lysis of primary CML tumor cells by Vδ1 T cells from healthy donors. Furthermore, we determined inherent cytotoxic reactivity by autologous patients’ Vδ1 T lymphocytes against primary CML tumor cells. Finally, the TCR clonality profiles showed in CML patients mostly polyclonal repertoires regardless of the TKI. Our results provide further evidence into γδ T cell antileukemia immunity in CML that might be beneficial for long-term disease control and treatment outcome.