High-plex imaging of RNA and proteins at subcellular resolution in fixed tissue by spatial molecular imaging.

Resolving the spatial distribution of RNA and protein in tissues at subcellular resolution is a challenge in the field of spatial biology. We describe spatial molecular imaging, a system that measures RNAs and proteins in intact biological samples at subcellular resolution by performing multiple cycles of nucleic acid hybridization of fluorescent molecular barcodes. We demonstrate that spatial molecular imaging has high sensitivity (one or two copies per cell) and very low error rate (0.0092 false calls per cell) and background (~0.04 counts per cell). The imaging system generates three-dimensional, super-resolution localization of analytes at ~2 million cells per sample. Cell segmentation is morphology based using antibodies, compatible with formalin-fixed, paraffin-embedded samples. We measured multiomic data (980 RNAs and 108 proteins) at subcellular resolution in formalin-fixed, paraffin-embedded tissues (nonsmall cell lung and breast cancer) and identified >18 distinct cell types, ten unique tumor microenvironments and 100 pairwise ligand-receptor interactions. Data on >800,000 single cells and ~260 million transcripts can be accessed at http://nanostring.com/CosMx-dataset .

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.

High-Plex Spatial RNA Profiling Reveals Cell Type‒Specific Biomarker Expression during Melanoma Development.

Early diagnosis of melanoma is critical for improved survival. However, the biomarkers of early melanoma evolution and their origin within the tumor and its microenvironment, including the keratinocytes, are poorly defined. To address this, we used spatial transcript profiling that maintains the morphological tumor context to measure the expression of >1,000 RNAs in situ in patient-derived formalin-fixed, paraffin-embedded tissue sections in primary melanoma and melanocytic nevi. We profiled 134 regions of interest (each 200 µm in diameter) enriched in melanocytes, neighboring keratinocytes, or immune cells. This approach captured distinct expression patterns across cell types and tumor types during melanoma development. Unexpectedly, we discovered that S100A8 is expressed by keratinocytes within the tumor microenvironment during melanoma growth. Immunohistochemistry of 252 tumors showed prominent keratinocyte-derived S100A8 expression in melanoma but not in benign tumors and confirmed the same pattern for S100A8's binding partner S100A9, suggesting that injury to the epidermis may be an early and readily detectable indicator of melanoma development. Together, our results establish a framework for high-plex, spatial, and cell type‒specific resolution of gene expression in archival tissue applicable to the development of biomarkers and characterization of tumor microenvironment interactions in tumor evolution.

Spatial proteomic characterization of HER2-positive breast tumors through neoadjuvant therapy predicts response.

The addition of HER2-targeted agents to neoadjuvant chemotherapy has dramatically improved pathological complete response (pCR) rates in early-stage, HER2-positive breast cancer. Nonetheless, up to 50% of patients have residual disease after treatment, while others are likely overtreated. Here, we performed multiplex spatial proteomic characterization of 122 samples from 57 HER2-positive breast tumors from the neoadjuvant TRIO-US B07 clinical trial sampled pre-treatment, after 14-21 d of HER2-targeted therapy and at surgery. We demonstrated that proteomic changes after a single cycle of HER2-targeted therapy aids the identification of tumors that ultimately undergo pCR, outperforming pre-treatment measures or transcriptomic changes. We further developed and validated a classifier that robustly predicted pCR using a single marker, CD45, measured on treatment, and showed that CD45-positive cell counts measured via conventional immunohistochemistry perform comparably. These results demonstrate robust biomarkers that can be used to enable the stratification of sensitive tumors early during neoadjuvant HER2-targeted therapy, with implications for tailoring subsequent therapy.

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.

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.

Inter- and intra-tumor heterogeneity of metastatic prostate cancer determined by digital spatial gene expression profiling.

Metastatic prostate cancer (mPC) comprises a spectrum of diverse phenotypes. However, the extent of inter- and intra-tumor heterogeneity is not established. Here we use digital spatial profiling (DSP) technology to quantitate transcript and protein abundance in spatially-distinct regions of mPCs. By assessing multiple discrete areas across multiple metastases, we find a high level of intra-patient homogeneity with respect to tumor phenotype. However, there are notable exceptions including tumors comprised of regions with high and low androgen receptor (AR) and neuroendocrine activity. While the vast majority of metastases examined are devoid of significant inflammatory infiltrates and lack PD1, PD-L1 and CTLA4, the B7-H3/CD276 immune checkpoint protein is highly expressed, particularly in mPCs with high AR activity. Our results demonstrate the utility of DSP for accurately classifying tumor phenotype, assessing tumor heterogeneity, and identifying aspects of tumor biology involving the immunological composition of metastases.

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.

