High-Throughput Mapping of Long-Range Neuronal Projection Using In Situ Sequencing.

Understanding neural circuits requires deciphering interactions among myriad cell types defined by spatial organization, connectivity, gene expression, and other properties. Resolving these cell types requires both single-neuron resolution and high throughput, a challenging combination with conventional methods. Here, we introduce barcoded anatomy resolved by sequencing (BARseq), a multiplexed method based on RNA barcoding for mapping projections of thousands of spatially resolved neurons in a single brain and relating those projections to other properties such as gene or Cre expression. Mapping the projections to 11 areas of 3,579 neurons in mouse auditory cortex using BARseq confirmed the laminar organization of the three top classes (intratelencephalic [IT], pyramidal tract-like [PT-like], and corticothalamic [CT]) of projection neurons. In depth analysis uncovered a projection type restricted almost exclusively to transcriptionally defined subtypes of IT neurons. By bridging anatomical and transcriptomic approaches at cellular resolution with high throughput, BARseq can potentially uncover the organizing principles underlying the structure and formation of neural circuits.

Clonal multi-omics reveals Bcor as a negative regulator of emergency dendritic cell development.

Despite advances in single-cell multi-omics, a single stem or progenitor cell can only be tested once. We developed clonal multi-omics, in which daughters of a clone act as surrogates of the founder, thereby allowing multiple independent assays per clone. With SIS-seq, clonal siblings in parallel "sister" assays are examined either for gene expression by RNA sequencing (RNA-seq) or for fate in culture. We identified, and then validated using CRISPR, genes that controlled fate bias for different dendritic cell (DC) subtypes. This included Bcor as a suppressor of plasmacytoid DC (pDC) and conventional DC type 2 (cDC2) numbers during Flt3 ligand-mediated emergency DC development. We then developed SIS-skew to examine development of wild-type and Bcor-deficient siblings of the same clone in parallel. We found Bcor restricted clonal expansion, especially for cDC2s, and suppressed clonal fate potential, especially for pDCs. Therefore, SIS-seq and SIS-skew can reveal the molecular and cellular mechanisms governing clonal fate.

Synthetic Lethal and Resistance Interactions with BET Bromodomain Inhibitors in Triple-Negative Breast Cancer.

BET bromodomain inhibitors (BBDIs) are candidate therapeutic agents for triple-negative breast cancer (TNBC) and other cancer types, but inherent and acquired resistance to BBDIs limits their potential clinical use. Using CRISPR and small-molecule inhibitor screens combined with comprehensive molecular profiling of BBDI response and resistance, we identified synthetic lethal interactions with BBDIs and genes that, when deleted, confer resistance. We observed synergy with regulators of cell cycle progression, YAP, AXL, and SRC signaling, and chemotherapeutic agents. We also uncovered functional similarities and differences among BRD2, BRD4, and BRD7. Although deletion of BRD2 enhances sensitivity to BBDIs, BRD7 loss leads to gain of TEAD-YAP chromatin binding and luminal features associated with BBDI resistance. Single-cell RNA-seq, ATAC-seq, and cellular barcoding analysis of BBDI responses in sensitive and resistant cell lines highlight significant heterogeneity among samples and demonstrate that BBDI resistance can be pre-existing or acquired.

Clonal tracking in cancer and metastasis.

The eradication of many cancers has proven challenging due to the presence of functionally and genetically heterogeneous clones maintained by rare cancer stem cells (CSCs), which contribute to disease progression, treatment refractoriness, and late relapse. The characterization of functional CSC activity has necessitated the development of modern clonal tracking strategies. This review describes viral-based and CRISPR-Cas9-based cellular barcoding, lineage tracing, and imaging-based approaches. DNA-based cellular barcoding technology is emerging as a powerful and robust strategy that has been widely applied to in vitro and in vivo model systems, including patient-derived xenograft models. This review also highlights the potential of these methods for use in the clinical and drug discovery contexts and discusses the important insights gained from such approaches.

Modeling glioblastoma heterogeneity as a dynamic network of cell states.

