A workflow for genome-wide mapping of archaeal transcription factors with ChIP-seq

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.

Unveiling the repressive mechanism of a PPS-like regulator (PspR) in polyhydroxyalkanoates biosynthesis network

Poly(3-hydroxybutyrate-co-3-hydroxyvalerate) (PHBV) is a type of polyhydroxyalkanoates (PHA) that exhibits numerous outstanding properties and is naturally synthesized and elaborately regulated in various microorganisms. However, the regulatory mechanism involving the specific regulator PhaR in Haloferax mediterranei, a major PHBV production model among Haloarchaea, is not well understood. In our previous study, we showed that deletion of the phosphoenolpyruvate (PEP) synthetase-like (pps-like) gene activates the cryptic phaC genes in H. mediterranei, resulting in enhanced PHBV accumulation. In this study, we demonstrated the specific function of the PPS-like protein as a negative regulator of phaR gene expression and PHBV synthesis. Chromatin immunoprecipitation (ChIP), in situ fluorescence reporting system, and in vitro electrophoretic mobility shift assay (EMSA) showed that the PPS-like protein can bind to the promoter region of phaRP. Computational modeling revealed a high structural similarity between the rifampin phosphotransferase (RPH) protein and the PPS-like protein, which has a conserved ATP-binding domain, a His domain, and a predicted DNA-binding domain. Key residues within this unique DNA-binding domain were subsequently validated through point mutation and functional evaluations. Based on these findings, we concluded that PPS-like protein, which we now renamed as PspR, has evolved into a repressor capable of regulating the key regulator PhaR, and thereby modulating PHBV synthesis. This regulatory network (PspR-PhaR) for PHA biosynthesis is likely widespread among haloarchaea, providing a novel approach to manipulate haloarchaea as a production platform for high-yielding PHA. KEY POINTS: • The repressive mechanism of a novel inhibitor PspR in the PHBV biosynthesis was demonstrated • PspR is widespread among the PHA accumulating haloarchaea • It is the first report of functional conversion from an enzyme to a trans-acting regulator in haloarchaea.

TroR is the primary regulator of the iron homeostasis transcription network in the halophilic archaeon Haloferax volcanii

Maintaining the intracellular iron concentration within the homeostatic range is vital to meet cellular metabolic needs and reduce oxidative stress. Previous research revealed that the haloarchaeon Halobacterium salinarum encodes four diphtheria toxin repressor (DtxR) family transcription factors (TFs) that together regulate the iron response through an interconnected transcriptional regulatory network (TRN). However, the conservation of the TRN and the metal specificity of DtxR TFs remained poorly understood. Here we identified and characterized the TRN of Haloferax volcanii for comparison. Genetic analysis demonstrated that Hfx. volcanii relies on three DtxR transcriptional regulators (Idr, SirR, and TroR), with TroR as the primary regulator of iron homeostasis. Bioinformatics and molecular approaches revealed that TroR binds a conserved cis-regulatory motif located ∼100 nt upstream of the start codon of iron-related target genes. Transcriptomics analysis demonstrated that, under conditions of iron sufficiency, TroR repressed iron uptake and induced iron storage mechanisms. TroR repressed the expression of one other DtxR TF, Idr. This reduced DtxR TRN complexity relative to that of Hbt. salinarum appeared correlated with natural variations in iron availability. Based on these data, we hypothesize that variable environmental conditions such as iron availability appear to select for increasing TRN complexity.

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.

MITF Contributes to the Body Color Differentiation of Sea Cucumbers Apostichopus japonicus through Expression Differences and Regulation of Downstream Genes

Melanin, which is a pigment produced in melanocytes, is an important contributor to sea cucumber body color. MITF is one of the most critical genes in melanocyte development and melanin synthesis pathways. However, how MITF regulates body color and differentiation in sea cucumbers is poorly understood. In this study, we analyzed the expression level and location of MITF in white, purple, and green sea cucumbers and identified the genes regulated by MITF using chromatin immunoprecipitation followed by sequencing. The mRNA and protein expression levels of MITF were all highest in purple morphs and lowest in white morphs. In situ hybridization indicated that MITF mRNA were mainly expressed in the epidermis. We also identified 984, 732, and 1191 peaks of MITF binding in green, purple, and white sea cucumbers, which were associated with 727, 557, and 887 genes, respectively. Our findings suggested that MITF contributed to the body color differentiation of green, purple, and white sea cucumbers through expression differences and regulation of downstream genes. These results provided a basis for future studies to determine the mechanisms underlying body color formation and provided insights into gene regulation in sea cucumbers.

TrmB Family Transcription Factor as a Thiol-Based Regulator of Oxidative Stress Response

