Career trajectories of MD-PhD physician scientists: The loss of women investigators.

Advances in biomedical research require a robust physician scientist workforce. Despite being equally successful at securing early career awards from the NIH as men, women MD-PhD physician scientists are less likely to serve as principal investigators on mid- and later careers awards. Here, we discuss the causes of gender disparities in academic medicine, the implications of losing highly trained women physician scientists, and the institutional and systemic changes needed to sustain this pool of talented investigators.

Transfer Learning Reveals Cancer-Associated Fibroblasts Are Associated with Epithelial-Mesenchymal Transition and Inflammation in Cancer Cells in Pancreatic Ductal Adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive malignancy characterized by an immunosuppressive tumor microenvironment enriched with cancer-associated fibroblasts (CAF). This study used a convergence approach to identify tumor cell and CAF interactions through the integration of single-cell data from human tumors with human organoid coculture experiments. Analysis of a comprehensive atlas of PDAC single-cell RNA sequencing data indicated that CAF density is associated with increased inflammation and epithelial-mesenchymal transition (EMT) in epithelial cells. Transfer learning using transcriptional data from patient-derived organoid and CAF cocultures provided in silico validation of CAF induction of inflammatory and EMT epithelial cell states. Further experimental validation in cocultures demonstrated integrin beta 1 (ITGB1) and vascular endothelial factor A (VEGFA) interactions with neuropilin-1 mediating CAF-epithelial cell cross-talk. Together, this study introduces transfer learning from human single-cell data to organoid coculture analyses for experimental validation of discoveries of cell-cell cross-talk and identifies fibroblast-mediated regulation of EMT and inflammation. SIGNIFICANCE: Adaptation of transfer learning to relate human single-cell RNA sequencing data to organoid-CAF cocultures facilitates discovery of human pancreatic cancer intercellular interactions and uncovers cross-talk between CAFs and tumor cells through VEGFA and ITGB1.

Leveraging multi-omics data to empower quantitative systems pharmacology in immuno-oncology.

Understanding the intricate interactions of cancer cells with the tumor microenvironment (TME) is a pre-requisite for the optimization of immunotherapy. Mechanistic models such as quantitative systems pharmacology (QSP) provide insights into the TME dynamics and predict the efficacy of immunotherapy in virtual patient populations/digital twins but require vast amounts of multimodal data for parameterization. Large-scale datasets characterizing the TME are available due to recent advances in bioinformatics for multi-omics data. Here, we discuss the perspectives of leveraging omics-derived bioinformatics estimates to inform QSP models and circumvent the challenges of model calibration and validation in immuno-oncology.

Personalized neoantigen vaccine and pembrolizumab in advanced hepatocellular carcinoma: a phase 1/2 trial.

Programmed cell death protein 1 (PD-1) inhibitors have modest efficacy as a monotherapy in hepatocellular carcinoma (HCC). A personalized therapeutic cancer vaccine (PTCV) may enhance responses to PD-1 inhibitors through the induction of tumor-specific immunity. We present results from a single-arm, open-label, phase 1/2 study of a DNA plasmid PTCV (GNOS-PV02) encoding up to 40 neoantigens coadministered with plasmid-encoded interleukin-12 plus pembrolizumab in patients with advanced HCC previously treated with a multityrosine kinase inhibitor. Safety and immunogenicity were assessed as primary endpoints, and treatment efficacy and feasibility were evaluated as secondary endpoints. The most common treatment-related adverse events were injection-site reactions, observed in 15 of 36 (41.6%) patients. No dose-limiting toxicities or treatment-related grade ≥3 events were observed. The objective response rate (modified intention-to-treat) per Response Evaluation Criteria in Solid Tumors 1.1 was 30.6% (11 of 36 patients), with 8.3% (3 of 36) of patients achieving a complete response. Clinical responses were associated with the number of neoantigens encoded in the vaccine. Neoantigen-specific T cell responses were confirmed in 19 of 22 (86.4%) evaluable patients by enzyme-linked immunosorbent spot assays. Multiparametric cellular profiling revealed active, proliferative and cytolytic vaccine-specific CD4+ and CD8+ effector T cells. T cell receptor ß-chain (TCRß) bulk sequencing results demonstrated vaccination-enriched T cell clone expansion and tumor infiltration. Single-cell analysis revealed posttreatment T cell clonal expansion of cytotoxic T cell phenotypes. TCR complementarity-determining region cloning of expanded T cell clones in the tumors following vaccination confirmed reactivity against vaccine-encoded neoantigens. Our results support the PTCV's mechanism of action based on the induction of antitumor T cells and show that a PTCV plus pembrolizumab has clinical activity in advanced HCC. ClinicalTrials.gov identifier: NCT04251117 .