Tumor cell heterogeneity is a crucial characteristic of malignant brain tumors and underpins phenomena such as therapy resistance and tumor recurrence. Advances in single-cell analysis have enabled the delineation of distinct cellular states of brain tumor cells, but the time-dependent changes in such states remain poorly understood. Here, we construct quantitative models of the time-dependent transcriptional variation of patient-derived glioblastoma (GBM) cells. We build the models by sampling and profiling barcoded GBM cells and their progeny over the course of 3 weeks and by fitting a mathematical model to estimate changes in GBM cell states and their growth rates. Our model suggests a hierarchical yet plastic organization of GBM, where the rates and patterns of cell state switching are partly patient-specific. Therapeutic interventions produce complex dynamic effects, including inhibition of specific states and altered differentiation. Our method provides a general strategy to uncover time-dependent changes in cancer cells and offers a way to evaluate and predict how therapy affects cell state composition.

Toward defining a 'lineage'--The case for dendritic cells.

The immune system consists of a heterogeneous ensemble of cell types that immunologists have tried to classify and order for decades. This classification has relied on varying criteria, resulting in major debates in the immunology community. Discovered in the late 1970s [1], dendritic cells (DCs) are no exception, and their membership to a distinct immune lineage is still vividly debated [2-6]. Here, we review recent work on the origin of DCs and discuss the possible definition of a separate 'DC lineage'.

Reproducibility of Illumina platform deep sequencing errors allows accurate determination of DNA barcodes in cells.

BACKGROUND: Next generation sequencing (NGS) of amplified DNA is a powerful tool to describe genetic heterogeneity within cell populations that can both be used to investigate the clonal structure of cell populations and to perform genetic lineage tracing. For applications in which both abundant and rare sequences are biologically relevant, the relatively high error rate of NGS techniques complicates data analysis, as it is difficult to distinguish rare true sequences from spurious sequences that are generated by PCR or sequencing errors. This issue, for instance, applies to cellular barcoding strategies that aim to follow the amount and type of offspring of single cells, by supplying these with unique heritable DNA tags. RESULTS: Here, we use genetic barcoding data from the Illumina HiSeq platform to show that straightforward read threshold-based filtering of data is typically insufficient to filter out spurious barcodes. Importantly, we demonstrate that specific sequencing errors occur at an approximately constant rate across different samples that are sequenced in parallel. We exploit this observation by developing a novel approach to filter out spurious sequences. CONCLUSIONS: Application of our new method demonstrates its value in the identification of true sequences amongst spurious sequences in biological data sets.

Enhanced multiplexing in mass cytometry using osmium and ruthenium tetroxide species.

Mass cytometry facilitates high-dimensional, quantitative, single-cell analysis. The method for sample multiplexing in mass cytometry, called mass-tag cellular barcoding (MCB), relies on the covalent reaction of bifunctional metal chelators with intracellular proteins. Here, we describe the use of osmium and ruthenium tetroxides (OsO4 and RuO4 ) that bind covalently with fatty acids in the cellular membranes and aromatic amino acids in proteins. Both OsO4 and RuO4 rapidly reacted and allowed for MCB with live cells, crosslinked cells, and permeabilized cells. Given the covalent nature of the labeling reaction, isotope leaching was not observed. OsO4 and RuO4 were used in a 20-sample barcoding protocol together with palladium isotopes. As mass channels occupied by osmium and ruthenium are not used for antibody detection the number of masses effectively utilized in a single experiment is expanded. OsO4 and RuO4 can therefore be used as MCB reagents for a wide range of mass cytometry workflows. © 2016 International Society for Advancement of Cytometry.

Transient partial permeabilization with saponin enables cellular barcoding prior to surface marker staining.

Fluorescent cellular barcoding and mass-tag cellular barcoding are cytometric methods that enable high sample throughput, minimize inter-sample variation, and reduce reagent consumption. Previously employed barcoding protocols require that barcoding be performed after surface marker staining, complicating combining the technique with measurement of alcohol-sensitive surface epitopes. This report describes a method of barcoding fixed cells after a transient partial permeabilization with 0.02% saponin that results in efficient and consistent barcode staining with fluorescent or mass-tagged reagents while preserving surface marker staining. This approach simplifies barcoding protocols and allows direct comparison of surface marker staining of multiple samples without concern for variations in the antibody cocktail volume, antigen-antibody ratio, or machine sensitivity. Using this protocol, cellular barcoding can be used to reliably detect subtle differences in surface marker expression.

