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.

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.

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.

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.