High-Spatial-Resolution Multi-Omics Sequencing via Deterministic Barcoding in Tissue.

We present deterministic barcoding in tissue for spatial omics sequencing (DBiT-seq) for co-mapping of mRNAs and proteins in a formaldehyde-fixed tissue slide via next-generation sequencing (NGS). Parallel microfluidic channels were used to deliver DNA barcodes to the surface of a tissue slide, and crossflow of two sets of barcodes, A1-50 and B1-50, followed by ligation in situ, yielded a 2D mosaic of tissue pixels, each containing a unique full barcode AB. Application to mouse embryos revealed major tissue types in early organogenesis as well as fine features like microvasculature in a brain and pigmented epithelium in an eye field. Gene expression profiles in 10-µm pixels conformed into the clusters of single-cell transcriptomes, allowing for rapid identification of cell types and spatial distributions. DBiT-seq can be adopted by researchers with no experience in microfluidics and may find applications in a range of fields including developmental biology, cancer biology, neuroscience, and clinical pathology.

Massively Parallel Interrogation of the Effects of Gene Expression Levels on Fitness.

Data of gene expression levels across individuals, cell types, and disease states is expanding, yet our understanding of how expression levels impact phenotype is limited. Here, we present a massively parallel system for assaying the effect of gene expression levels on fitness in Saccharomyces cerevisiae by systematically altering the expression level of ∼100 genes at ∼100 distinct levels spanning a 500-fold range at high resolution. We show that the relationship between expression levels and growth is gene and environment specific and provides information on the function, stoichiometry, and interactions of genes. Wild-type expression levels in some conditions are not optimal for growth, and genes whose fitness is greatly affected by small changes in expression level tend to exhibit lower cell-to-cell variability in expression. Our study addresses a fundamental gap in understanding the functional significance of gene expression regulation and offers a framework for evaluating the phenotypic effects of expression variation.

Slide-tags enables single-nucleus barcoding for multimodal spatial genomics.

Recent technological innovations have enabled the high-throughput quantification of gene expression and epigenetic regulation within individual cells, transforming our understanding of how complex tissues are constructed1-6. However, missing from these measurements is the ability to routinely and easily spatially localize these profiled cells. We developed a strategy, Slide-tags, in which single nuclei within an intact tissue section are tagged with spatial barcode oligonucleotides derived from DNA-barcoded beads with known positions. These tagged nuclei can then be used as an input into a wide variety of single-nucleus profiling assays. Application of Slide-tags to the mouse hippocampus positioned nuclei at less than 10 µm spatial resolution and delivered whole-transcriptome data that are indistinguishable in quality from ordinary single-nucleus RNA-sequencing data. To demonstrate that Slide-tags can be applied to a wide variety of human tissues, we performed the assay on brain, tonsil and melanoma. We revealed cell-type-specific spatially varying gene expression across cortical layers and spatially contextualized receptor-ligand interactions driving B cell maturation in lymphoid tissue. A major benefit of Slide-tags is that it is easily adaptable to almost any single-cell measurement technology. As a proof of principle, we performed multiomic measurements of open chromatin, RNA and T cell receptor (TCR) sequences in the same cells from metastatic melanoma, identifying transcription factor motifs driving cancer cell state transitions in spatially distinct microenvironments. Slide-tags offers a universal platform for importing the compendium of established single-cell measurements into the spatial genomics repertoire.

The development of terrestrial ecosystems emerging after glacier retreat.

The global retreat of glaciers is dramatically altering mountain and high-latitude landscapes, with new ecosystems developing from apparently barren substrates1-4. The study of these emerging ecosystems is critical to understanding how climate change interacts with microhabitat and biotic communities and determines the future of ice-free terrains1,5. Here, using a comprehensive characterization of ecosystems (soil properties, microclimate, productivity and biodiversity by environmental DNA metabarcoding6) across 46 proglacial landscapes worldwide, we found that all the environmental properties change with time since glaciers retreated, and that temperature modulates the accumulation of soil nutrients. The richness of bacteria, fungi, plants and animals increases with time since deglaciation, but their temporal patterns differ. Microorganisms colonized most rapidly in the first decades after glacier retreat, whereas most macroorganisms took longer. Increased habitat suitability, growing complexity of biotic interactions and temporal colonization all contribute to the increase in biodiversity over time. These processes also modify community composition for all the groups of organisms. Plant communities show positive links with all other biodiversity components and have a key role in ecosystem development. These unifying patterns provide new insights into the early dynamics of deglaciated terrains and highlight the need for integrated surveillance of their multiple environmental properties5.