Loss of Heterozygosity Drives Adaptation in Hybrid Yeast.

Hybridization is often considered maladaptive, but sometimes hybrids can invade new ecological niches and adapt to novel or stressful environments better than their parents. The genomic changes that occur following hybridization that facilitate genome resolution and/or adaptation are not well understood. Here, we examine hybrid genome evolution using experimental evolution of de novo interspecific hybrid yeast Saccharomyces cerevisiae × Saccharomyces uvarum and their parentals. We evolved these strains in nutrient-limited conditions for hundreds of generations and sequenced the resulting cultures identifying numerous point mutations, copy number changes, and loss of heterozygosity (LOH) events, including species-biased amplification of nutrient transporters. We focused on a particularly interesting example, in which we saw repeated LOH at the high-affinity phosphate transporter gene PHO84 in both intra- and interspecific hybrids. Using allele replacement methods, we tested the fitness of different alleles in hybrid and S. cerevisiae strain backgrounds and found that the LOH is indeed the result of selection on one allele over the other in both S. cerevisiae and the hybrids. This is an example where hybrid genome resolution is driven by positive selection on existing heterozygosity and demonstrates that even infrequent outcrossing may have lasting impacts on adaptation.

Differential paralog divergence modulates genome evolution across yeast species.

Evolutionary outcomes depend not only on the selective forces acting upon a species, but also on the genetic background. However, large timescales and uncertain historical selection pressures can make it difficult to discern such important background differences between species. Experimental evolution is one tool to compare evolutionary potential of known genotypes in a controlled environment. Here we utilized a highly reproducible evolutionary adaptation in Saccharomyces cerevisiae to investigate whether experimental evolution of other yeast species would select for similar adaptive mutations. We evolved populations of S. cerevisiae, S. paradoxus, S. mikatae, S. uvarum, and interspecific hybrids between S. uvarum and S. cerevisiae for ~200-500 generations in sulfate-limited continuous culture. Wild-type S. cerevisiae cultures invariably amplify the high affinity sulfate transporter gene, SUL1. However, while amplification of the SUL1 locus was detected in S. paradoxus and S. mikatae populations, S. uvarum cultures instead selected for amplification of the paralog, SUL2. We measured the relative fitness of strains bearing deletions and amplifications of both SUL genes from different species, confirming that, converse to S. cerevisiae, S. uvarum SUL2 contributes more to fitness in sulfate limitation than S. uvarum SUL1. By measuring the fitness and gene expression of chimeric promoter-ORF constructs, we were able to delineate the cause of this differential fitness effect primarily to the promoter of S. uvarum SUL1. Our data show evidence of differential sub-functionalization among the sulfate transporters across Saccharomyces species through recent changes in noncoding sequence. Furthermore, these results show a clear example of how such background differences due to paralog divergence can drive changes in genome evolution.

Genome-wide quantification of rare somatic mutations in normal human tissues using massively parallel sequencing.

We present the bottleneck sequencing system (BotSeqS), a next-generation sequencing method that simultaneously quantifies rare somatic point mutations across the mitochondrial and nuclear genomes. BotSeqS combines molecular barcoding with a simple dilution step immediately before library amplification. We use BotSeqS to show age- and tissue-dependent accumulations of rare mutations and demonstrate that somatic mutational burden in normal human tissues can vary by several orders of magnitude, depending on biologic and environmental factors. We further show major differences between the mutational patterns of the mitochondrial and nuclear genomes in normal tissues. Lastly, the mutation spectra of normal tissues were different from each other, but similar to those of the cancers that arose in them. This technology can provide insights into the number and nature of genetic alterations in normal tissues and can be used to address a variety of fundamental questions about the genomes of diseased tissues.

Aristolochic Acid in the Etiology of Renal Cell Carcinoma.

