TbsP and TrmB jointly regulate gapII to influence cell development phenotypes in the archaeon Haloferax volcanii.

Microbial cells must continually adapt their physiology in the face of changing environmental conditions. Archaea living in extreme conditions, such as saturated salinity, represent important examples of such resilience. The model salt-loving organism Haloferax volcanii exhibits remarkable plasticity in its morphology, biofilm formation, and motility in response to variations in nutrients and cell density. However, the mechanisms regulating these lifestyle transitions remain unclear. In prior research, we showed that the transcriptional regulator, TrmB, maintains the rod shape in the related species Halobacterium salinarum by activating the expression of enzyme-coding genes in the gluconeogenesis metabolic pathway. In Hbt. salinarum, TrmB-dependent production of glucose moieties is required for cell surface glycoprotein biogenesis. Here, we use a combination of genetics and quantitative phenotyping assays to demonstrate that TrmB is essential for growth under gluconeogenic conditions in Hfx. volcanii. The ∆trmB strain rapidly accumulated suppressor mutations in a gene encoding a novel transcriptional regulator, which we name trmB suppressor, or TbsP (a.k.a. "tablespoon"). TbsP is required for adhesion to abiotic surfaces (i.e., biofilm formation) and maintains wild-type cell morphology and motility. We use functional genomics and promoter fusion assays to characterize the regulons controlled by each of TrmB and TbsP, including joint regulation of the glucose-dependent transcription of gapII, which encodes an important gluconeogenic enzyme. We conclude that TrmB and TbsP coregulate gluconeogenesis, with downstream impacts on lifestyle transitions in response to nutrients in Hfx. volcanii.

A conserved transcription factor controls gluconeogenesis via distinct targets in hypersaline-adapted archaea with diverse metabolic capabilities.

Timely regulation of carbon metabolic pathways is essential for cellular processes and to prevent futile cycling of intracellular metabolites. In Halobacterium salinarum, a hypersaline adapted archaeon, a sugar-sensing TrmB family protein controls gluconeogenesis and other biosynthetic pathways. Notably, Hbt. salinarum does not utilize carbohydrates for energy, uncommon among Haloarchaea. We characterized a TrmB-family transcriptional regulator in a saccharolytic generalist, Haloarcula hispanica, to investigate whether the targets and function of TrmB, or its regulon, is conserved in related species with distinct metabolic capabilities. In Har. hispanica, TrmB binds to 15 sites in the genome and induces the expression of genes primarily involved in gluconeogenesis and tryptophan biosynthesis. An important regulatory control point in Hbt. salinarum, activation of ppsA and repression of pykA, is absent in Har. hispanica. Contrary to its role in Hbt. salinarum and saccharolytic hyperthermophiles, TrmB does not act as a global regulator: it does not directly repress the expression of glycolytic enzymes, peripheral pathways such as cofactor biosynthesis, or catabolism of other carbon sources in Har. hispanica. Cumulatively, these findings suggest rewiring of the TrmB regulon alongside metabolic network evolution in Haloarchaea.

The Hypersaline Archaeal Histones HpyA and HstA Are DNA Binding Proteins That Defy Categorization According to Commonly Used Functional Criteria.