Machine-learning-optimized Cas12a barcoding enables the recovery of single-cell lineages and transcriptional profiles.

The development of CRISPR-based barcoding methods creates an exciting opportunity to understand cellular phylogenies. We present a compact, tunable, high-capacity Cas12a barcoding system called dual acting inverted site array (DAISY). We combined high-throughput screening and machine learning to predict and optimize the 60-bp DAISY barcode sequences. After optimization, top-performing barcodes had ∼10-fold increased capacity relative to the best random-screened designs and performed reliably across diverse cell types. DAISY barcode arrays generated ∼12 bits of entropy and ∼66,000 unique barcodes. Thus, DAISY barcodes-at a fraction of the size of Cas9 barcodes-achieved high-capacity barcoding. We coupled DAISY barcoding with single-cell RNA-seq to recover lineages and gene expression profiles from ∼47,000 human melanoma cells. A single DAISY barcode recovered up to ∼700 lineages from one parental cell. This analysis revealed heritable single-cell gene expression and potential epigenetic modulation of memory gene transcription. Overall, Cas12a DAISY barcoding is an efficient tool for investigating cell-state dynamics.

Diverse clonal fates emerge upon drug treatment of homogeneous cancer cells.

Even among genetically identical cancer cells, resistance to therapy frequently emerges from a small subset of those cells1-7. Molecular differences in rare individual cells in the initial population enable certain cells to become resistant to therapy7-9; however, comparatively little is known about the variability in the resistance outcomes. Here we develop and apply FateMap, a framework that combines DNA barcoding with single-cell RNA sequencing, to reveal the fates of hundreds of thousands of clones exposed to anti-cancer therapies. We show that resistant clones emerging from single-cell-derived cancer cells adopt molecularly, morphologically and functionally distinct resistant types. These resistant types are largely predetermined by molecular differences between cells before drug addition and not by extrinsic factors. Changes in the dose and type of drug can switch the resistant type of an initial cell, resulting in the generation and elimination of certain resistant types. Samples from patients show evidence for the existence of these resistant types in a clinical context. We observed diversity in resistant types across several single-cell-derived cancer cell lines and cell types treated with a variety of drugs. The diversity of resistant types as a result of the variability in intrinsic cell states may be a generic feature of responses to external cues.

Ångström-resolution fluorescence microscopy.

Fluorescence microscopy, with its molecular specificity, is one of the major characterization methods used in the life sciences to understand complex biological systems. Super-resolution approaches1-6 can achieve resolution in cells in the range of 15 to 20 nm, but interactions between individual biomolecules occur at length scales below 10 nm and characterization of intramolecular structure requires Ångström resolution. State-of-the-art super-resolution implementations7-14 have demonstrated spatial resolutions down to 5 nm and localization precisions of 1 nm under certain in vitro conditions. However, such resolutions do not directly translate to experiments in cells, and Ångström resolution has not been demonstrated to date. Here we introdue a DNA-barcoding method, resolution enhancement by sequential imaging (RESI), that improves the resolution of fluorescence microscopy down to the Ångström scale using off-the-shelf fluorescence microscopy hardware and reagents. By sequentially imaging sparse target subsets at moderate spatial resolutions of >15 nm, we demonstrate that single-protein resolution can be achieved for biomolecules in whole intact cells. Furthermore, we experimentally resolve the DNA backbone distance of single bases in DNA origami with Ångström resolution. We use our method in a proof-of-principle demonstration to map the molecular arrangement of the immunotherapy target CD20 in situ in untreated and drug-treated cells, which opens possibilities for assessing the molecular mechanisms of targeted immunotherapy. These observations demonstrate that, by enabling intramolecular imaging under ambient conditions in whole intact cells, RESI closes the gap between super-resolution microscopy and structural biology studies and thus delivers information key to understanding complex biological systems.

Cycling cancer persister cells arise from lineages with distinct programs.