Retrospective identification of cell-intrinsic factors that mark pluripotency potential in rare somatic cells.

Pluripotency can be induced in somatic cells by the expression of OCT4, KLF4, SOX2, and MYC. Usually only a rare subset of cells reprogram, and the molecular characteristics of this subset remain unknown. We apply retrospective clone tracing to identify and characterize the rare human fibroblasts primed for reprogramming. These fibroblasts showed markers of increased cell cycle speed and decreased fibroblast activation. Knockdown of a fibroblast activation factor identified by our analysis increased the reprogramming efficiency. We provide evidence for a unified model in which cells can move into and out of the primed state over time, explaining how reprogramming appears deterministic at short timescales and stochastic at long timescales. Furthermore, inhibiting the activity of LSD1 enlarged the pool of cells that were primed for reprogramming. Thus, even homogeneous cell populations can exhibit heritable molecular variability that can dictate whether individual rare cells will reprogram or not.

Hybrid Macromolecular Constructs as a Platform for Spectral Nanoprobes for Advanced Cellular Barcoding in Flow Cytometry.

Herein, an advanced bioconjugation technique to synthesize hybrid polymer-antibody nanoprobes tailored for fluorescent cell barcoding in flow cytometry-based immunophenotyping of leukocytes is applied. A novel approach of attachment combining two fluorescent dyes on the copolymer precursor and its conjugation to antibody is employed to synthesize barcoded nanoprobes of antibody polymer dyes allowing up to six nanoprobes to be resolved in two-dimensional cytometry analysis. The major advantage of these nanoprobes is the construct design in which the selected antibody is labeled with an advanced copolymer bearing two types of fluorophores in different molar ratios. The cells after antibody recognition and binding to the target antigen have a characteristic double fluorescence signal for each nanoprobe providing a unique position on the dot plot, thus allowing antibody-based barcoding of cellular samples in flow cytometry assays. This technique is valuable for cellular assays that require low intersample variability and is demonstrated by the live cell barcoding of clinical samples with B cell abnormalities. In total, the samples from six various donors were successfully barcoded using only two detection channels. This barcoding of clinical samples enables sample preparation and measurement in a single tube.

Achieving single-cell-resolution lineage tracing in zebrafish by continuous barcoding mutations during embryogenesis.

Unraveling the lineage relationships of all descendants from a zygote is fundamental to advancing our understanding of developmental and stem cell biology. However, existing cell barcoding technologies in zebrafish lack the resolution to capture the majority of cell divisions during embryogenesis. A recently developed method, a substitution mutation-aided lineage-tracing system (SMALT), successfully reconstructed high-resolution cell phylogenetic trees for Drosophila melanogaster. Here, we implement the SMALT system in zebrafish, recording a median of 14 substitution mutations on a one-kilobase-pair barcoding sequence for one-day post-fertilization embryos. Leveraging this system, we reconstruct four cell lineage trees for zebrafish fin cells, encompassing both original and regenerated fin. Each tree consists of hundreds of internal nodes with a median bootstrap support of 99%. Analysis of the obtained cell lineage trees reveals that regenerated fin cells mainly originate from cells in the same part of the fins. Through multiple times sampling germ cells from the same individual, we show the stability of the germ cell pool and the early separation of germ cell and somatic cell progenitors. Our system offers the potential for reconstructing high-quality cell phylogenies across diverse tissues, providing valuable insights into development and disease in zebrafish.

Analysis of synthetic cellular barcodes in the genome and transcriptome with BARtab and bartools.

Cellular barcoding is a lineage-tracing methodology that couples heritable synthetic barcodes to high-throughput sequencing, enabling the accurate tracing of cell lineages across a range of biological contexts. Recent studies have extended these methods by incorporating lineage information into single-cell or spatial transcriptomics readouts. Leveraging the rich biological information within these datasets requires dedicated computational tools for dataset pre-processing and analysis. Here, we present BARtab, a portable and scalable Nextflow pipeline, and bartools, an open-source R package, designed to provide an integrated end-to-end cellular barcoding analysis toolkit. BARtab and bartools contain methods to simplify the extraction, quality control, analysis, and visualization of lineage barcodes from population-level, single-cell, and spatial transcriptomics experiments. We showcase the utility of our integrated BARtab and bartools workflow via the analysis of exemplar bulk, single-cell, and spatial transcriptomics experiments containing cellular barcoding information.