Histone proteins are found across diverse lineages of Archaea, many of which package DNA and form chromatin. However, previous research has led to the hypothesis that the histone-like proteins of high-salt-adapted archaea, or halophiles, function differently. The sole histone protein encoded by the model halophilic species Halobacterium salinarum, HpyA, is nonessential and expressed at levels too low to enable genome-wide DNA packaging. Instead, HpyA mediates the transcriptional response to salt stress. Here we compare the features of genome-wide binding of HpyA to those of HstA, the sole histone of another model halophile, Haloferax volcanii. hstA, like hpyA, is a nonessential gene. To better understand HpyA and HstA functions, protein-DNA binding data (chromatin immunoprecipitation sequencing [ChIP-seq]) of these halophilic histones are compared to publicly available ChIP-seq data from DNA binding proteins across all domains of life, including transcription factors (TFs), nucleoid-associated proteins (NAPs), and histones. These analyses demonstrate that HpyA and HstA bind the genome infrequently in discrete regions, which is similar to TFs but unlike NAPs, which bind a much larger genomic fraction. However, unlike TFs that typically bind in intergenic regions, HpyA and HstA binding sites are located in both coding and intergenic regions. The genome-wide dinucleotide periodicity known to facilitate histone binding was undetectable in the genomes of both species. Instead, TF-like and histone-like binding sequence preferences were detected for HstA and HpyA, respectively. Taken together, these data suggest that halophilic archaeal histones are unlikely to facilitate genome-wide chromatin formation and that their function defies categorization as a TF, NAP, or histone. IMPORTANCE Most cells in eukaryotic species-from yeast to humans-possess histone proteins that pack and unpack DNA in response to environmental cues. These essential proteins regulate genes necessary for important cellular processes, including development and stress protection. Although the histone fold domain originated in the domain of life Archaea, the function of archaeal histone-like proteins is not well understood relative to those of eukaryotes. We recently discovered that, unlike histones of eukaryotes, histones in hypersaline-adapted archaeal species do not package DNA and can act as transcription factors (TFs) to regulate stress response gene expression. However, the function of histones across species of hypersaline-adapted archaea still remains unclear. Here, we compare hypersaline histone function to a variety of DNA binding proteins across the tree of life, revealing histone-like behavior in some respects and specific transcriptional regulatory function in others.

A Bayesian non-parametric mixed-effects model of microbial growth curves.

Substantive changes in gene expression, metabolism, and the proteome are manifested in overall changes in microbial population growth. Quantifying how microbes grow is therefore fundamental to areas such as genetics, bioengineering, and food safety. Traditional parametric growth curve models capture the population growth behavior through a set of summarizing parameters. However, estimation of these parameters from data is confounded by random effects such as experimental variability, batch effects or differences in experimental material. A systematic statistical method to identify and correct for such confounding effects in population growth data is not currently available. Further, our previous work has demonstrated that parametric models are insufficient to explain and predict microbial response under non-standard growth conditions. Here we develop a hierarchical Bayesian non-parametric model of population growth that identifies the latent growth behavior and response to perturbation, while simultaneously correcting for random effects in the data. This model enables more accurate estimates of the biological effect of interest, while better accounting for the uncertainty due to technical variation. Additionally, modeling hierarchical variation provides estimates of the relative impact of various confounding effects on measured population growth.

The Ribbon-Helix-Helix Domain Protein CdrS Regulates the Tubulin Homolog ftsZ2 To Control Cell Division in Archaea.

Precise control of the cell cycle is central to the physiology of all cells. In prior work we demonstrated that archaeal cells maintain a constant size; however, the regulatory mechanisms underlying the cell cycle remain unexplored in this domain of life. Here, we use genetics, functional genomics, and quantitative imaging to identify and characterize the novel CdrSL gene regulatory network in a model species of archaea. We demonstrate the central role of these ribbon-helix-helix family transcription factors in the regulation of cell division through specific transcriptional control of the gene encoding FtsZ2, a putative tubulin homolog. Using time-lapse fluorescence microscopy in live cells cultivated in microfluidics devices, we further demonstrate that FtsZ2 is required for cell division but not elongation. The cdrS-ftsZ2 locus is highly conserved throughout the archaeal domain, and the central function of CdrS in regulating cell division is conserved across hypersaline adapted archaea. We propose that the CdrSL-FtsZ2 transcriptional network coordinates cell division timing with cell growth in archaea.IMPORTANCE Healthy cell growth and division are critical for individual organism survival and species long-term viability. However, it remains unknown how cells of the domain Archaea maintain a healthy cell cycle. Understanding the archaeal cell cycle is of paramount evolutionary importance given that an archaeal cell was the host of the endosymbiotic event that gave rise to eukaryotes. Here, we identify and characterize novel molecular players needed for regulating cell division in archaea. These molecules dictate the timing of cell septation but are dispensable for growth between divisions. Timing is accomplished through transcriptional control of the cell division ring. Our results shed light on mechanisms underlying the archaeal cell cycle, which has thus far remained elusive.