Oxidative stress causes cellular damage, including DNA mutations, protein dysfunction, and loss of membrane integrity. Here, we discovered that a TrmB (transcription regulator of mal operon) family protein (Pfam PF01978) composed of a single winged-helix DNA binding domain (InterPro IPR002831) can function as thiol-based transcriptional regulator of oxidative stress response. Using the archaeon Haloferax volcanii as a model system, we demonstrate that the TrmB-like OxsR is important for recovery of cells from hypochlorite stress. OxsR is shown to bind specific regions of genomic DNA, particularly during hypochlorite stress. OxsR-bound intergenic regions were found proximal to oxidative stress operons, including genes associated with thiol relay and low molecular weight thiol biosynthesis. Further analysis of a subset of these sites revealed OxsR to function during hypochlorite stress as a transcriptional activator and repressor. OxsR was shown to require a conserved cysteine (C24) for function and to use a CG-rich motif upstream of conserved BRE/TATA box promoter elements for transcriptional activation. Protein modeling suggested the C24 is located at a homodimer interface formed by antiparallel α helices, and that oxidation of this cysteine would result in the formation of an intersubunit disulfide bond. This covalent linkage may promote stabilization of an OxsR homodimer with the enhanced DNA binding properties observed in the presence of hypochlorite stress. The phylogenetic distribution TrmB family proteins, like OxsR, that have a single winged-helix DNA binding domain and conserved cysteine residue suggests this type of redox signaling mechanism is widespread in Archaea.

ChIP-Seq Occupancy Mapping of the Archaeal Transcription Machinery

Genome-wide occupancy studies for RNA polymerases and their basal transcription factors deliver information about transcription dynamics and the recruitment of transcription elongation and termination factors in eukaryotes and prokaryotes. The primary method to determine genome-wide occupancies is chromatin immunoprecipitation combined with deep sequencing (ChIP-seq). Archaea possess a transcription machinery that is evolutionarily closer related to its eukaryotic counterpart but it operates in a prokaryotic cellular context. Studies on archaeal transcription brought insight into the evolution of transcription machineries and the universality of transcription mechanisms. Because of the limited resolution of ChIP-seq, the close spacing of promoters and transcription units found in archaeal genomes pose a challenge for ChIP-seq and the ensuing data analysis. The extreme growth temperature of many established archaeal model organisms necessitates further adaptations. This chapter describes a version of ChIP-seq adapted for the basal transcription machinery of thermophilic archaea and some modifications to the data analysis.

Identifying Components of a Halobacterium salinarum N-Glycosylation Pathway

Whereas N-glycosylation is a seemingly universal process in Archaea, pathways of N-glycosylation have only been experimentally verified in a mere handful of species. Toward expanding the number of delineated archaeal N-glycosylation pathways, the involvement of the putative Halobacterium salinarum glycosyltransferases VNG1067G, VNG1066C, and VNG1062G in the assembly of an N-linked tetrasaccharide decorating glycoproteins in this species was addressed. Following deletion of each encoding gene, the impact on N-glycosylation of the S-layer glycoprotein and archaellins, major glycoproteins in this organism, was assessed by mass spectrometry. Likewise, the pool of dolichol phosphate, the lipid upon which this glycan is assembled, was also considered in each deletion strain. Finally, the impacts of such deletions were characterized in a series of biochemical, structural and physiological assays. The results revealed that VNG1067G, VNG1066C, and VNG1062G, renamed Agl25, Agl26, and Agl27 according to the nomenclature used for archaeal N-glycosylation pathway components, are responsible for adding the second, third and fourth sugars of the N-linked tetrasaccharide decorating Hbt. salinarum glycoproteins. Moreover, this study demonstrated how compromised N-glycosylation affects various facets of Hbt. salinarum cell behavior, including the transcription of archaellin-encoding genes.

An archaeal histone-like protein regulates gene expression in response to salt stress

Histones, ubiquitous in eukaryotes as DNA-packing proteins, find their evolutionary origins in archaea. Unlike the characterized histone proteins of a number of methanogenic and themophilic archaea, previous research indicated that HpyA, the sole histone encoded in the model halophile Halobacterium salinarum, is not involved in DNA packaging. Instead, it was found to have widespread but subtle effects on gene expression and to maintain wild type cell morphology. However, the precise function of halophilic histone-like proteins remain unclear. Here we use quantitative phenotyping, genetics, and functional genomics to investigate HpyA function. These experiments revealed that HpyA is important for growth and rod-shaped morphology in reduced salinity. HpyA preferentially binds DNA at discrete genomic sites under low salt to regulate expression of ion uptake, particularly iron. HpyA also globally but indirectly activates other ion uptake and nucleotide biosynthesis pathways in a salt-dependent manner. Taken together, these results demonstrate an alternative function for an archaeal histone-like protein as a transcriptional regulator, with its function tuned to the physiological stressors of the hypersaline environment.

The conserved ribonuclease aCPSF1 triggers genome-wide transcription termination of Archaea via a 3'-end cleavage mode

