Mapping single-cell developmental potential in health and disease with interpretable deep learning.

Single-cell RNA sequencing (scRNA-seq) has transformed our understanding of cell fate in developmental systems. However, identifying the molecular hallmarks of potency - the capacity of a cell to differentiate into other cell types - has remained challenging. Here, we introduce CytoTRACE 2, an interpretable deep learning framework for characterizing potency and differentiation states on an absolute scale from scRNA-seq data. Across 31 human and mouse scRNA-seq datasets encompassing 28 tissue types, CytoTRACE 2 outperformed existing methods for recovering experimentally determined potency levels and differentiation states covering the entire range of cellular ontogeny. Moreover, it reconstructed the temporal hierarchy of mouse embryogenesis across 62 timepoints; identified pan-tissue expression programs that discriminate major potency levels; and facilitated discovery of cellular phenotypes in cancer linked to survival and immunotherapy resistance. Our results illuminate a fundamental feature of cell biology and provide a broadly applicable platform for delineating single-cell differentiation landscapes in health and disease.

Racial disparities in early-stage NSCLC treatment: A Call for Action.

Non-small cell lung cancer (NSCLC) is associated with a 5-year survival rate of only 28%, however, when caught at an early stage, it can be cured with surgery or stereotactic body radiotherapy (SBRT). Unfortunately, racial disparities may result in limited access to care for some patients. Liu et al. analyzed 64,999 cases of early-stage NSCLC treated between 2005 and 2017 from the Florida cancer registry and showed that Black patients had 36% lower odds of receiving curative-intent surgery compared to their White counterparts. This study highlights significant racial disparities in treatment patterns that must be addressed urgently.

Early detection of malignant and pre-malignant peripheral nerve tumors using cell-free DNA fragmentomics.

Early detection of neurofibromatosis type 1 (NF1) associated peripheral nerve sheath tumors (PNST) informs clinical decision-making, potentially averting deadly outcomes. Here, we describe a cell-free DNA (cfDNA) fragmentomic approach which distinguishes non-malignant, pre-malignant and malignant forms of NF1 PNST. Using plasma samples from a novel cohort of 101 NF1 patients and 21 healthy controls, we validated that our previous cfDNA copy number alteration (CNA)-based approach identifies malignant peripheral nerve sheath tumor (MPNST) but cannot distinguish among benign and premalignant states. We therefore investigated the ability of fragment-based cfDNA features to differentiate NF1-associated tumors including binned genome-wide fragment length ratios, end motif analysis, and non-negative matrix factorization deconvolution of fragment lengths. Fragmentomic methods were able to differentiate pre-malignant states including atypical neurofibromas (AN). Fragmentomics also adjudicated AN cases suspicious for MPNST, correctly diagnosing samples noninvasively, which could have informed clinical management. Overall, this study pioneers the early detection of malignant and premalignant peripheral nerve sheath tumors in NF1 patients using plasma cfDNA fragmentomics. In addition to screening applications, this novel approach distinguishes atypical neurofibromas from benign plexiform neurofibromas and malignant peripheral nerve sheath tumors, enabling more precise clinical diagnosis and management.

Sensitive MRD Detection from Lymphatic Fluid after Surgery in HPV-Associated Oropharyngeal Cancer.

PURPOSE: Our goal was to demonstrate that lymphatic drainage fluid (lymph) has improved sensitivity in quantifying postoperative minimal residual disease (MRD) in locally advanced human papillomavirus (HPV)-associated oropharyngeal squamous cell carcinoma (OPSCC) compared with plasma, and leverage this novel biofluid for patient risk stratification. EXPERIMENTAL DESIGN: We prospectively collected lymph samples from neck drains of 106 patients with HPV (+) OPSCC, along with 67 matched plasma samples, 24 hours after surgery. PCR and next-generation sequencing were used to quantify cancer-associated cell-free HPV (cf-HPV) and tumor-informed variants in lymph and plasma. Next, lymph cf-HPV and variants were compared with TNM stage, extranodal extension (ENE), and composite definitions of high-risk pathology. We then created a machine learning model, informed by lymph MRD and clinicopathologic features, to compare with progression-free survival (PFS). RESULTS: Postoperative lymph was enriched with cf-HPV compared with plasma (P < 0.0001) and correlated with pN2 stage (P = 0.003), ENE (P < 0.0001), and trial-defined pathologic risk criteria (mean AUC = 0.78). In addition, the lymph mutation number and variant allele frequency were higher in pN2 ENE (+) necks than in pN1 ENE (+) (P = 0.03, P = 0.02) or pN0-N1 ENE (-) (P = 0.04, P = 0.03, respectively). The lymph MRD-informed risk model demonstrated inferior PFS in high-risk patients (AUC = 0.96, P < 0.0001). CONCLUSIONS: Variant and cf-HPV quantification, performed in 24-hour postoperative lymph samples, reflects single- and multifeature high-risk pathologic criteria. Incorporating lymphatic MRD and clinicopathologic feature analysis can stratify PFS early after surgery in patients with HPV (+) head and neck cancer. See related commentary by Shannon and Iyer, p. 1223.

