Mutations in mediator of RNA polymerase II transcription subunit 12 homolog (MED12, OMIM 300188) cause X-linked intellectual disability (XLID) disorders including FG, Lujan, and Ohdo syndromes. The Gli3-dependent Sonic Hedgehog (SHH) signaling pathway has been implicated in the original FG syndrome and Lujan syndrome. How are SHH-signaling defects related to the complex clinical phenotype of MED12-associated XLID syndromes are not fully understood.

Head and neck squamous cell carcinoma (HNSCC) is largely divided into two groups based on their etiology, human papillomavirus (HPV)-positive and -negative. Global DNA methylation changes are known to drive oncogene and tumor suppressor expression in primary HNSCC of both types. However, significant heterogeneity in DNA methylation within the groups results in different transcriptional profiles and clinical outcomes. We applied a meta-pathway analysis to link gene expression changes to DNA methylation in distinguishing HNSCC subtypes. This approach isolated specific epigenetic changes controlling expression in HPV- HNSCC that distinguish it from HPV+ HNSCC. Analysis of genes identified Hedgehog pathway activation specific to HPV- HNSCC. We confirmed that GLI1, the primary Hedgehog target, showed higher expression in tumors compared to normal samples with HPV- tumors having the highest GLI1 expression, suggesting that increased expression of GLI1 is a potential driver in HPV- HNSCC. Our algorithm for integration of DNA methylation and gene expression can infer biologically significant molecular pathways that may be exploited as therapeutics targets. Our results suggest that therapeutics targeting the Hedgehog pathway may be of benefit in HPV- HNSCC. Similar integrative analysis of high-throughput coupled DNA methylation and expression datasets may yield novel insights into deregulated pathways in other cancers.

Prior studies have demonstrated tumor-specific alternative splicing events in various solid tumor types. The role of alternative splicing in the development and progression of head and neck squamous cell carcinoma (HNSCC) is unclear. Our study queried exon-level expression to implicate splice variants in HNSCC tumors.

Genomics features with similar genome-wide distributions are generally hypothesized to be functionally related, for example, colocalization of histones and transcription start sites indicate chromatin regulation of transcription factor activity. Therefore, statistical algorithms to perform spatial, genome-wide correlation among genomic features are required.