Transcription termination defines accurate transcript 3′-ends and ensures programmed transcriptomes, making it critical to life. However, transcription termination mechanisms remain largely unknown in Archaea. Here, we reported the physiological significance of the newly identified general transcription termination factor of Archaea, the ribonuclease aCPSF1, and elucidated its 3′-end cleavage triggered termination mechanism. The depletion of Mmp-aCPSF1 in Methanococcus maripaludis caused a genome-wide transcription termination defect and disordered transcriptome. Transcript-3′end-sequencing revealed that transcriptions primarily terminate downstream of a uridine-rich motif where Mmp-aCPSF1 performed an endoribonucleolytic cleavage, and the endoribonuclease activity was determined to be essential to the in vivo transcription termination. Co-immunoprecipitation and chromatin-immunoprecipitation detected interactions of Mmp-aCPSF1 with RNA polymerase and chromosome. Phylogenetic analysis revealed that the aCPSF1 orthologs are ubiquitously distributed among the archaeal phyla, and two aCPSF1 orthologs from Lokiarchaeota and Thaumarchaeota could replace Mmp-aCPSF1 to terminate transcription of M. maripaludis. Therefore, the aCPSF1 dependent termination mechanism could be widely employed in Archaea, including Lokiarchaeota belonging to Asgard Archaea, the postulated archaeal ancestor of Eukaryotes. Strikingly, aCPSF1-dependent archaeal transcription termination reported here exposes a similar 3′-cleavage mode as the eukaryotic RNA polymerase II termination, thus would shed lights on understanding the evolutionary linking between archaeal and eukaryotic termination machineries.

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.

PvrA is a novel regulator that contributes to Pseudomonas aeruginosa pathogenesis by controlling bacterial utilization of long chain fatty acids

During infection of a host, Pseudomonas aeruginosa orchestrates global gene expression to adapt to the host environment and counter the immune attacks. P. aeruginosa harbours hundreds of regulatory genes that play essential roles in controlling gene expression. However, their contributions to the bacterial pathogenesis remain largely unknown. In this study, we analysed the transcriptomic profile of P. aeruginosa cells isolated from lungs of infected mice and examined the roles of upregulated regulatory genes in bacterial virulence. Mutation of a novel regulatory gene pvrA (PA2957) attenuated the bacterial virulence in an acute pneumonia model. Chromatin immunoprecipitation (ChIP)-Seq and genetic analyses revealed that PvrA directly regulates genes involved in phosphatidylcholine utilization and fatty acid catabolism. Mutation of the pvrA resulted in defective bacterial growth when phosphatidylcholine or palmitic acid was used as the sole carbon source. We further demonstrated that palmitoyl coenzyme A is a ligand for the PvrA, enhancing the binding affinity of PvrA to its target promoters. An arginine residue at position 136 was found to be essential for PvrA to bind palmitoyl coenzyme A. Overall, our results revealed a novel regulatory pathway that controls genes involved in phosphatidylcholine and fatty acid utilization and contributes to the bacterial virulence.

DNA mapping and kinetic modeling of the HrdB regulon in Streptomyces coelicolor

HrdB in streptomycetes is a principal sigma factor whose deletion is lethal. This is also the reason why its regulon has not been investigated so far. To overcome experimental obstacles, for investigating the HrdB regulon, we constructed a strain whose HrdB protein was tagged by an HA epitope. ChIP-seq experiment, done in 3 repeats, identified 2137 protein-coding genes organized in 337 operons, 75 small RNAs, 62 tRNAs, 6 rRNAs and 3 miscellaneous RNAs. Subsequent kinetic modeling of regulation of protein-coding genes with HrdB alone and with a complex of HrdB and a transcriptional cofactor RbpA, using gene expression time series, identified 1694 genes that were under their direct control. When using the HrdB-RbpA complex in the model, an increase of the model fidelity was found for 322 genes. Functional analysis revealed that HrdB controls the majority of gene groups essential for the primary metabolism and the vegetative growth. Particularly, almost all ribosomal protein-coding genes were found in the HrdB regulon. Analysis of promoter binding sites revealed binding motif at the -10 region and suggested the possible role of mono- or di-nucleotides upstream of the -10 element.

Key Concepts and Challenges in Archaeal Transcription

Transcription is enabled by RNA polymerase and general factors that allow its progress through the transcription cycle by facilitating initiation, elongation and termination. The transitions between specific stages of the transcription cycle provide opportunities for the global and gene-specific regulation of gene expression. The exact mechanisms and the extent to which the different steps of transcription are exploited for regulation vary between the domains of life, individual species and transcription units. However, a surprising degree of conservation is apparent. Similar key steps in the transcription cycle can be targeted by homologous or unrelated factors providing insights into the mechanisms of RNAP and the evolution of the transcription machinery. Archaea are bona fide prokaryotes but employ a eukaryote-like transcription system to express the information of bacteria-like genomes. Thus, archaea provide the means not only to study transcription mechanisms of interesting model systems but also to test key concepts of regulation in this arena. In this review, we discuss key principles of archaeal transcription, new questions that still await experimental investigation, and how novel integrative approaches hold great promise to fill this gap in our knowledge.

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.

Internal RNAs overlapping coding sequences can drive the production of alternative proteins in archaea

