Background: The Human Cell Atlas is a large international collaborative effort to map all cell types of the human body. Single-cell RNA sequencing can generate high-quality data for the delivery of such an atlas. However, delays between fresh sample collection and processing may lead to poor data and difficulties in experimental design. Results: This study assesses the effect of cold storage on fresh healthy spleen, esophagus, and lung from ≥ 5 donors over 72 h. We collect 240,000 high-quality single-cell transcriptomes with detailed cell type annotations and whole genome sequences of donors, enabling future eQTL studies. Our data provide a valuable resource for the study of these 3 organs and will allow cross-organ comparison of cell types. We see little effect of cold ischemic time on cell yield, total number of reads per cell, and other quality control metrics in any of the tissues within the first 24 h. However, we observe a decrease in the proportions of lung T cells at 72 h, higher percentage of mitochondrial reads, and increased contamination by background ambient RNA reads in the 72-h samples in the spleen, which is cell type specific. Conclusions: In conclusion, we present robust protocols for tissue preservation for up to 24 h prior to scRNA-seq analysis. This greatly facilitates the logistics of sample collection for Human Cell Atlas or clinical studies since it increases the time frames for sample processing.
Identifying the origin of cancer: Many cancers are classified on the basis of the organ or tissue from which they originated. However, identifying the specific cells and conditions that precede tumorigenesis can help us understand and better treat the resulting disease. Nowicki-Osuch et al . used a single-cell approach to investigate the cell of origin for Barrett’s esophagus (BE) and the mechanisms leading to the development of esophageal adenocarcinoma (EAC) (see the Perspective by Geboes and Hoorens). Analyses of healthy human esophageal tissues, mutational lineage tracing, and organoid models revealed that BE originates from the gastric cardia and that EAC arises from undifferentiated BE cells. This analysis provides a map of the transcriptional landscape of the healthy esophagus that can be compared with mouse models of disease. —LMZ
Intestinal metaplasia in the esophagus (Barrett's esophagus IM, or BE-IM) and stomach (GIM) are considered precursors for esophageal and gastric adenocarcinoma, respectively. We hypothesize that BE-IM and GIM follow parallel developmental trajectories in response to differing inflammatory insults. Here, we construct a single-cell RNA-sequencing atlas, supported by protein expression studies, of the entire gastrointestinal tract spanning physiologically normal and pathologic states including gastric metaplasia in the esophagus (E-GM), BE-IM, atrophic gastritis, and GIM. We demonstrate that BE-IM and GIM share molecular features, and individual cells simultaneously possess transcriptional properties of gastric and intestinal epithelia, suggesting phenotypic mosaicism. Transcriptionally E-GM resembles atrophic gastritis; genetically, it is clonal and has a lower mutational burden than BE-IM. Finally, we show that GIM and BE-IM acquire a protumorigenic, activated fibroblast microenvironment. These findings suggest that BE-IM and GIM can be considered molecularly similar entities in adjacent organs, opening the path for shared detection and treatment strategies. Significance: Our data capture the gradual molecular and phenotypic transition from a gastric to intestinal phenotype (IM) in the esophagus and stomach. Because BE-IM and GIM can predispose to cancer, this new understanding of a common developmental trajectory could pave the way for a more unified approach to detection and treatment. See related commentary by Stachler, p. 1291. This article is highlighted in the In This Issue feature, p. 1275
Barrett’s esophagus is a pre-malignant lesion that can progress to esophageal adenocarcinoma. We perform a multi-omic analysis of pre-cancer samples from 146 patients with a range of outcomes, comprising 642 person years of follow-up. Whole genome sequencing reveals complex structural variants and LINE-1 retrotransposons, as well as known copy number changes, occurring even prior to dysplasia. The structural variant burden captures the most variance across the cohort and genomic profiles do not always match consensus clinical pathology dysplasia grades. Increasing structural variant burden is associated with: high levels of chromothripsis and breakage-fusion-bridge events; increased expression of genes related to cell cycle checkpoint, DNA repair and chromosomal instability; and epigenetic silencing of Wnt signalling and cell cycle genes. Timing analysis reveals molecular events triggering genomic instability with more clonal expansion in dysplastic samples. Overall genomic complexity occurs early in the Barrett’s natural history and may inform the potential for cancer beyond the clinically discernible phenotype.
Esophageal adenocarcinoma is a prominent example of cancer characterized by frequent amplifications in oncogenes. However, the mechanisms leading to amplicons that involve breakage-fusion-bridge cycles and extrachromosomal DNA are poorly understood. Here, we use 710 esophageal adenocarcinoma cases with matched samples and patient-derived organoids to disentangle complex amplicons and their associated mechanisms. Short-read sequencing identifies ERBB2, MYC, MDM2, and HMGA2 as the most frequent oncogenes amplified in extrachromosomal DNAs. We resolve complex extrachromosomal DNA and breakage-fusion-bridge cycles amplicons by integrating of de-novo assemblies and DNA methylation in nine long-read sequenced cases. Complex amplicons shared between precancerous biopsy and late-stage tumor, an enrichment of putative enhancer elements and mobile element insertions are potential drivers of complex amplicons’ origin. We find that patient-derived organoids recapitulate extrachromosomal DNA observed in the primary tumors and single-cell DNA sequencing capture extrachromosomal DNA-driven clonal dynamics across passages. Prospectively, long-read and single-cell DNA sequencing technologies can lead to better prediction of clonal evolution in esophageal adenocarcinoma.
Protein aggregation is the abnormal association of misfolded proteins into larger, often insoluble structures that can be toxic during ageing and in protein aggregation-associated diseases. Previous research has established a role for the cytosolic Tsa1 peroxiredoxin in responding to protein misfolding stress. Tsa1 is also known to downregulate the cAMP/protein kinase A (PKA) pathway as part of the response to hydrogen peroxide stress. However, whether the cAMP/PKA pathway is involved in protein misfolding stress is not known. Using transcriptomics, we examined the response to protein misfolding stress and found up-regulation of numerous stress gene functions and down regulation of many genes related to protein synthesis and other growth-related processes consistent with the well characterized environmental stress response. The scope of the transcriptional response is largely similar in wild-type and tsa1 mutant strains, but the magnitude is dampened in the strain lacking Tsa1. We identified a direct protein interaction between Tsa1 and the Bcy1 regulatory subunit of PKA that is present under normal growth conditions and explains the observed differences in gene expression profiles. This interaction is increased in a redox-dependent manner in response to nascent protein misfolding, via Tsa1-mediated oxidation of Bcy1. Oxidation of Bcy1 causes a reduction in cAMP binding by Bcy1 which dampens PKA pathway activity, leading to a targeted reprogramming of gene expression. Redox regulation of the regulatory subunit of PKA provides a mechanism to mitigate the toxic consequences of protein misfolding stress that is distinct to stress caused by exogenous sources of reactive oxygen species.