BACKGROUND: Aristolochia species used in the practice of traditional herbal medicine contains aristolochic acid (AA), an established human carcinogen contributing to urothelial carcinomas of the upper urinary tract. AA binds covalently to genomic DNA, forming aristolactam (AL)-DNA adducts. Here we investigated whether AA is also an etiologic factor in clear cell renal cell carcinoma (ccRCC). METHODS: We conducted a population-based case-control study to investigate the linkage between Aristolochia prescription history, cumulative AA consumption, and ccRCC incidence in Taiwan (5,709 cases and 22,836 matched controls). The presence and level of mutagenic dA-AL-I adducts were determined in the kidney DNA of 51 Taiwanese ccRCC patients. The whole-exome sequences of ccRCC tumors from 10 Taiwanese ccRCC patients with prior exposure to AA were determined. RESULTS: Cumulative ingestion of more than 250 mg of AA increased risk of ccRCC (OR, 1.25), and we detected dA-AL-I adducts in 76% of Taiwanese ccRCC patients. Furthermore, the distinctive AA mutational signature was evident in six of 10 sequenced ccRCC exomes from Taiwanese patients. CONCLUSIONS: This study strongly suggests that AA contributes to the etiology of certain RCCs. IMPACT: The current study offers compelling evidence implicating AA in a significant fraction of the RCC arising in Taiwan and illustrates the power of integrating epidemiologic, molecular, and genetic data in the investigation of cancer etiology. Cancer Epidemiol Biomarkers Prev; 25(12); 1600-8. ©2016 AACR.

Mutational signature of aristolochic acid exposure as revealed by whole-exome sequencing.

In humans, exposure to aristolochic acid (AA) is associated with urothelial carcinoma of the upper urinary tract (UTUC). Exome sequencing of UTUCs from 19 individuals with documented exposure to AA revealed a remarkably large number of somatic mutations and an unusual mutational signature attributable to AA. Most of the mutations (72%) in these tumors were A:T-to-T:A transversions, located predominantly on the nontranscribed strand, with a strong preference for deoxyadenosine in a consensus sequence (T/CAG). This trinucleotide motif overlaps the canonical splice acceptor site, possibly accounting for the excess of splice site mutations observed in these tumors. The AA mutational fingerprint was found frequently in oncogenes and tumor suppressor genes in AA-associated UTUC. The AA mutational signature was observed in one patient's tumor from a UTUC cohort without previous indication of AA exposure. Together, these results directly link an established environmental mutagen to cancer through genome-wide sequencing and highlight its power to reveal individual exposure to carcinogens.

DNA resection at chromosome breaks promotes genome stability by constraining non-allelic homologous recombination.

DNA double-strand breaks impact genome stability by triggering many of the large-scale genome rearrangements associated with evolution and cancer. One of the first steps in repairing this damage is 5'→3' resection beginning at the break site. Recently, tools have become available to study the consequences of not extensively resecting double-strand breaks. Here we examine the role of Sgs1- and Exo1-dependent resection on genome stability using a non-selective assay that we previously developed using diploid yeast. We find that Saccharomyces cerevisiae lacking Sgs1 and Exo1 retains a very efficient repair process that is highly mutagenic to genome structure. Specifically, 51% of cells lacking Sgs1 and Exo1 repair a double-strand break using repetitive sequences 12-48 kb distal from the initial break site, thereby generating a genome rearrangement. These Sgs1- and Exo1-independent rearrangements depend partially upon a Rad51-mediated homologous recombination pathway. Furthermore, without resection a robust cell cycle arrest is not activated, allowing a cell with a single double-strand break to divide before repair, potentially yielding multiple progeny each with a different rearrangement. This profusion of rearranged genomes suggests that cells tolerate any dangers associated with extensive resection to inhibit mutagenic pathways such as break-distal recombination. The activation of break-distal recipient repeats and amplification of broken chromosomes when resection is limited raise the possibility that genome regions that are difficult to resect may be hotspots for rearrangements. These results may also explain why mutations in resection machinery are associated with cancer.

Tracking chromosome dynamics in live yeast cells: coordinated movement of rDNA homologs and anaphase disassembly of the nucleolus during meiosis.