Non-genetic mechanisms have recently emerged as important drivers of cancer therapy failure1, where some cancer cells can enter a reversible drug-tolerant persister state in response to treatment2. Although most cancer persisters remain arrested in the presence of the drug, a rare subset can re-enter the cell cycle under constitutive drug treatment. Little is known about the non-genetic mechanisms that enable cancer persisters to maintain proliferative capacity in the presence of drugs. To study this rare, transiently resistant, proliferative persister population, we developed Watermelon, a high-complexity expressed barcode lentiviral library for simultaneous tracing of each cell's clonal origin and proliferative and transcriptional states. Here we show that cycling and non-cycling persisters arise from different cell lineages with distinct transcriptional and metabolic programs. Upregulation of antioxidant gene programs and a metabolic shift to fatty acid oxidation are associated with persister proliferative capacity across multiple cancer types. Impeding oxidative stress or metabolic reprogramming alters the fraction of cycling persisters. In human tumours, programs associated with cycling persisters are induced in minimal residual disease in response to multiple targeted therapies. The Watermelon system enabled the identification of rare persister lineages that are preferentially poised to proliferate under drug pressure, thus exposing new vulnerabilities that can be targeted to delay or even prevent disease recurrence.

Comparative Analysis of Droplet-Based Ultra-High-Throughput Single-Cell RNA-Seq Systems.

Since its establishment in 2009, single-cell RNA sequencing (RNA-seq) has been a major driver behind progress in biomedical research. In developmental biology and stem cell studies, the ability to profile single cells confers particular benefits. Although most studies still focus on individual tissues or organs, the recent development of ultra-high-throughput single-cell RNA-seq has demonstrated potential power in characterizing more complex systems or even the entire body. However, although multiple ultra-high-throughput single-cell RNA-seq systems have attracted attention, no systematic comparison of these systems has been performed. Here, with the same cell line and bioinformatics pipeline, we developed directly comparable datasets for each of three widely used droplet-based ultra-high-throughput single-cell RNA-seq systems, inDrop, Drop-seq, and 10X Genomics Chromium. Although each system is capable of profiling single-cell transcriptomes, their detailed comparison revealed the distinguishing features and suitable applications for each system.

Combining environmental DNA and visual surveys can inform conservation planning for coral reefs.

Environmental DNA (eDNA) metabarcoding has the potential to revolutionize conservation planning by providing spatially and taxonomically comprehensive data on biodiversity and ecosystem conditions, but its utility to inform the design of protected areas remains untested. Here, we quantify whether and how identifying conservation priority areas within coral reef ecosystems differs when biodiversity information is collected via eDNA analyses or traditional visual census records. We focus on 147 coral reefs in Indonesia's hyper-diverse Wallacea region and show large discrepancies in the allocation and spatial design of conservation priority areas when coral reef species were surveyed with underwater visual techniques (fishes, corals, and algae) or eDNA metabarcoding (eukaryotes and metazoans). Specifically, incidental protection occurred for 55% of eDNA species when targets were set for species detected by visual surveys and 71% vice versa. This finding is supported by generally low overlap in detection between visual census and eDNA methods at species level, with more overlap at higher taxonomic ranks. Incomplete taxonomic reference databases for the highly diverse Wallacea reefs, and the complementary detection of species by the two methods, underscore the current need to combine different biodiversity data sources to maximize species representation in conservation planning.

Future of DNA-based insect monitoring.

Insects are crucial for ecosystem health but climate change and pesticide use are driving massive insect decline. To mitigate this loss, we need new and effective monitoring techniques. Over the past decade there has been a shift to DNA-based techniques. We describe key emerging techniques for sample collection. We suggest that the selection of tools should be broadened, and that DNA-based insect monitoring data need to be integrated more rapidly into policymaking. We argue that there are four key areas for advancement, including the generation of more complete DNA barcode databases to interpret molecular data, standardisation of molecular methods, scaling up of monitoring efforts, and integrating molecular tools with other technologies that allow continuous, passive monitoring based on images and/or laser imaging, detection, and ranging (LIDAR).

Arrayed in vivo barcoding for multiplexed sequence verification of plasmid DNA and demultiplexing of pooled libraries.

Sequence verification of plasmid DNA is critical for many cloning and molecular biology workflows. To leverage high-throughput sequencing, several methods have been developed that add a unique DNA barcode to individual samples prior to pooling and sequencing. However, these methods require an individual plasmid extraction and/or in vitro barcoding reaction for each sample processed, limiting throughput and adding cost. Here, we develop an arrayed in vivo plasmid barcoding platform that enables pooled plasmid extraction and library preparation for Oxford Nanopore sequencing. This method has a high accuracy and recovery rate, and greatly increases throughput and reduces cost relative to other plasmid barcoding methods or Sanger sequencing. We use in vivo barcoding to sequence verify >45 000 plasmids and show that the method can be used to transform error-containing dispersed plasmid pools into sequence-perfect arrays or well-balanced pools. In vivo barcoding does not require any specialized equipment beyond a low-overhead Oxford Nanopore sequencer, enabling most labs to flexibly process hundreds to thousands of plasmids in parallel.

