Reproducible evaluation of transposable element detectors with McClintock 2 guides accurate inference of Ty insertion patterns in yeast.

BACKGROUND: Many computational methods have been developed to detect non-reference transposable element (TE) insertions using short-read whole genome sequencing data. The diversity and complexity of such methods often present challenges to new users seeking to reproducibly install, execute, or evaluate multiple TE insertion detectors. RESULTS: We previously developed the McClintock meta-pipeline to facilitate the installation, execution, and evaluation of six first-generation short-read TE detectors. Here, we report a completely re-implemented version of McClintock written in Python using Snakemake and Conda that improves its installation, error handling, speed, stability, and extensibility. McClintock 2 now includes 12 short-read TE detectors, auxiliary pre-processing and analysis modules, interactive HTML reports, and a simulation framework to reproducibly evaluate the accuracy of component TE detectors. When applied to the model microbial eukaryote Saccharomyces cerevisiae, we find substantial variation in the ability of McClintock 2 components to identify the precise locations of non-reference TE insertions, with RelocaTE2 showing the highest recall and precision in simulated data. We find that RelocaTE2, TEMP, TEMP2 and TEBreak provide consistent estimates of [Formula: see text]50 non-reference TE insertions per strain and that Ty2 has the highest number of non-reference TE insertions in a species-wide panel of [Formula: see text]1000 yeast genomes. Finally, we show that best-in-class predictors for yeast applied to resequencing data have sufficient resolution to reveal a dyad pattern of integration in nucleosome-bound regions upstream of yeast tRNA genes for Ty1, Ty2, and Ty4, allowing us to extend knowledge about fine-scale target preferences revealed previously for experimentally-induced Ty1 insertions to spontaneous insertions for other copia-superfamily retrotransposons in yeast. CONCLUSION: McClintock ( https://github.com/bergmanlab/mcclintock/ ) provides a user-friendly pipeline for the identification of TEs in short-read WGS data using multiple TE detectors, which should benefit researchers studying TE insertion variation in a wide range of different organisms. Application of the improved McClintock system to simulated and empirical yeast genome data reveals best-in-class methods and novel biological insights for one of the most widely-studied model eukaryotes and provides a paradigm for evaluating and selecting non-reference TE detectors in other species.

Reproducible evaluation of transposable element detectors with McClintock 2 guides accurate inference of Ty insertion patterns in yeast.

BACKGROUND: Many computational methods have been developed to detect non-reference transposable element (TE) insertions using short-read whole genome sequencing data. The diversity and complexity of such methods often present challenges to new users seeking to reproducibly install, execute, or evaluate multiple TE insertion detectors. RESULTS: We previously developed the McClintock meta-pipeline to facilitate the installation, execution, and evaluation of six first-generation short-read TE detectors. Here, we report a completely re-implemented version of McClintock written in Python using Snakemake and Conda that improves its installation, error handling, speed, stability, and extensibility. McClintock 2 now includes 12 short-read TE detectors, auxiliary pre-processing and analysis modules, interactive HTML reports, and a simulation framework to reproducibly evaluate the accuracy of component TE detectors. When applied to the model microbial eukaryote Saccharomyces cerevisiae, we find substantial variation in the ability of McClintock 2 components to identify the precise locations of non-reference TE insertions, with RelocaTE2 showing the highest recall and precision in simulated data. We find that RelocaTE2, TEMP, TEMP2 and TEBreak provide a consistent and biologically meaningful view of non-reference TE insertions in a species-wide panel of ∻1000 yeast genomes, as evaluated by coverage-based abundance estimates and expected patterns of tRNA promoter targeting. Finally, we show that best-in-class predictors for yeast have sufficient resolution to reveal a dyad pattern of integration in nucleosome-bound regions upstream of yeast tRNA genes for Ty1, Ty2, and Ty4, allowing us to extend knowledge about fine-scale target preferences first revealed experimentally for Ty1 to natural insertions and related copia-superfamily retrotransposons in yeast. CONCLUSION: McClintock (https://github.com/bergmanlab/mcclintock/) provides a user-friendly pipeline for the identification of TEs in short-read WGS data using multiple TE detectors, which should benefit researchers studying TE insertion variation in a wide range of different organisms. Application of the improved McClintock system to simulated and empirical yeast genome data reveals best-in-class methods and novel biological insights for one of the most widely-studied model eukaryotes and provides a paradigm for evaluating and selecting non-reference TE detectors for other species.

Local assembly of long reads enables phylogenomics of transposable elements in a polyploid cell line.

Animal cell lines often undergo extreme genome restructuring events, including polyploidy and segmental aneuploidy that can impede de novo whole-genome assembly (WGA). In some species like Drosophila, cell lines also exhibit massive proliferation of transposable elements (TEs). To better understand the role of transposition during animal cell culture, we sequenced the genome of the tetraploid Drosophila S2R+ cell line using long-read and linked-read technologies. WGAs for S2R+ were highly fragmented and generated variable estimates of TE content across sequencing and assembly technologies. We therefore developed a novel WGA-independent bioinformatics method called TELR that identifies, locally assembles, and estimates allele frequency of TEs from long-read sequence data (https://github.com/bergmanlab/telr). Application of TELR to a ∼130x PacBio dataset for S2R+ revealed many haplotype-specific TE insertions that arose by transposition after initial cell line establishment and subsequent tetraploidization. Local assemblies from TELR also allowed phylogenetic analysis of paralogous TEs, which revealed that proliferation of TE families in vitro can be driven by single or multiple source lineages. Our work provides a model for the analysis of TEs in complex heterozygous or polyploid genomes that are recalcitrant to WGA and yields new insights into the mechanisms of genome evolution in animal cell culture.

Ongoing transposition in cell culture reveals the phylogeny of diverse Drosophila S2 sublines.

Cultured cells are widely used in molecular biology despite poor understanding of how cell line genomes change in vitro over time. Previous work has shown that Drosophila cultured cells have a higher transposable element content than whole flies, but whether this increase in transposable element content resulted from an initial burst of transposition during cell line establishment or ongoing transposition in cell culture remains unclear. Here, we sequenced the genomes of 25 sublines of Drosophila S2 cells and show that transposable element insertions provide abundant markers for the phylogenetic reconstruction of diverse sublines in a model animal cell culture system. DNA copy number evolution across S2 sublines revealed dramatically different patterns of genome organization that support the overall evolutionary history reconstructed using transposable element insertions. Analysis of transposable element insertion site occupancy and ancestral states support a model of ongoing transposition dominated by episodic activity of a small number of retrotransposon families. Our work demonstrates that substantial genome evolution occurs during long-term Drosophila cell culture, which may impact the reproducibility of experiments that do not control for subline identity.

Transposable element profiles reveal cell line identity and loss of heterozygosity in Drosophila cell culture.

Cell culture systems allow key insights into biological mechanisms yet suffer from irreproducible outcomes in part because of cross-contamination or mislabeling of cell lines. Cell line misidentification can be mitigated by the use of genotyping protocols, which have been developed for human cell lines but are lacking for many important model species. Here, we leverage the classical observation that transposable elements (TEs) proliferate in cultured Drosophila cells to demonstrate that genome-wide TE insertion profiles can reveal the identity and provenance of Drosophila cell lines. We identify multiple cases where TE profiles clarify the origin of Drosophila cell lines (Sg4, mbn2, and OSS_E) relative to published reports, and also provide evidence that insertions from only a subset of long-terminal repeat retrotransposon families are necessary to mark Drosophila cell line identity. We also develop a new bioinformatics approach to detect TE insertions and estimate intra-sample allele frequencies in legacy whole-genome sequencing data (called ngs_te_mapper2), which revealed loss of heterozygosity as a mechanism shaping the unique TE profiles that identify Drosophila cell lines. Our work contributes to the general understanding of the forces impacting metazoan genomes as they evolve in cell culture and paves the way for high-throughput protocols that use TE insertions to authenticate cell lines in Drosophila and other organisms.

A Meta-Analysis of Wolbachia Transcriptomics Reveals a Stage-Specific Wolbachia Transcriptional Response Shared Across Different Hosts.

Wolbachia is a genus containing obligate, intracellular endosymbionts with arthropod and nematode hosts. Numerous studies have identified differentially expressed transcripts in Wolbachia endosymbionts that potentially inform the biological interplay between these endosymbionts and their hosts, albeit with discordant results. Here, we re-analyze previously published Wolbachia RNA-Seq transcriptomics data sets using a single workflow consisting of the most up-to-date algorithms and techniques, with the aim of identifying trends or patterns in the pan-Wolbachia transcriptional response. We find that data from one of the early studies in filarial nematodes did not allow for robust conclusions about Wolbachia differential expression with these methods, suggesting the original interpretations should be reconsidered. Across datasets analyzed with this unified workflow, there is a general lack of global gene regulation with the exception of a weak transcriptional response resulting in the upregulation of ribosomal proteins in early larval stages. This weak response is observed across diverse Wolbachia strains from both nematode and insect hosts suggesting a potential pan-Wolbachia transcriptional response during host development that diverged more than 700 million years ago.

Complete Genome Assemblies for Three Variants of the Wolbachia Endosymbiont of Drosophila melanogaster.

Here, we report genome assemblies for three strains of Wolbachia pipientis, assembled from unenriched, unfiltered long-read shotgun sequencing data of geographically distinct strains of Drosophila melanogaster Our simple methodology can be applied to long-read data sets of other Wolbachia-infected species with limited Wolbachia-host lateral gene transfers to produce complete assemblies for this important model symbiont.

Inverted Regulation of Multidrug Efflux Pumps, Acid Resistance, and Porins in Benzoate-Evolved Escherichia coli K-12.

Benzoic acid, a partial uncoupler of the proton motive force (PMF), selects for sensitivity to chloramphenicol and tetracycline during the experimental evolution of Escherichia coli K-12. Transcriptomes of E. coli isolates evolved with benzoate showed the reversal of benzoate-dependent regulation, including the downregulation of multidrug efflux pump genes, the gene for the Gad acid resistance regulon, the nitrate reductase genes narHJ, and the gene for the acid-consuming hydrogenase Hyd-3. However, the benzoate-evolved strains had increased expression of OmpF and other large-hole porins that admit fermentable substrates and antibiotics. Candidate genes identified from benzoate-evolved strains were tested for their roles in benzoate tolerance and in chloramphenicol sensitivity. Benzoate or salicylate tolerance was increased by deletion of the Gad activator ariR or of the acid fitness island from slp to the end of the gadX gene encoding Gad regulators and the multidrug pump genes mdtEF Benzoate tolerance was also increased by deletion of multidrug component gene emrA, RpoS posttranscriptional regulator gene cspC, adenosine deaminase gene add, hydrogenase gene hyc (Hyd-3), and the RNA chaperone/DNA-binding regulator gene hfq Chloramphenicol resistance was decreased by mutations in genes for global regulators, such as RNA polymerase alpha subunit gene rpoA, the Mar activator gene rob, and hfq Deletion of lipopolysaccharide biosynthetic kinase gene rfaY decreased the rate of growth in chloramphenicol. Isolates from experimental evolution with benzoate had many mutations affecting aromatic biosynthesis and catabolism, such as aroF (encoding tyrosine biosynthesis) and apt (encoding adenine phosphoribosyltransferase). Overall, benzoate or salicylate exposure selects for the loss of multidrug efflux pumps and of hydrogenases that generate a futile cycle of PMF and upregulates porins that admit fermentable nutrients and antibiotics.IMPORTANCE Benzoic acid is a common food preservative, and salicylic acid (2-hydroxybenzoic acid) is the active form of aspirin. At high concentrations, benzoic acid conducts a proton across the membrane, depleting the proton motive force. In the absence of antibiotics, benzoate exposure selects against proton-driven multidrug efflux pumps and upregulates porins that admit fermentable substrates but that also allow the entry of antibiotics. Thus, evolution with benzoate and related molecules, such as salicylates, requires a trade-off for antibiotic sensitivity, a trade-off that could help define a stable gut microbiome. Benzoate and salicylate are naturally occurring plant signal molecules that may modulate the microbiomes of plants and animal digestive tracts so as to favor fermenters and exclude drug-resistant pathogens.

Experimental Evolution of Escherichia coli K-12 in the Presence of Proton Motive Force (PMF) Uncoupler Carbonyl Cyanide m-Chlorophenylhydrazone Selects for Mutations Affecting PMF-Driven Drug Efflux Pumps.

Experimental evolution of Escherichia coli K-12 with benzoate, a partial uncoupler of the proton motive force (PMF), selects for mutations that decrease antibiotic resistance. We conducted experimental evolution in the presence of carbonyl cyanide m-chlorophenylhydrazone (CCCP), a strong uncoupler. Cultures were serially diluted daily 1:100 in LBK medium containing 20 to 150 µM CCCP buffered at pH 6.5 or at pH 8.0. After 1,000 generations, the populations tolerated up to 150 µM CCCP. Sequenced isolates had mutations in mprA (emrR), which downregulates the EmrAB-TolC pump that exports CCCP. A mprA::kanR deletion conferred growth at 60 µM CCCP, though not at the higher levels resisted by evolved strains (150 µM). Some mprA mutant strains also had point mutations affecting emrA, but deletion of emrA abolished the CCCP resistance. Thus, CCCP-evolved isolates contained additional adaptations. One isolate lacked emrA or mprA mutations but had mutations in cecR (ybiH), whose product upregulates drug pumps YbhG and YbhFSR, and in gadE, which upregulates the multidrug pump MdtEF. A cecR::kanR deletion conferred partial resistance to CCCP. Other multidrug efflux genes that had mutations included ybhR and acrAB The acrB isolate was sensitive to the AcrAB substrates chloramphenicol and tetracycline. Other mutant genes in CCCP-evolved strains include rng (RNase G) and cyaA (adenylate cyclase). Overall, experimental evolution revealed a CCCP-dependent fitness advantage for mutations increasing CCCP efflux via EmrA and for mutations that may deactivate proton-driven pumps for drugs not present (cecR, gadE, acrAB, and ybhR). These results are consistent with our previous report of drug sensitivity associated with evolved benzoate tolerance.IMPORTANCE The genetic responses of bacteria to depletion of proton motive force (PMF), and their effects on drug resistance, are poorly understood. PMF drives export of many antibiotics, but the energy cost may decrease fitness when antibiotics are absent. Our evolution experiment reveals genetic mechanisms of adaptation to the PMF uncoupler CCCP, including selection for increased CCCP efflux but also against the expression of PMF-driven pumps for drugs not present. The results have implications for our understanding of the gut microbiome, which experiences high levels of organic acids that decrease PMF.

Experimental Evolution of Escherichia coli K-12 at High pH and with RpoS Induction.

Experimental evolution of Escherichia coli K-12 W3110 by serial dilutions for 2,200 generations at high pH extended the range of sustained growth from pH 9.0 to pH 9.3. pH 9.3-adapted isolates showed mutations in DNA-binding regulators and envelope proteins. One population showed an IS1 knockout of phoB (encoding the positive regulator of the phosphate regulon). A phoB::kanR knockout increased growth at high pH. phoB mutants are known to increase production of fermentation acids, which could enhance fitness at high pH. Mutations in pcnB [poly(A) polymerase] also increased growth at high pH. Three out of four populations showed deletions of torI, an inhibitor of TorR, which activates expression of torCAD (trimethylamine N-oxide respiration) at high pH. All populations showed point mutations affecting the stationary-phase sigma factor RpoS, either in the coding gene or in genes for regulators of RpoS expression. RpoS is required for survival at extremely high pH. In our microplate assay, rpoS deletion slightly decreased growth at pH 9.1. RpoS protein accumulated faster at pH 9 than at pH 7. The RpoS accumulation at high pH required the presence of one or more antiadaptors that block degradation (IraM, IraD, and IraP). Other genes with mutations after high-pH evolution encode regulators, such as those encoded by yobG (mgrB) (PhoPQ regulator), rpoN (nitrogen starvation sigma factor), malI, and purR, as well as envelope proteins, such as those encoded by ompT and yahO Overall, E. coli evolution at high pH selects for mutations in key transcriptional regulators, including phoB and the stationary-phase sigma factor RpoS.IMPORTANCEEscherichia coli in its native habitat encounters high-pH stress such as that of pancreatic secretions. Experimental evolution over 2,000 generations showed selection for mutations in regulatory factors, such as deletion of the phosphate regulator PhoB and mutations that alter the function of the global stress regulator RpoS. RpoS is induced at high pH via multiple mechanisms.

Acid Evolution of Escherichia coli K-12 Eliminates Amino Acid Decarboxylases and Reregulates Catabolism.

Acid-adapted strains of Escherichia coli K-12 W3110 were obtained by serial culture in medium buffered at pH 4.6 (M. M. Harden, A. He, K. Creamer, M. W. Clark, I. Hamdallah, K. A. Martinez, R. L. Kresslein, S. P. Bush, and J. L. Slonczewski, Appl Environ Microbiol 81:1932-1941, 2015, https://doi.org/10.1128/AEM.03494-14). Revised genomic analysis of these strains revealed insertion sequence (IS)-driven insertions and deletions that knocked out regulators CadC (acid induction of lysine decarboxylase), GadX (acid induction of glutamate decarboxylase), and FNR (anaerobic regulator). Each acid-evolved strain showed loss of one or more amino acid decarboxylase systems, which normally help neutralize external acid (pH 5 to 6) and increase survival in extreme acid (pH 2). Strains from populations B11, H9, and F11 had an IS5 insertion or IS-mediated deletion in cadC, while population B11 had a point mutation affecting the arginine activator adiY The cadC and adiY mutants failed to neutralize acid in the presence of exogenous lysine or arginine. In strain B11-1, reversion of an rpoC (RNA polymerase) mutation partly restored arginine-dependent neutralization. All eight strains showed deletion or downregulation of the Gad acid fitness island. Strains with the Gad deletion lost the ability to produce GABA (gamma-aminobutyric acid) and failed to survive extreme acid. Transcriptome sequencing (RNA-seq) of strain B11-1 showed upregulated genes for catabolism of diverse substrates but downregulated acid stress genes (the biofilm regulator ariR, yhiM, and Gad). Other strains showed downregulation of H2 consumption mediated by hydrogenases (hya and hyb) which release acid. Strains F9-2 and F9-3 had a deletion of fnr and showed downregulation of FNR-dependent genes (dmsABC, frdABCD, hybABO, nikABCDE, and nrfAC). Overall, strains that had evolved in buffered acid showed loss or downregulation of systems that neutralize unbuffered acid and showed altered regulation of catabolism.IMPORTANCE Experimental evolution of an enteric bacterium under a narrow buffered range of acid pH leads to loss of genes that enhance fitness above or below the buffered pH range, including loss of enzymes that may raise external pH in the absence of buffer. Prominent modes of evolutionary change involve IS-mediated insertions and deletions that knock out key regulators. Over generations of acid stress, catabolism undergoes reregulation in ways that differ for each evolving strain.

Benzoate- and Salicylate-Tolerant Strains of Escherichia coli K-12 Lose Antibiotic Resistance during Laboratory Evolution.

Escherichia coli K-12 W3110 grows in the presence of membrane-permeant organic acids that can depress cytoplasmic pH and accumulate in the cytoplasm. We conducted experimental evolution by daily diluting cultures in increasing concentrations of benzoic acid (up to 20 mM) buffered at external pH 6.5, a pH at which permeant acids concentrate in the cytoplasm. By 2,000 generations, clones isolated from evolving populations showed increasing tolerance to benzoate but were sensitive to chloramphenicol and tetracycline. Sixteen clones grew to stationary phase in 20 mM benzoate, whereas the ancestral strain W3110 peaked and declined. Similar growth occurred in 10 mM salicylate. Benzoate-evolved strains grew like W3110 in the absence of benzoate, in media buffered at pH 4.8, pH 7.0, or pH 9.0, or in 20 mM acetate or sorbate at pH 6.5. Genomes of 16 strains revealed over 100 mutations, including single-nucleotide polymorphisms (SNPs), large deletions, and insertion knockouts. Most strains acquired deletions in the benzoate-induced multiple antibiotic resistance (Mar) regulon or in associated regulators such as rob and cpxA, as well as the multidrug resistance (MDR) efflux pumps emrA, emrY, and mdtA Strains also lost or downregulated the Gad acid fitness regulon. In 5 mM benzoate or in 2 mM salicylate (2-hydroxybenzoate), most strains showed increased sensitivity to the antibiotics chloramphenicol and tetracycline; some strains were more sensitive than a marA knockout strain. Thus, our benzoate-evolved strains may reveal additional unknown drug resistance components. Benzoate or salicylate selection pressure may cause general loss of MDR genes and regulators. IMPORTANCE: Benzoate is a common food preservative, and salicylate is the primary active metabolite of aspirin. In the gut microbiome, genetic adaptation to salicylate may involve loss or downregulation of inducible multidrug resistance systems. This discovery implies that aspirin therapy may modulate the human gut microbiome to favor salicylate tolerance at the expense of drug resistance. Similar aspirin-associated loss of drug resistance might occur in bacterial pathogens found in arterial plaques.

Periplasmic Acid Stress Increases Cell Division Asymmetry (Polar Aging) of Escherichia coli.

Under certain kinds of cytoplasmic stress, Escherichia coli selectively reproduce by distributing the newer cytoplasmic components to new-pole cells while sequestering older, damaged components in cells inheriting the old pole. This phenomenon is termed polar aging or cell division asymmetry. It is unknown whether cell division asymmetry can arise from a periplasmic stress, such as the stress of extracellular acid, which is mediated by the periplasm. We tested the effect of periplasmic acid stress on growth and division of adherent single cells. We tracked individual cell lineages over five or more generations, using fluorescence microscopy with ratiometric pHluorin to measure cytoplasmic pH. Adherent colonies were perfused continually with LBK medium buffered at pH 6.00 or at pH 7.50; the external pH determines periplasmic pH. In each experiment, cell lineages were mapped to correlate division time, pole age and cell generation number. In colonies perfused at pH 6.0, the cells inheriting the oldest pole divided significantly more slowly than the cells inheriting the newest pole. In colonies perfused at pH 7.50 (near or above cytoplasmic pH), no significant cell division asymmetry was observed. Under both conditions (periplasmic pH 6.0 or pH 7.5) the cells maintained cytoplasmic pH values at 7.2-7.3. No evidence of cytoplasmic protein aggregation was seen. Thus, periplasmic acid stress leads to cell division asymmetry with minimal cytoplasmic stress.