N-Glycosylation Is Important for Halobacterium salinarum Archaellin Expression, Archaellum Assembly and Cell Motility.

Halobacterium salinarum are halophilic archaea that display directional swimming in response to various environmental signals, including light, chemicals and oxygen. In Hbt. salinarum, the building blocks (archaellins) of the archaeal swimming apparatus (the archaellum) are N-glycosylated. However, the physiological importance of archaellin N-glycosylation remains unclear. Here, a tetrasaccharide comprising a hexose and three hexuronic acids decorating the five archaellins was characterized by mass spectrometry. Such analysis failed to detect sulfation of the hexuronic acids, in contrast to earlier reports. To better understand the physiological significance of Hbt. salinarum archaellin N-glycosylation, a strain deleted of aglB, encoding the archaeal oligosaccharyltransferase, was generated. In this ΔaglB strain, archaella were not detected and only low levels of archaellins were released into the medium, in contrast to what occurs with the parent strain. Mass spectrometry analysis of the archaellins in ΔaglB cultures did not detect N-glycosylation. ΔaglB cells also showed a slight growth defect and were impaired for motility. Quantitative real-time PCR analysis revealed dramatically reduced transcript levels of archaellin-encoding genes in the mutant strain, suggesting that N-glycosylation is important for archaellin transcription, with downstream effects on archaellum assembly and function. Control of AglB-dependent post-translational modification of archaellins could thus reflect a previously unrecognized route for regulating Hbt. salinarum motility.

Copy number variation is associated with gene expression change in archaea.

Genomic instability, although frequently deleterious, is also an important mechanism for microbial adaptation to environmental change. Although widely studied in bacteria, in archaea the effect of genomic instability on organism phenotypes and fitness remains unclear. Here we use DNA segmentation methods to detect and quantify genome-wide copy number variation (CNV) in large compendia of high-throughput datasets in a model archaeal species, Halobacterium salinarum. CNV hotspots were identified throughout the genome. Some hotspots were strongly associated with changes in gene expression, suggesting a mechanism for phenotypic innovation. In contrast, CNV hotspots in other genomic loci left expression unchanged, suggesting buffering of certain phenotypes. The correspondence of CNVs with gene expression was validated with strain- and condition-matched transcriptomics and DNA quantification experiments at specific loci. Significant correlation of CNV hotspot locations with the positions of known insertion sequence (IS) elements suggested a mechanism for generating genomic instability. Given the efficient recombination capabilities in H. salinarum despite stability at the single nucleotide level, these results suggest that genomic plasticity mediated by IS element activity can provide a source of phenotypic innovation in extreme environments.

Transcriptional Regulation in Archaea: From Individual Genes to Global Regulatory Networks.

Archaea are major contributors to biogeochemical cycles, possess unique metabolic capabilities, and resist extreme stress. To regulate the expression of genes encoding these unique programs, archaeal cells use gene regulatory networks (GRNs) composed of transcription factor proteins and their target genes. Recent developments in genetics, genomics, and computational methods used with archaeal model organisms have enabled the mapping and prediction of global GRN structures. Experimental tests of these predictions have revealed the dynamical function of GRNs in response to environmental variation. Here, we review recent progress made in this area, from investigating the mechanisms of transcriptional regulation of individual genes to small-scale subnetworks and genome-wide global networks. At each level, archaeal GRNs consist of a hybrid of bacterial, eukaryotic, and uniquely archaeal mechanisms. We discuss this theme from the perspective of the role of individual transcription factors in genome-wide regulation, how these proteins interact to compile GRN topological structures, and how these topologies lead to emergent, high-level GRN functions. We conclude by discussing how systems biology approaches are a fruitful avenue for addressing remaining challenges, such as discovering gene function and the evolution of GRNs.

Systematic Discovery of Archaeal Transcription Factor Functions in Regulatory Networks through Quantitative Phenotyping Analysis.