Scientists have unprecedented access to a wide variety of high-quality datasets. These datasets, which are often independently curated, commonly use unstructured spreadsheets to store their data. Standardized annotations are essential to perform synthesis studies across investigators, but are often not used in practice. Therefore, accurately combining records in spreadsheets from differing studies requires tedious and error-prone human curation. These efforts result in a significant time and cost barrier to synthesis research. We propose an information retrieval inspired algorithm, Synthesize, that merges unstructured data automatically based on both column labels and values. Application of the Synthesize algorithm to cancer and ecological datasets had high accuracy (on the order of 85-100%). We further implement Synthesize in an open source web application, Synthesizer (https://github.com/lisagandy/synthesizer). The software accepts input as spreadsheets in comma separated value (CSV) format, visualizes the merged data, and outputs the results as a new spreadsheet. Synthesizer includes an easy to use graphical user interface, which enables the user to finish combining data and obtain perfect accuracy. Future work will allow detection of units to automatically merge continuous data and application of the algorithm to other data formats, including databases.

Screening of patients for cancer-driving mutations is now used for cancer prognosis, remission scoring and treatment selection. Although recently emerged targeted next-generation sequencing-based approaches offer promising diagnostic capabilities, there are still limitations. There is a pressing clinical need for a well-validated, rapid, cost-effective mutation profiling system in patient specimens. Given their speed and cost-effectiveness, quantitative PCR mutation detection techniques are well suited for the clinical environment. The qBiomarker mutation PCR array has high sensitivity and shorter turnaround times compared with other methods. However, a direct comparison with existing viable alternatives are required to assess its true potential and limitations.

Identifying potential resistance mechanisms while tumour cells still respond to therapy is critical to delay acquired resistance.

Hepatocyte growth factor (HGF) signaling through its receptor Met has been implicated in hepatocellular carcinoma tumorigenesis and progression. Met interaction with integrins is shown to modulate the downstream signaling to Akt and ERK (extracellular-regulated kinase). In this study, we developed a mechanistically detailed systems biology model of HGF/Met signaling pathway that incorporated specific interactions with integrins to investigate the efficacy of integrin-binding peptide, AXT050, as monotherapy and in combination with other therapeutics targeting this pathway. Here we report that the modeled dynamics of the response to AXT050 revealed that receptor trafficking is sufficient to explain the effect of Met-integrin interactions on HGF signaling. Furthermore, the model predicted patient-specific synergy and antagonism of efficacy and potency for combination of AXT050 with sorafenib, cabozantinib, and rilotumumab. Overall, the model provides a valuable framework for studying the efficacy of drugs targeting receptor tyrosine kinase interaction with integrins, and identification of synergistic drug combinations for the patients.

Cancer cells disseminate and seed in distant organs, where they can remain dormant for many years before forming clinically detectable metastases. Here we studied how disseminated tumor cells sense and remodel the extracellular matrix (ECM) to sustain dormancy. ECM proteomics revealed that dormant cancer cells assemble a type III collagen-enriched ECM niche. Tumor-derived type III collagen is required to sustain tumor dormancy, as its disruption restores tumor cell proliferation through DDR1-mediated STAT1 signaling. Second-harmonic generation two-photon microscopy further revealed that the dormancy-to-reactivation transition is accompanied by changes in type III collagen architecture and abundance. Analysis of clinical samples revealed that type III collagen levels were increased in tumors from patients with lymph node-negative head and neck squamous cell carcinoma compared to patients who were positive for lymph node colonization. Our data support the idea that the manipulation of these mechanisms could serve as a barrier to metastasis through disseminated tumor cell dormancy induction.

Pancreatic ductal adenocarcinoma (PDAC) is projected to be the second leading cause of cancer death in the USA by 2030. Immune checkpoint inhibitors fail to control most PDAC tumors because of PDAC's extensive immunosuppressive microenvironment and poor immune infiltration, a phenotype also seen in other non-inflamed (ie, 'cold') tumors. Identifying novel ways to enhance immunotherapy efficacy in PDAC is critical. Dipeptidyl peptidase (DPP) inhibition can enhance immunotherapy efficacy in other cancer types; however, the impact of DPP inhibition on PDAC tumors remains unexplored.

Novel immunotherapy combination therapies have improved outcomes for patients with hepatocellular carcinoma (HCC), but responses are limited to a subset of patients and recurrence can also occur. 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. We applied spatial transcriptomics (ST) profiling to characterize the tumor microenvironment in HCC resection specimens from a 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. 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 markers expression suggesting strong activity of these cells. Cancer-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 tumor. The ECM remodeling is consistent with proliferative fibrosis in association with immune-mediated tumor regression. Among the patients with major pathologic response, 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 cancer-CAF interactions and expression of cancer stem cell markers, potentially mediating early tumor immune escape and recurrence in this patient. 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.

The development of single-cell RNA sequencing (scRNA-seq) has allowed high-resolution analysis of cell-type diversity and transcriptional networks controlling cell-fate specification. To identify the transcriptional networks governing human retinal development, we performed scRNA-seq analysis on 16 time points from developing retina as well as four early stages of retinal organoid differentiation. We identified evolutionarily conserved patterns of gene expression during retinal progenitor maturation and specification of all seven major retinal cell types. Furthermore, we identified gene-expression differences between developing macula and periphery and between distinct populations of horizontal cells. We also identified species-specific patterns of gene expression during human and mouse retinal development. Finally, we identified an unexpected role for ATOH7 expression in regulation of photoreceptor specification during late retinogenesis. These results provide a roadmap to future studies of human retinal development and may help guide the design of cell-based therapies for treating retinal dystrophies.

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.

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.

Pancreatic ductal adenocarcinoma (PDAC) frequently presents with metastasis, but the molecular programs in human PDAC cells that drive invasion are not well understood. Using an experimental pipeline enabling PDAC organoid isolation and collection based on invasive phenotype, we assessed the transcriptomic programs associated with invasion in our organoid model. We identified differentially expressed genes in invasive organoids compared with matched noninvasive organoids from the same patients, and we confirmed that the encoded proteins were enhanced in organoid invasive protrusions. We identified 3 distinct transcriptomic groups in invasive organoids, 2 of which correlated directly with the morphological invasion patterns and were characterized by distinct upregulated pathways. Leveraging publicly available single-cell RNA-sequencing data, we mapped our transcriptomic groups onto human PDAC tissue samples, highlighting differences in the tumor microenvironment between transcriptomic groups and suggesting that non-neoplastic cells in the tumor microenvironment can modulate tumor cell invasion. To further address this possibility, we performed computational ligand-receptor analysis and validated the impact of multiple ligands (TGF-β1, IL-6, CXCL12, MMP9) on invasion and gene expression in an independent cohort of fresh human PDAC organoids. Our results identify molecular programs driving morphologically defined invasion patterns and highlight the tumor microenvironment as a potential modulator of these programs.

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.

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.

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.

The epidermal growth factor receptor (EGFR) monoclonal IgG1 antibody cetuximab is approved for first-line treatment of recurrent and metastatic (R/M) HNSCC as a part of the standard of care EXTREME regimen (platinum/5-fluorouracil/cetuximab). This regimen has relatively high response and disease control rates but is generally not curative and many patients will experience recurrent disease and/or metastasis. Therefore, there is a great need to identify predictive biomarkers for recurrence and disease progression in cetuximab-treated HNSCC patients to facilitate patient management and allow for treatment modification. The goal of this work is to assess the potential of activating interleukin-1 (IL-1) ligands (IL-1 alpha [IL-1α], IL-1 beta [IL-1β]) as predictive biomarkers of survival outcomes in HNSCC patients treated with cetuximab-based chemotherapy.

Breast cancer metastasis involves lymphatic dissemination in addition to hematogenous spreading. Although stromal lymphatic vessels (LVs) serve as initial metastatic routes, roles of organ-residing LVs are underinvestigated. Here we show that lymphatic endothelial cells (LECs), a component of LVs within pre-metastatic niches, are conditioned by triple-negative breast cancer (TNBC) cells to accelerate metastasis. LECs within the lungs and lymph nodes, conditioned by tumour-secreted factors, express CCL5 that is not expressed either in normal LECs or in cancer cells, and direct tumour dissemination into these tissues. Moreover, tumour-conditioned LECs promote angiogenesis in these organs, allowing tumour extravasation and colonization. Mechanistically, tumour cell-secreted IL6 causes Stat3 phosphorylation in LECs. This pStat3 induces HIF-1α and VEGF, and a pStat3-pc-Jun-pATF-2 ternary complex induces CCL5 expression in LECs. This study demonstrates anti-metastatic activities of multiple repurposed drugs, blocking a self-reinforcing paracrine loop between breast cancer cells and LECs.