Prokaryotic genomes show a high level of information compaction often with different molecules transcribed from the same locus. Although antisense RNAs have been relatively well studied, RNAs in the same strand, internal RNAs (intraRNAs), are still poorly understood. The question of how common is the translation of overlapping reading frames remains open. We address this question in the model archaeon Halobacterium salinarum. In the present work we used differential RNA-seq (dRNA-seq) in H. salinarum NRC-1 to locate intraRNA signals in subsets of internal transcription start sites (iTSS) and establish the open reading frames associated to them (intraORFs). Using C-terminally flagged proteins, we experimentally observed isoforms accurately predicted by intraRNA translation for kef1, acs3 and orc4 genes. We also recovered from the literature and mass spectrometry databases several instances of protein isoforms consistent with intraRNA translation such as the gas vesicle protein gene gvpC1. We found evidence for intraRNAs in horizontally transferred genes such as the chaperone dnaK and the aerobic respiration related cydA in both H. salinarum and Escherichia coli. Also, intraRNA translation evidence in H. salinarum, E. coli and yeast of a universal elongation factor (aEF-2, fusA and eEF-2) suggests that this is an ancient phenomenon present in all domains of life.

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.

The Transcriptional Regulator TFB-RF1 Activates Transcription of a Putative ABC Transporter in Pyrococcus furiosus

Transcription factor B recruiting factor 1 (TFB-RF1; PF1088) is a transcription regulator which activates transcription on archaeal promoters containing weak TFB recognition elements (BRE) by recruiting TFB to the promoter. The mechanism of activation is described in detail, but nothing is known about the biological function of this protein in Pyrococcus furiosus. The protein is located in an operon structure together with the hypothetical gene pf1089 and western blot as well as end-point RT-PCR experiments revealed an extremely low expression rate of both proteins. Furthermore, conditions to induce the expression of the operon are not known. By introducing an additional copy of tfb-RF1 using a Pyrococcus shuttle vector we could circumvent the lacking expression of both proteins under standard growth conditions as indicated by western blot as well as end-point RT-PCR experiments. A ChIP-seq experiment revealed an additional binding site of TFB-RF1 in the upstream region of the pf1011/1012 operon, beside the expected target of the pf1089/tfb-RF1 region. This operon codes for a putative ABC transporter which is most-related to a multidrug export system and in vitro analysis using gel shift assays, DNase I footprinting and in vitro transcription confirmed the activator function of TFB-RF1 on the corresponding promoter. These findings are also in agreement with in vivo data, as RT-qPCR experiments also indicate transcriptional activation of both operons. Taken together, the overexpression strategy of tfb-RF1 enabled the identification of an additional operon of the TFB-RF1 regulon which indicates a transport-related function and provides a promising starting position to decipher the physiological function of the TFB-RF1 gene regulatory network in P. furiosus.

Genome-Wide Profiling of DNA Methyltransferases in Mammalian Cells

Chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) is currently the method of choice to determine binding sites of chromatin-associated factors in a genome-wide manner. Here, we describe a method to investigate the binding preferences of mammalian DNA methyltransferases (DNMT) based on ChIP-seq using biotin-tagging. Stringent ChIP of DNMT proteins based on the strong interaction between biotin and avidin circumvents limitations arising from low antibody specificity and ensures reproducible enrichment. DNMT-bound DNA fragments are ligated to sequencing adaptors, amplified and sequenced on a high-throughput sequencing instrument. Bioinformatic analysis gives valuable information about the binding preferences of DNMTs genome-wide and around promoter regions. This method is unconventional due to the use of genetically engineered cells; however, it allows specific and reliable determination of DNMT binding.

A transcription network of interlocking positive feedback loops maintains intracellular iron balance in archaea

Iron is required for key metabolic processes but is toxic in excess. This circumstance forces organisms across the tree of life to tightly regulate iron homeostasis. In hypersaline lakes dominated by archaeal species, iron levels are extremely low and subject to environmental change; however, mechanisms regulating iron homeostasis in archaea remain unclear. In previous work, we demonstrated that two transcription factors (TFs), Idr1 and Idr2, collaboratively regulate aspects of iron homeostasis in the model species Halobacterium salinarum. Here we show that Idr1 and Idr2 are part of an extended regulatory network of four TFs of the bacterial DtxR family that maintains intracellular iron balance. We demonstrate that each TF directly regulates at least one of the other DtxR TFs at the level of transcription. Dynamical modeling revealed interlocking positive feedback loop architecture, which exhibits bistable or oscillatory network dynamics depending on iron availability. TF knockout mutant phenotypes are consistent with model predictions. Together, our results support that this network regulates iron homeostasis despite variation in extracellular iron levels, consistent with dynamical properties of interlocking feedback architecture in eukaryotes. These results suggest that archaea use bacterial-type TFs in a eukaryotic regulatory network topology to adapt to harsh environments.

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

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.

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.

In vitro site-specific recombination mediated by the tyrosine recombinase XerA of Thermoplasma acidophilum