MLL1 regulates cytokine-driven cell migration and metastasis.

Cell migration is a critical contributor to metastasis. Cytokine production and its role in cancer cell migration have been traditionally associated with immune cells. We find that the histone methyltransferase Mixed-Lineage Leukemia 1 (MLL1) controls 3D cell migration via cytokines, IL-6, IL-8, and TGF-ß1, secreted by the cancer cells themselves. MLL1, with its scaffold protein Menin, controls actin filament assembly via the IL-6/8/pSTAT3/Arp3 axis and myosin contractility via the TGF-ß1/Gli2/ROCK1/2/pMLC2 axis, which together regulate dynamic protrusion generation and 3D cell migration. MLL1 also regulates cell proliferation via mitosis-based and cell cycle-related pathways. Mice bearing orthotopic MLL1-depleted tumors exhibit decreased lung metastatic burden and longer survival. MLL1 depletion leads to lower metastatic burden even when controlling for the difference in primary tumor growth rates. Combining MLL1-Menin inhibitor with paclitaxel abrogates tumor growth and metastasis, including preexistent metastasis. These results establish MLL1 as a potent regulator of cell migration and highlight the potential of targeting MLL1 in patients with metastatic disease.

Entinostat, nivolumab and ipilimumab for women with advanced HER2-negative breast cancer: a phase Ib trial.

We report the results of 24 women, 50% (N = 12) with hormone receptor-positive breast cancer and 50% (N = 12) with advanced triple-negative breast cancer, treated with entinostat + nivolumab + ipilimumab from the dose escalation (N = 6) and expansion cohort (N = 18) of ETCTN-9844 ( NCT02453620 ). The primary endpoint was safety. Secondary endpoints were overall response rate, clinical benefit rate, progression-free survival and change in tumor CD8:FoxP3 ratio. There were no dose-limiting toxicities. Among evaluable participants (N = 20), the overall response rate was 25% (N = 5), with 40% (N = 4) in triple-negative breast cancer and 10% (N = 1) in hormone receptor-positive breast cancer. The clinical benefit rate was 40% (N = 8), and progression-free survival at 6 months was 50%. Exploratory analyses revealed that changes in myeloid cells may contribute to responses; however, no correlation was noted between changes in CD8:FoxP3 ratio, PD-L1 status and tumor mutational burden and response. These findings support further investigation of this treatment in a phase II trial.

CD4 T cell-activating neoantigens enhance personalized cancer vaccine efficacy.

Personalized cancer vaccines aim to activate and expand cytotoxic antitumor CD8+ T cells to recognize and kill tumor cells. However, the role of CD4+ T cell activation in the clinical benefit of these vaccines is not well defined. We previously established a personalized neoantigen vaccine (PancVAX) for the pancreatic cancer cell line Panc02, which activates tumor-specific CD8+ T cells but required combinatorial checkpoint modulators to achieve therapeutic efficacy. To determine the effects of neoantigen-specific CD4+ T cell activation, we generated a vaccine (PancVAX2) targeting both major histocompatibility complex class I- (MHCI-) and MHCII-specific neoantigens. Tumor-bearing mice vaccinated with PancVAX2 had significantly improved control of tumor growth and long-term survival benefit without concurrent administration of checkpoint inhibitors. PancVAX2 significantly enhanced priming and recruitment of neoantigen-specific CD8+ T cells into the tumor with lower PD-1 expression after reactivation compared with the CD8+ vaccine alone. Vaccine-induced neoantigen-specific Th1 CD4+ T cells in the tumor were associated with decreased Tregs. Consistent with this, PancVAX2 was associated with more proimmune myeloid-derived suppressor cells and M1-like macrophages in the tumor, demonstrating a less immunosuppressive tumor microenvironment. This study demonstrates the biological importance of prioritizing and including CD4+ T cell-specific neoantigens for personalized cancer vaccine modalities.

Systems immunology spanning tumors, lymph nodes, and periphery.

