Enricherator: A Bayesian Method for Inferring Regularized Genome-wide Enrichments from Sequencing Count Data.

A pervasive question in biological research studying gene regulation, chromatin structure, or genomics is where, and to what extent, does a signal of interest arise genome-wide? This question is addressed using a variety of methods relying on high-throughput sequencing data as their final output, including ChIP-seq for protein-DNA interactions,1 GapR-seq for measuring supercoiling,2 and HBD-seq or DRIP-seq for R-loop positioning.3,4 Current computational methods to calculate genome-wide enrichment of the signal of interest usually do not properly handle the count-based nature of sequencing data, they often do not make use of the local correlation structure of sequencing data, and they do not apply any regularization of enrichment estimates. This can result in unrealistic estimates of the true underlying biological enrichment of interest, unrealistically low estimates of confidence in point estimates of enrichment (or no estimates of confidence at all), unrealistic gyrations in enrichment estimates at very close (<10 bp) genomic loci due to noise inherent in sequencing data, and in a multiple-hypothesis testing problem during interpretation of genome-wide enrichment estimates. We developed a tool called Enricherator to infer genome-wide enrichments from sequencing count data. Enricherator uses the variational Bayes algorithm to fit a generalized linear model to sequencing count data and to sample from the approximate posterior distribution of enrichment estimates (https://github.com/jwschroeder3/enricherator). Enrichments inferred by Enricherator more precisely identify known binding sites in cases where low coverage between binding sites leads to false-positive peak calls in these noisy regions of the genome; these benefits extend to published datasets.

Tools and methods for high-throughput single-cell imaging with the mother machine.

Despite much progress, image processing remains a significant bottleneck for high-throughput analysis of microscopy data. One popular platform for single-cell time-lapse imaging is the mother machine, which enables long-term tracking of microbial cells under precisely controlled growth conditions. While several mother machine image analysis pipelines have been developed in the past several years, adoption by a non-expert audience remains a challenge. To fill this gap, we implemented our own software, MM3, as a plugin for the multidimensional image viewer napari. napari-MM3 is a complete and modular image analysis pipeline for mother machine data, which takes advantage of the high-level interactivity of napari. Here, we give an overview of napari-MM3 and test it against several well-designed and widely-used image analysis pipelines, including BACMMAN and DeLTA. Researchers often analyze mother machine data with custom scripts using varied image analysis methods, but a quantitative comparison of the output of different pipelines has been lacking. To this end, we show that key single-cell physiological parameter correlations and distributions are robust to the choice of analysis method. However, we also find that small changes in thresholding parameters can systematically alter parameters extracted from single-cell imaging experiments. Moreover, we explicitly show that in deep learning based segmentation, "what you put is what you get" (WYPIWYG) - i.e., pixel-level variation in training data for cell segmentation can propagate to the model output and bias spatial and temporal measurements. Finally, while the primary purpose of this work is to introduce the image analysis software that we have developed over the last decade in our lab, we also provide information for those who want to implement mother-machine-based high-throughput imaging and analysis methods in their research.

Deep mutational scanning highlights a role for cytosolic regions in Hrd1 function.

Misfolded endoplasmic reticulum (ER) proteins are degraded through a process called ER-associated degradation (ERAD). Soluble, lumenal ERAD targets are recognized, retrotranslocated across the ER membrane, ubiquitinated, extracted from the membrane, and degraded by the proteasome using an ERAD pathway containing a ubiquitin ligase called Hrd1. To determine how Hrd1 mediates these processes, we developed a deep mutational scanning approach to identify residues involved in Hrd1 function, including those exclusively required for lumenal degradation. We identify several regions required for different Hrd1 functions. Most surprisingly, we find two cytosolic regions of Hrd1 required for lumenal ERAD substrate degradation. Using in vivo and in vitro approaches, we define roles for disordered regions between structural elements that are required for Hrd1 autoubiquitination and substrate interaction. Our results demonstrate that disordered cytosolic regions promote substrate retrotranslocation by controlling Hrd1 activation and establishing directionality of retrotranslocation for lumenal substrate across the ER membrane.

RNase H genes cause distinct impacts on RNA:DNA hybrid formation and mutagenesis genome wide.

RNA:DNA hybrids compromise replication fork progression and genome integrity in all cells. The overall impacts of naturally occurring RNA:DNA hybrids on genome integrity, and the relative contributions of ribonucleases H to mitigating the negative effects of hybrids, remain unknown. Here, we investigate the contributions of RNases HII (RnhB) and HIII (RnhC) to hybrid removal, DNA replication, and mutagenesis genome wide. Deletion of either rnhB or rnhC triggers RNA:DNA hybrid accumulation but with distinct patterns of mutagenesis and hybrid accumulation. Across all cells, hybrids accumulate strongly in noncoding RNAs and 5'-UTRs of coding sequences. For ΔrnhB, hybrids accumulate preferentially in untranslated regions and early in coding sequences. We show that hybrid accumulation is particularly sensitive to gene expression in ΔrnhC cells. DNA replication in ΔrnhC cells is disrupted, leading to transversions and structural variation. Our results resolve the outstanding question of how hybrids in native genomic contexts cause mutagenesis and shape genome organization.

RNase H genes cause distinct impacts on RNA:DNA hybrid formation and mutagenesis genome-wide.

RNA:DNA hybrids such as R-loops affect genome integrity and DNA replication fork progression. The overall impacts of naturally occurring RNA:DNA hybrids on genome integrity, and the relative contributions of ribonucleases H to mitigating the negative effects of hybrids, remain unknown. Here, we investigate the contributions of RNases HII (RnhB) and HIII (RnhC) to hybrid removal, DNA replication, and mutagenesis genome-wide. Deletion of either rnhB or rnhC triggers RNA:DNA hybrid accumulation, but with distinct patterns of mutagenesis and hybrid accumulation. Across all cells, hybrids accumulate most strongly in non-coding RNAs and 5'-UTRs of coding sequences. For Δ rnhB , hybrids accumulate preferentially in untranslated regions and early in coding sequences. Hybrid accumulation is particularly sensitive to gene expression in Δ rnhC ; in cells lacking RnhC, DNA replication is disrupted leading to transversions and structural variation. Our results resolve the outstanding question of how hybrids in native genomic contexts interact with replication to cause mutagenesis and shape genome organization.

The nucleotide messenger (p)ppGpp is an anti-inducer of the purine synthesis transcription regulator PurR in Bacillus.

The nucleotide messenger (p)ppGpp allows bacteria to adapt to fluctuating environments by reprogramming the transcriptome. Despite its well-recognized role in gene regulation, (p)ppGpp is only known to directly affect transcription in Proteobacteria by binding to the RNA polymerase. Here, we reveal a different mechanism of gene regulation by (p)ppGpp in Firmicutes: (p)ppGpp directly binds to the transcription factor PurR to downregulate purine biosynthesis gene expression upon amino acid starvation. We first identified PurR as a receptor of (p)ppGpp in Bacillus anthracis. A co-structure with Bacillus subtilis PurR reveals that (p)ppGpp binds to a PurR pocket reminiscent of the active site of phosphoribosyltransferase enzymes that has been repurposed to serve a purely regulatory role, where the effectors (p)ppGpp and PRPP compete to allosterically control transcription. PRPP inhibits PurR DNA binding to induce transcription of purine synthesis genes, whereas (p)ppGpp antagonizes PRPP to enhance PurR DNA binding and repress transcription. A (p)ppGpp-refractory purR mutant in B. subtilis fails to downregulate purine synthesis genes upon amino acid starvation. Our work establishes the precedent of (p)ppGpp as an effector of a classical transcription repressor and reveals the key function of (p)ppGpp in regulating nucleotide synthesis through gene regulation, from soil bacteria to pathogens.

The roles of replication-transcription conflict in mutagenesis and evolution of genome organization.

Replication-transcription conflicts promote mutagenesis and give rise to evolutionary signatures, with fundamental importance to genome stability ranging from bacteria to metastatic cancer cells. This review focuses on the interplay between replication-transcription conflicts and the evolution of gene directionality. In most bacteria, the majority of genes are encoded on the leading strand of replication such that their transcription is co-directional with the direction of DNA replication fork movement. This gene strand bias arises primarily due to negative selection against deleterious consequences of head-on replication-transcription conflict. However, many genes remain head-on. Can head-on orientation provide some benefit? We combine insights from both mechanistic and evolutionary studies, review published work, and analyze gene expression data to evaluate an emerging model that head-on genes are temporal targets for adaptive mutagenesis during stress. We highlight the alternative explanation that genes in the head-on orientation may simply be the result of genomic inversions and relaxed selection acting on nonessential genes. We seek to clarify how the mechanisms of replication-transcription conflict, in concert with other mutagenic mechanisms, balanced by natural selection, have shaped bacterial genome evolution.

Methyltransferase DnmA is responsible for genome-wide N6-methyladenosine modifications at non-palindromic recognition sites in Bacillus subtilis.

The genomes of organisms from all three domains of life harbor endogenous base modifications in the form of DNA methylation. In bacterial genomes, methylation occurs on adenosine and cytidine residues to include N6-methyladenine (m6A), 5-methylcytosine (m5C), and N4-methylcytosine (m4C). Bacterial DNA methylation has been well characterized in the context of restriction-modification (RM) systems, where methylation regulates DNA incision by the cognate restriction endonuclease. Relative to RM systems less is known about how m6A contributes to the epigenetic regulation of cellular functions in Gram-positive bacteria. Here, we characterize site-specific m6A modifications in the non-palindromic sequence GACGmAG within the genomes of Bacillus subtilis strains. We demonstrate that the yeeA gene is a methyltransferase responsible for the presence of m6A modifications. We show that methylation from YeeA does not function to limit DNA uptake during natural transformation. Instead, we identify a subset of promoters that contain the methylation consensus sequence and show that loss of methylation within promoter regions causes a decrease in reporter expression. Further, we identify a transcriptional repressor that preferentially binds an unmethylated promoter used in the reporter assays. With these results we suggest that m6A modifications in B. subtilis function to promote gene expression.

Metabolic Remodeling during Biofilm Development of Bacillus subtilis.

Biofilms are structured communities of tightly associated cells that constitute the predominant state of bacterial growth in natural and human-made environments. Although the core genetic circuitry that controls biofilm formation in model bacteria such as Bacillus subtilis has been well characterized, little is known about the role that metabolism plays in this complex developmental process. Here, we performed a time-resolved analysis of the metabolic changes associated with pellicle biofilm formation and development in B. subtilis by combining metabolomic, transcriptomic, and proteomic analyses. We report surprisingly widespread and dynamic remodeling of metabolism affecting central carbon metabolism, primary biosynthetic pathways, fermentation pathways, and secondary metabolism. Most of these metabolic alterations were hitherto unrecognized as biofilm associated. For example, we observed increased activity of the tricarboxylic acid (TCA) cycle during early biofilm growth, a shift from fatty acid biosynthesis to fatty acid degradation, reorganization of iron metabolism and transport, and a switch from acetate to acetoin fermentation. Close agreement between metabolomic, transcriptomic, and proteomic measurements indicated that remodeling of metabolism during biofilm development was largely controlled at the transcriptional level. Our results also provide insights into the transcription factors and regulatory networks involved in this complex metabolic remodeling. Following upon these results, we demonstrated that acetoin production via acetolactate synthase is essential for robust biofilm growth and has the dual role of conserving redox balance and maintaining extracellular pH. This report represents a comprehensive systems-level investigation of the metabolic remodeling occurring during B. subtilis biofilm development that will serve as a useful road map for future studies on biofilm physiology.IMPORTANCE Bacterial biofilms are ubiquitous in natural environments and play an important role in many clinical, industrial, and ecological settings. Although much is known about the transcriptional regulatory networks that control biofilm formation in model bacteria such as Bacillus subtilis, very little is known about the role of metabolism in this complex developmental process. To address this important knowledge gap, we performed a time-resolved analysis of the metabolic changes associated with bacterial biofilm development in B. subtilis by combining metabolomic, transcriptomic, and proteomic analyses. Here, we report a widespread and dynamic remodeling of metabolism affecting central carbon metabolism, primary biosynthetic pathways, fermentation pathways, and secondary metabolism. This report serves as a unique hypothesis-generating resource for future studies on bacterial biofilm physiology. Outside the biofilm research area, this work should also prove relevant to any investigators interested in microbial physiology and metabolism.

Discovery of a dual protease mechanism that promotes DNA damage checkpoint recovery.

The DNA damage response is a signaling pathway found throughout biology. In many bacteria the DNA damage checkpoint is enforced by inducing expression of a small, membrane bound inhibitor that delays cell division providing time to repair damaged chromosomes. How cells promote checkpoint recovery after sensing successful repair is unknown. By using a high-throughput, forward genetic screen, we identified two unrelated proteases, YlbL and CtpA, that promote DNA damage checkpoint recovery in Bacillus subtilis. Deletion of both proteases leads to accumulation of the checkpoint protein YneA. We show that DNA damage sensitivity and increased cell elongation in protease mutants depends on yneA. Further, expression of YneA in protease mutants was sufficient to inhibit cell proliferation. Finally, we show that both proteases interact with YneA and that one of the two proteases, CtpA, directly cleaves YneA in vitro. With these results, we report the mechanism for DNA damage checkpoint recovery in bacteria that use membrane bound cell division inhibitors.

Visualizing bacterial DNA replication and repair with molecular resolution.

Although DNA replication and repair in bacteria have been extensively studied for many decades, in recent years the development of single-molecule microscopy has provided a new perspective on these fundamental processes. Because single-molecule imaging super-resolves the nanometer-scale dynamics of molecules, and because single-molecule imaging is sensitive to heterogeneities within a sample, this nanoscopic microscopy technique measures the motions, localizations, and interactions of proteins in real time without averaging ensemble observations, both in vitro and in vivo. In this Review, we provide an overview of several recent single-molecule fluorescence microscopy studies on DNA replication and repair. These experiments have shown that, in both Escherichia coli and Bacillus subtilis the DNA replication proteins are highly dynamic. In particular, even highly processive replicative DNA polymerases exchange to and from the replication fork on the scale of a few seconds. Furthermore, single-molecule investigations of the DNA mismatch repair (MMR) pathway have measured the complex interactions between MMR proteins, replication proteins, and DNA. Single-molecule imaging will continue to improve our understanding of fundamental processes in bacteria including DNA replication and repair.

Sources of spontaneous mutagenesis in bacteria.

Mutations in an organism's genome can arise spontaneously, that is, in the absence of exogenous stress and prior to selection. Mutations are often neutral or deleterious to individual fitness but can also provide genetic diversity driving evolution. Mutagenesis in bacteria contributes to the already serious and growing problem of antibiotic resistance. However, the negative impacts of spontaneous mutagenesis on human health are not limited to bacterial antibiotic resistance. Spontaneous mutations also underlie tumorigenesis and evolution of drug resistance. To better understand the causes of genetic change and how they may be manipulated in order to curb antibiotic resistance or the development of cancer, we must acquire a mechanistic understanding of the major sources of mutagenesis. Bacterial systems are particularly well-suited to studying mutagenesis because of their fast growth rate and the panoply of available experimental tools, but efforts to understand mutagenic mechanisms can be complicated by the experimental system employed. Here, we review our current understanding of mutagenic mechanisms in bacteria and describe the methods used to study mutagenesis in bacterial systems.

Mutagenic cost of ribonucleotides in bacterial DNA.

Replicative DNA polymerases misincorporate ribonucleoside triphosphates (rNTPs) into DNA approximately once every 2,000 base pairs synthesized. Ribonucleotide excision repair (RER) removes ribonucleoside monophosphates (rNMPs) from genomic DNA, replacing the error with the appropriate deoxyribonucleoside triphosphate (dNTP). Ribonucleotides represent a major threat to genome integrity with the potential to cause strand breaks. Furthermore, it has been shown in the bacterium Bacillus subtilis that loss of RER increases spontaneous mutagenesis. Despite the high rNTP error rate and the effect on genome integrity, the mechanism underlying mutagenesis in RER-deficient bacterial cells remains unknown. We performed mutation accumulation lines and genome-wide mutational profiling of B. subtilis lacking RNase HII, the enzyme that incises at single rNMP residues initiating RER. We show that loss of RER in B. subtilis causes strand- and sequence-context-dependent GC → AT transitions. Using purified proteins, we show that the replicative polymerase DnaE is mutagenic within the sequence context identified in RER-deficient cells. We also found that DnaE does not perform strand displacement synthesis. Given the use of nucleotide excision repair (NER) as a backup pathway for RER in RNase HII-deficient cells and the known mutagenic profile of DnaE, we propose that misincorporated ribonucleotides are removed by NER followed by error-prone resynthesis with DnaE.

Pseudomonas aeruginosa defense systems against microbicidal oxidants.

The most abundant oxidants controlling bacterial colonization on mucosal barrier epithelia are hypochlorous acid (HOCl), hypobromous acid (HOBr) and hypothiocyanous acid (HOSCN). All three oxidants are highly antimicrobial but little is known about their relative efficacies, their respective cellular targets, or what specific responses they elicit in bacteria. To address these important questions, we directly tested the individual oxidants on the virulent Pseudomonas aeruginosa strain PA14. We discovered that HOCl and HOBr work almost interchangeably, impacting non-growing bacterial cultures more significantly than actively growing bacteria, and eliciting similar stress responses, including the heat shock response. HOSCN treatment is distinctly different, affecting primarily actively growing PA14 and evoking stress responses suggestive of membrane damage. What all three oxidants have in common, however, is their ability to cause substantial protein aggregation. This effect became particularly obvious in strains lacking polyphosphate, a newly recognized chemical chaperone. Treatment of PA14 with the FDA-approved anti-inflammatory drug mesalamine, which has recently been shown to attenuate polyP production in a wide range of bacteria, effectively decreased the resistance of PA14 toward all three oxidants, suggesting that we have discovered a novel, targetable defense system in P. aeruginosa.

Complete Genome Sequence of Undomesticated Bacillus subtilis Strain NCIB 3610.

Bacillus subtilis is a Gram-positive bacterium that serves as an important experimental system. B. subtilis NCIB 3610 is an undomesticated strain that exhibits phenotypes lost from the more common domesticated laboratory strains. Here, we announce the complete genome sequence of DK1042, a genetically competent derivative of NCIB 3610.

Implementation and Data Analysis of Tn-seq, Whole-Genome Resequencing, and Single-Molecule Real-Time Sequencing for Bacterial Genetics.

Few discoveries have been more transformative to the biological sciences than the development of DNA sequencing technologies. The rapid advancement of sequencing and bioinformatics tools has revolutionized bacterial genetics, deepening our understanding of model and clinically relevant organisms. Although application of newer sequencing technologies to studies in bacterial genetics is increasing, the implementation of DNA sequencing technologies and development of the bioinformatics tools required for analyzing the large data sets generated remain a challenge for many. In this minireview, we have chosen to summarize three sequencing approaches that are particularly useful for bacterial genetics. We provide resources for scientists new to and interested in their application. Here, we discuss the analysis of data from transposon mutagenesis followed by deep sequencing (Tn-seq) to determine gene disruptions differentially represented in a mutant population and Illumina sequencing for identification of suppressor or other mutations, and we summarize single-molecule real-time (SMRT) sequencing for de novo genome assembly and the use of the output data for detection of DNA base modifications.

Single-Molecule DNA Polymerase Dynamics at a Bacterial Replisome in Live Cells.

PolC is one of two essential replicative DNA polymerases found in the Gram-positive bacterium Bacillus subtilis. The B. subtilis replisome is eukaryotic-like in that it relies on a two DNA polymerase system for chromosomal replication. To quantitatively image how the replicative DNA polymerase PolC functions in B. subtilis, we applied photobleaching-assisted microscopy, three-dimensional superresolution imaging, and single-particle tracking to examine the in vivo behavior of PolC at single-molecule resolution. We report the stoichiometry of PolC proteins within each cell and within each replisome, we elucidate the diffusion characteristics of individual PolC molecules, and we quantify the exchange dynamics for PolC engaged in lagging strand synthesis. We show that PolC is highly dynamic: this DNA polymerase is constantly recruited to and released from a centrally located replisome, providing, to our knowledge, new insight into the organization and dynamics of the replisome in bacterial cells.

Complex Haploinsufficiency-Based Genetic Analysis of the NDR/Lats Kinase Cbk1 Provides Insight into Its Multiple Functions in Candida albicans.

Although the analysis of genetic interactions and networks is a powerful approach to understanding biology, it has not been applied widely to the pathogenic yeast Candida albicans Here, we describe the use of both screening and directed genetic interaction studies based on complex haploinsufficiency to probe the function of the R: egulation of A: ce2 and M: orphogenesis (RAM) pathway in C. albicans A library of 5200 Tn7-mutagenized derivatives of a parental strain heterozygous at CBK1, the key kinase in the RAM pathway, was screened for alterations in serum-induced filamentation. Following confirmation of phenotypes and identification of insertion sites by sequencing, a set of 36 unique double heterozygous strains showing complex haploinsufficiency was obtained. In addition to a large set of genes regulated by the RAM transcription factor Ace2, genes related to cell wall biosynthesis, cell cycle, polarity, oxidative stress, and nitrogen utilization were identified. Follow-up analysis led to the first demonstration that the RAM pathway is required for oxidative stress tolerance in a manner related to the two-component-regulated kinase Chk1 and revealed a potential direct connection between the RAM pathway and the essential Mps1 spindle pole-related kinase. In addition, genetic interactions with CDC42-related genes MSB1, a putative scaffold protein, and RGD3, a putative Rho GTPase-activating protein (GAP) were identified. We also provide evidence that Rgd3 is a GAP for Cdc42 and show that its localization and phosphorylation are dependent on Cbk1.

The Effect of Local Sequence Context on Mutational Bias of Genes Encoded on the Leading and Lagging Strands.

All organisms must replicate their genetic information accurately to ensure its faithful transmission. DNA polymerase errors provide an important source of genetic variation that can drive evolution. Understanding the origins of genetic variation will inform our understanding of evolution and the development of genetic diseases. A number of factors have been proposed to influence mutagenesis [1-10]. Here, we used mutation accumulation lines, whole-genome sequencing, and whole-transcriptome analysis to study the locations and rate at which mutations arise in bacteria with as little selection bias as possible [11, 12]. Our analysis of greater than 7,000 replication errors in over 180 sequenced lines that underwent a total of more than 370,000 generations has provided new insights into how DNA polymerase errors sculpt genetic variation and drive evolution. Homopolymer run enrichment outside of genes causes insertions and deletions in these regions. Genes encoded in the lagging strand are transcribed such that RNA polymerase and DNA polymerase collide head-on. Head-on genes have been proposed to mutate at a higher rate than genes transcribed codirectionally with DNA polymerase progression due to conflicts between transcription and DNA replication [6, 10]. We did not detect associations between the number of base pair substitutions in genes and their orientation or expression. Strikingly, any higher mutation rate for head-on genes can be explained by differing sequence composition between the leading and lagging strands and the error bias for DNA polymerase in specific sequence contexts. Therefore, we find local sequence context is the major determinant of mutagenesis in bacteria.

Single-molecule motions and interactions in live cells reveal target search dynamics in mismatch repair.

MutS is responsible for initiating the correction of DNA replication errors. To understand how MutS searches for and identifies rare base-pair mismatches, we characterized the dynamic movement of MutS and the replisome in real time using superresolution microscopy and single-molecule tracking in living cells. We report that MutS dynamics are heterogeneous in cells, with one MutS population exploring the nucleoid rapidly, while another MutS population moves to and transiently dwells at the replisome region, even in the absence of appreciable mismatch formation. Analysis of MutS motion shows that the speed of MutS is correlated with its separation distance from the replisome and that MutS motion slows when it enters the replisome region. We also show that mismatch detection increases MutS speed, supporting the model for MutS sliding clamp formation after mismatch recognition. Using variants of MutS and the replication processivity clamp to impair mismatch repair, we find that MutS dynamically moves to and from the replisome before mismatch binding to scan for errors. Furthermore, a block to DNA synthesis shows that MutS is only capable of binding mismatches near the replisome. It is well-established that MutS engages in an ATPase cycle, which is necessary for signaling downstream events. We show that a variant of MutS with a nucleotide binding defect is no longer capable of dynamic movement to and from the replisome, showing that proper nucleotide binding is critical for MutS to localize to the replisome in vivo. Our results provide mechanistic insight into the trafficking and movement of MutS in live cells as it searches for mismatches.