Gene regulatory networks (GRNs) are critical for dynamic transcriptional responses to environmental stress. However, the mechanisms by which GRN regulation adjusts physiology to enable stress survival remain unclear. Here we investigate the functions of transcription factors (TFs) within the global GRN of the stress-tolerant archaeal microorganism Halobacterium salinarum. We measured growth phenotypes of a panel of TF deletion mutants in high temporal resolution under heat shock, oxidative stress, and low-salinity conditions. To quantitate the noncanonical functional forms of the growth trajectories observed for these mutants, we developed a novel modeling framework based on Gaussian process regression and functional analysis of variance (FANOVA). We employ unique statistical tests to determine the significance of differential growth relative to the growth of the control strain. This analysis recapitulated known TF functions, revealed novel functions, and identified surprising secondary functions for characterized TFs. Strikingly, we observed that the majority of the TFs studied were required for growth under multiple stress conditions, pinpointing regulatory connections between the conditions tested. Correlations between quantitative phenotype trajectories of mutants are predictive of TF-TF connections within the GRN. These phenotypes are strongly concordant with predictions from statistical GRN models inferred from gene expression data alone. With genome-wide and targeted data sets, we provide detailed functional validation of novel TFs required for extreme oxidative stress and heat shock survival. Together, results presented in this study suggest that many TFs function under multiple conditions, thereby revealing high interconnectivity within the GRN and identifying the specific TFs required for communication between networks responding to disparate stressors. IMPORTANCE To ensure survival in the face of stress, microorganisms employ inducible damage repair pathways regulated by extensive and complex gene networks. Many archaea, microorganisms of the third domain of life, persist under extremes of temperature, salinity, and pH and under other conditions. In order to understand the cause-effect relationships between the dynamic function of the stress network and ultimate physiological consequences, this study characterized the physiological role of nearly one-third of all regulatory proteins known as transcription factors (TFs) in an archaeal organism. Using a unique quantitative phenotyping approach, we discovered functions for many novel TFs and revealed important secondary functions for known TFs. Surprisingly, many TFs are required for resisting multiple stressors, suggesting cross-regulation of stress responses. Through extensive validation experiments, we map the physiological roles of these novel TFs in stress response back to their position in the regulatory network wiring. This study advances understanding of the mechanisms underlying how microorganisms resist extreme stress. Given the generality of the methods employed, we expect that this study will enable future studies on how regulatory networks adjust cellular physiology in a diversity of organisms.

Detecting differential growth of microbial populations with Gaussian process regression.

Microbial growth curves are used to study differential effects of media, genetics, and stress on microbial population growth. Consequently, many modeling frameworks exist to capture microbial population growth measurements. However, current models are designed to quantify growth under conditions for which growth has a specific functional form. Extensions to these models are required to quantify the effects of perturbations, which often exhibit nonstandard growth curves. Rather than assume specific functional forms for experimental perturbations, we developed a general and robust model of microbial population growth curves using Gaussian process (GP) regression. GP regression modeling of high-resolution time-series growth data enables accurate quantification of population growth and allows explicit control of effects from other covariates such as genetic background. This framework substantially outperforms commonly used microbial population growth models, particularly when modeling growth data from environmentally stressed populations. We apply the GP growth model and develop statistical tests to quantify the differential effects of environmental perturbations on microbial growth across a large compendium of genotypes in archaea and yeast. This method accurately identifies known transcriptional regulators and implicates novel regulators of growth under standard and stress conditions in the model archaeal organism Halobacterium salinarum For yeast, our method correctly identifies known phenotypes for a diversity of genetic backgrounds under cyclohexamide stress and also detects previously unidentified oxidative stress sensitivity across a subset of strains. Together, these results demonstrate that the GP models are interpretable, recapitulating biological knowledge of growth response while providing new insights into the relevant parameters affecting microbial population growth.

Systems biology approaches to defining transcription regulatory networks in halophilic archaea.

