Beyond genome sequencing: Lineage tracking with barcodes to study the dynamics of evolution, infection, and cancer

Optimizing Nucleic Acid Delivery Systems through Barcode Technology

Conventional biological experiments often focus on in vitro assays because of the inherent limitations when handling multiple variables in vivo, including labor-intensive and time-consuming procedures. Often only a subset of samples demonstrating significant efficacy in the in vitro assays can be evaluated in vivo. Nonetheless, because of the low correlation between the in vitro and in vivo tests, evaluation of the variables under examination in vivo and not solely in vitro is critical. An emerging approach to achieve high-throughput in vivo tests involves using a barcode system consisting of various nucleotide combinations. Unique barcodes for each variant enable the simultaneous testing of multiple entities, eliminating the need for separate individual tests. Subsequently, to identify crucial parameters, samples were collected and analyzed using barcode sequencing. This review explores the development of barcode design and its applications, including the evaluation of nucleic acid delivery systems and the optimization of gene expression in vivo.

Extracting, filtering and simulating cellular barcodes using CellBarcode tools

Identifying true DNA cellular barcodes among polymerase chain reaction and sequencing errors is challenging. Current tools are restricted in the diversity of barcode types supported or the analysis strategies implemented. As such, there is a need for more versatile and efficient tools for barcode extraction, as well as for tools to investigate which factors impact barcode detection and which filtering strategies to best apply. Here we introduce the package CellBarcode and its barcode simulation kit, CellBarcodeSim, that allows efficient and versatile barcode extraction and filtering for a range of barcode types from bulk or single-cell sequencing data using a variety of filtering strategies. Using the barcode simulation kit and biological data, we explore the technical and biological factors influencing barcode identification and provide a decision tree on how to optimize barcode identification for different barcode settings. We believe that CellBarcode and CellBarcodeSim have the capability to enhance the reproducibility and interpretation of barcode results across studies.

Quantitative systems-based prediction of antimicrobial resistance evolution

Predicting evolution is a fundamental problem in biology with practical implications for treating antimicrobial resistance, which is a complex system-level phenomenon. In this perspective article, we explore the limits of predicting antimicrobial resistance evolution, quantitatively define the predictability and repeatability of microevolutionary processes, and speculate on how these quantities vary across temporal, biological, and complexity scales. The opportunities and challenges for predicting antimicrobial resistance in the context of systems biology are also discussed. Based on recent research, we conclude that the evolution of antimicrobial resistance can be predicted using a systems biology approach integrating quantitative models with multiscale data from microbial evolution experiments.

High-throughput library transgenesis in Caenorhabditis elegans via Transgenic Arrays Resulting in Diversity of Integrated Sequences (TARDIS)

High-throughput transgenesis using synthetic DNA libraries is a powerful method for systematically exploring genetic function. Diverse synthesized libraries have been used for protein engineering, identification of protein-protein interactions, characterization of promoter libraries, developmental and evolutionary lineage tracking, and various other exploratory assays. However, the need for library transgenesis has effectively restricted these approaches to single-cell models. Here, we present Transgenic Arrays Resulting in Diversity of Integrated Sequences (TARDIS), a simple yet powerful approach to large-scale transgenesis that overcomes typical limitations encountered in multicellular systems. TARDIS splits the transgenesis process into a two-step process: creation of individuals carrying experimentally introduced sequence libraries, followed by inducible extraction and integration of individual sequences/library components from the larger library cassette into engineered genomic sites. Thus, transformation of a single individual, followed by lineage expansion and functional transgenesis, gives rise to thousands of genetically unique transgenic individuals. We demonstrate the power of this system using engineered, split selectable TARDIS sites in Caenorhabditis elegans to generate (1) a large set of individually barcoded lineages and (2) transcriptional reporter lines from predefined promoter libraries. We find that this approach increases transformation yields up to approximately 1000-fold over current single-step methods. While we demonstrate the utility of TARDIS using C. elegans, in principle the process is adaptable to any system where experimentally generated genomic loci landing pads and diverse, heritable DNA elements can be generated.

gUMI-BEAR, a modular, unsupervised population barcoding method to track variants and evolution at high resolution

Cellular lineage tracking provides a means to observe population makeup at the clonal level, allowing exploration of heterogeneity, evolutionary and developmental processes and individual clones' relative fitness. It has thus contributed significantly to understanding microbial evolution, organ differentiation and cancer heterogeneity, among others. Its use, however, is limited because existing methods are highly specific, expensive, labour-intensive, and, critically, do not allow the repetition of experiments. To address these issues, we developed gUMI-BEAR (genomic Unique Molecular Identifier Barcoded Enriched Associated Regions), a modular, cost-effective method for tracking populations at high resolution. We first demonstrate the system's application and resolution by applying it to track tens of thousands of Saccharomyces cerevisiae lineages growing together under varying environmental conditions applied across multiple generations, revealing fitness differences and lineage-specific adaptations. Then, we demonstrate how gUMI-BEAR can be used to perform parallel screening of a huge number of randomly generated variants of the Hsp82 gene. We further show how our method allows isolation of variants, even if their frequency in the population is low, thus enabling unsupervised identification of modifications that lead to a behaviour of interest.

Barcoding Populations of Pseudomonas fluorescens SBW25

In recent years, evolutionary biologists have developed an increasing interest in the use of barcoding strategies to study eco-evolutionary dynamics of lineages within evolving populations and communities. Although barcoded populations can deliver unprecedented insight into evolutionary change, barcoding microbes presents specific technical challenges. Here, strategies are described for barcoding populations of the model bacterium Pseudomonas fluorescens SBW25, including the design and cloning of barcoded regions, preparation of libraries for amplicon sequencing, and quantification of resulting barcoded lineages. In so doing, we hope to aid the design and implementation of barcoding methodologies in a broad range of model and non-model organisms.

Best Practices in Designing, Sequencing, and Identifying Random DNA Barcodes