Identification of collateral sensitivity and evolutionary landscape of chemotherapy-induced drug resistance using cellular barcoding technology.

Background: One of the most significant challenges impeding cancer treatment effectiveness is drug resistance. Combining evolutionary understanding with drug resistance can pave the way for the identification of second-line drug options that can overcome drug resistance. Although capecitabine and irinotecan are commonly used therapeutic agents in the treatment of CRC patients, resistance to these agents is common. The underlying clonal dynamics of resistance to these agents using high-resolution barcode technology and identification of effective second-line drugs in this context remain unclear. Methods and materials: Caco-2 and HT-29 cell lines were barcoded, and then capecitabine and irinotecan resistant derivatives of these cell lines were established. The frequencies of barcodes from resistant cell lines and harvested medium, longitudinally, were determined. Collateral drug sensitivity testing was carried out on resistant Caco-2 and HT-29 cell lines using single agents or drug combinations. The SyngeryFinder tool was used to analyse drug combination testing. Results: In Caco-2 and HT-29 cell lines, barcode frequency measurements revealed clonal dynamics of capecitabine and irinotecan formed by both pre-existing and de novo barcodes, indicating the presence of polyclonal drug resistance. The temporal dynamics of clonal evolution in Caco-2 and HT-29 cell lines were demonstrated by longitudinal analysis of pre-existing and de novo barcodes from harvested medium. In Caco-2 and HT-29 cell lines, collateral drug sensitivity revealed a number of drugs that were effective alone and in combination. Conclusion: The use of barcoding technology reveals the clonal dynamics of chemotherapy-induced drug resistance not only from harvested cell populations, but also from longitudinal sampling throughout the course of clonal evolution. Second-line drugs that sensitize drug-resistant CRC cell lines are identified through collateral drug testing.

How colonization bottlenecks, tissue niches, and transmission strategies shape protozoan infections.

Protozoan pathogens such as Plasmodium spp., Leishmania spp., Toxoplasma gondii, and Trypanosoma spp. are often associated with high-mortality, acute and chronic diseases of global health concern. For transmission and immune evasion, protozoans have evolved diverse strategies to interact with a range of host tissue environments. These interactions are linked to disease pathology, yet our understanding of the association between parasite colonization and host homeostatic disruption is limited. Recently developed techniques for cellular barcoding have the potential to uncover the biology regulating parasite transmission, dissemination, and the stability of infection. Understanding bottlenecks to infection and the in vivo tissue niches that facilitate chronic infection and spread has the potential to reveal new aspects of parasite biology.

Mechanisms of Acquired Resistance and Tolerance to EGFR Targeted Therapy in Non-Small Cell Lung Cancer.

Non-small cell lung cancers (NSCLC) harboring activating mutations of the epidermal growth factor receptor (EGFR) are treated with specific tyrosine kinase inhibitors (EGFR-TKIs) of this receptor, resulting in clinically responses that can generally last several months. Unfortunately, EGFR-targeted therapy also favors the emergence of drug tolerant or resistant cells, ultimately resulting in tumor relapse. Recently, cellular barcoding strategies have arisen as a powerful tool to investigate the clonal evolution of these subpopulations in response to anti-cancer drugs. In this review, we provide an overview of the currently available treatment options for NSCLC, focusing on EGFR targeted therapy, and discuss the common mechanisms of resistance to EGFR-TKIs. We also review the characteristics of drug-tolerant persister (DTP) cells and the mechanistic basis of drug tolerance in EGFR-mutant NSCLC. Lastly, we address how cellular barcoding can be applied to investigate the response and the behavior of DTP cells upon EGFR-TKI treatment.

Characterization of the minimal residual disease state reveals distinct evolutionary trajectories of human glioblastoma.

