Human lung cancer harbors spatially organized stem-immunity hubs associated with response to immunotherapy.

The organization of immune cells in human tumors is not well understood. Immunogenic tumors harbor spatially localized multicellular 'immunity hubs' defined by expression of the T cell-attracting chemokines CXCL10/CXCL11 and abundant T cells. Here, we examined immunity hubs in human pre-immunotherapy lung cancer specimens and found an association with beneficial response to PD-1 blockade. Critically, we discovered the stem-immunity hub, a subtype of immunity hub strongly associated with favorable PD-1-blockade outcome. This hub is distinct from mature tertiary lymphoid structures and is enriched for stem-like TCF7+PD-1+CD8+ T cells, activated CCR7+LAMP3+ dendritic cells and CCL19+ fibroblasts as well as chemokines that organize these cells. Within the stem-immunity hub, we find preferential interactions between CXCL10+ macrophages and TCF7-CD8+ T cells as well as between mature regulatory dendritic cells and TCF7+CD4+ and regulatory T cells. These results provide a picture of the spatial organization of the human intratumoral immune response and its relevance to patient immunotherapy outcomes.

Therapy-associated remodeling of pancreatic cancer revealed by single-cell spatial transcriptomics and optimal transport analysis.

In combination with cell intrinsic properties, interactions in the tumor microenvironment modulate therapeutic response. We leveraged high-plex single-cell spatial transcriptomics to dissect the remodeling of multicellular neighborhoods and cell-cell interactions in human pancreatic cancer associated with specific malignant subtypes and neoadjuvant chemotherapy/radiotherapy. We developed Spatially Constrained Optimal Transport Interaction Analysis (SCOTIA), an optimal transport model with a cost function that includes both spatial distance and ligand-receptor gene expression. Our results uncovered a marked change in ligand-receptor interactions between cancer-associated fibroblasts and malignant cells in response to treatment, which was supported by orthogonal datasets, including an ex vivo tumoroid co-culture system. Overall, this study demonstrates that characterization of the tumor microenvironment using high-plex single-cell spatial transcriptomics allows for identification of molecular interactions that may play a role in the emergence of chemoresistance and establishes a translational spatial biology paradigm that can be broadly applied to other malignancies, diseases, and treatments.

Spatial analysis of human lung cancer reveals organized immune hubs enriched for stem-like CD8 T cells and associated with immunotherapy response.

The organization of immune cells in human tumors is not well understood. Immunogenic tumors harbor spatially-localized multicellular 'immunity hubs' defined by expression of the T cell-attracting chemokines CXCL10/CXCL11 and abundant T cells. Here, we examined immunity hubs in human pre-immunotherapy lung cancer specimens, and found that they were associated with beneficial responses to PD-1-blockade. Immunity hubs were enriched for many interferon-stimulated genes, T cells in multiple differentiation states, and CXCL9/10/11 + macrophages that preferentially interact with CD8 T cells. Critically, we discovered the stem-immunity hub, a subtype of immunity hub strongly associated with favorable PD-1-blockade outcomes, distinct from mature tertiary lymphoid structures, and enriched for stem-like TCF7+PD-1+ CD8 T cells and activated CCR7 + LAMP3 + dendritic cells, as well as chemokines that organize these cells. These results elucidate the spatial organization of the human intratumoral immune response and its relevance to patient immunotherapy outcomes.

Single-nucleus and spatial transcriptome profiling of pancreatic cancer identifies multicellular dynamics associated with neoadjuvant treatment.

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal and treatment-refractory cancer. Molecular stratification in pancreatic cancer remains rudimentary and does not yet inform clinical management or therapeutic development. Here, we construct a high-resolution molecular landscape of the cellular subtypes and spatial communities that compose PDAC using single-nucleus RNA sequencing and whole-transcriptome digital spatial profiling (DSP) of 43 primary PDAC tumor specimens that either received neoadjuvant therapy or were treatment naive. We uncovered recurrent expression programs across malignant cells and fibroblasts, including a newly identified neural-like progenitor malignant cell program that was enriched after chemotherapy and radiotherapy and associated with poor prognosis in independent cohorts. Integrating spatial and cellular profiles revealed three multicellular communities with distinct contributions from malignant, fibroblast and immune subtypes: classical, squamoid-basaloid and treatment enriched. Our refined molecular and cellular taxonomy can provide a framework for stratification in clinical trials and serve as a roadmap for therapeutic targeting of specific cellular phenotypes and multicellular interactions.

