Monitoring of Leukemia Clones in B-cell Acute Lymphoblastic Leukemia at Diagnosis and During Treatment by Single-cell DNA Amplicon Sequencing.

Acute lymphoblastic leukemia (ALL) is characterized by the presence of chromosomal changes, including numerical changes, translocations, and deletions, which are often associated with additional single-nucleotide mutations. In this study, we used single cell-targeted DNA sequencing to evaluate the clonal heterogeneity of B-ALL at diagnosis and during chemotherapy treatment. We designed a custom DNA amplicon library targeting mutational hotspot regions (in 110 genes) present in ALL, and we measured the presence of mutations and small insertions/deletions (indels) in bone marrow or blood samples from 12 B-ALL patients, with a median of 7973 cells per sample. Nine of the 12 cases showed at least 1 subclonal mutation, of which cases with PAX5 alterations or high hyperdiploidy (with intermediate to good prognosis) showed a high number of subclones (1 to 7) at diagnosis, defined by a variety of mutations in the JAK/STAT, RAS, or FLT3 signaling pathways. Cases with RAS pathway mutations had multiple mutations in FLT3, NRAS, KRAS, or BRAF in various clones. For those cases where we detected multiple mutational clones at diagnosis, we also studied blood samples during the first weeks of chemotherapy treatment. The leukemia clones disappeared during treatment with various kinetics, and few cells with mutations were easily detectable, even at low frequency (<0.1%). Our data illustrate that about half of the B-ALL cases show >2 subclones at diagnosis and that even very rare mutant cells can be detected at diagnosis or during treatment by single cell-targeted DNA sequencing.

Evolutionary predictability of genetic versus nongenetic resistance to anticancer drugs in melanoma.

Therapy resistance arises from heterogeneous drug-tolerant persister cells or minimal residual disease (MRD) through genetic and nongenetic mechanisms. A key question is whether specific molecular features of the MRD ecosystem determine which of these two distinct trajectories will eventually prevail. We show that, in melanoma exposed to mitogen-activated protein kinase therapeutics, emergence of a transient neural crest stem cell (NCSC) population in MRD concurs with the development of nongenetic resistance. This increase relies on a glial cell line-derived neurotrophic factor-dependent signaling cascade, which activates the AKT survival pathway in a focal adhesion kinase (FAK)-dependent manner. Ablation of the NCSC population through FAK inhibition delays relapse in patient-derived tumor xenografts. Strikingly, all tumors that ultimately escape this treatment exhibit resistance-conferring genetic alterations and increased sensitivity to extracellular signal-regulated kinase inhibition. These findings identify an approach that abrogates the nongenetic resistance trajectory in melanoma and demonstrate that the cellular composition of MRD deterministically imposes distinct drug resistance evolutionary paths.

Single-cell DNA amplicon sequencing reveals clonal heterogeneity and evolution in T-cell acute lymphoblastic leukemia.

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive leukemia that is most frequent in children and is characterized by the presence of few chromosomal rearrangements and 10 to 20 somatic mutations in protein-coding regions at diagnosis. The majority of T-ALL cases harbor activating mutations in NOTCH1 together with mutations in genes implicated in kinase signaling, transcriptional regulation, or protein translation. To obtain more insight in the level of clonal heterogeneity at diagnosis and during treatment, we used single-cell targeted DNA sequencing with the Tapestri platform. We designed a custom ALL panel and obtained accurate single-nucleotide variant and small insertion-deletion mutation calling for 305 amplicons covering 110 genes in about 4400 cells per sample and time point. A total of 108 188 cells were analyzed for 25 samples of 8 T-ALL patients. We typically observed a major clone at diagnosis (>35% of the cells) accompanied by several minor clones of which some were less than 1% of the total number of cells. Four patients had >2 NOTCH1 mutations, some of which present in minor clones, indicating a strong pressure to acquire NOTCH1 mutations in developing T-ALL cells. By analyzing longitudinal samples, we detected the presence and clonal nature of residual leukemic cells and clones with a minor presence at diagnosis that evolved to clinically relevant major clones at later disease stages. Therefore, single-cell DNA amplicon sequencing is a sensitive assay to detect clonal architecture and evolution in T-ALL.

Ethanol exposure increases mutation rate through error-prone polymerases.