TTSBBC: triplex target site biomarkers and barcodes in cancer.

The technology of triplex-forming oligonucleotides (TFOs) provides an approach to manipulate genes at the DNA level. TFOs bind to specific sites on genomic DNA, creating a unique intermolecular triple-helix DNA structure through Hoogsteen hydrogen bonding. This targeting by TFOs is site-specific and the locations TFOs bind are referred to as TFO target sites (TTS). Triplexes have been observed to selectively influence gene expression, homologous recombination, mutations, protein binding, and DNA damage. These sites typically feature a poly-purine sequence in duplex DNA, and the characteristics of these TTS sequences greatly influence the formation of the triplex. We introduce TTSBBC, a novel analysis and visualization platform designed to explore features of TTS sequences to enable users to design and validate TTSs. The web server can be freely accessed at https://kowalski-labapps.dellmed.utexas.edu/TTSBBC/.

Diversity of plant DNA in stool is linked to dietary quality, age, and household income.

Eating a varied diet is a central tenet of good nutrition. Here, we develop a molecular tool to quantify human dietary plant diversity by applying DNA metabarcoding with the chloroplast trnL-P6 marker to 1,029 fecal samples from 324 participants across two interventional feeding studies and three observational cohorts. The number of plant taxa per sample (plant metabarcoding richness or pMR) correlated with recorded intakes in interventional diets and with indices calculated from a food frequency questionnaire in typical diets (ρ = 0.40 to 0.63). In adolescents unable to collect validated dietary survey data, trnL metabarcoding detected 111 plant taxa, with 86 consumed by more than one individual and four (wheat, chocolate, corn, and potato family) consumed by >70% of individuals. Adolescent pMR was associated with age and household income, replicating prior epidemiologic findings. Overall, trnL metabarcoding promises an objective and accurate measure of the number and types of plants consumed that is applicable to diverse human populations.

Improving the Accuracy of Bulk Fitness Assays by Correcting Barcode Processing Biases.

Measuring the fitnesses of genetic variants is a fundamental objective in evolutionary biology. A standard approach for measuring microbial fitnesses in bulk involves labeling a library of genetic variants with unique sequence barcodes, competing the labeled strains in batch culture, and using deep sequencing to track changes in the barcode abundances over time. However, idiosyncratic properties of barcodes can induce nonuniform amplification or uneven sequencing coverage that causes some barcodes to be over- or under-represented in samples. This systematic bias can result in erroneous read count trajectories and misestimates of fitness. Here, we develop a computational method, named REBAR (Removing the Effects of Bias through Analysis of Residuals), for inferring the effects of barcode processing bias by leveraging the structure of systematic deviations in the data. We illustrate this approach by applying it to two independent data sets, and demonstrate that this method estimates and corrects for bias more accurately than standard proxies, such as GC-based corrections. REBAR mitigates bias and improves fitness estimates in high-throughput assays without introducing additional complexity to the experimental protocols, with potential applications in a range of experimental evolution and mutation screening contexts.

phyloBARCODER: A Web Tool for Phylogenetic Classification of Eukaryote Metabarcodes Using Custom Reference Databases.