Visualizing in deceased COVID-19 patients how SARS-CoV-2 attacks the respiratory and olfactory mucosae but spares the olfactory bulb.

Anosmia, the loss of smell, is a common and often the sole symptom of COVID-19. The onset of the sequence of pathobiological events leading to olfactory dysfunction remains obscure. Here, we have developed a postmortem bedside surgical procedure to harvest endoscopically samples of respiratory and olfactory mucosae and whole olfactory bulbs. Our cohort of 85 cases included COVID-19 patients who died a few days after infection with SARS-CoV-2, enabling us to catch the virus while it was still replicating. We found that sustentacular cells are the major target cell type in the olfactory mucosa. We failed to find evidence for infection of olfactory sensory neurons, and the parenchyma of the olfactory bulb is spared as well. Thus, SARS-CoV-2 does not appear to be a neurotropic virus. We postulate that transient insufficient support from sustentacular cells triggers transient olfactory dysfunction in COVID-19. Olfactory sensory neurons would become affected without getting infected.

Best Practices for Spatial Profiling for Breast Cancer Research with the GeoMx® Digital Spatial Profiler.

Breast cancer is a heterogenous disease with variability in tumor cells and in the surrounding tumor microenvironment (TME). Understanding the molecular diversity in breast cancer is critical for improving prediction of therapeutic response and prognostication. High-plex spatial profiling of tumors enables characterization of heterogeneity in the breast TME, which can holistically illuminate the biology of tumor growth, dissemination and, ultimately, response to therapy. The GeoMx Digital Spatial Profiler (DSP) enables researchers to spatially resolve and quantify proteins and RNA transcripts from tissue sections. The platform is compatible with both formalin-fixed paraffin-embedded and frozen tissues. RNA profiling was developed at the whole transcriptome level for human and mouse samples and protein profiling of 100-plex for human samples. Tissue can be optically segmented for analysis of regions of interest or cell populations to study biology-directed tissue characterization. The GeoMx Breast Cancer Consortium (GBCC) is composed of breast cancer researchers who are developing innovative approaches for spatial profiling to accelerate biomarker discovery. Here, the GBCC presents best practices for GeoMx profiling to promote the collection of high-quality data, optimization of data analysis and integration of datasets to advance collaboration and meta-analyses. Although the capabilities of the platform are presented in the context of breast cancer research, they can be generalized to a variety of other tumor types that are characterized by high heterogeneity.

The CD155/TIGIT axis promotes and maintains immune evasion in neoantigen-expressing pancreatic cancer.

The CD155/TIGIT axis can be co-opted during immune evasion in chronic viral infections and cancer. Pancreatic adenocarcinoma (PDAC) is a highly lethal malignancy, and immune-based strategies to combat this disease have been largely unsuccessful to date. We corroborate prior reports that a substantial portion of PDAC harbors predicted high-affinity MHC class I-restricted neoepitopes and extend these findings to advanced/metastatic disease. Using multiple preclinical models of neoantigen-expressing PDAC, we demonstrate that intratumoral neoantigen-specific CD8+ T cells adopt multiple states of dysfunction, resembling those in tumor-infiltrating lymphocytes of PDAC patients. Mechanistically, genetic and/or pharmacologic modulation of the CD155/TIGIT axis was sufficient to promote immune evasion in autochthonous neoantigen-expressing PDAC. Finally, we demonstrate that the CD155/TIGIT axis is critical in maintaining immune evasion in PDAC and uncover a combination immunotherapy (TIGIT/PD-1 co-blockade plus CD40 agonism) that elicits profound anti-tumor responses in preclinical models, now poised for clinical evaluation.

Spatially organized multicellular immune hubs in human colorectal cancer.