Ethanol is a ubiquitous environmental stressor that is toxic to all lifeforms. Here, we use the model eukaryote Saccharomyces cerevisiae to show that exposure to sublethal ethanol concentrations causes DNA replication stress and an increased mutation rate. Specifically, we find that ethanol slows down replication and affects localization of Mrc1, a conserved protein that helps stabilize the replisome. In addition, ethanol exposure also results in the recruitment of error-prone DNA polymerases to the replication fork. Interestingly, preventing this recruitment through mutagenesis of the PCNA/Pol30 polymerase clamp or deleting specific error-prone polymerases abolishes the mutagenic effect of ethanol. Taken together, this suggests that the mutagenic effect depends on a complex mechanism, where dysfunctional replication forks lead to recruitment of error-prone polymerases. Apart from providing a general mechanistic framework for the mutagenic effect of ethanol, our findings may also provide a route to better understand and prevent ethanol-associated carcinogenesis in higher eukaryotes.

IAMBEE: a web-service for the identification of adaptive pathways from parallel evolved clonal populations.

IAMBEE is a web server designed for the Identification of Adaptive Mutations in Bacterial Evolution Experiments (IAMBEE). Input data consist of genotype information obtained from independently evolved clonal populations or strains that show the same adapted behavior (phenotype). To distinguish adaptive from passenger mutations, IAMBEE searches for neighborhoods in an organism-specific interaction network that are recurrently mutated in the adapted populations. This search for recurrently mutated network neighborhoods, as proxies for pathways is driven by additional information on the functional impact of the observed genetic changes and their dynamics during adaptive evolution. In addition, the search explicitly accounts for the differences in mutation rate between the independently evolved populations. Using this approach, IAMBEE allows exploiting parallel evolution to identify adaptive pathways. The web-server is freely available at http://bioinformatics.intec.ugent.be/iambee/ with no login requirement.

Experimental Design, Population Dynamics, and Diversity in Microbial Experimental Evolution.

In experimental evolution, laboratory-controlled conditions select for the adaptation of species, which can be monitored in real time. Despite the current popularity of such experiments, nature's most pervasive biological force was long believed to be observable only on time scales that transcend a researcher's life-span, and studying evolution by natural selection was therefore carried out solely by comparative means. Eventually, microorganisms' propensity for fast evolutionary changes proved us wrong, displaying strong evolutionary adaptations over a limited time, nowadays massively exploited in laboratory evolution experiments. Here, we formulate a guide to experimental evolution with microorganisms, explaining experimental design and discussing evolutionary dynamics and outcomes and how it is used to assess ecoevolutionary theories, improve industrially important traits, and untangle complex phenotypes. Specifically, we give a comprehensive overview of the setups used in experimental evolution. Additionally, we address population dynamics and genetic or phenotypic diversity during evolution experiments and expand upon contributing factors, such as epistasis and the consequences of (a)sexual reproduction. Dynamics and outcomes of evolution are most profoundly affected by the spatiotemporal nature of the selective environment, where changing environments might lead to generalists and structured environments could foster diversity, aided by, for example, clonal interference and negative frequency-dependent selection. We conclude with future perspectives, with an emphasis on possibilities offered by fast-paced technological progress. This work is meant to serve as an introduction to those new to the field of experimental evolution, as a guide to the budding experimentalist, and as a reference work to the seasoned expert.

CRISPR-FRT targets shared sites in a knock-out collection for off-the-shelf genome editing.

CRISPR advances genome engineering by directing endonuclease sequence specificity with a guide RNA molecule (gRNA). For precisely targeting a gene for modification, each genetic construct requires a unique gRNA. By generating a gRNA against the flippase recognition target (FRT) site, a common genetic element shared by multiple genetic collections, CRISPR-FRT circumvents this design constraint to provide a broad platform for fast, scarless, off-the-shelf genome engineering.

Hitting with a BAM: Selective Killing by Lectin-Like Bacteriocins.