Random DNA barcodes are a versatile tool for tracking cell lineages, with applications ranging from development to cancer to evolution. Here, we review and critically evaluate barcode designs as well as methods of barcode sequencing and initial processing of barcode data. We first demonstrate how various barcode design decisions affect data quality and propose a new design that balances all considerations that we are currently aware of. We then discuss various options for the preparation of barcode sequencing libraries, including inline indices and Unique Molecular Identifiers (UMIs). Finally, we test the performance of several established and new bioinformatic pipelines for the extraction of barcodes from raw sequencing reads and for error correction. We find that both alignment and regular expression-based approaches work well for barcode extraction, and that error-correction pipelines designed specifically for barcode data are superior to generic ones. Overall, this review will help researchers to approach their barcoding experiments in a deliberate and systematic way.

Accurate estimation of molecular counts from amplicon sequence data with unique molecular identifiers

MOTIVATION

Amplicon sequencing is widely applied to explore heterogeneity and rare variants in genetic populations. Resolving true biological variants and quantifying their abundance is crucial for downstream analyses, but measured abundances are distorted by stochasticity and bias in amplification, plus errors during polymerase chain reaction (PCR) and sequencing. One solution attaches unique molecular identifiers (UMIs) to sample sequences before amplification. Counting UMIs instead of sequences provides unbiased estimates of abundance. While modern methods improve over naïve counting by UMI identity, most do not account for UMI reuse or collision, and they do not adequately model PCR and sequencing errors in the UMIs and sample sequences.

RESULTS

We introduce Deduplication and Abundance estimation with UMIs (DAUMI), a probabilistic framework to detect true biological amplicon sequences and accurately estimate their deduplicated abundance. DAUMI recognizes UMI collision, even on highly similar sequences, and detects and corrects most PCR and sequencing errors in the UMI and sampled sequences. DAUMI performs better on simulated and real data compared to other UMI-aware clustering methods.

AVAILABILITY AND IMPLEMENTATION

Source code is available at https://github.com/DormanLab/AmpliCI.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Parasite potpourri: Toxoplasma gondii lineage tracing in vivo

Using a CRISPR/Cas9-based method, Wincott et al. generated a stable, complex Toxoplasma gondii population composed of 96 barcoded clonal lineages. By tracking the population structure in vivo, they determine that - contrary to expectations - the pathway to infecting the brain is widely permissive for T. gondii.

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.

Shepherd: Accurate Clustering for Correcting DNA Barcode Errors

Motivation

DNA barcodes are short, random nucleotide sequences introduced into cell populations to track the relative counts of hundreds of thousands of individual lineages over time. Lineage tracking is widely applied, e.g. to understand evolutionary dynamics in microbial populations and the progression of breast cancer in humans. Barcode sequences are unknown upon insertion and must be identified using next-generation sequencing technology, which is error prone. In this study, we frame the barcode error correction task as a clustering problem with the aim to identify true barcode sequences from noisy sequencing data. We present Shepherd, a novel clustering method that is based on an indexing system of barcode sequences using k-mers, and a Bayesian statistical test incorporating a substitution error rate to distinguish true from error sequences.

Results

When benchmarking with synthetic data, Shepherd provides barcode count estimates that are significantly more accurate than state-of-the-art methods, producing 10–150 times fewer spurious lineages. For empirical data, Shepherd produces results that are consistent with the improvements seen on synthetic data. These improvements enable higher resolution lineage tracking and more accurate estimates of biologically relevant quantities, e.g. the detection of small effect mutations.

Availability and implementation

A Python implementation of Shepherd is freely available at: https://www.github.com/Nik-Tavakolian/Shepherd.

Supplementary information

Supplementary data are available at Bioinformatics online.

A Cell Double-Barcoding System for Quantitative Evaluation of Primary Tumors and Metastasis in Animals That Uncovers Clonal-Specific Anti-Cancer Drug Effects

Simple Summary

The main problem in treating advanced cancers is a metastatic spread when individual cancer cells leave the primary tumor and colonize to distant organs. In drug development, it is important to quantitatively assess effects of novel drug candidates on both primary tumors and metastasis. Unfortunately, current methods of monitoring metastasis in mouse models have low sensitivity and are not quantitative. Here, we developed a methodology to monitor drug effects on metastasis that is quantitative and has a very high sensitivity and resolution. In fact, it allows monitoring effects of drugs on individual cancer cells in animals.

Abstract

Imaging in monitoring metastasis in mouse models has low sensitivity and is not quantitative. Cell DNA barcoding, demonstrating high sensitivity and resolution, allows monitoring effects of drugs on the number of tumor and metastatic clones. However, this technology is not suitable for comparison of sizes of metastatic clones in different animals, for example, drug treated and untreated, due to high biological and technical variability upon tumor and metastatic growth and isolation of barcodes from tissue DNA. However, both numbers of clones and their sizes are critical parameters for analysis of drug effects. Here we developed a modification of the barcoding approach for monitoring drug effects on tumors and metastasis that is quantitative, highly sensitive and highly reproducible. This novel cell double-barcoding system allows simultaneously following the fate of two or more cell variants or cell lines in xenograft models in vivo, and also following the fates of individual clones within each of these populations. This system allows comparing effects of drugs on different cell populations and thus normalizing drug effects by drug-resistant lines, which corrects for both biological and technical variabilities and significantly increases the reproducibility of results. Using this barcoding system, we uncovered that effects of a novel DYRK1B kinase inhibitor FX9847 on primary tumors and metastasis is clone-dependent, while a distinct drug osimertinib demonstrated clone-independent effects on cancer cell populations. Overall, a cell double-barcoding approach can significantly enrich our understanding of drug effects in basic research and preclinical studies.

Bacterial age distribution in soil - Generational gaps in adjacent hot and cold spots

Resource patchiness and aqueous phase fragmentation in soil may induce large differences local growth conditions at submillimeter scales. These are translated to vast differences in bacterial age from cells dividing every thirty minutes in close proximity to plant roots to very old cells experiencing negligible growth in adjacent nutrient poor patches. In this study, we link bacterial population demographics with localized soil and hydration conditions to predict emerging generation time distributions and estimate mean bacterial cell ages using mechanistic and heuristic models of bacterial life in soil. Results show heavy-tailed distributions of generation times that resemble a power law for certain conditions, suggesting that we may find bacterial cells of vastly different ages living side by side within small soil volumes. Our results imply that individual bacteria may exist concurrently with all of their ancestors, resulting in an archive of bacterial cells with traits that have been gained (and lost) throughout time-a feature unique to microbial life. This reservoir of bacterial strains and the potential for the reemergence of rare strains with specific functions may be critical for ecosystem stability and function.