Immune responses to cancer are highly variable, with mismatch repair-deficient (MMRd) tumors exhibiting more anti-tumor immunity than mismatch repair-proficient (MMRp) tumors. To understand the rules governing these varied responses, we transcriptionally profiled 371,223 cells from colorectal tumors and adjacent normal tissues of 28 MMRp and 34 MMRd individuals. Analysis of 88 cell subsets and their 204 associated gene expression programs revealed extensive transcriptional and spatial remodeling across tumors. To discover hubs of interacting malignant and immune cells, we identified expression programs in different cell types that co-varied across tumors from affected individuals and used spatial profiling to localize coordinated programs. We discovered a myeloid cell-attracting hub at the tumor-luminal interface associated with tissue damage and an MMRd-enriched immune hub within the tumor, with activated T cells together with malignant and myeloid cells expressing T cell-attracting chemokines. By identifying interacting cellular programs, we reveal the logic underlying spatially organized immune-malignant cell networks.

Shotgun transcriptome, spatial omics, and isothermal profiling of SARS-CoV-2 infection reveals unique host responses, viral diversification, and drug interactions.

In less than nine months, the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) killed over a million people, including >25,000 in New York City (NYC) alone. The COVID-19 pandemic caused by SARS-CoV-2 highlights clinical needs to detect infection, track strain evolution, and identify biomarkers of disease course. To address these challenges, we designed a fast (30-minute) colorimetric test (LAMP) for SARS-CoV-2 infection from naso/oropharyngeal swabs and a large-scale shotgun metatranscriptomics platform (total-RNA-seq) for host, viral, and microbial profiling. We applied these methods to clinical specimens gathered from 669 patients in New York City during the first two months of the outbreak, yielding a broad molecular portrait of the emerging COVID-19 disease. We find significant enrichment of a NYC-distinctive clade of the virus (20C), as well as host responses in interferon, ACE, hematological, and olfaction pathways. In addition, we use 50,821 patient records to find that renin-angiotensin-aldosterone system inhibitors have a protective effect for severe COVID-19 outcomes, unlike similar drugs. Finally, spatial transcriptomic data from COVID-19 patient autopsy tissues reveal distinct ACE2 expression loci, with macrophage and neutrophil infiltration in the lungs. These findings can inform public health and may help develop and drive SARS-CoV-2 diagnostic, prevention, and treatment strategies.

Temporal and spatial heterogeneity of host response to SARS-CoV-2 pulmonary infection.

The relationship of SARS-CoV-2 pulmonary infection and severity of disease is not fully understood. Here we show analysis of autopsy specimens from 24 patients who succumbed to SARS-CoV-2 infection using a combination of different RNA and protein analytical platforms to characterize inter-patient and intra-patient heterogeneity of pulmonary virus infection. There is a spectrum of high and low virus cases associated with duration of disease. High viral cases have high activation of interferon pathway genes and a predominant M1-like macrophage infiltrate. Low viral cases are more heterogeneous likely reflecting inherent patient differences in the evolution of host response, but there is consistent indication of pulmonary epithelial cell recovery based on napsin A immunohistochemistry and RNA expression of surfactant and mucin genes. Using a digital spatial profiling platform, we find the virus corresponds to distinct spatial expression of interferon response genes demonstrating the intra-pulmonary heterogeneity of SARS-CoV-2 infection.

Comprehensive analysis of cutaneous and uveal melanoma liver metastases.

BACKGROUND: The profound disparity in response to immune checkpoint blockade (ICB) by cutaneous melanoma (CM) and uveal melanoma (UM) patients is not well understood. Therefore, we characterized metastases of CM and UM from the same metastatic site (liver), in order to dissect the potential underlying mechanism in differential response on ICB. METHODS: Tumor liver samples from CM (n=38) and UM (n=28) patients were analyzed at the genomic (whole exome sequencing), transcriptional (RNA sequencing) and protein (immunohistochemistry and GeoMx Digital Spatial Profiling) level. RESULTS: Comparison of CM and UM metastases from the same metastatic site revealed that, although originating from the same melanocyte lineage, CM and UM differed in somatic mutation profile, copy number profile, tumor mutational burden (TMB) and consequently predicted neoantigens. A higher melanin content and higher expression of the melanoma differentiation antigen MelanA was observed in liver metastases of UM patients. No difference in B2M and human leukocyte antigen-DR (HLA-DR) expression was observed. A higher expression of programmed cell death ligand 1 (PD-L1) was found in CM compared with UM liver metastases, although the majority of CM and UM liver metastases lacked PD-L1 expression. There was no difference in the extent of immune infiltration observed between CM and UM metastases, with the exception of a higher expression of CD163 (p<0.0001) in CM liver samples. While the extent of immune infiltration was similar for CM and UM metastases, the ratio of exhausted CD8 T cells to cytotoxic T cells, to total CD8 T cells and to Th1 cells, was significantly higher in UM metastases. CONCLUSIONS: While TMB was different between CM and UM metastases, tumor immune infiltration was similar. The greater dependency on PD-L1 as an immune checkpoint in CM and the identification of higher exhaustion ratios in UM may both serve as explanations for the difference in response to ICB. Consequently, in order to improve current treatment for metastatic UM, reversal of T cell exhaustion beyond programmed cell death 1 blockade should be considered.