To survive complex and changing environmental conditions, microorganisms use gene regulatory networks (GRNs) composed of interacting regulatory transcription factors (TFs) to control the timing and magnitude of gene expression. Genome-wide datasets; such as transcriptomics and protein-DNA interactions; and experiments such as high throughput growth curves; facilitate the construction of GRNs and provide insight into TF interactions occurring under stress. Systems biology approaches integrate these datasets into models of GRN architecture as well as statistical and/or dynamical models to understand the function of networks occurring in cells. Previously, these types of studies have focused on traditional model organisms (e.g. Escherichia coli, yeast). However, recent advances in archaeal genetics and other tools have enabled a systems approach to understanding GRNs in these relatively less studied archaeal model organisms. In this report, we outline a systems biology workflow for generating and integrating data focusing on the TF regulator. We discuss experimental design, outline the process of data collection, and provide the tools required to produce high confidence regulons for the TFs of interest. We provide a case study as an example of this workflow, describing the construction of a GRN centered on multi-TF coordinate control of gene expression governing the oxidative stress response in the hypersaline-adapted archaeon Halobacterium salinarum.

Chemosensory regulation of a HEAT-repeat protein couples aggregation and sporulation in Myxococcus xanthus.

Chemosensory systems are complex, highly modified two-component systems (TCS) used by bacteria to control various biological functions ranging from motility to sporulation. Chemosensory systems and TCS both modulate phosphorelays comprised of histidine kinases and response regulators, some of which are single-domain response regulators (SD-RRs) such as CheY. In this study, we have identified and characterized the Che7 chemosensory system of Myxococcus xanthus, a common soil bacterium which displays multicellular development in response to stress. Both genetic and biochemical analyses indicate that the Che7 system regulates development via a direct interaction between the SD-RR CheY7 and a HEAT repeat domain-containing protein, Cpc7. Phosphorylation of the SD-RR affects the interaction with its target, and residues within the α4-ß5-α5 fold of the REC domain govern this interaction. The identification of the Cpc7 interaction with CheY7 extends the diversity of known targets for SD-RRs in biological systems.

Draft Genome of a Type 4 Pilus Defective Myxococcus xanthus Strain, DZF1.

Myxococcus xanthus is a member of the Myxococcales order within the deltaproteobacterial subdivision. Here, we report the whole-genome shotgun sequence of the type IV pilus (T4P) defective strain DZF1, which includes many genes found in strain DZ2 but absent from strain DK1622.

Draft Genome Sequence of Myxococcus xanthus Wild-Type Strain DZ2, a Model Organism for Predation and Development.

Myxococcus xanthus is a member of the Myxococcales order within the Deltaproteobacteria subdivision. The myxobacteria reside in soil, have relatively large genomes, and display complex life cycles. Here, we report the whole-genome shotgun sequence of strain DZ2, which includes unique genes not found previously in strain DK1622.

Inactivation of a novel response regulator is necessary for biofilm formation and host colonization by Vibrio fischeri.

The marine bacterium Vibrio fischeri uses a biofilm to promote colonization of its eukaryotic host Euprymna scolopes. This biofilm depends on the symbiosis polysaccharide (syp) locus, which is transcriptionally regulated by the RscS-SypG two-component regulatory system. An additional response regulator (RR), SypE, exerts both positive and negative control over biofilm formation. SypE is a novel RR protein, with its three putative domains arranged in a unique configuration: a central phosphorylation receiver (REC) domain flanked by two effector domains with putative enzymatic activities (serine kinase and serine phosphatase). To determine how SypE regulates biofilm formation and host colonization, we generated a library of SypE domain mutants. Our results indicate that the N-terminus inhibits biofilm formation, while the C-terminus plays a positive role. The phosphorylation state of SypE appears to regulate these opposing activities, as disruption of the putative site of phosphorylation results in a protein that constitutively inhibits biofilm formation. Furthermore, SypE restricts host colonization: (i) sypE mutants with constitutive inhibitory activity fail to efficiently initiate host colonization and (ii) loss of sypE partially alleviates the colonization defect of an rscS mutant. We conclude that SypE must be inactivated to promote symbiotic colonization by V. fischeri.

The putative hybrid sensor kinase SypF coordinates biofilm formation in Vibrio fischeri by acting upstream of two response regulators, SypG and VpsR.