Lectin-like bacteriocins (LlpAs) are secreted by proteobacteria and selectively kill strains of their own or related species, and they are composed of two B-lectin domains with divergent sequences. In Pseudomonas spp., initial binding of these antibacterial proteins to cells is mediated by the carboxy-terminal domain through d-rhamnose residues present in the common polysaccharide antigen of their lipopolysaccharide, whereas the amino-terminal domain accounts for strain selectivity of killing. Here, we show that spontaneous LlpA-resistant mutants carry mutations in one of three surface-exposed moieties of the essential ß-barrel outer membrane protein insertase BamA, the core component of the BAM complex. Polymorphism of this loop in different Pseudomonas groups is linked to LlpA susceptibility, and targeted cells all share the same signature motif in this loop. Since heterologous expression of such a bamA gene confers LlpA susceptibility upon a resistant strain, BamA represents the primary bacteriocin selectivity determinant in pseudomonads. Contrary to modular bacteriocins that require uptake via the Tol or Ton system, parasitism of BamA as an LlpA receptor advocates a novel bacteriocin killing mechanism initiated by impairment of the BAM machinery.IMPORTANCE Bacteria secrete a variety of molecules to eliminate microbial rivals. Bacteriocins are a pivotal group of peptides and proteins that assist in this fight, specifically killing related bacteria. In Gram-negative bacteria, these antibacterial proteins often comprise distinct domains for initial binding to a target cell's surface and subsequent killing via enzymatic or pore-forming activity. Here, we show that lectin-like bacteriocins, a family of bacteriocins that lack the prototypical modular toxin architecture, also stand out by parasitizing BamA, the core component of the outer membrane protein assembly machinery. A particular surface-exposed loop of BamA, critical for its function, serves as a key discriminant for cellular recognition, and polymorphisms in this loop determine whether a strain is susceptible or immune to a particular bacteriocin. These findings suggest a novel mechanism of contact-dependent killing that does not require cellular uptake. The evolutionary advantage of piracy of an essential cellular compound is highlighted by the observation that contact-dependent growth inhibition, a distinct antagonistic system, can equally take advantage of this receptor.

1-((2,4-Dichlorophenethyl)Amino)-3-Phenoxypropan-2-ol Kills Pseudomonas aeruginosa through Extensive Membrane Damage.

The ever increasing multidrug-resistance of clinically important pathogens and the lack of novel antibiotics have resulted in a true antibiotic crisis where many antibiotics are no longer effective. Further complicating the treatment of bacterial infections are antibiotic-tolerant persister cells. Besides being responsible for the recalcitrant nature of chronic infections, persister cells greatly contribute to the observed antibiotic tolerance in biofilms and even facilitate the emergence of antibiotic resistance. Evidently, eradication of these persister cells could greatly improve patient outcomes and targeting persistence may provide an alternative approach in combatting chronic infections. We recently characterized 1-((2,4-dichlorophenethyl)amino)-3-phenoxypropan-2-ol (SPI009), a novel anti-persister molecule capable of directly killing persisters from both Gram-negative and Gram-positive pathogens. SPI009 potentiates antibiotic activity in several in vitro and in vivo infection models and possesses promising anti-biofilm activity. Strikingly, SPI009 restores antibiotic sensitivity even in resistant strains. In this study, we investigated the mode of action of this novel compound using several parallel approaches. Genetic analyses and a macromolecular synthesis assays suggest that SPI009 acts by causing extensive membrane damage. This hypothesis was confirmed by liposome leakage assay and membrane permeability studies, demonstrating that SPI009 rapidly impairs the bacterial outer and inner membranes. Evaluation of SPI009-resistant mutants, which only could be generated under severe selection pressure, suggested a possible role for the MexCD-OprJ efflux pump. Overall, our results demonstrate the extensive membrane-damaging activity of SPI009 and confirm its clinical potential in the development of novel anti-persister therapies.

Network-Based Identification of Adaptive Pathways in Evolved Ethanol-Tolerant Bacterial Populations.

Efficient production of ethanol for use as a renewable fuel requires organisms with a high level of ethanol tolerance. However, this trait is complex and increased tolerance therefore requires mutations in multiple genes and pathways. Here, we use experimental evolution for a system-level analysis of adaptation of Escherichia coli to high ethanol stress. As adaptation to extreme stress often results in complex mutational data sets consisting of both causal and noncausal passenger mutations, identifying the true adaptive mutations in these settings is not trivial. Therefore, we developed a novel method named IAMBEE (Identification of Adaptive Mutations in Bacterial Evolution Experiments). IAMBEE exploits the temporal profile of the acquisition of mutations during evolution in combination with the functional implications of each mutation at the protein level. These data are mapped to a genome-wide interaction network to search for adaptive mutations at the level of pathways. The 16 evolved populations in our data set together harbored 2,286 mutated genes with 4,470 unique mutations. Analysis by IAMBEE significantly reduced this number and resulted in identification of 90 mutated genes and 345 unique mutations that are most likely to be adaptive. Moreover, IAMBEE not only enabled the identification of previously known pathways involved in ethanol tolerance, but also identified novel systems such as the AcrAB-TolC efflux pump and fatty acids biosynthesis and even allowed to gain insight into the temporal profile of adaptation to ethanol stress. Furthermore, this method offers a solid framework for identifying the molecular underpinnings of other complex traits as well.