Recurrence of solid tumors renders patients vulnerable to advanced, treatment-refractory disease state with mutational and oncogenic landscape distinctive from initial diagnosis. Improving outcomes for recurrent cancers requires a better understanding of cell populations that expand from the post-therapy, minimal residual disease (MRD) state. We profile barcoded tumor stem cell populations through therapy at tumor initiation, MRD, and recurrence in our therapy-adapted, patient-derived xenograft models of glioblastoma (GBM). Tumors show distinct patterns of recurrence in which clonal populations exhibit either a pre-existing fitness advantage or an equipotency fitness acquired through therapy. Characterization of the MRD state by single-cell and bulk RNA sequencing reveals a tumor-intrinsic immunomodulatory signature with prognostic significance at the transcriptomic level and in proteomic analysis of cerebrospinal fluid (CSF) collected from patients with GBM. Our results provide insight into the innate and therapy-driven dynamics of human GBM and the prognostic value of interrogating the MRD state in solid cancers.

Erythropoietin directly remodels the clonal composition of murine hematopoietic multipotent progenitor cells.

The cytokine erythropoietin (EPO) is a potent inducer of erythrocyte development and one of the most prescribed biopharmaceuticals. The action of EPO on erythroid progenitor cells is well established, but its direct action on hematopoietic stem and progenitor cells (HSPCs) is still debated. Here, using cellular barcoding, we traced the differentiation of hundreds of single murine HSPCs, after ex vivo EPO exposure and transplantation, in five different hematopoietic cell lineages, and observed the transient occurrence of high-output myeloid-erythroid-megakaryocyte-biased and myeloid-B-cell-dendritic cell-biased clones. Single-cell RNA sequencing analysis of ex vivo EPO-exposed HSPCs revealed that EPO induced the upregulation of erythroid associated genes in a subset of HSPCs, overlapping with multipotent progenitor (MPP) 1 and MPP2. Transplantation of barcoded EPO-exposed MPP2 confirmed their enrichment in myeloid-erythroid-biased clones. Collectively, our data show that EPO does act directly on MPP independent of the niche and modulates fate by remodeling the clonal composition of the MPP pool.

Cellular barcoding of protozoan pathogens reveals the within-host population dynamics of Toxoplasma gondii host colonization.

Cellular barcoding techniques are powerful tools to understand microbial pathogenesis. However, barcoding strategies have not been broadly applied to protozoan parasites, which have unique genomic structures and virulence strategies compared with viral and bacterial pathogens. Here, we present a CRISPR-based method to barcode protozoa, which we successfully apply to Toxoplasma gondii and Trypanosoma brucei. Using libraries of barcoded T. gondii, we evaluate shifts in the population structure from acute to chronic infection of mice. Contrary to expectation, most barcodes were present in the brain one month post-intraperitoneal infection in both inbred CBA/J and outbred Swiss mice. Although parasite cyst number and barcode diversity declined over time, barcodes representing a minor fraction of the inoculum could become a dominant population in the brain by three months post-infection. These data establish a cellular barcoding approach for protozoa and evidence that the blood-brain barrier is not a major bottleneck to colonization by T. gondii.

In Vivo Tracking of Hematopoietic Stem and Progenitor Cell Ontogeny by Cellular Barcoding.

Cellular barcoding is a powerful technique that allows for high-throughput mapping of the fate of single cells, notably hematopoietic stem and progenitor cells (HSPCs) after transplantation. Unique artificial DNA fragments, termed barcodes, are stably inserted into HSPCs using lentiviral transduction, making sure that each individual cell receives a single unique barcode. Barcoded HSPCs are transplanted into sublethally irradiated mice where they reconstitute the hematopoietic system through proliferation and differentiation. During this process, the barcode of each HSPC is inherited by all of its daughter cells and their subsequent mature hematopoietic cell progeny. After sorting mature hematopoietic cell subsets, their barcodes can be retrieved from genomic DNA through nested PCR and sequencing. Analysis of barcode sequencing results allows for determination of clonal relationships between the mature cells, that is, which cell types were produced by a single barcoded HSPC, as well as the heterogeneity of the initial HSPC population. Here, we give a detailed protocol of a complete HSPC cellular barcoding experiment, starting with barcode lentivirus production, isolation, transduction, and transplantation of HSPCs, isolation of target cells followed by PCR amplification and sequencing of DNA barcodes. Finally, we describe the basic filtering and analysis steps of barcode sequencing data to ensure high-quality results.