Functionalized Lineage Tracing for the Study and Manipulation of Heterogeneous Cell Populations

The ability to track and isolate unique cell lineages from large heterogeneous populations increases the resolution at which cellular processes can be understood under normal and pathogenic states beyond snapshots obtained from single-cell RNA sequencing (scRNA-seq). Here, we describe the Control of Lineages by Barcode Enabled Recombinant Transcription (COLBERT) method in which unique single guide RNA (sgRNA) barcodes are used as functional tags to identify and recall specific lineages of interest. An sgRNA barcode is stably integrated and actively transcribed, such that all cellular progeny will contain the parental barcode and produce a functional sgRNA. The sgRNA barcode has all the benefits of a DNA barcode and added functionalities. Once a barcode pertaining to a lineage of interest is identified, the lineage of interest can be isolated using an activator variant of Cas9 (such as dCas9-VPR) and a barcode-matched sequence upstream of a fluorescent reporter gene. CRISPR activation of the fluorescent reporter will only occur in cells producing the matched sgRNA barcode, allowing precise identification and isolation of lineages of interest from heterogeneous populations.

Facile Assembly of Combinatorial Mutagenesis Libraries Using Nicking Mutagenesis

Combinatorial mutagenesis is a method where multiple user-defined mutations are encoded at defined positions in a sequence. Combinatorial mutagenic libraries can be used in a variety of applications including evaluating fundamental questions about molecular evolution, directed evolution workflows for enzyme engineering, and in better understanding of biological processes like antibody affinity maturation. Here, we show a method of combinatorial mutagenesis utilizing the template-based nicking mutagenesis with several modifications. We show an example for generating a combinatorial library with 14 mutated positions, a total of 16,384 library variants, and a protocol for the generation of large, user-defined combinatorial libraries. The reader can use this protocol to create such libraries in 2 days.

Evolution of metastasis: new tools and insights

Metastasis is an evolutionary process occurring across multiple organs and timescales. Due to its continuous and dynamic nature, this multifaceted process has been challenging to investigate and remains incompletely understood, in part due to the lack of tools capable of probing genomic evolution at high enough resolution. However, technological advances in genetic sequencing and editing have provided new and powerful methods to refine our understanding of the complex series of events that lead to metastatic dissemination. In this review, we summarize the latest genetic and lineage-tracing approaches developed to unravel the genetic evolution of metastasis. The findings that have emerged have enhanced our comprehension of the mechanistic trajectories and timescales of metastasis and could provide new strategies for therapy.

Pathogen clonal expansion underlies multiorgan dissemination and organ-specific outcomes during murine systemic infection

The dissemination of pathogens through blood and their establishment within organs lead to severe clinical outcomes. However, the within-host dynamics that underlie pathogen spread to and clearance from systemic organs remain largely uncharacterized. In animal models of infection, the observed pathogen population results from the combined contributions of bacterial replication, persistence, death, and dissemination, each of which can vary across organs. Quantifying the contribution of each these processes is required to interpret and understand experimental phenotypes. Here, we leveraged STAMPR, a new barcoding framework, to investigate the population dynamics of extraintestinal pathogenic Escherichia coli, a common cause of bacteremia, during murine systemic infection. We show that while bacteria are largely cleared by most organs, organ-specific clearance failures are pervasive and result from dramatic expansions of clones representing less than 0.0001% of the inoculum. Clonal expansion underlies the variability in bacterial burden between animals, and stochastic dissemination of clones profoundly alters the pathogen population structure within organs. Despite variable pathogen expansion events, host bottlenecks are consistent yet highly sensitive to infection variables, including inoculum size and macrophage depletion. We adapted our barcoding methodology to facilitate multiplexed validation of bacterial fitness determinants identified with transposon mutagenesis and confirmed the importance of bacterial hexose metabolism and cell envelope homeostasis pathways for organ-specific pathogen survival. Collectively, our findings provide a comprehensive map of the population biology that underlies bacterial systemic infection and a framework for barcode-based high-resolution mapping of infection dynamics.

Precise measurement of the fitness effects of spontaneous mutations by droplet digital PCR in Burkholderia cenocepacia

Understanding how mutations affect survivability is a key component to knowing how organisms and complex traits evolve. However, most mutations have a minor effect on fitness and these effects are difficult to resolve using traditional molecular techniques. Therefore, there is a dire need for more accurate and precise fitness measurements methods. Here, we measured the fitness effects in Burkholderia cenocepacia HI2424 mutation accumulation (MA) lines using droplet-digital polymerase chain reaction (ddPCR). Overall, the fitness measurements from ddPCR-MA are correlated positively with fitness measurements derived from traditional phenotypic marker assays (r = 0.297, P = 0.05), but showed some differences. First, ddPCR had significantly lower measurement variance in fitness (F = 3.78, P < 2.6 × 10-13) in control experiments. Second, the mean fitness from ddPCR-MA measurements were significantly lower than phenotypic marker assays (-0.0041 vs -0.0071, P = 0.006). Consistent with phenotypic marker assays, ddPCR-MA measurements observed multiple (27/43) lineages that significantly deviated from mean fitness, suggesting that a majority of the mutations are neutral or slightly deleterious and intermixed with a few mutations that have extremely large effects. Of these mutations, we found a significant excess of mutations within DNA excinuclease and Lys R transcriptional regulators that have extreme deleterious and beneficial effects, indicating that modifications to transcription and replication may have a strong effect on organismal fitness. This study demonstrates the power of ddPCR as a ubiquitous method for high-throughput fitness measurements in both DNA- and RNA-based organisms regardless of cell type or physiology.

High-throughput characterization of mutations in genes that drive clonal evolution using multiplex adaptome capture sequencing