Temporal and Spatial Heterogeneity of Host Response to SARS-CoV-2 Pulmonary Infection.

The relationship of SARS-CoV-2 lung infection and severity of pulmonary disease is not fully understood. We analyzed autopsy specimens from 24 patients who succumbed to SARS-CoV-2 infection using a combination of different RNA and protein analytical platforms to characterize inter- and intra- patient heterogeneity of pulmonary virus infection. There was a spectrum of high and low virus cases that was associated with duration of disease and activation of interferon pathway genes. Using a digital spatial profiling platform, the virus corresponded to distinct spatial expression of interferon response genes and immune checkpoint genes demonstrating the intra-pulmonary heterogeneity of SARS-CoV-2 infection.

Perceptions of Personalized Medicine in an Academic Health System: Educational Findings.

OBJECTIVE: Prior reports demonstrate that personalized medicine implementation in clinical care is lacking. Given the program focus at Duke University on personalized medicine, we assessed health care providers' perspectives on their preparation and educational needs to effectively integrate personalized medicine tools and applications into their clinical practices. METHODS: Data from 78 health care providers who participated in a larger study of personalized and precision medicine at Duke University were analyzed using Qualtrics (descriptive statistics). Individuals age 18 years and older were recruited for the larger study through broad email contacts across the university and health system. All participants completed an online 35-question survey that was developed, pilot-tested, and administered by a team of interdisciplinary researchers and clinicians at the Center for Applied Genomics and Precision Medicine. RESULTS: Overall, providers reported being ill-equipped to implement personalized medicine in clinical practice. Many respondents identified educational resources as critical for strengthening personalized medicine implementation in both research and clinical practice. Responses did not differ significantly between specialists and primary providers or by years since completion of the medical degree. CONCLUSIONS: Survey findings support prior calls for provider and patient education in personalized medicine. Respondents identified focus areas in training, education, and research for improving personalized medicine uptake. Given respondents' emphasis on educational needs, now may be an ideal time to address these needs in clinical training and public education programs.

Perspectives on genetic and genomic technologies in an academic medical center: the duke experience.

UNLABELLED: In this age of personalized medicine, genetic and genomic testing is expected to become instrumental in health care delivery, but little is known about its actual implementation in clinical practice. METHODS: We surveyed Duke faculty and healthcare providers to examine the extent of genetic and genomic testing adoption. We assessed providers' use of genetic and genomic testing options and indications in clinical practice, providers' awareness of pharmacogenetic applications, and providers' opinions on returning research-generated genetic test results to participants. Most clinician respondents currently use family history routinely in their clinical practice, but only 18 percent of clinicians use pharmacogenetics. Only two respondents correctly identified the number of drug package inserts with pharmacogenetic indications. We also found strong support for the return of genetic research results to participants. Our results demonstrate that while Duke healthcare providers are enthusiastic about genomic technologies, use of genomic tools outside of research has been limited. Respondents favor return of research-based genetic results to participants, but clinicians lack knowledge about pharmacogenetic applications. We identified challenges faced by this institution when implementing genetic and genomic testing into patient care that should inform a policy and education agenda to improve provider support and clinician-researcher partnerships.

Histone modifications within the human X centromere region.