Genomic and epigenomic analysis of plasma cell-free DNA identifies stemness features associated with worse survival in AR -altered lethal prostate cancer.

Metastatic castration-resistant prostate cancer (mCRPC) resistant to androgen receptor (AR)-targeted agents is often lethal. Unfortunately, biomarkers for this deadly disease remain under investigation, and underpinning mechanisms are ill-understood. Here, we applied deep sequencing to ∼100 mCRPC patients prior to the initiation of first-line AR-targeted therapy, which detected AR /enhancer alterations in over a third of patients, which correlated with lethality. To delve into the mechanism underlying why these patients with cell-free AR /enhancer alterations developed more lethal prostate cancer, we next performed genome-wide cell-free DNA epigenomics. Strikingly, we found that binding sites for transcription factors associated with developmental stemness were nucleosomally more accessible. These results were corroborated using cell-free DNA methylation data, as well as tumor RNA sequencing from a held-out cohort of mCRPC patients. Thus, we validated the importance of AR /enhancer alterations as a prognostic biomarker in lethal mCRPC, and showed that the underlying mechanism for lethality involves reprogramming developmental states toward increased stemness.

High-dimensional deconstruction of pancreatic cancer identifies tumor microenvironmental and developmental stemness features that predict survival.

Numerous cell states are known to comprise the pancreatic ductal adenocarcinoma (PDAC) tumor microenvironment (TME). However, the developmental stemness and co-occurrence of these cell states remain poorly defined. Here, we performed single-cell RNA sequencing (scRNA-seq) on a cohort of treatment-naive PDAC time-of-diagnosis endoscopic ultrasound-guided fine needle biopsy (EUS-FNB) samples (n = 25). We then combined these samples with surgical resection (n = 6) and publicly available samples to increase statistical power (n = 80). Following annotation into 25 distinct cell states, cells were scored for developmental stemness, and a customized version of the Ecotyper tool was used to identify communities of co-occurring cell states in bulk RNA-seq samples (n = 268). We discovered a tumor microenvironmental community comprised of aggressive basal-like malignant cells, tumor-promoting SPP1+ macrophages, and myofibroblastic cancer-associated fibroblasts associated with especially poor prognosis. We also found a developmental stemness continuum with implications for survival that is present in both malignant cells and cancer-associated fibroblasts (CAFs). We further demonstrated that high-dimensional analyses predictive of survival are feasible using standard-of-care, time-of-diagnosis EUS-FNB specimens. In summary, we identified tumor microenvironmental and developmental stemness characteristics from a high-dimensional gene expression analysis of PDAC using human tissue specimens, including time-of-diagnosis EUS-FNB samples. These reveal new connections between tumor microenvironmental composition, CAF and malignant cell stemness, and patient survival that could lead to better upfront risk stratification and more personalized upfront clinical decision-making.

Pre-radiotherapy ctDNA liquid biopsy for risk stratification of oligometastatic non-small cell lung cancer.

The optimal treatment paradigm for patients with oligometastatic non-small cell lung cancer (NSCLC) remains unclear. Some patients with oligometastatic disease experience prolonged remission after locally consolidative radiation therapy (RT), while others harbor micrometastatic disease (below limits of detection by imaging) and benefit from systemic therapy. To risk-stratify and identify the patients most likely to benefit from locally consolidative RT, we performed a multi-institutional cohort study of 1487 patients with oligometastatic NSCLC undergoing liquid biopsy analysis of circulating tumor DNA (ctDNA). In total, 1880 liquid biopsies were performed and approximately 20% of patients (n = 309) had ctDNA measured prior to RT and after their diagnosis of oligometastatic disease. Patients with undetectable ctDNA (pathogenic or likely pathogenic variants in plasma using the Tempus xF assay) before RT had significantly improved progression-free survival (PFS) (P = 0.004) and overall survival (OS) (P = 0.030). ctDNA maximum variant allele frequency (VAF) pre-RT and ctDNA mutational burden pre-RT were both significantly inversely correlated with PFS (maximum VAF P = 0.008, mutational burden P = 0.003) and OS (maximum VAF P = 0.007, mutational burden P = 0.045). These findings were corroborated by multivariate Cox proportional hazards models that included eight additional clinical and genomic parameters. Overall, these data suggest that in patients with oligometastatic NSCLC, pre-RT ctDNA can potentially identify the patients most likely to benefit from locally consolidative RT and experience prolonged PFS and OS. Similarly, ctDNA may be useful to identify undiagnosed micrometastatic disease where it may be appropriate to prioritize systemic therapies.