Colonization of the Hawaiian squid Euprymna scolopes by the marine bacterium Vibrio fischeri requires the symbiosis polysaccharide (syp) gene cluster, which contributes to symbiotic initiation by promoting biofilm formation on the surface of the symbiotic organ. We previously described roles for the syp-encoded response regulator SypG and an unlinked gene encoding the sensor kinase RscS in controlling syp transcription and inducing syp-dependent cell-cell aggregation phenotypes. Here, we report the involvement of an additional syp-encoded regulator, the putative sensor kinase SypF, in promoting biofilm formation. Through the isolation of an increased activity allele, sypF1, we determined that SypF can function to induce syp transcription as well as a variety of biofilm phenotypes, including wrinkled colony formation, adherence to glass, and pellicle formation. SypF1-mediated transcription of the syp cluster was entirely dependent on SypG. However, the biofilm phenotypes were reduced, not eliminated, in the sypG mutant. These phenotypes were also reduced in a mutant deleted for sypE, another syp-encoded response regulator. However, SypF1 still induced phenotypes in a sypG sypE double mutant, suggesting that SypF1 might activate another regulator(s). Our subsequent work revealed that the residual SypF1-induced biofilm formation depended on VpsR, a putative response regulator, and cellulose biosynthesis. These data support a model in which a network of regulators and at least two polysaccharide loci contribute to biofilm formation in V. fischeri.

RscS functions upstream of SypG to control the syp locus and biofilm formation in Vibrio fischeri.

Two-component signal transduction systems, composed of sensor kinase (SK) and response regulator (RR) proteins, allow bacterial cells to adapt to changes such as environmental flux or the presence of a host. RscS is an SK required for Vibrio fischeri to initiate a symbiotic partnership with the Hawaiian squid Euprymna scolopes, likely due to its role in controlling the symbiosis polysaccharide (syp) genes and thus biofilm formation. To determine which RR(s) functions downstream of RscS, we performed epistasis experiments with a library of 35 RR mutants. We found that several RRs contributed to RscS-mediated biofilm formation in V. fischeri. However, only the syp-encoded symbiosis regulator SypG was required for both biofilm phenotypes and syp transcription induced by RscS. These data support the hypothesis that RscS functions upstream of SypG to induce biofilm formation. In addition, this work also revealed a role for the syp-encoded RR SypE in biofilm formation. To our knowledge, no other study has used a large-scale epistasis approach to elucidate two-component signaling pathways. Therefore, this work both contributes to our understanding of regulatory pathways important for symbiotic colonization by V. fischeri and establishes a paradigm for evaluating two-component pathways in the genomics era.

Two-component response regulators of Vibrio fischeri: identification, mutagenesis, and characterization.

Two-component signal transduction systems are utilized by prokaryotic and eukaryotic cells to sense and respond to environmental stimuli, both to maintain homeostasis and to rapidly adapt to changing conditions. Studies have begun to emerge that utilize a large-scale mutagenesis approach to analyzing these systems in prokaryotic organisms. Due to the recent availability of its genome sequence, such a global approach is now possible for the marine bioluminescent bacterium Vibrio fischeri, which exists either in a free-living state or as a mutualistic symbiont within a host organism such as the Hawaiian squid species Euprymna scolopes. In this work, we identified 40 putative two-component response regulators encoded within the V. fischeri genome. Based on the type of effector domain present, we classified six as NarL type, 13 as OmpR type, and six as NtrC type; the remaining 15 lacked a predicted DNA-binding domain. We subsequently mutated 35 of these genes via a vector integration approach and analyzed the resulting mutants for roles in bioluminescence, motility, and competitive colonization of squid. Through these assays, we identified three novel regulators of V. fischeri luminescence and seven regulators that altered motility. Furthermore, we found 11 regulators with a previously undescribed effect on competitive colonization of the host squid. Interestingly, five of the newly characterized regulators each affected two or more of the phenotypes examined, strongly suggesting interconnectivity among systems. This work represents the first large-scale mutagenesis of a class of genes in V. fischeri using a genomic approach and emphasizes the importance of two-component signal transduction in bacterium-host interactions.