The immune system defines a complex network of tissues and cell types that orchestrate responses across the body in a dynamic manner. The local and systemic interactions between immune and cancer cells contribute to disease progression. Lymphocytes are activated in lymph nodes, traffic through the periphery, and impact cancer progression through their interactions with tumor cells. As a result, therapeutic response and resistance are mediated across tissues, and a comprehensive understanding of lymphocyte dynamics requires a systems-level approach. In this review, we highlight experimental and computational methods that can leverage the study of leukocyte trafficking through an immunomics lens and reveal how adaptive immunity shapes cancer.

Age-associated Senescent - T Cell Signaling Promotes Type 3 Immunity that Inhibits the Biomaterial Regenerative Response.

Aging is associated with immunological changes that compromise response to infections and vaccines, exacerbate inflammatory diseases and can potentially mitigate tissue repair. Even so, age-related changes to the immune response to tissue damage and regenerative medicine therapies remain unknown. Here, it is characterized how aging induces changes in immunological signatures that inhibit tissue repair and therapeutic response to a clinical regenerative biological scaffold derived from extracellular matrix. Signatures of inflammation and interleukin (IL)-17 signaling increased with injury and treatment both locally and regionally in aged animals, and computational analysis uncovered age-associated senescent-T cell communication that promotes type 3 immunity in T cells. Local inhibition of type 3 immune activation using IL17-neutralizing antibodies improves healing and restores therapeutic response to the regenerative biomaterial, promoting muscle repair in older animals. These results provide insights into tissue immune dysregulation that occurs with aging that can be targeted to rejuvenate repair.

Inferring cellular and molecular processes in single-cell data with non-negative matrix factorization using Python, R and GenePattern Notebook implementations of CoGAPS.

Non-negative matrix factorization (NMF) is an unsupervised learning method well suited to high-throughput biology. However, inferring biological processes from an NMF result still requires additional post hoc statistics and annotation for interpretation of learned features. Here, we introduce a suite of computational tools that implement NMF and provide methods for accurate and clear biological interpretation and analysis. A generalized discussion of NMF covering its benefits, limitations and open questions is followed by four procedures for the Bayesian NMF algorithm Coordinated Gene Activity across Pattern Subsets (CoGAPS). Each procedure will demonstrate NMF analysis to quantify cell state transitions in a public domain single-cell RNA-sequencing dataset. The first demonstrates PyCoGAPS, our new Python implementation that enhances runtime for large datasets, and the second allows its deployment in Docker. The third procedure steps through the same single-cell NMF analysis using our R CoGAPS interface. The fourth introduces a beginner-friendly CoGAPS platform using GenePattern Notebook, aimed at users with a working conceptual knowledge of data analysis but without a basic proficiency in the R or Python programming language. We also constructed a user-facing website to serve as a central repository for information and instructional materials about CoGAPS and its application programming interfaces. The expected timing to setup the packages and conduct a test run is around 15 min, and an additional 30 min to conduct analyses on a precomputed result. The expected runtime on the user's desired dataset can vary from hours to days depending on factors such as dataset size or input parameters.

Antibody dependent cell-mediated cytotoxicity selection pressure induces diverse mechanisms of resistance.

Targeted monoclonal antibody therapy has emerged as a powerful therapeutic strategy for cancer. However, only a minority of patients have durable responses and the development of resistance remains a major clinical obstacle. Antibody-dependent cell-mediated cytotoxicity (ADCC) represents a crucial therapeutic mechanism of action; however, few studies have explored ADCC resistance. Using multiple in vitro models of ADCC selection pressure, we have uncovered both shared and distinct resistance mechanisms. Persistent ADCC selection pressure yielded ADCC-resistant cells that are characterized by a loss of NK cell conjugation and this shared resistance phenotype is associated with cell-line dependent modulation of cell surface proteins that contribute to immune synapse formation and NK cell function. We employed single-cell RNA sequencing and proteomic screens to interrogate molecular mechanisms of resistance. We demonstrate that ADCC resistance involves upregulation of interferon/STAT1 and DNA damage response signaling as well as activation of the immunoproteasome. Here, we identify pathways that modulate ADCC sensitivity and report strategies to enhance ADCC-mediated elimination of cancer cells. ADCC resistance could not be reversed with combinatorial treatment approaches. Hence, our findings indicate that tumor cells utilize multiple strategies to inhibit NK cell mediated-ADCC. Future research and development of NK cell-based immunotherapies must incorporate plans to address or potentially prevent the induction of resistance.

Immune landscape of tertiary lymphoid structures in hepatocellular carcinoma (HCC) treated with neoadjuvant immune checkpoint blockade.