Understanding how cells are likely to evolve can guide medical interventions and bioengineering efforts that must contend with unwanted mutations. The adaptome of a cell-the neighborhood of genetic changes that are most likely to drive adaptation in a given environment-can be mapped by tracking rare beneficial variants during the early stages of clonal evolution. We used multiplex adaptome capture sequencing (mAdCap-seq), a procedure that combines unique molecular identifiers and hybridization-based enrichment, to characterize mutations in eight Escherichia coli genes known to be under selection in a laboratory environment. We tracked 301 mutations at frequencies as low as 0.01% and inferred the fitness effects of 240 of these mutations. There were distinct molecular signatures of selection on protein structure and function for the three genes with the most beneficial mutations. Our results demonstrate how mAdCap-seq can be used to deeply profile a targeted portion of a cell's adaptome.

Refined Quantification of Infection Bottlenecks and Pathogen Dissemination with STAMPR

Pathogen population dynamics during infection are critical determinants of infection susceptibility and define patterns of dissemination. However, deciphering these dynamics, particularly founding population sizes in host organs and patterns of dissemination between organs, is difficult because measuring bacterial burden alone is insufficient to observe these patterns. Introduction of allelic diversity into otherwise identical bacteria using DNA barcodes enables sequencing-based measurements of these parameters, in a method known as STAMP (Sequence Tag-based Analysis of Microbial Populations). However, bacteria often undergo unequal expansion within host organs, resulting in marked differences in the frequencies of barcodes in input and output libraries. Here, we show that these differences confound STAMP-based analyses of founding population sizes and dissemination patterns. We present STAMPR, a successor to STAMP, which accounts for such population expansions. Using data from systemic infection of barcoded extraintestinal pathogenic E. coli, we show that this new framework, along with the metrics it yields, enhances the fidelity of measurements of bottlenecks and dissemination patterns. STAMPR was also validated on an independent barcoded Pseudomonas aeruginosa data set, uncovering new patterns of dissemination within the data. This framework (available at https://github.com/hullahalli/stampr_rtisan), when coupled with barcoded data sets, enables a more complete assessment of within-host bacterial population dynamics. IMPORTANCE Barcoded bacteria are often employed to monitor pathogen population dynamics during infection. The accuracy of these measurements is diminished by unequal bacterial expansion rates. Here, we develop computational tools to circumvent this limitation and establish additional metrics that collectively enhance the fidelity of measuring within-host pathogen founding population sizes and dissemination patterns. These new tools will benefit future studies of the dynamics of pathogens and symbionts within their respective hosts and may have additional barcode-based applications beyond host-microbe interactions.

Does transcriptional heterogeneity facilitate the development of genetic drug resistance?

Non-genetic forms of antimicrobial (drug) resistance can result from cell-to-cell variability that is not encoded in the genetic material. Data from recent studies also suggest that non-genetic mechanisms can facilitate the development of genetic drug resistance. We speculate on how the interplay between non-genetic and genetic mechanisms may affect microbial adaptation and evolution during drug treatment. We argue that cellular heterogeneity arising from fluctuations in gene expression, epigenetic modifications, as well as genetic changes contribute to drug resistance at different timescales, and that the interplay between these mechanisms enhance pathogen resistance. Accordingly, developing a better understanding of the role of non-genetic mechanisms in drug resistance and how they interact with genetic mechanisms will enhance our ability to combat antimicrobial resistance. Also see the video abstract here: https://youtu.be/aefGpdh-bgU.

Optimization of multi-site nicking mutagenesis for generation of large, user-defined combinatorial libraries

Generating combinatorial libraries of specific sets of mutations are essential for addressing protein engineering questions involving contingency in molecular evolution, epistatic relationships between mutations, as well as functional antibody and enzyme engineering. Here we present optimization of a combinatorial mutagenesis method involving template-based nicking mutagenesis, which allows for the generation of libraries with >99% coverage for tens of thousands of user-defined variants. The non-optimized method resulted in low library coverage, which could be rationalized by a model of oligonucleotide annealing bias resulting from the nucleotide mismatch free-energy difference between mutagenic oligo and template. The optimized method mitigated this thermodynamic bias using longer primer sets and faster annealing conditions. Our updated method, applied to two antibody fragments, delivered between 99.0% (32451/32768 library members) to >99.9% coverage (32757/32768) for our desired libraries in 2 days and at an approximate 140-fold sequencing depth of coverage.

Sexual Dimorphism through the Lens of Genome Manipulation, Forward Genetics, and Spatiotemporal Sequencing

Sexual reproduction often leads to selection that favors the evolution of sex-limited traits or sex-specific variation for shared traits. These sexual dimorphisms manifest due to sex-specific genetic architectures and sex-biased gene expression across development, yet the molecular mechanisms underlying these patterns are largely unknown. The first step is to understand how sexual dimorphisms arise across the genotype-phenotype-fitness map. The emergence of "4D genome technologies" allows for efficient, high-throughput, and cost-effective manipulation and observations of this process. Studies of sexual dimorphism will benefit from combining these technological advances (e.g., precision genome editing, inducible transgenic systems, and single-cell RNA sequencing) with clever experiments inspired by classic designs (e.g., bulked segregant analysis, experimental evolution, and pedigree tracing). This perspective poses a synthetic view of how manipulative approaches coupled with cutting-edge observational methods and evolutionary theory are poised to uncover the molecular genetic basis of sexual dimorphism with unprecedented resolution. We outline hypothesis-driven experimental paradigms for identifying genetic mechanisms of sexual dimorphism among tissues, across development, and over evolutionary time.

Multiplexed Competition in a Synthetic Squid Light Organ Microbiome Using Barcode-Tagged Gene Deletions

Beneficial microbes play essential roles in the health and development of their hosts. However, the complexity of animal microbiomes and general genetic intractability of their symbionts have made it difficult to study the coevolved mechanisms for establishing and maintaining specificity at the microbe-animal host interface.

Beneficial symbioses between microbes and their eukaryotic hosts are ubiquitous and have widespread impacts on host health and development. The binary symbiosis between the bioluminescent bacterium Vibrio fischeri and its squid host Euprymna scolopes serves as a model system to study molecular mechanisms at the microbe-animal interface. To identify colonization factors in this system, our lab previously conducted a global transposon insertion sequencing (INSeq) screen and identified over 300 putative novel squid colonization factors in V. fischeri. To pursue mechanistic studies on these candidate genes, we present an approach to quickly generate barcode-tagged gene deletions and perform high-throughput squid competition experiments with detection of the proportion of each strain in the mixture by barcode sequencing (BarSeq). Our deletion approach improves on previous techniques based on splicing by overlap extension PCR (SOE-PCR) and tfoX-based natural transformation by incorporating a randomized barcode that results in unique DNA sequences within each deletion scar. Amplicon sequencing of the pool of barcoded strains before and after colonization faithfully reports on known colonization factors and provides increased sensitivity over colony counting methods. BarSeq enables rapid and sensitive characterization of the molecular factors involved in establishing the Vibrio-squid symbiosis and provides a valuable tool to interrogate the molecular dialogue at microbe-animal host interfaces.

IMPORTANCE Beneficial microbes play essential roles in the health and development of their hosts. However, the complexity of animal microbiomes and general genetic intractability of their symbionts have made it difficult to study the coevolved mechanisms for establishing and maintaining specificity at the microbe-animal host interface. Model symbioses are therefore invaluable for studying the mechanisms of beneficial microbe-host interactions. Here, we present a combined barcode-tagged deletion and BarSeq approach to interrogate the molecular dialogue that ensures specific and reproducible colonization of the Hawaiian bobtail squid by Vibrio fischeri. The ability to precisely manipulate the bacterial genome, combined with multiplex colonization assays, will accelerate the use of this valuable model system for mechanistic studies of how environmental microbes—both beneficial and pathogenic—colonize specific animal hosts.

Can mobile genetic elements rescue genes from extinction?

Bacteria and other prokaryotes evolve primarily through rapid changes in their gene content by quickly losing and gaining genes whenever an ecological opportunity emerges. As gene loss and horizontal gene transfer (HGT) appear to be the most common events across the prokaryotic tree of life, we need to think beyond gradual sequence evolution if we wish to understand the microbial world. Especially genes that reside on mobile genetic elements (MGEs) may spread much more rapidly through a microbial population than genes that reside on the bacterial chromosome. This raises the question: why are some genes associated with MGEs, while others are not? Here, I briefly review a recently proposed class of genes for which we have coined the term "rescuable genes". The fitness effect of carrying these genes is so small, either constantly or on average, that they are prone to be lost from a microbial population. I argue that HGT, even when costly to the individual cells, may play an important role in maintaining these rescuable genes in microbial communities.

High-throughput analysis of adaptation using barcoded strains of Saccharomyces cerevisiae

BACKGROUND

Experimental evolution of microbes can be used to empirically address a wide range of questions about evolution and is increasingly employed to study complex phenomena ranging from genetic evolution to evolutionary rescue. Regardless of experimental aims, fitness assays are a central component of this type of research, and low-throughput often limits the scope and complexity of experimental evolution studies. We created an experimental evolution system in Saccharomyces cerevisiae that utilizes genetic barcoding to overcome this challenge.

RESULTS

We first confirm that barcode insertions do not alter fitness and that barcode sequencing can be used to efficiently detect fitness differences via pooled competition-based fitness assays. Next, we examine the effects of ploidy, chemical stress, and population bottleneck size on the evolutionary dynamics and fitness gains (adaptation) in a total of 76 experimentally evolving, asexual populations by conducting 1,216 fitness assays and analyzing 532 longitudinal-evolutionary samples collected from the evolving populations. In our analysis of these data we describe the strengths of this experimental evolution system and explore sources of error in our measurements of fitness and evolutionary dynamics.

CONCLUSIONS

Our experimental treatments generated distinct fitness effects and evolutionary dynamics, respectively quantified via multiplexed fitness assays and barcode lineage tracking. These findings demonstrate the utility of this new resource for designing and improving high-throughput studies of experimental evolution. The approach described here provides a framework for future studies employing experimental designs that require high-throughput multiplexed fitness measurements.

Contingent evolution of alternative metabolic network topologies determines whether cross-feeding evolves

Metabolic exchange is widespread in natural microbial communities and an important driver of ecosystem structure and diversity, yet it remains unclear what determines whether microbes evolve division of labor or maintain metabolic autonomy. Here we use a mechanistic model to study how metabolic strategies evolve in a constant, one resource environment, when metabolic networks are allowed to freely evolve. We find that initially identical ancestral communities of digital organisms follow different evolutionary trajectories, as some communities become dominated by a single, autonomous lineage, while others are formed by stably coexisting lineages that cross-feed on essential building blocks. Our results show how without presupposed cellular trade-offs or external drivers such as temporal niches, diverse metabolic strategies spontaneously emerge from the interplay between ecology, spatial structure, and metabolic constraints that arise during the evolution of metabolic networks. Thus, in the long term, whether microbes remain autonomous or evolve metabolic division of labour is an evolutionary contingency.

Slightly beneficial genes are retained by bacteria evolving DNA uptake despite selfish elements

Horizontal gene transfer (HGT) and gene loss result in rapid changes in the gene content of bacteria. While HGT aids bacteria to adapt to new environments, it also carries risks such as selfish genetic elements (SGEs). Here, we use modelling to study how HGT of slightly beneficial genes impacts growth rates of bacterial populations, and if bacterial collectives can evolve to take up DNA despite selfish elements. We find four classes of slightly beneficial genes: indispensable, enrichable, rescuable, and unrescuable genes. Rescuable genes - genes with small fitness benefits that are lost from the population without HGT - can be collectively retained by a community that engages in costly HGT. While this 'gene-sharing' cannot evolve in well-mixed cultures, it does evolve in a spatial population like a biofilm. Despite enabling infection by harmful SGEs, the uptake of foreign DNA is evolutionarily maintained by the hosts, explaining the coexistence of bacteria and SGEs.

Diversity in biology: definitions, quantification and models

Diversity indices are useful single-number metrics for characterizing a complex distribution of a set of attributes across a population of interest. The utility of these different metrics or sets of metrics depends on the context and application, and whether a predictive mechanistic model exists. In this topical review, we first summarize the relevant mathematical principles underlying heterogeneity in a large population, before outlining the various definitions of 'diversity' and providing examples of scientific topics in which its quantification plays an important role. We then review how diversity has been a ubiquitous concept across multiple fields, including ecology, immunology, cellular barcoding experiments, and socioeconomic studies. Since many of these applications involve sampling of populations, we also review how diversity in small samples is related to the diversity in the entire population. Features that arise in each of these applications are highlighted.

Chromosomal barcoding of E. coli populations reveals lineage diversity dynamics at high resolution

Evolutionary dynamics in large asexual populations is strongly influenced by multiple competing beneficial lineages, most of which segregate at very low frequencies. However, technical barriers to tracking a large number of these rare lineages in bacterial populations have so far prevented a detailed elucidation of evolutionary dynamics. Here, we overcome this hurdle by developing a chromosomal-barcoding technique that allows simultaneous tracking of approximately 450,000 distinct lineages in Escherichia coli, which we use to test the effect of sub-inhibitory concentrations of common antibiotics on the evolutionary dynamics of low-frequency lineages. We find that populations lose lineage diversity at distinct rates that correspond to their antibiotic regimen. We also determine that some lineages have similar fates across independent experiments. By analysing the trajectory dynamics, we attribute the reproducible fates of these lineages to the presence of pre-existing beneficial mutations, and we demonstrate how the relative contribution of pre-existing and de novo mutations varies across drug regimens. Finally, we reproduce the observed lineage dynamics by simulations. Altogether, our results provide a valuable methodology for studying bacterial evolution as well as insights into evolution under sub-inhibitory antibiotic levels.

Cellular Barcoding Identifies Clonal Substitution as a Hallmark of Local Recurrence in a Surgical Model of Head and Neck Squamous Cell Carcinoma

Local recurrence after surgery for head and neck squamous cell carcinoma (HNSCC) remains a common event associated with a dismal prognosis. Improving this outcome requires a better understanding of cancer cell populations that expand from postsurgical minimal residual disease (MRD). Therefore, we assessed clonal dynamics in a surgical model of barcoded HNSCC growing in the submental region of immunodeficient mice. Clonal substitution and massive reduction of clonal heterogeneity emerged as hallmarks of local recurrence, as the clones dominating in less heterogeneous recurrences were scarce in their matched primary tumors. These lineages were selected by their ability to persist after surgery and competitively expand from MRD. Clones enriched in recurrences exhibited both private and shared genetic features and likely originated from ancestors shared with clones dominating in primary tumors. They demonstrated high invasiveness and epithelial-to-mesenchymal transition, eventually providing an attractive target for obtaining better local control for these tumors.

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 population^1-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 adaptation^6-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 theory^11-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.

Trusting the hand that feeds: microbes evolve to anticipate a serial transfer protocol as individuals or collectives

BACKGROUND

Experimental evolution of microbes often involves a serial transfer protocol, where microbes are repeatedly diluted by transfer to a fresh medium, starting a new growth cycle. This has revealed that evolution can be remarkably reproducible, where microbes show parallel adaptations both on the level of the phenotype as well as the genotype. However, these studies also reveal a strong potential for divergent evolution, leading to diversity both between and within replicate populations. We here study how in silico evolved Virtual Microbe "wild types" (WTs) adapt to a serial transfer protocol to investigate generic evolutionary adaptations, and how these adaptations can be manifested by a variety of different mechanisms.

RESULTS

We show that all WTs evolve to anticipate the regularity of the serial transfer protocol by adopting a fine-tuned balance of growth and survival. This anticipation is done by evolving either a high yield mode, or a high growth rate mode. We find that both modes of anticipation can be achieved by individual lineages and by collectives of microbes. Moreover, these different outcomes can be achieved with or without regulation, although the individual-based anticipation without regulation is less well adapted in the high growth rate mode.

CONCLUSIONS

All our in silico WTs evolve to trust the hand that feeds by evolving to anticipate the periodicity of a serial transfer protocol, but can do so by evolving two distinct growth strategies. Furthermore, both these growth strategies can be accomplished by gene regulation, a variety of different polymorphisms, and combinations thereof. Our work reveals that, even under controlled conditions like those in the lab, it may not be possible to predict individual evolutionary trajectories, but repeated experiments may well result in only a limited number of possible outcomes.

Key challenges facing data-driven multicellular systems biology

Increasingly sophisticated experiments, coupled with large-scale computational models, have the potential to systematically test biological hypotheses to drive our understanding of multicellular systems. In this short review, we explore key challenges that must be overcome to achieve robust, repeatable data-driven multicellular systems biology. If these challenges can be solved, we can grow beyond the current state of isolated tools and datasets to a community-driven ecosystem of interoperable data, software utilities, and computational modeling platforms. Progress is within our grasp, but it will take community (and financial) commitment.

Cancer stem cells in radiation response: current views and future perspectives in radiation oncology

Purpose: Despite technological improvement and advances in biology-driven patient stratification, many patients still fail radiotherapy resulting in loco-regional and distant recurrence. Tumor heterogeneity remains a key challenge to effective cancer treatment, and reliable stratification of cancer patients for prediction of outcomes is highly important. Intratumoral heterogeneity is manifested at the different levels, including different tumorigenic properties of cancer cells. Since John Dick et al. isolated leukemia initiating cells in 1990, the populations of tumor initiating or cancer stem cells (CSCs) were identified and characterized also for a broad spectrum of solid tumor types. The properties of CSCs are of considerable clinical relevance: CSCs have self-renewal and tumor initiating potential, and the metastases are initiated by the CSC clones with the ability to disseminate from the primary tumor site. Conclusion: Evidence from both, experimental and clinical studies demonstrates that the probability of achieving local tumor control by radiation therapy depends on the complete eradication of CSC populations. The number, properties and molecular signature of CSCs are highly predictive for clinical outcome of radiotherapy, whereas targeted therapies against CSCs combined with conventional treatment are expected to provide an improved clinical response and prevent tumor relapse. In this review, we discuss the modern methods to study CSCs in radiation biology, the role of CSCs in personalized cancer therapy as well as future directions for CSC research in translational radiooncology.

SodaPop: a forward simulation suite for the evolutionary dynamics of asexual populations on protein fitness landscapes

Motivation

Protein evolution is determined by forces at multiple levels of biological organization. Random mutations have an immediate effect on the biophysical properties, structure and function of proteins. These same mutations also affect the fitness of the organism. However, the evolutionary fate of mutations, whether they succeed to fixation or are purged, also depends on population size and dynamics. There is an emerging interest, both theoretically and experimentally, to integrate these two factors in protein evolution. Although there are several tools available for simulating protein evolution, most of them focus on either the biophysical or the population-level determinants, but not both. Hence, there is a need for a publicly available computational tool to explore both the effects of protein biophysics and population dynamics on protein evolution.

Results

To address this need, we developed SodaPop, a computational suite to simulate protein evolution in the context of the population dynamics of asexual populations. SodaPop accepts as input several fitness landscapes based on protein biochemistry or other user-defined fitness functions. The user can also provide as input experimental fitness landscapes derived from deep mutational scanning approaches or theoretical landscapes derived from physical force field estimates. Here, we demonstrate the broad utility of SodaPop with different applications describing the interplay of selection for protein properties and population dynamics. SodaPop is designed such that population geneticists can explore the influence of protein biochemistry on patterns of genetic variation, and that biochemists and biophysicists can explore the role of population size and demography on protein evolution.

Availability and implementation

Source code and binaries are freely available at https://github.com/louisgt/SodaPop under the GNU GPLv3 license. The software is implemented in C++ and supported on Linux, Mac OS/X and Windows.

Supplementary information

Supplementary data are available at Bioinformatics online.

Exploiting Molecular Barcodes in High-Throughput Cellular Assays

Multiplexing strategies, which greatly increase the number of simultaneously measured parameters in single experiments, are now being widely implemented by both the pharmaceutical industry and academic researchers. Color has long been used to identify biological signals and, when combined with molecular barcodes, has substantially enhanced the depth of multiplexed sample characterization. Moreover, the recent advent of DNA barcodes has led to an explosion of innovative cell sequencing approaches. Novel barcoding strategies also show great promise for encoding spatial information in transcriptomic studies, and for precise assessment of molecular abundance. Both color- and DNA-based barcodes can be conveniently analyzed with either a microscope or a cytometer, or via DNA sequencing. Here we review the basic principles of several technologies used to create barcodes and detail the type of samples that can be identified with such tags.

Experimental evolution: its principles and applications in developing stress-tolerant yeasts

Stress tolerance and resistance in industrial yeast strains are important attributes for cost-effective bioprocessing. The source of stress-tolerant yeasts ranges from extremophilic environments to laboratory engineered strains. However, industrial stress-tolerant yeasts are very rare in nature as the natural environment forces them to evolve traits that optimize survival and reproduction and not the ability to withstand harsh habitat-irrelevant industrial conditions. Experimental evolution is a frequent method used to uncover the mechanisms of evolution and microbial adaption towards environmental stresses. It optimizes biological systems by means of adaptation to environmental stresses and thus has immense power of development of robust stress-tolerant yeasts. This mini-review briefly outlines the basics and implications of evolution experiments and their applications to industrial biotechnology. This work is meant to serve as an introduction to those new to the field of experimental evolution, and as a guide to biologists working in the field of yeast stress response. Future perspectives of experimental evolution for potential biotechnological applications have also been elucidated.

Bartender: a fast and accurate clustering algorithm to count barcode reads

Motivation

Barcode sequencing (bar-seq) is a high-throughput, and cost effective method to assay large numbers of cell lineages or genotypes in complex cell pools. Because of its advantages, applications for bar-seq are quickly growing-from using neutral random barcodes to study the evolution of microbes or cancer, to using pseudo-barcodes, such as shRNAs or sgRNAs to simultaneously screen large numbers of cell perturbations. However, the computational pipelines for bar-seq clustering are not well developed. Available methods often yield a high frequency of under-clustering artifacts that result in spurious barcodes, or over-clustering artifacts that group distinct barcodes together. Here, we developed Bartender, an accurate clustering algorithm to detect barcodes and their abundances from raw next-generation sequencing data.

Results

In contrast with existing methods that cluster based on sequence similarity alone, Bartender uses a modified two-sample proportion test that also considers cluster size. This modification results in higher accuracy and lower rates of under- and over-clustering artifacts. Additionally, Bartender includes unique molecular identifier handling and a 'multiple time point' mode that matches barcode clusters between different clustering runs for seamless handling of time course data. Bartender is a set of simple-to-use command line tools that can be performed on a laptop at comparable run times to existing methods.

Availability and implementation

Bartender is available at no charge for non-commercial use at https://github.com/LaoZZZZZ/bartender-1.1.

Contact

sasha.levy@stonybrook.edu or song.wu@stonybrook.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Control of Lineage-Specific Gene Expression by Functionalized gRNA Barcodes

Lineage tracking delivers essential quantitative insight into dynamic, probabilistic cellular processes, such as somatic tumor evolution and differentiation. Methods for high diversity lineage quantitation rely on sequencing a population of DNA barcodes. However, manipulation of specific individual lineages is not possible with this approach. To address this challenge, we developed a functionalized lineage tracing tool, Control of Lineages by Barcode Enabled Recombinant Transcription (COLBERT), that enables high diversity lineage tracing and lineage-specific manipulation of gene expression. This modular platform utilizes expressed barcode gRNAs to both track cell lineages and direct lineage-specific gene expression.

Experimental Studies of Evolutionary Dynamics in Microbes

Evolutionary dynamics in laboratory microbial evolution experiments can be surprisingly complex. In the past two decades, observations of these dynamics have challenged simple models of adaptation and have shown that clonal interference, hitchhiking, ecological diversification, and contingency are widespread. In recent years, advances in high-throughput strain maintenance and phenotypic assays, the dramatically reduced cost of genome sequencing, and emerging methods for lineage barcoding have made it possible to observe evolutionary dynamics at unprecedented resolution. These new methods can now begin to provide detailed measurements of key aspects of fitness landscapes and of evolutionary outcomes across a range of systems. These measurements can highlight challenges to existing theoretical models and guide new theoretical work towards the complications that are most widely important.

Chromosomal barcoding as a tool for multiplexed phenotypic characterization of laboratory evolved lineages

Adaptive laboratory evolution is an important tool to evolve organisms to increased tolerance towards different physical and chemical stress. It is applied to study the evolution of antibiotic resistance as well as genetic mechanisms underlying improvements in production strains. Adaptive evolution experiments can be automated in a high-throughput fashion. However, the characterization of the resulting lineages can become a time consuming task, when the performance of each lineage is evaluated individually. Here, we present a novel method for the markerless insertion of randomized genetic barcodes into the genome of Escherichia coli using a novel dual-auxotrophic selection approach. The barcoded E. coli library allows multiplexed phenotyping of evolved strains in pooled competition experiments. We use the barcoded library in an adaptive evolution experiment; evolving resistance towards three common antibiotics. Comparing this multiplexed phenotyping with conventional susceptibility testing and growth-rate measurements we can show a significant positive correlation between the two approaches. Use of barcoded bacterial strain libraries for individual adaptive evolution experiments drastically reduces the workload of characterizing the resulting phenotypes and enables prioritization of lineages for in-depth characterization. In addition, barcoded clones open up new ways to profile community dynamics or to track lineages in vivo or situ.

DNA Barcoding Reveals Habitual Clonal Dominance of Myeloma Plasma Cells in the Bone Marrow Microenvironment

Multiple myeloma (MM) is a hematological malignancy resulting from the uncontrolled proliferation of antibody-producing plasma cells in the bone marrow. At diagnosis, independent plasma cell tumors are found throughout the skeleton. The recirculation of mutant plasma cells from the initial lesion and their recolonization of distant marrow sites are thought to occur by a process similar to solid tumor metastasis. However, the efficiency of this bone marrow homing process and the proportion of disseminated cells that actively divide and contribute to new tumor growth in MM are both unknown. We used the C57BL/KaLwRij mouse model of myeloma, lentiviral-mediated DNA barcoding of 5TGM1 myeloma cells, and next-generation sequencing to investigate the relative efficiency of plasma cell migration to, and growth within, the bone marrow. This approach revealed three major findings: firstly, establishment of metastasis within the bone marrow was extremely inefficient, with approximately 0.01% of circulating myeloma cells becoming resident long term in the bone marrow of each long bone; secondly, the individual cells of each metastasis exhibited marked differences in their proliferative fates, with the majority of final tumor burden within a bone being attributable to the progeny of between 1 and 8 cells; and, thirdly, the proliferative fate of individual clonal plasma cells differed at each bone marrow site in which the cells “landed.” These findings suggest that individual myeloma plasma cells are subjected to vastly different selection pressures within the bone marrow microenvironment, highlighting the importance of niche-driven factors, which determine the disease course and outcome.

Variability in fitness effects can preclude selection of the fittest

Evolutionary biologists often predict the outcome of natural selection on an allele by measuring its effects on lifetime survival and reproduction of individual carriers. However, alleles affecting traits like sex, evolvability, and cooperation can cause fitness effects that depend heavily on differences in the environmental, social, and genetic context of individuals carrying the allele. This variability makes it difficult to summarize the evolutionary fate of an allele based solely on its effects on any one individual. Attempts to average over this variability can sometimes salvage the concept of fitness. In other cases evolutionary outcomes can only be predicted by considering the entire genealogy of an allele, thus limiting the utility of individual fitness altogether. We describe a number of intriguing new evolutionary phenomena that have emerged in studies that explicitly model long-term lineage dynamics and discuss implications for the evolution of infectious diseases.

Digitally Barcoding Mycobacterium tuberculosis Reveals In Vivo Infection Dynamics in the Macaque Model of Tuberculosis

Infection with Mycobacterium tuberculosis causes a spectrum of outcomes; the majority of individuals contain but do not eliminate the infection, while a small subset present with primary active tuberculosis (TB) disease. This variability in infection outcomes is recapitulated at the granuloma level within each host, such that some sites of infection can be fully cleared while others progress. Understanding the spectrum of TB outcomes requires new tools to deconstruct the mechanisms underlying differences in granuloma fate. Here, we use novel genome-encoded barcodes to uniquely tag individual M. tuberculosis bacilli, enabling us to quantitatively track the trajectory of each infecting bacterium in a macaque model of TB. We also introduce a robust bioinformatics pipeline capable of identifying and counting barcode sequences within complex mixtures and at various read depths. By coupling this tagging strategy with serial positron emission tomography coregistered with computed tomography (PET/CT) imaging of lung pathology in macaques, we define a lesional map of M. tuberculosis infection dynamics. We find that there is no significant infection bottleneck, but there are significant constraints on productive bacterial trafficking out of primary granulomas. Our findings validate our barcoding approach and demonstrate its utility in probing lesion-specific biology and dissemination. This novel technology has the potential to greatly enhance our understanding of local dynamics in tuberculosis.

Classically, M. tuberculosis infection was thought to result in either latent infection or active disease. More recently, the field has recognized that there is a spectrum of M. tuberculosis infection clinical outcomes. Within a single host, this spectrum is recapitulated at the granuloma level, where there can simultaneously be lesional sterilization and poorly contained disease. To better understand the lesional biology of TB infection, we digitally barcoded M. tuberculosis to quantitatively track the fate of each infecting bacterium. By combining this technology with serial PET-CT imaging, we can dynamically track both bacterial populations and granuloma trajectories. We demonstrate that there is little constraint on the bacterial population at the time of infection. However, the granuloma imposes a strong bottleneck on dissemination, and the subset of granulomas at risk of dissemination can be distinguished by physical features.