Machine Learning Links T-cell Function and Spatial Localization to Neoadjuvant Immunotherapy and Clinical Outcome in Pancreatic Cancer.

Tumor molecular data sets are becoming increasingly complex, making it nearly impossible for humans alone to effectively analyze them. Here, we demonstrate the power of using machine learning (ML) to analyze a single-cell, spatial, and highly multiplexed proteomic data set from human pancreatic cancer and reveal underlying biological mechanisms that may contribute to clinical outcomes. We designed a multiplex immunohistochemistry antibody panel to compare T-cell functionality and spatial localization in resected tumors from treatment-naïve patients with localized pancreatic ductal adenocarcinoma (PDAC) with resected tumors from a second cohort of patients treated with neoadjuvant agonistic CD40 (anti-CD40) monoclonal antibody therapy. In total, nearly 2.5 million cells from 306 tissue regions collected from 29 patients across both cohorts were assayed, and over 1,000 tumor microenvironment (TME) features were quantified. We then trained ML models to accurately predict anti-CD40 treatment status and disease-free survival (DFS) following anti-CD40 therapy based on TME features. Through downstream interpretation of the ML models' predictions, we found anti-CD40 therapy reduced canonical aspects of T-cell exhaustion within the TME, as compared with treatment-naïve TMEs. Using automated clustering approaches, we found improved DFS following anti-CD40 therapy correlated with an increased presence of CD44+CD4+ Th1 cells located specifically within cellular neighborhoods characterized by increased T-cell proliferation, antigen experience, and cytotoxicity in immune aggregates. Overall, our results demonstrate the utility of ML in molecular cancer immunology applications, highlight the impact of anti-CD40 therapy on T cells within the TME, and identify potential candidate biomarkers of DFS for anti-CD40-treated patients with PDAC.

Imputing Single-Cell Protein Abundance in Multiplex Tissue Imaging.

Multiplex tissue imaging are a collection of increasingly popular single-cell spatial proteomics and transcriptomics assays for characterizing biological tissues both compositionally and spatially. However, several technical issues limit the utility of multiplex tissue imaging, including the limited number of RNAs and proteins that can be assayed, tissue loss, and protein probe failure. In this work, we demonstrate how machine learning methods can address these limitations by imputing protein abundance at the single-cell level using multiplex tissue imaging datasets from a breast cancer cohort. We first compared machine learning methods' strengths and weaknesses for imputing single-cell protein abundance. Machine learning methods used in this work include regularized linear regression, gradient-boosted regression trees, and deep learning autoencoders. We also incorporated cellular spatial information to improve imputation performance. Using machine learning, single-cell protein expression can be imputed with mean absolute error ranging between 0.05-0.3 on a [0,1] scale. Our results demonstrate (1) the feasibility of imputing single-cell abundance levels for many proteins using machine learning to overcome the technical constraints of multiplex tissue imaging and (2) how including cellular spatial information can substantially enhance imputation results.

Machine learning links T cell function and spatial localization to neoadjuvant immunotherapy and clinical outcome in pancreatic cancer.

Tumor molecular datasets are becoming increasingly complex, making it nearly impossible for humans alone to effectively analyze them. Here, we demonstrate the power of using machine learning to analyze a single-cell, spatial, and highly multiplexed proteomic dataset from human pancreatic cancer and reveal underlying biological mechanisms that may contribute to clinical outcome. A novel multiplex immunohistochemistry antibody panel was used to audit T cell functionality and spatial localization in resected tumors from treatment-naive patients with localized pancreatic ductal adenocarcinoma (PDAC) compared to a second cohort of patients treated with neoadjuvant agonistic CD40 (αCD40) monoclonal antibody therapy. In total, nearly 2.5 million cells from 306 tissue regions collected from 29 patients across both treatment cohorts were assayed, and more than 1,000 tumor microenvironment (TME) features were quantified. We then trained machine learning models to accurately predict αCD40 treatment status and disease-free survival (DFS) following αCD40 therapy based upon TME features. Through downstream interpretation of the machine learning models' predictions, we found αCD40 therapy to reduce canonical aspects of T cell exhaustion within the TME, as compared to treatment-naive TMEs. Using automated clustering approaches, we found improved DFS following αCD40 therapy to correlate with the increased presence of CD44+ CD4+ Th1 cells located specifically within cellular spatial neighborhoods characterized by increased T cell proliferation, antigen-experience, and cytotoxicity in immune aggregates. Overall, our results demonstrate the utility of machine learning in molecular cancer immunology applications, highlight the impact of αCD40 therapy on T cells within the TME, and identify potential candidate biomarkers of DFS for αCD40-treated patients with PDAC.