Circulating tumor DNA based minimal residual disease detection and adjuvant treatment decision-making for muscle-invasive bladder cancer guided by modern clinical trials.

Up to 430,000 cases of bladder cancer are diagnosed each year worldwide. A proposed method for non-invasive monitoring has been to utilize a "liquid biopsy." Liquid biopsy has been proposed as a non-invasive method of testing biomarkers in bodily fluids in order to detect and survey cancer. The liquid biopsy could be utilized to obtain information regarding circulating tumor cells, circulating cell-free tumor DNA, circulating cell-free tumor RNA, and more. It is currently being investigated to help guide adjuvant therapy and improve oncological outcomes. We highlight an array of exciting past and ongoing clinical trials regarding ctDNA and adjuvant therapy in regard to urothelial carcinoma which we believe to be amongst the leaders in the field.

First-in-human prospective trial of sonobiopsy in high-grade glioma patients using neuronavigation-guided focused ultrasound.

Sonobiopsy is an emerging technology that combines focused ultrasound (FUS) with microbubbles to enrich circulating brain disease-specific biomarkers for noninvasive molecular diagnosis of brain diseases. Here, we report the first-in-human prospective trial of sonobiopsy in high-grade glioma patients to evaluate its feasibility and safety in enriching plasma circulating tumor biomarkers. A nimble FUS device integrated with a clinical neuronavigation system was used to perform sonobiopsy following an established clinical workflow for neuronavigation. Analysis of blood samples collected before and after FUS sonication showed that sonobiopsy enriched plasma circulating tumor DNA (ctDNA), including a maximum increase of 1.6-fold for the mononucleosome cell-free DNA (cfDNA) fragments (120-280 bp), 1.9-fold for the patient-specific tumor variant ctDNA level, and 5.6-fold for the TERT mutation ctDNA level. Histological analysis of surgically resected tumors confirmed the safety of the procedure. Transcriptome analysis of sonicated and nonsonicated tumor tissues found that FUS sonication modulated cell physical structure-related genes. Only 2 out of 17,982 total detected genes related to the immune pathways were upregulated. These feasibility and safety data support the continued investigation of sonobiopsy for noninvasive molecular diagnosis of brain diseases.

PACT: a pipeline for analysis of circulating tumor DNA.

MOTIVATION: Detection of genomic alterations in circulating tumor DNA (ctDNA) is currently used for active clinical monitoring of cancer progression and treatment response. While methods for analysis of small mutations are more developed, strategies for detecting structural variants (SVs) in ctDNA are limited. Additionally, reproducibly calling small-scale mutations, copy number alterations, and SVs in ctDNA is challenging due to the lack to unified tools for these different classes of variants. RESULTS: We developed a unified pipeline for the analysis of ctDNA [Pipeline for the Analysis of ctDNA (PACT)] that accurately detects SVs and consistently outperformed similar tools when applied to simulated, cell line, and clinical data. We provide PACT in the form of a Common Workflow Language pipeline which can be run by popular workflow management systems in high-performance computing environments. AVAILABILITY AND IMPLEMENTATION: PACT is freely available at https://github.com/ChrisMaherLab/PACT.

CD4 T cells and toxicity from immune checkpoint blockade.

Immune-related toxicities, otherwise known as immune-related adverse events (irAEs), occur in a substantial fraction of cancer patients treated with immune checkpoint inhibitors (ICIs). Ranging from asymptomatic to life-threatening, ICI-induced irAEs can result in hospital admission, high-dose corticosteroid treatment, ICI discontinuation, and in some cases, death. A deeper understanding of the factors underpinning severe irAE development will be essential for improved irAE prediction and prevention, toward maximizing the benefits and safety profiles of ICIs. In recent work, we applied mass cytometry, single-cell RNA sequencing, single-cell V(D)J sequencing, bulk RNA sequencing, and bulk T-cell receptor (TCR) sequencing to identify pretreatment determinants of severe irAE development in patients with advanced melanoma. Across 71 patients separated into three cohorts, we found that two baseline features in circulation-elevated activated CD4 effector memory T-cell abundance and TCR diversity-are associated with severe irAE development, independent of the affected organ system within 3 months of ICI treatment initiation. Here, we provide an extended perspective on this work, synthesize and discuss related literature, and summarize practical considerations for clinical translation. Collectively, these findings lay a foundation for data-driven and mechanistic insights into irAE development, with the potential to reduce ICI morbidity and mortality in the future.

