Search

We present Light-Seq, an approach for multiplexed spatial indexing of intact biological samples using light-directed DNA barcoding in fixed cells and tissues followed by ex situ sequencing. Light-Seq combines spatially targeted, rapid photocrosslinking of DNA barcodes onto complementary DNAs in situ with a one-step DNA stitching reaction to create pooled, spatially indexed sequencing libraries. This light-directed barcoding enables in situ selection of multiple cell populations in intact fixed tissue samples for full-transcriptome sequencing based on location, morphology or protein stains, without cellular dissociation. Applying Light-Seq to mouse retinal sections, we recovered thousands of differentially enriched transcripts from three cellular layers and discovered biomarkers for a very rare neuronal subtype, dopaminergic amacrine cells, from only four to eight individual cells per section. Light-Seq provides an accessible workflow to combine in situ imaging and protein staining with next generation sequencing of the same cells, leaving the sample intact for further analysis post-sequencing.

Orthogonal sequence library design is an essential task in bioengineering. Typical design approaches scale quadratically in the size of the candidate sequence space. As such, exhaustive searches of sequence space to maximize library size are computationally intractable with existing methods. Here, we present SeqWalk, a time and memory efficient method for designing maximally-sized orthogonal sequence libraries using the sequence symmetry minimization heuristic. SeqWalk encodes sequence design constraints in a de Bruijn graph representation of sequence space, enabling the application of efficient graph traversal techniques to the problem of orthogonal DNA sequence design. We demonstrate the scalability of SeqWalk by designing a provably maximal set of > 106 orthogonal 25nt sequences in less than 20 seconds on a single standard CPU core. We additionally derive fundamental bounds on orthogonal sequence library size under a variety of design constraints.

Circulating extracellular vesicles (EVs)-biological nanomaterials shed from most mammalian cells-have emerged as promising biomarkers, drug delivery vesicles, and treatment modulators. While different types of vesicles are being explored for these applications, it is becoming clear that human EVs are quite heterogeneous even in homogeneous or monoclonal cell populations. Since it is the surface EV protein composition that will largely dictate their biological behavior, high-throughput single EV profiling methods are needed to better define EV subpopulations. Here, we present an antibody-based immunosequencing method that allows multiplexed measurement of protein molecules from individual nanometer-sized EVs. We use droplet microfluidics to compartmentalize and barcode individual EVs. The barcodes/antibody-DNA are then sequenced to determine protein composition. Using this highly sensitive technology, we detected specific proteins at the single EV level. We expect that this technology can be further adapted for multiplexed protein analysis of any nanoparticle.