In organisms with circular chromosomes, such as bacteria and archaea, an odd number of homologous recombination events can generate a chromosome dimer. Such chromosome dimers cannot be segregated unless they are converted to monomers before cell division. In Escherichia coli, dimer-to-monomer conversion is mediated by the paralogous XerC and XerD recombinases at a specific dif site in the replication termination region. Dimer resolution requires the highly conserved cell division protein/chromosome translocase FtsK, and this site-specific chromosome resolution system is present or predicted in most bacteria. However, most archaea have only XerA, a homologue of the bacterial XerC/D proteins, but no homologues of FtsK. In addition, the molecular mechanism of XerA-mediated chromosome resolution in archaea has been less thoroughly elucidated than those of the corresponding bacterial systems. In this study, we identified two XerA-binding sites (dif1 and dif2) in the Thermoplasma acidophilum chromosome. In vitro site-specific recombination assays showed that dif2, but not dif1, serves as a target site for XerA-mediated chromosome resolution. Mutational analysis indicated that not only the core consensus sequence of dif2, but also its flanking regions play important roles in the recognition and recombination reactions mediated by XerA.

The genome-scale DNA-binding profile of BarR, a β-alanine responsive transcription factor in the archaeon Sulfolobus acidocaldarius

Background

The Leucine-responsive Regulatory Protein (Lrp) family is a widespread family of regulatory transcription factors in prokaryotes. BarR is an Lrp-like transcription factor in the model archaeon Sulfolobus acidocaldarius that activates the expression of a β-alanine aminotransferase gene, which is involved in β-alanine degradation. In contrast to classical Lrp-like transcription factors, BarR is not responsive to any of the α-amino acids but interacts specifically with β-alanine. Besides the juxtaposed β-alanine aminotransferase gene, other regulatory targets of BarR have not yet been identified although β-alanine is the precursor of coenzyme A and thus an important central metabolite. The aim of this study is to extend the knowledge of the DNA-binding characteristics of BarR and of its corresponding regulon from a local to a genome-wide perspective.

Results

We characterized the genome-wide binding profile of BarR using chromatin immunoprecipation combined with high-throughput sequencing (ChIP-seq). This revealed 21 genomic binding loci. High-enrichment binding regions were validated to interact with purified BarR protein in vitro using electrophoretic mobility shift assays and almost all targets were also shown to harbour a conserved semi-palindromic binding motif. Only a small subset of enriched genomic sites are located in intergenic regions at a relative short distance to a promoter, and qRT-PCR analysis demonstrated that only one additional operon is under activation of BarR, namely the glutamine synthase operon. The latter is also a target of other Lrp-like transcription factors. Detailed inspection of the BarR ChIP-seq profile at the β-alanine aminotransferase promoter region in combination with binding motif predictions indicate that the operator structure is more complicated than previously anticipated, consisting of multiple (major and auxiliary) operators.

Conclusions

BarR has a limited regulon, and includes also glutamine synthase genes besides the previously characterized β-alanine aminotransferase. Regulation of glutamine synthase is suggestive of a link between β-alanine and α-amino acid metabolism in S. acidocaldarius. Furthermore, this work reveals that the BarR regulon overlaps with that of other Lrp-like regulators.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2890-0) contains supplementary material, which is available to authorized users.

Genome-wide binding analysis of the transcriptional regulator TrmBL1 in Pyrococcus furiosus

Background

Several in vitro studies document the function of the transcriptional regulator TrmBL1 of Pyrococcus furiosus. These data indicate that the protein can act as repressor or activator and is mainly involved in transcriptional control of sugar uptake and in the switch between glycolysis and gluconeogenesis. The aim of this study was to complement the in vitro data with an in vivo analysis using ChIP-seq to explore the genome-wide binding profile of TrmBL1 under glycolytic and gluconeogenic growth conditions.

Results

The ChIP-seq analysis revealed under gluconeogenic growth conditions 28 TrmBL1 binding sites where the TGM is located upstream of coding regions and no binding sites under glycolytic conditions. The experimental confirmation of the binding sites using qPCR, EMSA, DNase I footprinting and in vitro transcription experiments validated the in vivo identified TrmBL1 binding sites. Furthermore, this study provides evidence that TrmBL1 is also involved in transcriptional regulation of additional cellular processes e.g. amino acid metabolism, transcriptional control or metabolic pathways. In the initial setup we were interested to include the binding analysis of TrmB, an additional member of the TrmB family, but western blot experiments and the ChIP-seq data indicated that the corresponding gene is deleted in our Pyrococcus strain. A detailed analysis of a new type strain demonstrated that a 16 kb fragment containing the trmb gene is almost completely deleted after the first re-cultivation.

Conclusions

The identified binding sites in the P. furiosus genome classified TrmBL1 as a more global regulator as hitherto known. Furthermore, the high resolution of the mapped binding positions enabled reliable predictions, if TrmBL1 activates (binding site upstream of the promoter) or represses transcription (binding site downstream) of the corresponding genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2360-0) contains supplementary material, which is available to authorized users.

Identification of an archaeal mercury regulon by chromatin immunoprecipitation