Emerging Roles of Circulating Tumor DNA for Increased Precision and Personalization in Radiation Oncology.

Recent breakthroughs in circulating tumor DNA (ctDNA) technologies present a compelling opportunity to combine this emerging liquid biopsy approach with the field of radiogenomics, the study of how tumor genomics correlate with radiotherapy response and radiotoxicity. Canonically, ctDNA levels reflect metastatic tumor burden, although newer ultrasensitive technologies can be used after curative-intent radiotherapy of localized disease to assess ctDNA for minimal residual disease (MRD) detection or for post-treatment surveillance. Furthermore, several studies have demonstrated the potential utility of ctDNA analysis across various cancer types managed with radiotherapy or chemoradiotherapy, including sarcoma and cancers of the head and neck, lung, colon, rectum, bladder, and prostate . Additionally, because peripheral blood mononuclear cells are routinely collected alongside ctDNA to filter out mutations associated with clonal hematopoiesis, these cells are also available for single nucleotide polymorphism analysis and could potentially be used to detect patients at high risk for radiotoxicity. Lastly, future ctDNA assays will be utilized to better assess locoregional MRD in order to more precisely guide adjuvant radiotherapy after surgery in cases of localized disease, and guide ablative radiotherapy in cases of oligometastatic disease.

Genomic approaches to cancer and minimal residual disease detection using circulating tumor DNA.

Liquid biopsies using cell-free circulating tumor DNA (ctDNA) are being used frequently in both research and clinical settings. ctDNA can be used to identify actionable mutations to personalize systemic therapy, detect post-treatment minimal residual disease (MRD), and predict responses to immunotherapy. ctDNA can also be isolated from a range of different biofluids, with the possibility of detecting locoregional MRD with increased sensitivity if sampling more proximally than blood plasma. However, ctDNA detection remains challenging in early-stage and post-treatment MRD settings where ctDNA levels are minuscule giving a high risk for false negative results, which is balanced with the risk of false positive results from clonal hematopoiesis. To address these challenges, researchers have developed ever-more elegant approaches to lower the limit of detection (LOD) of ctDNA assays toward the part-per-million range and boost assay sensitivity and specificity by reducing sources of low-level technical and biological noise, and by harnessing specific genomic and epigenomic features of ctDNA. In this review, we highlight a range of modern assays for ctDNA analysis, including advancements made to improve the signal-to-noise ratio. We further highlight the challenge of detecting ultra-rare tumor-associated variants, overcoming which will improve the sensitivity of post-treatment MRD detection and open a new frontier of personalized adjuvant treatment decision-making.

ctDNA Monitoring for Small Cell Lung Cancer: Ready for Prime Time?

Small cell lung cancer (SCLC) is the deadliest form of lung cancer and has few precision medicine approaches available. A recent study analyzed circulating tumor DNA (ctDNA) in 33 patients with extensive-stage SCLC and showed that ctDNA levels and response patterns correlate strongly with clinical response and survival outcomes. See related article by Sivapalan et al., p. 2310.

Epigenetic and transcriptomic characterization reveals progression markers and essential pathways in clear cell renal cell carcinoma.

Identifying tumor-cell-specific markers and elucidating their epigenetic regulation and spatial heterogeneity provides mechanistic insights into cancer etiology. Here, we perform snRNA-seq and snATAC-seq in 34 and 28 human clear cell renal cell carcinoma (ccRCC) specimens, respectively, with matched bulk proteogenomics data. By identifying 20 tumor-specific markers through a multi-omics tiered approach, we reveal an association between higher ceruloplasmin (CP) expression and reduced survival. CP knockdown, combined with spatial transcriptomics, suggests a role for CP in regulating hyalinized stroma and tumor-stroma interactions in ccRCC. Intratumoral heterogeneity analysis portrays tumor cell-intrinsic inflammation and epithelial-mesenchymal transition (EMT) as two distinguishing features of tumor subpopulations. Finally, BAP1 mutations are associated with widespread reduction of chromatin accessibility, while PBRM1 mutations generally increase accessibility, with the former affecting five times more accessible peaks than the latter. These integrated analyses reveal the cellular architecture of ccRCC, providing insights into key markers and pathways in ccRCC tumorigenesis.

First-in-human prospective trial of sonobiopsy in glioblastoma patients using neuronavigation-guided focused ultrasound.