Human centromeres are multi-megabase regions of highly ordered arrays of alpha satellite DNA that are separated from chromosome arms by unordered alpha satellite monomers and other repetitive elements. Complexities in assembling such large repetitive regions have limited detailed studies of centromeric chromatin organization. However, a genomic map of the human X centromere has provided new opportunities to explore genomic architecture of a complex locus. We used ChIP to examine the distribution of modified histones within centromere regions of multiple X chromosomes. Methylation of H3 at lysine 4 coincided with DXZ1 higher order alpha satellite, the site of CENP-A localization. Heterochromatic histone modifications were distributed across the 400-500 kb pericentromeric regions. The large arrays of alpha satellite and gamma satellite DNA were enriched for both euchromatic and heterochromatic modifications, implying that some pericentromeric repeats have multiple chromatin characteristics. Partial truncation of the X centromere resulted in reduction in the size of the CENP-A/Cenp-A domain and increased heterochromatic modifications in the flanking pericentromere. Although the deletion removed approximately 1/3 of centromeric DNA, the ratio of CENP-A to alpha satellite array size was maintained in the same proportion, suggesting that a limited, but defined linear region of the centromeric DNA is necessary for kinetochore assembly. Our results indicate that the human X centromere contains multiple types of chromatin, is organized similarly to smaller eukaryotic centromeres, and responds to structural changes by expanding or contracting domains.

Oligonucleotide transformation of yeast reveals mismatch repair complexes to be differentially active on DNA replication strands.

Transformation of both prokaryotes and eukaryotes with single-stranded oligonucleotides can transfer sequence information from the oligonucleotide to the chromosome. We have studied this process using oligonucleotides that correct a -1 frameshift mutation in the LYS2 gene of Saccharomyces cerevisiae. We demonstrate that transformation by oligonucleotides occurs preferentially on the lagging strand of replication and is strongly inhibited by the mismatch-repair system. These results are consistent with a mechanism in which oligonucleotides anneal to single-stranded regions of DNA at a replication fork and serve as primers for DNA synthesis. Because the mispairs the primers create are efficiently removed by the mismatch-repair system, single-stranded oligonucleotides can be used to probe mismatch-repair function in a chromosomal context. Removal of mispairs created by annealing of the single-stranded oligonucleotides to the chromosomal DNA is as expected, with 7-nt loops being recognized solely by MutS beta and 1-nt loops being recognized by both MutS alpha and MutS beta. We also find evidence for Mlh1-independent repair of 7-nt, but not 1-nt, loops. Unexpectedly, we find a strand asymmetry of mismatch-repair function; transformation is blocked more efficiently by MutS alpha on the lagging strand of replication, whereas MutS beta does not show a significant strand bias. These results suggest an inherent strand-related difference in how the yeast MutS alpha and MutS beta complexes access and/or repair mismatches that arise in the context of DNA replication.

The effects of mismatch repair and RAD1 genes on interchromosomal crossover recombination in Saccharomyces cerevisiae.

We have previously shown that recombination between 400-bp substrates containing only 4-bp differences, when present in an inverted repeat orientation, is suppressed by >20-fold in wild-type strains of S. cerevisiae. Among the genes involved in this suppression were three genes involved in mismatch repair--MSH2, MSH3, and MSH6--and one in nucleotide excision repair, RAD1. We now report the involvement of these genes in interchromosomal recombination occurring via crossovers using these same short substrates. In these experiments, recombination was stimulated by a double-strand break generated by the HO endonuclease and can occur between completely identical (homologous) substrates or between nonidentical (homeologous) substrates. In addition, a unique feature of this system is that recombining DNA strands can be given a choice of either type of substrate. We find that interchromosomal crossover recombination with these short substrates is severely inhibited in the absence of MSH2, MSH3, or RAD1 and is relatively insensitive to the presence of mismatches. We propose that crossover recombination with these short substrates requires the products of MSH2, MSH3, and RAD1 and that these proteins have functions in recombination in addition to the removal of terminal nonhomology. We further propose that the observed insensitivity to homeology is a result of the difference in recombinational mechanism and/or the timing of the observed recombination events. These results are in contrast with those obtained using longer substrates and may be particularly relevant to recombination events between the abundant short repeated sequences that characterize the genomes of higher eukaryotes.

A new reversion assay for measuring all possible base pair substitutions in Saccharomyces cerevisiae.

A TRP5-based reversion system that allows the rates of all possible base pair substitutions to be measured when the TRP5 locus is in both orientations relative to a defined origin of replication has been developed. This system should be useful for a wide variety of mutation and repair studies in yeast.