Mercury is a heavy metal and toxic to all forms of life. Metal exposure can invoke a response to improve survival. In archaea, several components of a mercury response system have been identified, but it is not known whether metal transport is a member of this system. To identify such missing components, a peptide-tagged MerR transcription factor was used to localize enriched chromosome regions by chromosome immunoprecipitation combined with DNA sequence analysis. Such regions could serve as secondary regulatory binding sites to control the expression of additional genes associated with mercury detoxification. Among the 31 highly enriched loci, a subset of five was pursued as potential candidates based on their current annotations. Quantitative reverse transcription-PCR analysis of these regions with and without mercury treatment in WT and mutant strains lacking merR indicated significant regulatory responses under these conditions. Of these, a Family 5 extracellular solute-binding protein and the MarR transcription factor shown previously to control responses to oxidation were most strongly affected. Inactivation of the solute-binding protein by gene disruption increased the resistance of mutant cells to mercury challenge. Inductively coupled plasma-MS analysis of the mutant cell line following metal challenge indicated there was less intracellular mercury compared with the isogenic WT strain. Together, these regulated genes comprise new members of the archaeal MerR regulon and reveal a cascade of transcriptional control not previously demonstrated in this model organism.

Growth-Phase-Specific Modulation of Cell Morphology and Gene Expression by an Archaeal Histone Protein

In all three domains of life, organisms use nonspecific DNA-binding proteins to compact and organize the genome as well as to regulate transcription on a global scale. Histone is the primary eukaryotic nucleoprotein, and its evolutionary roots can be traced to the archaea. However, not all archaea use this protein as the primary DNA-packaging component, raising questions regarding the role of histones in archaeal chromatin function. Here, quantitative phenotyping, transcriptomic, and proteomic assays were performed on deletion and overexpression mutants of the sole histone protein of the hypersaline-adapted haloarchaeal model organism Halobacterium salinarum. This protein is highly conserved among all sequenced haloarchaeal species and maintains hallmark residues required for eukaryotic histone functions. Surprisingly, despite this conservation at the sequence level, unlike in other archaea or eukaryotes, H. salinarum histone is required to regulate cell shape but is not necessary for survival. Genome-wide expression changes in histone deletion strains were global, significant but subtle in terms of fold change, bidirectional, and growth phase dependent. Mass spectrometric proteomic identification of proteins from chromatin enrichments yielded levels of histone and putative nucleoid-associated proteins similar to those of transcription factors, consistent with an open and transcriptionally active genome. Taken together, these data suggest that histone in H. salinarum plays a minor role in DNA compaction but important roles in growth-phase-dependent gene expression and regulation of cell shape. Histone function in haloarchaea more closely resembles a regulator of gene expression than a chromatin-organizing protein like canonical eukaryotic histone.

Histones comprise the major protein component of eukaryotic chromatin and are required for both genome packaging and global regulation of expression. The current paradigm maintains that archaea whose genes encode histone also use these proteins to package DNA. In contrast, here we demonstrate that the sole histone encoded in the genome of the salt-adapted archaeon Halobacterium salinarum is both unessential and unlikely to be involved in DNA compaction despite conservation of residues important for eukaryotic histones. Rather, H. salinarum histone is required for global regulation of gene expression and cell shape. These data are consistent with the hypothesis that H. salinarum histone, strongly conserved across all other known salt-adapted archaea, serves a novel role in gene regulation and cell shape maintenance. Given that archaea possess the ancestral form of eukaryotic histone, this study has important implications for understanding the evolution of histone function.

Development of New Modular Genetic Tools for Engineering the Halophilic Archaeon Halobacterium salinarum

Our ability to genetically manipulate living organisms is usually constrained by the efficiency of the genetic tools available for the system of interest. In this report, we present the design, construction and characterization of a set of four new modular vectors, the pHsal series, for engineering Halobacterium salinarum, a model halophilic archaeon widely used in systems biology studies. The pHsal shuttle vectors are organized in four modules: (i) the E. coli’s specific part, containing a ColE1 origin of replication and an ampicillin resistance marker, (ii) the resistance marker and (iii) the replication origin, which are specific to H. salinarum and (iv) the cargo, which will carry a sequence of interest cloned in a multiple cloning site, flanked by universal M13 primers. Each module was constructed using only minimal functional elements that were sequence edited to eliminate redundant restriction sites useful for cloning. This optimization process allowed the construction of vectors with reduced sizes compared to currently available platforms and expanded multiple cloning sites. Additionally, the strong constitutive promoter of the fer2 gene was sequence optimized and incorporated into the platform to allow high-level expression of heterologous genes in H. salinarum. The system also includes a new minimal suicide vector for the generation of knockouts and/or the incorporation of chromosomal tags, as well as a vector for promoter probing using a GFP gene as reporter. This new set of optimized vectors should strongly facilitate the engineering of H. salinarum and similar strategies could be implemented for other archaea.

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.

A novel DNA-binding protein, PhaR, plays a central role in the regulation of polyhydroxyalkanoate accumulation and granule formation in the haloarchaeon Haloferax mediterranei