Cancer-induced immune responses affect tumor progression and therapeutic response. In multiple murine models and clinical datasets, we identified large variations of neutrophils and macrophages, which define “immune subtypes” of triple negative breast cancer (TNBC) including neutrophil-enriched (NES) and macrophage-enriched subtypes (MES). Different tumor-intrinsic pathways and mutual regulation between macrophages/monocytes and neutrophils contribute to the development of dichotomous myeloid compartment. MES contains predominantly macrophages that are CCR2-dependent and exhibit variable responses to immune checkpoint blockade (ICB). NES exhibits systemic and local accumulation of immunosuppressive neutrophils (or granulocytic myeloid-derived suppressor cells (gMDSCs), is resistant to ICB, and contains a minority of macrophages that appear to be unaffected by CCR2 knockout. A MES-to-NES conversion mediated acquired ICB resistance of initially sensitive MES models. Our results demonstrate diverse myeloid cell frequencies, functionality, and potential roles in immunotherapies, and highlight the need to better understand the inter-patient heterogeneity of the myeloid compartment.

We report the single-strand Recombinase Polymerase Amplification (ssRPA) method, which merges the fast, isothermal amplification of RPA with subsequent rapid conversion of the double-strand DNA amplicon to single strands, and hence enables facile hybridization-based, high-specificity readout. We demonstrate the utility of ssRPA for sensitive and rapid (4 copies per 50 µL reaction within 10 min, or 8 copies within 8 min) visual detection of SARS-CoV-2 RNA spiked samples, as well as clinical saliva and nasopharyngeal swabs in VTM or water, on lateral flow devices. The ssRPA method promises rapid, sensitive, and accessible RNA detection to facilitate mass testing in the COVID-19 pandemic.

Droplet‐based single cell sequencing technologies, such as inDrop, Drop‐seq, and 10X Genomics, are catalyzing a revolution in the understanding of biology. Barcoding beads are key components for these technologies. What is limiting today are barcoding beads that are easy to fabricate, can efficiently deliver primers into drops, and thus achieve high detection efficiency. Here, this work reports an approach to fabricate dissolvable polyacrylamide beads, by crosslinking acrylamide with disulfide bridges that can be cleaved with dithiothreitol. The beads can be rapidly dissolved in drops and release DNA barcode primers. The dissolvable beads are easy to synthesize, and the primer cost for the beads is significantly lower than that for the previous barcoding beads. Furthermore, the dissolvable beads can be loaded into drops with >95% loading efficiency of a single bead per drop and the dissolution of beads does not influence reverse transcription or the polymerase chain reaction (PCR) in drops. Based on this approach, the dissolvable beads are used for single cell RNA and protein analysis., This work reports an approach to fabricate dissolvable polyacrylamide beads, by crosslinking acrylamide with disulfide bridges that can be cleaved with dithiothreitol. The beads can be rapidly dissolved in drops and release DNA barcode primers for droplet DNA barcoding, which demonstrate great potential for high‐throughput single cell sequencing.

ARID1A (AT-rich interacting domain-containing protein 1A) is one of the frequently mutated epigenetic regulators in a wide spectrum of cancers. Recent studies have shown that ARID1A deficiency induces global changes in the epigenetic landscape of enhancers and gene regulations. These complex changes make it challenging to identify the driving mechanisms of ARID1A deficiency in promoting cancer progression. Using a mouse model, we identified the anti-senescence effect of ARID1A deficiency, which promotes the progression of pancreatic intraepithelial neoplasia (PanIN). To identify the underlying mechanism and the players, we profiled the transcriptomes of individual PanIN lesions. As a result, we observed not only the anti-senescence effect but also the complex trans-differentiation of neoplastic cells. Furthermore, we discovered that significantly upregulated aldehyde dehydrogenase (ALDH) help to inhibit Kras-induced reactive oxygen species (ROS) production and lead to the reduction of the cellular senescence, which results in the accelerated PanIN progression.

Cancer cells function as primary architects of the tumor microenvironment. Yet, the molecular features of cancer cells that govern stromal cell phenotypes remain unclear. Here, we show that cancer-associated fibroblasts (CAFs) are distinguishable on the basis of gene expression signatures they acquire in co-culture with epithelial- or mesenchymal-like lung adenocarcinoma (LUAD) cells. High expression of the EMT activator ZEB1 endows LUAD cells with the capacity to activate a soluble factor exchange that leads to CAF reprogramming, to generate CAF-led invasive projections in multicellular aggregates, and to respond to pro-metastatic signals from CAFs in mice. Thus, ZEB1-expressing LUAD cells are positioned at the apex of a signaling hierarchy in the tumor microenvironment.

Enterococci are a part of human microbiota and a leading cause of multidrug resistant infections. Here, we identify a family of Enterococcus pore-forming toxins (Epxs) in E. faecalis, E. faecium, and E. hirae strains isolated across the globe. Structural studies reveal that Epxs form a branch of β-barrel pore-forming toxins with a β-barrel protrusion (designated the top domain) sitting atop the cap domain. Through a genome-wide CRISPR-Cas9 screen, we identify human leukocyte antigen class I (HLA-I) complex as a receptor for two members (Epx2 and Epx3), which preferentially recognize human HLA-I and homologous MHC-I of equine, bovine, and porcine, but not murine, origin. Interferon exposure, which stimulates MHC-I expression, sensitizes human cells and intestinal organoids to Epx2 and Epx3 toxicity. Co-culture with Epx2-harboring E. faecium damages human peripheral blood mononuclear cells and intestinal organoids, and this toxicity is neutralized by an Epx2 antibody, demonstrating the toxin-mediated virulence of Epx-carrying Enterococcus.

In the studies of single-cell genomics, the large endeavor has been focused on the detection of the permanent changes in the genome. On the other hand, spontaneous DNA damage frequently occurs and results in transient single-stranded changes to the genome until they are repaired. So far, successful profiling of these dynamic changes has not been demonstrated by single-cell whole-genome amplification methods. Here we reported a novel single-cell WGA method: Linearly Produced Semiamplicon based Split Amplification Reaction (LPSSAR), which allows, for the first time, the genome-wide detection of the DNA damage associated single nucleotide variants (dSNVs) in single human cells. The sequence-based detection of dSNVs allows the direct characterization of the major damage signature that occurred in human cells. In the analysis of the abundance of dSNVs along the genome, we observed two modules of dSNV abundance, instead of a homogeneous abundance of dSNVs. Interestingly, we found that the two modules are associated with the A/B topological compartments of the genome. This result suggests that the genome topology directly influences genome stability. Furthermore, with the detection of a large number of dSNVs in single cells, we showed that only under a stringent filtering condition, can we distinguish the de novo mutations from the dSNVs and achieve a reliable estimation of the total level of de novo mutations in a single cell.

Single-cell technologies have emerged as advanced tools to study various biological processes that demand the single cell resolution. To detect subtle heterogeneity in the transcriptome, high accuracy and sensitivity are still desired for single-cell RNA-seq. We describe here multiple annealing and dC-tailing-based quantitative single-cell RNA-seq (MATQ-seq) with ~90% capture efficiency. In addition, MATQ-seq is a total RNA assay allowing for detection of nonpolyadenylated transcripts.

The mammalian cochlea loses its ability to regenerate new hair cells prior to the onset of hearing. In contrast, the adult vestibular system can produce new hair cells in response to damage, or by reprogramming of supporting cells with the hair cell transcription factor Atoh1. We used RNA-seq and ATAC-seq to probe the transcriptional and epigenetic responses of utricle supporting cells to damage and Atoh1 transduction. We show that the regenerative response of the utricle correlates with a more accessible chromatin structure in utricle supporting cells compared to their cochlear counterparts. We also provide evidence that Atoh1 transduction of supporting cells is able to promote increased transcriptional accessibility of some hair cell genes. Our study offers a possible explanation for regenerative differences between sensory organs of the inner ear, but shows that additional factors to Atoh1 may be required for optimal reprogramming of hair cell fate.

Single-cell technologies have emerged as advanced tools to study various biological processes that demand the single cell resolution. To detect subtle heterogeneity in the transcriptome, high accuracy and sensitivity are still desired for single-cell RNA-seq. We describe here multiple annealing and dC-tailing-based quantitative single-cell RNA-seq (MATQ-seq) with ~90% capture efficiency. In addition, MATQ-seq is a total RNA assay allowing for detection of nonpolyadenylated transcripts.

SUMMARY Cancer cells function as primary architects of the tumor microenvironment. However, the molecular features of cancer cells that govern stromal cell phenotypes remain unclear. Here, we show that cancer-associated fibroblast (CAF) heterogeneity is driven by lung adenocarcinoma (LUAD) cells at either end of the epithelial-to-mesenchymal transition (EMT) spectrum. LUAD cells that have high expression of the EMT-activating transcription factor ZEB1 reprogram CAFs through a ZEB1-dependent secretory program and direct CAFs to the tips of invasive projections through a ZEB1-driven CAF repulsion process. The EMT, in turn, sensitizes LUAD cells to pro-metastatic signals from CAFs. Thus, CAFs respond to contextual cues from LUAD cells to promote metastasis., In brief Bota-Rabassedas et al. show that EMT in lung adenocarcinoma cells activates a secretory process that governs CAF heterogeneity and, in turn, sensitizes lung adenocarcinoma cells to pro-metastatic signals from CAFs. Thus, EMT positions lung adenocarcinoma cells at the apex of a signaling hierarchy in the tumor microenvironment., Graphical Abstract

Decoding the molecular composition of individual Ngn3 + endocrine progenitors (EPs) during pancreatic morphogenesis could provide insight into the mechanisms regulating hormonal cell fate. Here, we identify population markers and extensive cellular diversity including four EP subtypes reflecting EP maturation using high-resolution single-cell RNA-sequencing of the e14.5 and e16.5 mouse pancreas. While e14.5 and e16.5 EPs are constantly born and share select genes, these EPs are overall transcriptionally distinct concomitant with changes in the underlying epithelium. As a consequence, e16.5 EPs are not the same as e14.5 EPs: e16.5 EPs have a higher propensity to form beta cells. Analysis of e14.5 and e16.5 EP chromatin states reveals temporal shifts, with enrichment of beta cell motifs in accessible regions at later stages. Finally, we provide transcriptional maps outlining the route progenitors take as they make cell fate decisions, which can be applied to advance the in vitro generation of beta cells., Endocrine progenitors form early in pancreatic development but the diversity of this cell population is unclear. Here, the authors use single cell RNA sequencing of the mouse pancreas at e14.5 and e16.5 to show that endocrine progenitors are temporally distinct and those formed later are more likely to become beta cells

The quantification of transcriptional variation in single cells, particularly within the same cell population, is currently limited by the low sensitivity and high technical noise of single-cell RNA-seq assays. We report multiple annealing and dC-tailing-based quantitative single-cell RNA-seq (MATQ-seq), a highly sensitive and quantitative method for single-cell sequencing of total RNA. By systematically determining technical noise, we show that MATQ-seq captures genuine biological variation between whole transcriptomes of single cells.

Blockade of angiogenesis can retard tumour growth, but may also paradoxically increase metastasis. This paradox may be resolved by vessel normalization, which involves increased pericyte coverage, improved tumour vessel perfusion, reduced vascular permeability, and consequently mitigated hypoxia. Although these processes alter tumour progression, their regulation is poorly understood. Here we show that type 1 T helper (TH1) cells play a crucial role in vessel normalization. Bioinformatic analyses revealed that gene expression features related to vessel normalization correlate with immunostimulatory pathways, especially T lymphocyte infiltration or activity. To delineate the causal relationship, we used various mouse models with vessel normalization or T lymphocyte deficiencies. Although disruption of vessel normalization reduced T lymphocyte infiltration as expected, reciprocal depletion or inactivation of CD4+ T lymphocytes decreased vessel normalization, indicating a mutually regulatory loop. In addition, activation of CD4+ T lymphocytes by immune checkpoint blockade increased vessel normalization. TH1 cells that secrete interferon-γ are a major population of cells associated with vessel normalization. Patient-derived xenograft tumours growing in immunodeficient mice exhibited enhanced hypoxia compared to the original tumours in immunocompetent humans, and hypoxia was reduced by adoptive TH1 transfer. Our findings elucidate an unexpected role of TH1 cells in vasculature and immune reprogramming. TH1 cells may be a marker and a determinant of both immune checkpoint blockade and anti-angiogenesis efficacy.

SummaryHematopoietic stem cells (HSCs) possess unique gene expression programs that enforce their identity and regulate lineage commitment. Long non-coding RNAs (lncRNAs) have emerged as important regulators of gene expression and cell fate decisions, although their functions in HSCs are unclear. Here we profiled the transcriptome of purified HSCs by deep sequencing and identified 323 unannotated lncRNAs. Comparing their expression in differentiated lineages revealed 159 lncRNAs enriched in HSCs, some of which are likely HSC specific (LncHSCs). These lncRNA genes share epigenetic features with protein-coding genes, including regulated expression via DNA methylation, and knocking down two LncHSCs revealed distinct effects on HSC self-renewal and lineage commitment. We mapped the genomic binding sites of one of these candidates and found enrichment for key hematopoietic transcription factor binding sites, especially E2A. Together, these results demonstrate that lncRNAs play important roles in regulating HSCs, providing an additional layer to the genetic circuitry controlling HSC function.