Neoadjuvant immunotherapy is thought to produce long-term remissions through induction of antitumor immune responses before removal of the primary tumor. Tertiary lymphoid structures (TLS), germinal center-like structures that can arise within tumors, may contribute to the establishment of immunological memory in this setting, but understanding of their role remains limited. Here, we investigated the contribution of TLS to antitumor immunity in hepatocellular carcinoma (HCC) treated with neoadjuvant immunotherapy. We found that neoadjuvant immunotherapy induced the formation of TLS, which were associated with superior pathologic response, improved relapse free survival, and expansion of the intratumoral T and B cell repertoire. While TLS in viable tumor displayed a highly active mature morphology, in areas of tumor regression we identified an involuted TLS morphology, which was characterized by dispersion of the B cell follicle and persistence of a T cell zone enriched for ongoing antigen presentation and T cell-mature dendritic cell interactions. Involuted TLS showed increased expression of T cell memory markers and expansion of CD8+ cytotoxic and tissue resident memory clonotypes. Collectively, these data reveal the circumstances of TLS dissolution and suggest a functional role for late-stage TLS as sites of T cell memory formation after elimination of viable tumor.

Digitize your Biology! Modeling multicellular systems through interpretable cell behavior.

Cells are fundamental units of life, constantly interacting and evolving as dynamical systems. While recent spatial multi-omics can quantitate individual cells' characteristics and regulatory programs, forecasting their evolution ultimately requires mathematical modeling. We develop a conceptual framework-a cell behavior hypothesis grammar-that uses natural language statements (cell rules) to create mathematical models. This allows us to systematically integrate biological knowledge and multi-omics data to make them computable. We can then perform virtual "thought experiments" that challenge and extend our understanding of multicellular systems, and ultimately generate new testable hypotheses. In this paper, we motivate and describe the grammar, provide a reference implementation, and demonstrate its potential through a series of examples in tumor biology and immunotherapy. Altogether, this approach provides a bridge between biological, clinical, and systems biology researchers for mathematical modeling of biological systems at scale, allowing the community to extrapolate from single-cell characterization to emergent multicellular behavior.

Barcoding intracellular reverse transcription enables high-throughput phenotype-coupled T cell receptor analyses.

Assays linking cellular phenotypes with T cell or B cell antigen receptor sequences are crucial for characterizing adaptive immune responses. Existing methodologies are limited by low sample throughput and high cost. Here, we present INtraCEllular Reverse Transcription with Sorting and sequencing (INCERTS), an approach that combines molecular indexing of receptor repertoires within intact cells and fluorescence-activated cell sorting (FACS). We demonstrate that INCERTS enables efficient processing of millions of cells from pooled human peripheral blood mononuclear cell (PBMC) samples while retaining robust association between T cell receptor (TCR) sequences and cellular phenotypes. We used INCERTS to discover antigen-specific TCRs from patients with cancer immunized with a novel mutant KRAS peptide vaccine. After ex vivo stimulation, 28 uniquely barcoded samples were pooled prior to FACS into peptide-reactive and non-reactive CD4+ and CD8+ populations. Combining complementary patient-matched single-cell RNA sequencing (scRNA-seq) data enabled retrieval of full-length, paired TCR alpha and beta chain sequences for future validation of therapeutic utility.

Spatial transcriptomics analysis of neoadjuvant cabozantinib and nivolumab in advanced hepatocellular carcinoma identifies independent mechanisms of resistance and recurrence.

BACKGROUND: Novel immunotherapy combination therapies have improved outcomes for patients with hepatocellular carcinoma (HCC), but responses are limited to a subset of patients. Little is known about the inter- and intra-tumor heterogeneity in cellular signaling networks within the HCC tumor microenvironment (TME) that underlie responses to modern systemic therapy. METHODS: We applied spatial transcriptomics (ST) profiling to characterize the tumor microenvironment in HCC resection specimens from a prospective clinical trial of neoadjuvant cabozantinib, a multi-tyrosine kinase inhibitor that primarily blocks VEGF, and nivolumab, a PD-1 inhibitor in which 5 out of 15 patients were found to have a pathologic response at the time of resection. RESULTS: ST profiling demonstrated that the TME of responding tumors was enriched for immune cells and cancer-associated fibroblasts (CAF) with pro-inflammatory signaling relative to the non-responders. The enriched cancer-immune interactions in responding tumors are characterized by activation of the PAX5 module, a known regulator of B cell maturation, which colocalized with spots with increased B cell marker expression suggesting strong activity of these cells. HCC-CAF interactions were also enriched in the responding tumors and were associated with extracellular matrix (ECM) remodeling as there was high activation of FOS and JUN in CAFs adjacent to the tumor. The ECM remodeling is consistent with proliferative fibrosis in association with immune-mediated tumor regression. Among the patients with major pathologic responses, a single patient experienced early HCC recurrence. ST analysis of this clinical outlier demonstrated marked tumor heterogeneity, with a distinctive immune-poor tumor region that resembles the non-responding TME across patients and was characterized by HCC-CAF interactions and expression of cancer stem cell markers, potentially mediating early tumor immune escape and recurrence in this patient. CONCLUSIONS: These data show that responses to modern systemic therapy in HCC are associated with distinctive molecular and cellular landscapes and provide new targets to enhance and prolong responses to systemic therapy in HCC.