Polyhydroxyalkanoates (PHAs) are synthesized and assembled as PHA granules that undergo well-regulated formation in many microorganisms. However, this regulation remains unclear in haloarchaea. In this study, we identified a PHA granule-associated regulator (PhaR) that negatively regulates the expression of both its own gene and the granule structural gene phaP in the same operon (phaRP) in Haloferax mediterranei. Chromatin immunoprecipitation-quantitative PCR (ChIP-qPCR) assays demonstrated a significant interaction between PhaR and the phaRP promoter in vivo. Scanning mutagenesis of the phaRP promoter revealed a specific cis-element as the possible binding position of the PhaR. The haloarchaeal homologs of the PhaR contain a novel conserved domain that belongs to a swapped-hairpin barrel fold family found in AbrB-like proteins. Amino acid substitution indicated that this AbrB-like domain is critical for the repression activity of PhaR. In addition, the phaRP promoter had a weaker activity in the PHA-negative strains, implying a function of the PHA granules in titration of the PhaR. Moreover, the H. mediterranei strain lacking phaR was deficient in PHA accumulation and produced granules with irregular shapes. Interestingly, the PhaR itself can promote PHA synthesis and granule formation in a PhaP-independent manner. Collectively, our results demonstrated that the haloarchaeal PhaR is a novel bifunctional protein that plays the central role in the regulation of PHA accumulation and granule formation in H. mediterranei.

Inference of expanded Lrp-like feast/famine transcription factor targets in a non-model organism using protein structure-based prediction

Widespread microbial genome sequencing presents an opportunity to understand the gene regulatory networks of non-model organisms. This requires knowledge of the binding sites for transcription factors whose DNA-binding properties are unknown or difficult to infer. We adapted a protein structure-based method to predict the specificities and putative regulons of homologous transcription factors across diverse species. As a proof-of-concept we predicted the specificities and transcriptional target genes of divergent archaeal feast/famine regulatory proteins, several of which are encoded in the genome of Halobacterium salinarum. This was validated by comparison to experimentally determined specificities for transcription factors in distantly related extremophiles, chromatin immunoprecipitation experiments, and cis-regulatory sequence conservation across eighteen related species of halobacteria. Through this analysis we were able to infer that Halobacterium salinarum employs a divergent local trans-regulatory strategy to regulate genes (carA and carB) involved in arginine and pyrimidine metabolism, whereas Escherichia coli employs an operon. The prediction of gene regulatory binding sites using structure-based methods is useful for the inference of gene regulatory relationships in new species that are otherwise difficult to infer.

Analysis of the transcriptional regulator GlpR, promoter elements, and posttranscriptional processing involved in fructose-induced activation of the phosphoenolpyruvate-dependent sugar phosphotransferase system in Haloferax mediterranei