Identification of 1-((2,4-Dichlorophenethyl)Amino)-3-Phenoxypropan-2-ol, a Novel Antibacterial Compound Active against Persisters of Pseudomonas aeruginosa.

Antibiotics typically fail to completely eradicate a bacterial population, leaving a small fraction of transiently antibiotic-tolerant persister cells intact. Persisters are therefore seen to be a major cause of treatment failure and greatly contribute to the recalcitrant nature of chronic infections. The current study focused on Pseudomonas aeruginosa, a Gram-negative pathogen belonging to the notorious ESKAPE group of pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species) and, due to increasing resistance against most conventional antibiotics, posing a serious threat to human health. Greatly contributing to the difficult treatment of P. aeruginosa infections is the presence of persister cells, and elimination of these cells would therefore significantly improve patient outcomes. In this study, a small-molecule library was screened for compounds that, in combination with the fluoroquinolone antibiotic ofloxacin, reduced the number of P. aeruginosa persisters compared to the number achieved with treatment with the antibiotic alone. Based on the early structure-activity relationship, 1-((2,4-dichlorophenethyl)amino)-3-phenoxypropan-2-ol (SPI009) was selected for further characterization. Combination of SPI009 with mechanistically distinct classes of antibiotics reduced the number of persisters up to 106-fold in both lab strains and clinical isolates of P. aeruginosa Further characterization of the compound revealed a direct and efficient killing of persister cells. SPI009 caused no erythrocyte damage and demonstrated minor cytotoxicity. In conclusion, we identified a novel antipersister compound active against P. aeruginosa with promising applications for the design of novel, case-specific combination therapies in the fight against chronic infections.

Adaptive tuning of mutation rates allows fast response to lethal stress in Escherichia coli.

While specific mutations allow organisms to adapt to stressful environments, most changes in an organism's DNA negatively impact fitness. The mutation rate is therefore strictly regulated and often considered a slowly-evolving parameter. In contrast, we demonstrate an unexpected flexibility in cellular mutation rates as a response to changes in selective pressure. We show that hypermutation independently evolves when different Escherichia coli cultures adapt to high ethanol stress. Furthermore, hypermutator states are transitory and repeatedly alternate with decreases in mutation rate. Specifically, population mutation rates rise when cells experience higher stress and decline again once cells are adapted. Interestingly, we identified cellular mortality as the major force driving the quick evolution of mutation rates. Together, these findings show how organisms balance robustness and evolvability and help explain the prevalence of hypermutation in various settings, ranging from emergence of antibiotic resistance in microbes to cancer relapses upon chemotherapy.

Elucidation of the Mode of Action of a New Antibacterial Compound Active against Staphylococcus aureus and Pseudomonas aeruginosa.

Nosocomial and community-acquired infections caused by multidrug resistant bacteria represent a major human health problem. Thus, there is an urgent need for the development of antibiotics with new modes of action. In this study, we investigated the antibacterial characteristics and mode of action of a new antimicrobial compound, SPI031 (N-alkylated 3, 6-dihalogenocarbazol 1-(sec-butylamino)-3-(3,6-dichloro-9H-carbazol-9-yl)propan-2-ol), which was previously identified in our group. This compound exhibits broad-spectrum antibacterial activity, including activity against the human pathogens Staphylococcus aureus and Pseudomonas aeruginosa. We found that SPI031 has rapid bactericidal activity (7-log reduction within 30 min at 4x MIC) and that the frequency of resistance development against SPI031 is low. To elucidate the mode of action of SPI031, we performed a macromolecular synthesis assay, which showed that SPI031 causes non-specific inhibition of macromolecular biosynthesis pathways. Liposome leakage and membrane permeability studies revealed that SPI031 rapidly exerts membrane damage, which is likely the primary cause of its antibacterial activity. These findings were supported by a mutational analysis of SPI031-resistant mutants, a transcriptome analysis and the identification of transposon mutants with altered sensitivity to the compound. In conclusion, our results show that SPI031 exerts its antimicrobial activity by causing membrane damage, making it an interesting starting point for the development of new antibacterial therapies.