A prerequisite for determination of chromosome dynamics in live cells is development of a method for staining or marking the chromosome of interest. We describe here a unique chromosome-tracking system that differentially marks two large chromosome segments from homologs in the budding yeast Saccharomyces cerevisiae. Using yeast genetics and the special features at the repetitive ribosomal RNA (rRNA) gene cluster, we incorporated arrays of the tet operator and the lac operator into each repeat of the two rDNA homologs by homologous recombination. Expression of tet repressor-fused green fluorescent protein and lac repressor-fused red fluorescent protein in engineered cells led to the differential labeling of rDNA homologs. Using live-cell three-dimensional fluorescence microscopy, we showed that homologs undergo contraction and expansion cycles in an actin-dependent manner during meiosis and that chromosome mobility appears to be correlated with nuclear positioning. Our observations further revealed that, in contrast to mitosis, in meiosis the yeast nucleolus, the site of rRNA processing, was disassembled upon anaphase onset, suggesting a differential regulation of the rDNA array during meiotic chromosome segregation. Because rRNA genes are highly conserved, a similar chromosome-engineering approach may be adaptable in other eukaryotes for functional assays of chromosome organization in live cells.

Competitive repair by naturally dispersed repetitive DNA during non-allelic homologous recombination.

Genome rearrangements often result from non-allelic homologous recombination (NAHR) between repetitive DNA elements dispersed throughout the genome. Here we systematically analyze NAHR between Ty retrotransposons using a genome-wide approach that exploits unique features of Saccharomyces cerevisiae purebred and Saccharomyces cerevisiae/Saccharomyces bayanus hybrid diploids. We find that DNA double-strand breaks (DSBs) induce NAHR-dependent rearrangements using Ty elements located 12 to 48 kilobases distal to the break site. This break-distal recombination (BDR) occurs frequently, even when allelic recombination can repair the break using the homolog. Robust BDR-dependent NAHR demonstrates that sequences very distal to DSBs can effectively compete with proximal sequences for repair of the break. In addition, our analysis of NAHR partner choice between Ty repeats shows that intrachromosomal Ty partners are preferred despite the abundance of potential interchromosomal Ty partners that share higher sequence identity. This competitive advantage of intrachromosomal Tys results from the relative efficiencies of different NAHR repair pathways. Finally, NAHR generates deleterious rearrangements more frequently when DSBs occur outside rather than within a Ty repeat. These findings yield insights into mechanisms of repeat-mediated genome rearrangements associated with evolution and cancer.

Persistent mechanical linkage between sister chromatids throughout anaphase.

In budding yeast, we have found that sister rDNA arrays marked with fluorescent probes can be visualized as two distinguishable strands during metaphase. Upon anaphase, these arm loci are drawn into the spindle, where they adopt a cruciform-like structure and stretch 2.5-fold as they migrate to the poles. Therefore, while sister rDNA arrays appear separated in metaphase, mechanical linkages between sister arm loci persist throughout anaphase in yeast, as shown in grasshopper spermatocytes (Paliulis and Nicklas 2004). These linkages are partially dependent on the protector of cohesin, SGO1. In anaphase, the spatially regulated dissolution of these mechanical linkages serves to prevent premature sister separation and restrain the rate of spindle elongation. Thus, sister separation is temporally controlled and linkages between sister chromatids contribute to the regulation of anaphase spindle elongation.

Structural changes in Mcm5 protein bypass Cdc7-Dbf4 function and reduce replication origin efficiency in Saccharomyces cerevisiae.

Eukaryotic chromosomal replication is a complicated process with many origins firing at different efficiencies and times during S phase. Prereplication complexes are assembled on all origins in G(1) phase, and yet only a subset of complexes is activated during S phase by DDK (for Dbf4-dependent kinase) (Cdc7-Dbf4). The yeast mcm5-bob1 (P83L) mutation bypasses DDK but results in reduced intrinsic firing efficiency at 11 endogenous origins and at origins located on minichromosomes. Origin efficiency may result from Mcm5 protein assuming an altered conformation, as predicted from the atomic structure of an archaeal MCM (for minichromosome maintenance) homologue. Similarly, an intragenic mutation in a residue predicted to interact with P83L suppresses the mcm5-bob1 bypass phenotype. We propose DDK phosphorylation of the MCM complex normally results in a single, highly active conformation of Mcm5, whereas the mcm5-bob1 mutation produces a number of conformations, only one of which is permissive for origin activation. Random adoption of these alternate states by the mcm5-bob1 protein can explain both how origin firing occurs independently of DDK and why origin efficiency is reduced. Because similar mutations in mcm2 and mcm4 cannot bypass DDK, Mcm5 protein may be a unique Mcm protein that is the final target of DDK regulation.