We developed phyloBARCODER (https://github.com/jun-inoue/phyloBARCODER), a new web tool that can identify short DNA sequences to the species level using metabarcoding. phyloBARCODER estimates phylogenetic trees based on the uploaded anonymous DNA sequences and reference sequences from databases. Without such phylogenetic contexts, alternative, similarity-based methods independently identify species names and anonymous sequences of the same group by pairwise comparisons between queries and database sequences, with the caveat that they must match exactly or very closely. By putting metabarcoding sequences into a phylogenetic context, phyloBARCODER accurately identifies (i) species or classification of query sequences and (ii) anonymous sequences associated with the same species or even with populations of query sequences, with clear and accurate explanations. Version 1 of phyloBARCODER stores a database comprising all eukaryotic mitochondrial gene sequences. Moreover, by uploading their own databases, phyloBARCODER users can conduct species identification specialized for sequences obtained from a local geographic region or those of nonmitochondrial genes, e.g. ITS or rbcL.

mBARq: a versatile and user-friendly framework for the analysis of DNA barcodes from transposon insertion libraries, knockout mutants, and isogenic strain populations.

MOTIVATION: DNA barcoding has become a powerful tool for assessing the fitness of strains in a variety of studies, including random transposon mutagenesis screens, attenuation of site-directed mutants, and population dynamics of isogenic strain pools. However, the statistical analysis, visualization, and contextualization of the data resulting from such experiments can be complex and require bioinformatic skills. RESULTS: Here, we developed mBARq, a user-friendly tool designed to simplify these steps for diverse experimental setups. The tool is seamlessly integrated with an intuitive web app for interactive data exploration via the STRING and KEGG databases to accelerate scientific discovery. AVAILABILITY AND IMPLEMENTATION: The tool is implemented in Python. The source code is freely available (https://github.com/MicrobiologyETHZ/mbarq) and the web app can be accessed at: https://microbiomics.io/tools/mbarq-app.

High-throughput protein characterization by complementation using DNA barcoded fragment libraries.

Our ability to predict, control, or design biological function is fundamentally limited by poorly annotated gene function. This can be particularly challenging in non-model systems. Accordingly, there is motivation for new high-throughput methods for accurate functional annotation. Here, we used complementation of auxotrophs and DNA barcode sequencing (Coaux-Seq) to enable high-throughput characterization of protein function. Fragment libraries from eleven genetically diverse bacteria were tested in twenty different auxotrophic strains of Escherichia coli to identify genes that complement missing biochemical activity. We recovered 41% of expected hits, with effectiveness ranging per source genome, and observed success even with distant E. coli relatives like Bacillus subtilis and Bacteroides thetaiotaomicron. Coaux-Seq provided the first experimental validation for 53 proteins, of which 11 are less than 40% identical to an experimentally characterized protein. Among the unexpected function identified was a sulfate uptake transporter, an O-succinylhomoserine sulfhydrylase for methionine synthesis, and an aminotransferase. We also identified instances of cross-feeding wherein protein overexpression and nearby non-auxotrophic strains enabled growth. Altogether, Coaux-Seq's utility is demonstrated, with future applications in ecology, health, and engineering.

High-resolution lineage tracking reveals travelling wave of adaptation in laboratory yeast.

In rapidly adapting asexual populations, including many microbial pathogens and viruses, numerous mutant lineages often compete for dominance within the population1-5. These complex evolutionary dynamics determine the outcomes of adaptation, but have been difficult to observe directly. Previous studies have used whole-genome sequencing to follow molecular adaptation6-10; however, these methods have limited resolution in microbial populations. Here we introduce a renewable barcoding system to observe evolutionary dynamics at high resolution in laboratory budding yeast. We find nested patterns of interference and hitchhiking even at low frequencies. These events are driven by the continuous appearance of new mutations that modify the fates of existing lineages before they reach substantial frequencies. We observe how the distribution of fitness within the population changes over time, and find a travelling wave of adaptation that has been predicted by theory11-17. We show that clonal competition creates a dynamical 'rich-get-richer' effect: fitness advantages that are acquired early in evolution drive clonal expansions, which increase the chances of acquiring future mutations. However, less-fit lineages also routinely leapfrog over strains of higher fitness. Our results demonstrate that this combination of factors, which is not accounted for in existing models of evolutionary dynamics, is critical in determining the rate, predictability and molecular basis of adaptation.

Squalene accumulation in cholesterol auxotrophic lymphomas prevents oxidative cell death.

Cholesterol is essential for cells to grow and proliferate. Normal mammalian cells meet their need for cholesterol through its uptake or de novo synthesis1, but the extent to which cancer cells rely on each of these pathways remains poorly understood. Here, using a competitive proliferation assay on a pooled collection of DNA-barcoded cell lines, we identify a subset of cancer cells that is auxotrophic for cholesterol and thus highly dependent on its uptake. Through metabolic gene expression analysis, we pinpoint the loss of squalene monooxygenase expression as a cause of cholesterol auxotrophy, particularly in ALK+ anaplastic large cell lymphoma (ALCL) cell lines and primary tumours. Squalene monooxygenase catalyses the oxidation of squalene to 2,3-oxidosqualene in the cholesterol synthesis pathway and its loss results in accumulation of the upstream metabolite squalene, which is normally undetectable. In ALK+ ALCLs, squalene alters the cellular lipid profile and protects cancer cells from ferroptotic cell death, providing a growth advantage under conditions of oxidative stress and in tumour xenografts. Finally, a CRISPR-based genetic screen identified cholesterol uptake by the low-density lipoprotein receptor as essential for the growth of ALCL cells in culture and as patient-derived xenografts. This work reveals that the cholesterol auxotrophy of ALCLs is a targetable liability and, more broadly, that systematic approaches can be used to identify nutrient dependencies unique to individual cancer types.