Sonobiopsy is an emerging technology that combines focused ultrasound (FUS) with microbubbles to enrich circulating brain disease-specific biomarkers for noninvasive molecular diagnosis of brain diseases. Here, we report the first-in-human prospective trial of sonobiopsy in glioblastoma patients to evaluate its feasibility and safety in enriching circulating tumor biomarkers. A nimble FUS device integrated with a clinical neuronavigation system was used to perform sonobiopsy following an established clinical workflow for neuronavigation. Analysis of blood samples collected before and after FUS sonication showed enhanced plasma circulating tumor biomarker levels. Histological analysis of surgically resected tumors confirmed the safety of the procedure. Transcriptome analysis of sonicated and unsonicated tumor tissues found that FUS sonication modulated cell physical structure-related genes but evoked minimal inflammatory response. These feasibility and safety data support the continued investigation of sonobiopsy for noninvasive molecular diagnosis of brain diseases.

Pre-radiotherapy ctDNA liquid biopsy for risk stratification of oligometastatic non-small cell lung cancer.

The optimal treatment for patients with oligometastatic non-small cell lung cancer (NSCLC) remains unclear. Some patients with oligometastatic disease can experience prolonged remission after locally consolidative radiation therapy (RT), while others harbor micrometastatic disease (below current limits of detection by imaging) that may benefit from further prioritization of systemic therapy. To better risk-stratify this population and identify the patients most likely to benefit from locally consolidative radiation therapy, we performed a multi-institutional cohort study of patients with oligometastatic NSCLC undergoing liquid biopsy analysis of circulating tumor DNA (ctDNA). Among this real-world cohort of 1,487 patients undergoing analysis (using the Tempus xF assay), a total of 1,880 ctDNA liquid biopsies along with paired clinical data were obtained across various timepoints. Approximately 20% (n=309) of patients had ctDNA obtained prior to RT and after their diagnosis of oligometastatic disease. Samples were de-identified and analyzed for mutational burden and variant frequencies of detectable deleterious (or likely deleterious) mutations in plasma. Patients with undetectable ctDNA before RT had significantly improved progression-free survival and overall survival compared to patients with detectable ctDNA prior to RT. In patients that received RT, 598 pathogenic (or likely deleterious) variants were identified. ctDNA mutational burden pre-RT and ctDNA maximum variant allele frequency (VAF) pre-RT were both significantly inversely correlated with both progression-free (P = 0.0031 for mutational burden, P = 0.0084 for maximum VAF) and overall survival (P = 0.045 for mutational burden, P = 0.0073 for maximum VAF). Patients without detectable ctDNA prior to RT had significantly improved progression-free survival (P = 0.004) and overall survival (P = 0.03) compared to patients with detectable ctDNA prior to RT. These data suggest that in patients with oligometastatic NSCLC, pre-radiotherapy ctDNA analysis can potentially identify the patients most likely to benefit from locally consolidative RT and experience prolonged progression-free and overall survival. Similarly, ctDNA may be useful to identify those patients with undiagnosed micrometastatic disease, in whom it may be appropriate to prioritize systemic therapy.

Urine cell-free DNA multi-omics to detect MRD and predict survival in bladder cancer patients.

Circulating tumor DNA (ctDNA) sensitivity remains subpar for molecular residual disease (MRD) detection in bladder cancer patients. To remedy this problem, we focused on the biofluid most proximal to the disease, urine, and analyzed urine tumor DNA in 74 localized bladder cancer patients. We integrated ultra-low-pass whole genome sequencing (ULP-WGS) with urine cancer personalized profiling by deep sequencing (uCAPP-Seq) to achieve sensitive MRD detection and predict overall survival. Variant allele frequency, inferred tumor mutational burden, and copy number-derived tumor fraction levels in urine cell-free DNA (cfDNA) significantly predicted pathologic complete response status, far better than plasma ctDNA was able to. A random forest model incorporating these urine cfDNA-derived factors with leave-one-out cross-validation was 87% sensitive for predicting residual disease in reference to gold-standard surgical pathology. Both progression-free survival (HR = 3.00, p = 0.01) and overall survival (HR = 4.81, p = 0.009) were dramatically worse by Kaplan-Meier analysis for patients predicted by the model to have MRD, which was corroborated by Cox regression analysis. Additional survival analyses performed on muscle-invasive, neoadjuvant chemotherapy, and held-out validation subgroups corroborated these findings. In summary, we profiled urine samples from 74 patients with localized bladder cancer and used urine cfDNA multi-omics to detect MRD sensitively and predict survival accurately.