Transcriptomic forecasting with neural ordinary differential equations.

Single-cell transcriptomics technologies can uncover changes in the molecular states that underlie cellular phenotypes. However, understanding the dynamic cellular processes requires extending from inferring trajectories from snapshots of cellular states to estimating temporal changes in cellular gene expression. To address this challenge, we have developed a neural ordinary differential-equation-based method, RNAForecaster, for predicting gene expression states in single cells for multiple future time steps in an embedding-independent manner. We demonstrate that RNAForecaster can accurately predict future expression states in simulated single-cell transcriptomic data with cellular tracking over time. We then show that by using metabolic labeling single-cell RNA sequencing (scRNA-seq) data from constitutively dividing cells, RNAForecaster accurately recapitulates many of the expected changes in gene expression during progression through the cell cycle over a 3-day period. Thus, RNAForecaster enables short-term estimation of future expression states in biological systems from high-throughput datasets with temporal information.

Informing virtual clinical trials of hepatocellular carcinoma with spatial multi-omics analysis of a human neoadjuvant immunotherapy clinical trial.

Human clinical trials are important tools to advance novel systemic therapies improve treatment outcomes for cancer patients. The few durable treatment options have led to a critical need to advance new therapeutics in hepatocellular carcinoma (HCC). Recent human clinical trials have shown that new combination immunotherapeutic regimens provide unprecedented clinical response in a subset of patients. Computational methods that can simulate tumors from mathematical equations describing cellular and molecular interactions are emerging as promising tools to simulate the impact of therapy entirely in silico. To facilitate designing dosing regimen and identifying potential biomarkers, we developed a new computational model to track tumor progression at organ scale while reflecting the spatial heterogeneity in the tumor at tissue scale in HCC. This computational model is called a spatial quantitative systems pharmacology (spQSP) platform and it is also designed to simulate the effects of combination immunotherapy. We then validate the results from the spQSP system by leveraging real-world spatial multi-omics data from a neoadjuvant HCC clinical trial combining anti-PD-1 immunotherapy and a multitargeted tyrosine kinase inhibitor (TKI) cabozantinib. The model output is compared with spatial data from Imaging Mass Cytometry (IMC). Both IMC data and simulation results suggest closer proximity between CD8 T cell and macrophages among non-responders while the reverse trend was observed for responders. The analyses also imply wider dispersion of immune cells and less scattered cancer cells in responders' samples. We also compared the model output with Visium spatial transcriptomics analyses of samples from post-treatment tumor resections in the original clinical trial. Both spatial transcriptomic data and simulation results identify the role of spatial patterns of tumor vasculature and TGFß in tumor and immune cell interactions. To our knowledge, this is the first spatial tumor model for virtual clinical trials at a molecular scale that is grounded in high-throughput spatial multi-omics data from a human clinical trial.

Boosting Single-Cell RNA Sequencing Analysis with Simple Neural Attention.

A limitation of current deep learning (DL) approaches for single-cell RNA sequencing (scRNAseq) analysis is the lack of interpretability. Moreover, existing pipelines are designed and trained for specific tasks used disjointly for different stages of analysis. We present scANNA, a novel interpretable DL model for scRNAseq studies that leverages neural attention to learn gene associations. After training, the learned gene importance (interpretability) is used to perform downstream analyses (e.g., global marker selection and cell-type classification) without retraining. ScANNA's performance is comparable to or better than state-of-the-art methods designed and trained for specific standard scRNAseq analyses even though scANNA was not trained for these tasks explicitly. ScANNA enables researchers to discover meaningful results without extensive prior knowledge or training separate task-specific models, saving time and enhancing scRNAseq analyses.