Draft Genome Sequence of Pseudomonas putida BW11M1, a Banana Rhizosphere Isolate with a Diversified Antimicrobial Armamentarium.

In this study, we report the draft genome ofPseudomonas putidaBW11M1, a banana rhizosphere isolate producing various antimicrobial compounds, including a lectin-like bacteriocin, an R-type tailocin, the cyclic lipopeptide xantholysin, and the fatty acid-derived pseudopyronine.

A study of SeqA subcellular localization in Escherichia coli using photo-activated localization microscopy.

Escherichia coli (E. coli) cells replicate their genome once per cell cycle to pass on genetic information to the daughter cells. The SeqA protein binds the origin of replication, oriC, after DNA replication initiation and sequesters it from new initiations in order to prevent overinitiation. Conventional fluorescence microscopy studies of SeqA localization in bacterial cells have shown that the protein is localized to discrete foci. In this study we have used photo-activated localization microscopy (PALM) to determine the localization of SeqA molecules, tagged with fluorescent proteins, with a localization precision of 20-30 nm with the aim to visualize the SeqA subcellular structures in more detail than previously possible. SeqA-PAmCherry was imaged in wild type E. coli, expressed from plasmid or genetically engineered into the bacterial genome, replacing the native seqA gene. Unsynchronized cells as well as cells with a synchronized cell cycle were imaged at various time points, in order to investigate the evolution of SeqA localization during the cell cycle. We found that SeqA indeed localized into discrete foci but these were not the only subcellular localizations of the protein. A significant amount of SeqA-PAmCherry molecules was localized outside the foci and in a fraction of cells we saw patterns indicating localization at the membrane. Using quantitative PALM, we counted protein copy numbers per cell, protein copy numbers per focus, the numbers of foci per cell and the sizes of the SeqA clusters. The data showed broad cell-to-cell variation and we did not observe a correlation between SeqA-PAmCherry protein numbers and the cell cycle under the experimental conditions of this study. The numbers of SeqA-PAmCherry molecules per focus as well as the foci sizes also showed broad distributions indicating that the foci are likely not characterized by a fixed number of molecules. We also imaged an E. coli strain devoid of the dam methylase (Δdam) and observed that SeqA-PAmCherry no longer formed foci, and was dispersed throughout the cell and localized to the plasma membrane more readily. We discuss our results in the context of the limitations of the technique.

Frequency-based haplotype reconstruction from deep sequencing data of bacterial populations.

Clonal populations accumulate mutations over time, resulting in different haplotypes. Deep sequencing of such a population in principle provides information to reconstruct these haplotypes and the frequency at which the haplotypes occur. However, this reconstruction is technically not trivial, especially not in clonal systems with a relatively low mutation frequency. The low number of segregating sites in those systems adds ambiguity to the haplotype phasing and thus obviates the reconstruction of genome-wide haplotypes based on sequence overlap information.Therefore, we present EVORhA, a haplotype reconstruction method that complements phasing information in the non-empty read overlap with the frequency estimations of inferred local haplotypes. As was shown with simulated data, as soon as read lengths and/or mutation rates become restrictive for state-of-the-art methods, the use of this additional frequency information allows EVORhA to still reliably reconstruct genome-wide haplotypes. On real data, we show the applicability of the method in reconstructing the population composition of evolved bacterial populations and in decomposing mixed bacterial infections from clinical samples.

Identification and characterization of an anti-pseudomonal dichlorocarbazol derivative displaying anti-biofilm activity.

Pseudomonas aeruginosa strains resistant towards all currently available antibiotics are increasingly encountered, raising the need for new anti-pseudomonal drugs. We therefore conducted a medium-throughput screen of a small-molecule collection resulting in the identification of the N-alkylated 3,6-dihalogenocarbazol 1-(sec-butylamino)-3-(3,6-dichloro-9H-carbazol-9-yl)propan-2-ol (MIC = 18.5 µg mL⁻¹). This compound, compound 1, is bacteriostatic towards a broad spectrum of Gram-positive and Gram-negative pathogens, including P. aeruginosa. Importantly, 1 also eradicates mature biofilms of P. aeruginosa. 1 displays no cytotoxicity against various human cell types, pointing to its potential for further development as a novel antibacterial drug.