Organ-specific immunity: A tissue analysis framework for investigating local immune responses to SARS-CoV-2.

Local immune activation at mucosal surfaces, mediated by mucosal lymphoid tissues, is vital for effective immune responses against pathogens. While pathogens like severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can spread to multiple organs, patients with coronavirus disease 2019 (COVID-19) primarily experience inflammation and damage in their lungs. To investigate this apparent organ-specific immune response, we develop an analytical framework that recognizes the significance of mucosal lymphoid tissues. This framework combines histology, immunofluorescence, spatial transcript profiling, and mathematical modeling to identify cellular and gene expression differences between the lymphoid tissues of the lung and the gut and predict the determinants of those differences. Our findings indicate that mucosal lymphoid tissues are pivotal in organ-specific immune response to SARS-CoV-2, mediating local inflammation and tissue damage and contributing to immune dysfunction. The framework developed here has potential utility in the study of long COVID and may streamline biomarker discovery and treatment design for diseases with differential pathologies at the organ level.

Spatially resolved whole transcriptome profiling in human and mouse tissue using Digital Spatial Profiling.

Emerging spatial profiling technology has enabled high-plex molecular profiling in biological tissues, preserving the spatial and morphological context of gene expression. Here, we describe expanding the chemistry for the Digital Spatial Profiling platform to quantify whole transcriptomes in human and mouse tissues using a wide range of spatial profiling strategies and sample types. We designed multiplexed in situ hybridization probes targeting the protein-coding genes of the human and mouse transcriptomes, referred to as the human or mouse Whole Transcriptome Atlas (WTA). Human and mouse WTAs were validated in cell lines for concordance with orthogonal gene expression profiling methods in regions ranging from ∼10-500 cells. By benchmarking against bulk RNA-seq and fluorescence in situ hybridization, we show robust transcript detection down to ∼100 transcripts per region. To assess the performance of WTA across tissue and sample types, we applied WTA to biological questions in cancer, molecular pathology, and developmental biology. Spatial profiling with WTA detected expected gene expression differences between tumor and tumor microenvironment, identified disease-specific gene expression heterogeneity in histological structures of the human kidney, and comprehensively mapped transcriptional programs in anatomical substructures of nine organs in the developing mouse embryo. Digital Spatial Profiling technology with the WTA assays provides a flexible method for spatial whole transcriptome profiling applicable to diverse tissue types and biological contexts.

Rate volatility and asymmetric segregation diversify mutation burden in cells with mutator alleles.

Mutations that compromise mismatch repair (MMR) or DNA polymerase ε or δ exonuclease domains produce mutator phenotypes capable of fueling cancer evolution. Here, we investigate how combined defects in these pathways expands genetic heterogeneity in cells of the budding yeast, Saccharomyces cerevisiae, using a single-cell resolution approach that tallies all mutations arising from individual divisions. The distribution of replication errors present in mother cells after the initial S-phase was broader than expected for a single uniform mutation rate across all cell divisions, consistent with volatility of the mutator phenotype. The number of mismatches that then segregated to the mother and daughter cells co-varied, suggesting that each division is governed by a different underlying genome-wide mutation rate. The distribution of mutations that individual cells inherit after the second S-phase is further broadened by the sequential actions of semiconservative replication and mitotic segregation of chromosomes. Modeling suggests that this asymmetric segregation may diversify mutation burden in mutator-driven tumors.

Multiplex digital spatial profiling of proteins and RNA in fixed tissue.

Digital Spatial Profiling (DSP) is a method for highly multiplex spatial profiling of proteins or RNAs suitable for use on formalin-fixed, paraffin-embedded (FFPE) samples. The approach relies on (1) multiplexed readout of proteins or RNAs using oligonucleotide tags; (2) oligonucleotide tags attached to affinity reagents (antibodies or RNA probes) through a photocleavable (PC) linker; and (3) photocleaving light projected onto the tissue sample to release PC oligonucleotides in any spatial pattern across a region of interest (ROI) covering 1 to ~5,000 cells. DSP is capable of single-cell sensitivity within an ROI using the antibody readout, with RNA detection feasible down to ~600 individual mRNA transcripts. We show spatial profiling of up to 44 proteins and 96 genes (928 RNA probes) in lymphoid, colorectal tumor and autoimmune tissues by using the nCounter system and 1,412 genes (4,998 RNA probes) by using next-generation sequencing (NGS). DSP may be used to profile not only proteins and RNAs in biobanked samples but also immune markers in patient samples, with potential prognostic and predictive potential for clinical decision-making.

Neutrophil content predicts lymphocyte depletion and anti-PD1 treatment failure in NSCLC.

Immune checkpoint inhibitor (ICI) treatment has recently become a first-line therapy for many non-small cell lung cancer (NSCLC) patients. Unfortunately, most NSCLC patients are refractory to ICI monotherapy, and initial attempts to address this issue with secondary therapeutics have proven unsuccessful. To identify entities precluding CD8+ T cell accumulation in this process, we performed unbiased analyses on flow cytometry, gene expression, and multiplexed immunohistochemical data from a NSCLC patient cohort. The results revealed the presence of a myeloid-rich subgroup, which was devoid of CD4+ and CD8+ T cells. Of all myeloid cell types assessed, neutrophils were the most highly associated with the myeloid phenotype. Additionally, the ratio of CD8+ T cells to neutrophils (CD8/PMN) within the tumor mass optimally distinguished between active and myeloid cases. This ratio was also capable of showing the separation of patients responsive to ICI therapy from those with stable or progressive disease in 2 independent cohorts. Tumor-bearing mice treated with a combination of anti-PD1 and SX-682 (CXCR1/2 inhibitor) displayed relocation of lymphocytes from the tumor periphery into a malignant tumor, which was associated with induction of IFN-γ-responsive genes. These results suggest that neutrophil antagonism may represent a viable secondary therapeutic strategy to enhance ICI treatment outcomes.

Normally lethal amino acid substitutions suppress an ultramutator DNA Polymerase δ variant.

In yeast, the pol3-01,L612M double mutant allele, which causes defects in DNA polymerase delta (Pol δ) proofreading (pol3-01) and nucleotide selectivity (pol3-L612M), confers an "ultramutator" phenotype that rapidly drives extinction of haploid and diploid MMR-proficient cells. Here, we investigate antimutator mutations that encode amino acid substitutions in Pol δ that suppress this lethal phenotype. We find that most of the antimutator mutations individually suppress the pol3-01 and pol3-L612M mutator phenotypes. The locations of many of the amino acid substitutions in Pol δ resemble those of previously identified antimutator substitutions; however, two novel mutations encode substitutions (R674G and Q697R) of amino acids in the fingers domain that coordinate the incoming dNTP. These mutations are lethal without pol3-L612M and markedly change the mutation spectra produced by the pol3-01,L612M mutator allele, suggesting that they alter nucleotide selection to offset the pol3-L612M mutator phenotype. Consistent with this hypothesis, mutations and drug treatments that perturb dNTP pool levels disproportionately influence the viability of pol3-L612M,R674G and pol3-L612M,Q697R cells. Taken together, our findings suggest that mutation rate can evolve through genetic changes that alter the balance of dNTP binding and dissociation from DNA polymerases.