Among all known archaeal strains, the phosphoenolpyruvate-dependent phosphotransferase system (PTS) for fructose utilization is used primarily by haloarchaea, which thrive in hypersaline environments, whereas the molecular details of the regulation of the archaeal PTS under fructose induction remain unclear. In this study, we present a comprehensive examination of the regulatory mechanism of the fructose PTS in the haloarchaeon Haloferax mediterranei. With gene knockout and complementation, microarray analysis, and chromatin immunoprecipitation-quantitative PCR (ChIP-qPCR), we revealed that GlpR is the indispensable activator, which specifically binds to the PTS promoter (PPTS) during fructose induction. Further promoter-scanning mutation indicated that three sites located upstream of the H. mediterranei PPTS, which are conserved in most haloarchaeal PPTSs, are involved in this induction. Interestingly, two PTS transcripts (named T8 and T17) with different lengths of 5' untranslated region (UTR) were observed, and promoter or 5' UTR swap experiments indicated that the shorter 5' UTR was most likely generated from the longer one. Notably, the translation efficiency of the transcript with this shorter 5' UTR was significantly higher and the ratio of T8 (with the shorter 5' UTR) to T17 increased during fructose induction, implying that a posttranscriptional mechanism is also involved in PTS activation. With these insights into the molecular regulation of the haloarchaeal PTS, we have proposed a working model for haloarchaea in response to environmental fructose.

Promoter element arising from the fusion of standard BioBrick parts

We characterize the appearance of a constitutive promoter element in the commonly used cI repressor-encoding BioBrick BBa_C0051. We have termed this promoter element pKAT. Full pKAT activity is created by the ordered assembly of sequences in BBa_C0051 downstream of the cI gene encoding the 11 amino acid LVA proteolytic degradation tag, a BioBrick standard double-TAA stop codon, a genetic barcode, and part of the RFC10 SpeI-XbaI BioBrick scar. Placing BBa_C0051 or other pKAT containing parts upstream of other functional RNA coding elements in a polycistronic context may therefore lead to the unintended transcription of the downstream elements. The frequent reuse of pKAT or pKAT-like containing basic parts in the Registry of Biological Parts has resulted in approximately 5% of registry parts encoding at least one instance of a predicted pKAT promoter located directly upstream of a ribosome binding site and ATG start codon. This example highlights that even seemingly simple modifications of a part's sequence (in this case addition of degradation tags and barcodes) may be sufficient to unexpectedly change the contextual behavior of a part and reaffirms the inherent challenge in carefully characterizing the behavior of standardized biological parts across a broad range of reasonable use scenarios.

Next-generation sequencing technologies and their impact on microbial genomics

Next-generation sequencing technologies have had a dramatic impact in the field of genomic research through the provision of a low cost, high-throughput alternative to traditional capillary sequencers. These new sequencing methods have surpassed their original scope and now provide a range of utility-based applications, which allow for a more comprehensive analysis of the structure and content of microbial genomes than was previously possible. With the commercialization of a third generation of sequencing technologies imminent, we discuss the applications of current next-generation sequencing methods and explore their impact on and contribution to microbial genome research.

Nanobody(R)-based chromatin immunoprecipitation/micro-array analysis for genome-wide identification of transcription factor DNA binding sites

Nanobodies® are single-domain antibody fragments derived from camelid heavy-chain antibodies. Because of their small size, straightforward production in Escherichia coli, easy tailoring, high affinity, specificity, stability and solubility, nanobodies® have been exploited in various biotechnological applications. A major challenge in the post-genomics and post-proteomics era is the identification of regulatory networks involving nucleic acid-protein and protein-protein interactions. Here, we apply a nanobody® in chromatin immunoprecipitation followed by DNA microarray hybridization (ChIP-chip) for genome-wide identification of DNA-protein interactions. The Lrp-like regulator Ss-LrpB, arguably one of the best-studied specific transcription factors of the hyperthermophilic archaeon Sulfolobus solfataricus, was chosen for this proof-of-principle nanobody®-assisted ChIP. Three distinct Ss-LrpB-specific nanobodies®, each interacting with a different epitope, were generated for ChIP. Genome-wide ChIP-chip with one of these nanobodies® identified the well-established Ss-LrpB binding sites and revealed several unknown target sequences. Furthermore, these ChIP-chip profiles revealed auxiliary operator sites in the open reading frame of Ss-lrpB. Our work introduces nanobodies® as a novel class of affinity reagents for ChIP. Taking into account the unique characteristics of nanobodies®, in particular, their short generation time, nanobody®-based ChIP is expected to further streamline ChIP-chip and ChIP-Seq experiments, especially in organisms with no (or limited) possibility of genetic manipulation.

A Monte Carlo-based framework enhances the discovery and interpretation of regulatory sequence motifs

BACKGROUND

Discovery of functionally significant short, statistically overrepresented subsequence patterns (motifs) in a set of sequences is a challenging problem in bioinformatics. Oftentimes, not all sequences in the set contain a motif. These non-motif-containing sequences complicate the algorithmic discovery of motifs. Filtering the non-motif-containing sequences from the larger set of sequences while simultaneously determining the identity of the motif is, therefore, desirable and a non-trivial problem in motif discovery research.

RESULTS

We describe MotifCatcher, a framework that extends the sensitivity of existing motif-finding tools by employing random sampling to effectively remove non-motif-containing sequences from the motif search. We developed two implementations of our algorithm; each built around a commonly used motif-finding tool, and applied our algorithm to three diverse chromatin immunoprecipitation (ChIP) data sets. In each case, the motif finder with the MotifCatcher extension demonstrated improved sensitivity over the motif finder alone. Our approach organizes candidate functionally significant discovered motifs into a tree, which allowed us to make additional insights. In all cases, we were able to support our findings with experimental work from the literature.

CONCLUSIONS

Our framework demonstrates that additional processing at the sequence entry level can significantly improve the performance of existing motif-finding tools. For each biological data set tested, we were able to propose novel biological hypotheses supported by experimental work from the literature. Specifically, in Escherichia coli, we suggested binding site motifs for 6 non-traditional LexA protein binding sites; in Saccharomyces cerevisiae, we hypothesize 2 disparate mechanisms for novel binding sites of the Cse4p protein; and in Halobacterium sp. NRC-1, we discoverd subtle differences in a general transcription factor (GTF) binding site motif across several data sets. We suggest that small differences in our discovered motif could confer specificity for one or more homologous GTF proteins. We offer a free implementation of the MotifCatcher software package at http://www.bme.ucdavis.edu/facciotti/resources_data/software/.

Overview of the genetic tools in the Archaea

This section provides an overview of the genetic systems developed in the Archaea. Genetic manipulation is possible in many members of the halophiles, methanogens, Sulfolobus, and Thermococcales. We describe the selection/counterselection principles utilized in each of these groups, which consist of antibiotics and their resistance markers, and auxotrophic host strains and complementary markers. The latter strategy utilizes techniques similar to those developed in yeast. However, Archaea are resistant to many of the antibiotics routinely used for selection in the Bacteria, and a number of strategies specific to the Archaea have been developed. In addition, examples utilizing the genetic systems developed for each group will be briefly described.