An omic and multidimensional spatial atlas from serial biopsies of an evolving metastatic breast cancer.

Mechanisms of therapeutic resistance and vulnerability evolve in metastatic cancers as tumor cells and extrinsic microenvironmental influences change during treatment. To support the development of methods for identifying these mechanisms in individual people, here we present an omic and multidimensional spatial (OMS) atlas generated from four serial biopsies of an individual with metastatic breast cancer during 3.5 years of therapy. This resource links detailed, longitudinal clinical metadata that includes treatment times and doses, anatomic imaging, and blood-based response measurements to clinical and exploratory analyses, which includes comprehensive DNA, RNA, and protein profiles; images of multiplexed immunostaining; and 2- and 3-dimensional scanning electron micrographs. These data report aspects of heterogeneity and evolution of the cancer genome, signaling pathways, immune microenvironment, cellular composition and organization, and ultrastructure. We present illustrative examples of how integrative analyses of these data reveal potential mechanisms of response and resistance and suggest novel therapeutic vulnerabilities.

Inverting the model of genomics data sharing with the NHGRI Genomic Data Science Analysis, Visualization, and Informatics Lab-space.

The NHGRI Genomic Data Science Analysis, Visualization, and Informatics Lab-space (AnVIL; https://anvilproject.org) was developed to address a widespread community need for a unified computing environment for genomics data storage, management, and analysis. In this perspective, we present AnVIL, describe its ecosystem and interoperability with other platforms, and highlight how this platform and associated initiatives contribute to improved genomic data sharing efforts. The AnVIL is a federated cloud platform designed to manage and store genomics and related data, enable population-scale analysis, and facilitate collaboration through the sharing of data, code, and analysis results. By inverting the traditional model of data sharing, the AnVIL eliminates the need for data movement while also adding security measures for active threat detection and monitoring and provides scalable, shared computing resources for any researcher. We describe the core data management and analysis components of the AnVIL, which currently consists of Terra, Gen3, Galaxy, RStudio/Bioconductor, Dockstore, and Jupyter, and describe several flagship genomics datasets available within the AnVIL. We continue to extend and innovate the AnVIL ecosystem by implementing new capabilities, including mechanisms for interoperability and responsible data sharing, while streamlining access management. The AnVIL opens many new opportunities for analysis, collaboration, and data sharing that are needed to drive research and to make discoveries through the joint analysis of hundreds of thousands to millions of genomes along with associated clinical and molecular data types.

Single-cell spatial architectures associated with clinical outcome in head and neck squamous cell carcinoma.

There is increasing evidence that the spatial organization of cells within the tumor-immune microenvironment (TiME) of solid tumors influences survival and response to therapy in numerous cancer types. Here, we report results and demonstrate the applicability of quantitative single-cell spatial proteomics analyses in the TiME of primary and recurrent human papillomavirus (HPV)-negative head and neck squamous cell carcinoma (HNSCC) tumors. Single-cell compositions of a nine patient, primary and recurrent (n = 18), HNSCC cohort is presented, followed by deeper investigation into the spatial architecture of the TiME and its relationship with clinical variables and progression free survival (PFS). Multiple spatial algorithms were used to quantify the spatial landscapes of immune cells within TiMEs and demonstrate that neoplastic tumor-immune cell spatial compartmentalization, rather than mixing, is associated with longer PFS. Mesenchymal (αSMA+) cellular neighborhoods describe distinct immune landscapes associated with neoplastic tumor-immune compartmentalization and improved patient outcomes. Results from this investigation are concordant with studies in other tumor types, suggesting that trends in TiME cellular heterogeneity and spatial organization may be shared across cancers and may provide prognostic value in multiple cancer types.

MCMICRO: a scalable, modular image-processing pipeline for multiplexed tissue imaging.

Highly multiplexed tissue imaging makes detailed molecular analysis of single cells possible in a preserved spatial context. However, reproducible analysis of large multichannel images poses a substantial computational challenge. Here, we describe a modular and open-source computational pipeline, MCMICRO, for performing the sequential steps needed to transform whole-slide images into single-cell data. We demonstrate the use of MCMICRO on tissue and tumor images acquired using multiple imaging platforms, thereby providing a solid foundation for the continued development of tissue imaging software.

Multiomics analysis of serial PARP inhibitor treated metastatic TNBC inform on rational combination therapies.

In a pilot study, we evaluated the feasibility of real-time deep analysis of serial tumor samples from triple negative breast cancer patients to identify mechanisms of resistance and treatment opportunities as they emerge under therapeutic stress engendered by poly-ADP-ribose polymerase (PARP) inhibitors (PARPi). In a BRCA-mutant basal breast cancer exceptional long-term survivor, a striking tumor destruction was accompanied by a marked infiltration of immune cells containing CD8 effector cells, consistent with pre-clinical evidence for association between STING mediated immune activation and benefit from PARPi and immunotherapy. Tumor cells in the exceptional responder underwent extensive protein network rewiring in response to PARP inhibition. In contrast, there were minimal changes in the ecosystem of a luminal androgen receptor rapid progressor, likely due to indifference to the effects of PARP inhibition. Together, identification of PARPi-induced emergent changes could be used to select patient specific combination therapies, based on tumor and immune state changes.

Galaxy-ML: An accessible, reproducible, and scalable machine learning toolkit for biomedicine.

Supervised machine learning is an essential but difficult to use approach in biomedical data analysis. The Galaxy-ML toolkit (https://galaxyproject.org/community/machine-learning/) makes supervised machine learning more accessible to biomedical scientists by enabling them to perform end-to-end reproducible machine learning analyses at large scale using only a web browser. Galaxy-ML extends Galaxy (https://galaxyproject.org), a biomedical computational workbench used by tens of thousands of scientists across the world, with a suite of tools for all aspects of supervised machine learning.

Leukocyte Heterogeneity in Pancreatic Ductal Adenocarcinoma: Phenotypic and Spatial Features Associated with Clinical Outcome.

Immunotherapies targeting aspects of T cell functionality are efficacious in many solid tumors, but pancreatic ductal adenocarcinoma (PDAC) remains refractory to these treatments. Deeper understanding of the PDAC immune ecosystem is needed to identify additional therapeutic targets and predictive biomarkers for therapeutic response and resistance monitoring. To address these needs, we quantitatively evaluated leukocyte contexture in 135 human PDACs at single-cell resolution by profiling density and spatial distribution of myeloid and lymphoid cells within histopathologically defined regions of surgical resections from treatment-naive and presurgically (neoadjuvant)-treated patients and biopsy specimens from metastatic PDAC. Resultant data establish an immune atlas of PDAC heterogeneity, identify leukocyte features correlating with clinical outcomes, and, through an in silico study, provide guidance for use of PDAC tissue microarrays to optimally measure intratumoral immune heterogeneity. Atlas data have direct applicability as a reference for evaluating immune responses to investigational neoadjuvant PDAC therapeutics where pretherapy baseline specimens are not available. SIGNIFICANCE: We provide a phenotypic and spatial immune atlas of human PDAC identifying leukocyte composition at steady state and following standard neoadjuvant therapies. These data have broad utility as a resource that can inform on leukocyte responses to emerging therapies where baseline tissues were not acquired.This article is highlighted in the In This Issue feature, p. 1861.

A long-read RNA-seq approach to identify novel transcripts of very large genes.

RNA-seq is widely used for studying gene expression, but commonly used sequencing platforms produce short reads that only span up to two exon junctions per read. This makes it difficult to accurately determine the composition and phasing of exons within transcripts. Although long-read sequencing improves this issue, it is not amenable to precise quantitation, which limits its utility for differential expression studies. We used long-read isoform sequencing combined with a novel analysis approach to compare alternative splicing of large, repetitive structural genes in muscles. Analysis of muscle structural genes that produce medium (Nrap: 5 kb), large (Neb: 22 kb), and very large (Ttn: 106 kb) transcripts in cardiac muscle, and fast and slow skeletal muscles identified unannotated exons for each of these ubiquitous muscle genes. This also identified differential exon usage and phasing for these genes between the different muscle types. By mapping the in-phase transcript structures to known annotations, we also identified and quantified previously unannotated transcripts. Results were confirmed by endpoint PCR and Sanger sequencing, which revealed muscle-type-specific differential expression of these novel transcripts. The improved transcript identification and quantification shown by our approach removes previous impediments to studies aimed at quantitative differential expression of ultralong transcripts.

The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2020 update.

Galaxy (https://galaxyproject.org) is a web-based computational workbench used by tens of thousands of scientists across the world to analyze large biomedical datasets. Since 2005, the Galaxy project has fostered a global community focused on achieving accessible, reproducible, and collaborative research. Together, this community develops the Galaxy software framework, integrates analysis tools and visualizations into the framework, runs public servers that make Galaxy available via a web browser, performs and publishes analyses using Galaxy, leads bioinformatics workshops that introduce and use Galaxy, and develops interactive training materials for Galaxy. Over the last two years, all aspects of the Galaxy project have grown: code contributions, tools integrated, users, and training materials. Key advances in Galaxy's user interface include enhancements for analyzing large dataset collections as well as interactive tools for exploratory data analysis. Extensions to Galaxy's framework include support for federated identity and access management and increased ability to distribute analysis jobs to remote resources. New community resources include large public servers in Europe and Australia, an increasing number of regional and local Galaxy communities, and substantial growth in the Galaxy Training Network.

G-OnRamp: Generating genome browsers to facilitate undergraduate-driven collaborative genome annotation.

Scientists are sequencing new genomes at an increasing rate with the goal of associating genome contents with phenotypic traits. After a new genome is sequenced and assembled, structural gene annotation is often the first step in analysis. Despite advances in computational gene prediction algorithms, most eukaryotic genomes still benefit from manual gene annotation. This requires access to good genome browsers to enable annotators to visualize and evaluate multiple lines of evidence (e.g., sequence similarity, RNA sequencing [RNA-Seq] results, gene predictions, repeats) and necessitates many volunteers to participate in the work. To address the technical barriers to creating genome browsers, the Genomics Education Partnership (GEP; https://gep.wustl.edu/) has partnered with the Galaxy Project (https://galaxyproject.org) to develop G-OnRamp (http://g-onramp.org), a web-based platform for creating UCSC Genome Browser Assembly Hubs and JBrowse genome browsers. G-OnRamp also converts a JBrowse instance into an Apollo instance for collaborative genome annotations in research and educational settings. The genome browsers produced can be transferred to the CyVerse Data Store for long-term access. G-OnRamp enables researchers to easily visualize their experimental results, educators to create Course-based Undergraduate Research Experiences (CUREs) centered on genome annotation, and students to participate in genomics research. In the process, students learn about genes/genomes and about how to utilize large datasets. Development of G-OnRamp was guided by extensive user feedback. Sixty-five researchers/educators from >40 institutions participated through in-person workshops, which produced >20 genome browsers now available for research and education. Genome browsers generated for four parasitoid wasp species have been used in a CURE engaging students at 15 colleges and universities. Our assessment results in the classroom demonstrate that the genome browsers produced by G-OnRamp are effective tools for engaging undergraduates in research and in enabling their contributions to the scientific literature in genomics. Expansion of such genomics research/education partnerships will be beneficial to researchers, faculty, and students alike.

The Human Tumor Atlas Network: Charting Tumor Transitions across Space and Time at Single-Cell Resolution.

Crucial transitions in cancer-including tumor initiation, local expansion, metastasis, and therapeutic resistance-involve complex interactions between cells within the dynamic tumor ecosystem. Transformative single-cell genomics technologies and spatial multiplex in situ methods now provide an opportunity to interrogate this complexity at unprecedented resolution. The Human Tumor Atlas Network (HTAN), part of the National Cancer Institute (NCI) Cancer Moonshot Initiative, will establish a clinical, experimental, computational, and organizational framework to generate informative and accessible three-dimensional atlases of cancer transitions for a diverse set of tumor types. This effort complements both ongoing efforts to map healthy organs and previous large-scale cancer genomics approaches focused on bulk sequencing at a single point in time. Generating single-cell, multiparametric, longitudinal atlases and integrating them with clinical outcomes should help identify novel predictive biomarkers and features as well as therapeutically relevant cell types, cell states, and cellular interactions across transitions. The resulting tumor atlases should have a profound impact on our understanding of cancer biology and have the potential to improve cancer detection, prevention, and therapeutic discovery for better precision-medicine treatments of cancer patients and those at risk for cancer.

How Machine Learning Will Transform Biomedicine.

This Perspective explores the application of machine learning toward improved diagnosis and treatment. We outline a vision for how machine learning can transform three broad areas of biomedicine: clinical diagnostics, precision treatments, and health monitoring, where the goal is to maintain health through a range of diseases and the normal aging process. For each area, early instances of successful machine learning applications are discussed, as well as opportunities and challenges for machine learning. When these challenges are met, machine learning promises a future of rigorous, outcomes-based medicine with detection, diagnosis, and treatment strategies that are continuously adapted to individual and environmental differences.

A harmonized meta-knowledgebase of clinical interpretations of somatic genomic variants in cancer.

Precision oncology relies on accurate discovery and interpretation of genomic variants, enabling individualized diagnosis, prognosis and therapy selection. We found that six prominent somatic cancer variant knowledgebases were highly disparate in content, structure and supporting primary literature, impeding consensus when evaluating variants and their relevance in a clinical setting. We developed a framework for harmonizing variant interpretations to produce a meta-knowledgebase of 12,856 aggregate interpretations. We demonstrated large gains in overlap between resources across variants, diseases and drugs as a result of this harmonization. We subsequently demonstrated improved matching between a patient cohort and harmonized interpretations of potential clinical significance, observing an increase from an average of 33% per individual knowledgebase to 57% in aggregate. Our analyses illuminate the need for open, interoperable sharing of variant interpretation data. We also provide a freely available web interface (search.cancervariants.org) for exploring the harmonized interpretations from these six knowledgebases.

High-dimensional multiplexed immunohistochemical characterization of immune contexture in human cancers.

Biomarker assessments of tumor specimens is widely used in cancer research to audit tumor cell intrinsic as well as tumor cell extrinsic features including the diversity of immune, stromal, and mesenchymal cells. To comprehensively and quantitatively audit the tumor-immune microenvironment (TiME), we developed a novel multiplex immunohistochemistry (mIHC) platform and computational image processing workflow using a single formalin-fixed paraffin-embedded (FFPE) tissue section. Herein, we validated this platform using nine matched primary newly diagnosed and recurrent head and neck squamous cell carcinoma (HNSCC) sections sequentially subjected to immunodetection with a panel of 29 antibodies identifying malignant tumor cells, and 17 distinct leukocyte lineages and their functional states. Image cytometric analysis was applied to interpret chromogenic signals from digitally scanned and coregistered light microscopy-based images enabling identification and quantification of individual tumor cells, structural features, immune cell phenotypes and their functional state. In agreement with our previous study via a 12-plex imaging mIHC platform, myeloid-inflamed status in newly diagnosed primary tumors associated with significantly short progression free survival, independent of lymphoid-inflamed status. Spatial distribution of tumor and immune cell lineages in TiME was also examined and revealed statistically significant CD8+ T cell exclusion from tumor nests, whereas regulatory T cells and myeloid cells, when present in close proximity to tumor cells, highly associated with rapid cancer recurrence. These findings indicate presence of differential immune-spatial profiles in newly diagnosed and recurrent HNSCC, and establish the robustness of the 29-plex mIHC platform and associated analytics for quantitative analysis of single tissue sections revealing longitudinal TiME changes.

Cloud bursting galaxy: federated identity and access management.

MOTIVATION: Large biomedical datasets, such as those from genomics and imaging, are increasingly being stored on commercial and institutional cloud computing platforms. This is because cloud-scale computing resources, from robust backup to high-speed data transfer to scalable compute and storage, are needed to make these large datasets usable. However, one challenge for large-scale biomedical data on the cloud is providing secure access, especially when datasets are distributed across platforms. While there are open Web protocols for secure authentication and authorization, these protocols are not in wide use in bioinformatics and are difficult to use for even technologically sophisticated users. RESULTS: We have developed a generic and extensible approach for securely accessing biomedical datasets distributed across cloud computing platforms. Our approach combines OpenID Connect and OAuth2, best-practice Web protocols for authentication and authorization, together with Galaxy (https://galaxyproject.org), a web-based computational workbench used by thousands of scientists across the world. With our enhanced version of Galaxy, users can access and analyze data distributed across multiple cloud computing providers without any special knowledge of access/authorization protocols. Our approach does not require users to share permanent credentials (e.g. username, password, API key), instead relying on automatically generated temporary tokens that refresh as needed. Our approach is generalizable to most identity providers and cloud computing platforms. To the best of our knowledge, Galaxy is the only computational workbench where users can access biomedical datasets across multiple cloud computing platforms using best-practice Web security approaches and thereby minimize risks of unauthorized data access and credential use. AVAILABILITY AND IMPLEMENTATION: Freely available for academic and commercial use under the open-source Academic Free License (https://opensource.org/licenses/AFL-3.0) from the following Github repositories: https://github.com/galaxyproject/galaxy and https://github.com/galaxyproject/cloudauthz.

Innate αß T Cells Mediate Antitumor Immunity by Orchestrating Immunogenic Macrophage Programming.

Unconventional T-lymphocyte populations are emerging as important regulators of tumor immunity. Despite this, the role of TCRαß+CD4-CD8-NK1.1- innate αß T cells (iαßT) in pancreatic ductal adenocarcinoma (PDA) has not been explored. We found that iαßTs represent ∼10% of T lymphocytes infiltrating PDA in mice and humans. Intratumoral iαßTs express a distinct T-cell receptor repertoire and profoundly immunogenic phenotype compared with their peripheral counterparts and conventional lymphocytes. iαßTs comprised ∼75% of the total intratumoral IL17+ cells. Moreover, iαßT-cell adoptive transfer is protective in both murine models of PDA and human organotypic systems. We show that iαßT cells induce a CCR5-dependent immunogenic macrophage reprogramming, thereby enabling marked CD4+ and CD8+ T-cell expansion/activation and tumor protection. Collectively, iαßTs govern fundamental intratumoral cross-talk between innate and adaptive immune populations and are attractive therapeutic targets. SIGNIFICANCE: We found that iαßTs are a profoundly activated T-cell subset in PDA that slow tumor growth in murine and human models of disease. iαßTs induce a CCR5-dependent immunogenic tumor-associated macrophage program, T-cell activation and expansion, and should be considered as novel targets for immunotherapy.See related commentary by Banerjee et al., p. 1164.This article is highlighted in the In This Issue feature, p. 1143.