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

Perturbed N-glycosylation of Halobacterium salinarum archaellum filaments leads to filament bundling and compromised cell motility

The swimming device of archaea-the archaellum-presents asparagine (N)-linked glycans. While N-glycosylation serves numerous roles in archaea, including enabling their survival in extreme environments, how this post-translational modification contributes to cell motility remains under-explored. Here, we report the cryo-EM structure of archaellum filaments from the haloarchaeon Halobacterium salinarum, where archaellins, the building blocks of the archaellum, are N-glycosylated, and the N-glycosylation pathway is well-resolved. We further determined structures of archaellum filaments from two N-glycosylation mutant strains that generate truncated glycans and analyzed their motility. While cells from the parent strain exhibited unidirectional motility, the N-glycosylation mutant strain cells swam in ever-changing directions within a limited area. Although these mutant strain cells presented archaellum filaments that were highly similar in architecture to those of the parent strain, N-linked glycan truncation greatly affected interactions between archaellum filaments, leading to dramatic clustering of both isolated and cell-attached filaments. We propose that the N-linked tetrasaccharides decorating archaellins act as physical spacers that minimize the archaellum filament aggregation that limits cell motility.

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.

Agl28 and Agl29 are key components of a Halobacterium salinarum N-glycosylation pathway

Although Halobacterim salinarum provided the first example of N-glycosylation outside the Eukarya, only recently has attention focused on delineating the pathway responsible for the assembly of the N-linked tetrasaccharide decorating selected proteins in this haloarchaeon. In the present report, the roles of VNG1053G and VNG1054G, two proteins encoded by genes clustered together with a set of genes demonstrated to encode N-glycosylation pathway components, were considered. Relying on both bioinformatics and gene deletion and subsequent mass spectrometry analysis of known N-glycosylated proteins, VNG1053G was determined to be the glycosyltransferase responsible for addition of the linking glucose, while VNG1054G was deemed to be the flippase that translocates the lipid-bound tetrasaccharide across the plasma membrane to face the cell exterior, or to contribute to such activity. As observed with Hbt. salinarum lacking other components of the N-glycosylation machinery, both cell growth and motility were compromised in the absence of VNG1053G or VNG1054G. Thus, given their demonstrated roles in Hbt. salinarum N-glycosylation, VNG1053G and VNG1054G were re-annotated as Agl28 and Agl29, according to the nomenclature used to define archaeal N-glycosylation pathway components.

An N-linked tetrasaccharide from Halobacterium salinarum presents a novel modification, sulfation of iduronic acid at the O-3 position

Halobacterium salinarum, a halophilic archaeon that grows at near-saturating salt concentrations, provided the first example of N-glycosylation outside Eukarya. Yet, almost 50 years later, numerous aspects of such post-translational protein processing in this microorganism remain to be determined, including the architecture of glycoprotein-bound glycans. In the present report, nuclear magnetic resonance spectroscopy was used to define a tetrasaccharide N-linked to both archaellins, building blocks of the archaeal swimming device (the archaellum), and the S-layer glycoprotein that comprises the protein shell surrounding the Hbt. salinarum cell as β-GlcA(2S)-(1 → 4)-α-IdoA(3S)-(1 → 4)-β-GlcA-(1 → 4)-β-Glc-Asn. The structure of this tetrasaccharide fills gaps remaining from previous studies, including confirmation of the first known inclusion of iduronic acid in an archaeal N-linked glycan. At the same time, the sulfation of this iduronic acid at the O-3 position has not, to the best of our knowledge, been previously seen. As such, this may represent yet another unique facet of N-glycosylation in Archaea.

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.

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.

AMiGA: Software for Automated Analysis of Microbial Growth Assays

The analysis of microbial growth is one of the central methods in the field of microbiology. Microbial growth dynamics can be characterized by meaningful parameters, including carrying capacity, exponential growth rate, and growth lag. However, microbial assays with clinical isolates, fastidious organisms, or microbes under stress often produce atypical growth shapes that do not follow the classical microbial growth pattern. Here, we introduce the analysis of microbial growth assays (AMiGA) software, which streamlines the analysis of growth curves without any assumptions about their shapes. AMiGA can pool replicates of growth curves and infer summary statistics for biologically meaningful growth parameters. In addition, AMiGA can quantify death phases and characterize diauxic shifts. It can also statistically test for differential growth under distinct experimental conditions. Altogether, AMiGA streamlines the organization, analysis, and visualization of microbial growth assays. IMPORTANCE Our current understanding of microbial physiology relies on the simple method of measuring microbial populations' sizes over time and under different conditions. Many advances have increased the throughput of those assays and enabled the study of nonlab-adapted microbes under diverse conditions that widely affect their growth dynamics. Our software provides an all-in-one tool for estimating the growth parameters of microbial cultures and testing for differential growth in a high-throughput and user-friendly fashion without any underlying assumptions about how microbes respond to their growth conditions.

CopR, a Global Regulator of Transcription to Maintain Copper Homeostasis in Pyrococcus furiosus

Although copper is in many cases an essential micronutrient for cellular life, higher concentrations are toxic. Therefore, all living cells have developed strategies to maintain copper homeostasis. In this manuscript, we have analyzed the transcriptome-wide response of Pyrococcus furiosus to increased copper concentrations and described the essential role of the putative copper-sensing metalloregulator CopR in the detoxification process. To this end, we employed biochemical and biophysical methods to characterize the role of CopR. Additionally, a copR knockout strain revealed an amplified sensitivity in comparison to the parental strain towards increased copper levels, which designates an essential role of CopR for copper homeostasis. To learn more about the CopR-regulated gene network, we performed differential gene expression and ChIP-seq analysis under normal and 20 μM copper-shock conditions. By integrating the transcriptome and genome-wide binding data, we found that CopR binds to the upstream regions of many copper-induced genes. Negative-stain transmission electron microscopy and 2D class averaging revealed an octameric assembly formed from a tetramer of dimers for CopR, similar to published crystal structures from the Lrp family. In conclusion, we propose a model for CopR-regulated transcription and highlight the regulatory network that enables Pyrococcus to respond to increased copper concentrations.

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.

Archaeal transcription

Increasingly sophisticated biochemical and genetic techniques are unraveling the regulatory factors and mechanisms that control gene expression in the Archaea. While some similarities in regulatory strategies are universal, archaeal-specific regulatory strategies are emerging to complement a complex patchwork of shared archaeal-bacterial and archaeal-eukaryotic regulatory mechanisms employed in the archaeal domain. The prokaryotic archaea encode core transcription components with homology to the eukaryotic transcription apparatus and also share a simplified eukaryotic-like initiation mechanism, but also deploy tactics common to bacterial systems to regulate promoter usage and influence elongation-termination decisions. We review the recently established complete archaeal transcription cycle, highlight recent findings of the archaeal transcription community and detail the expanding post-initiation regulation imposed on archaeal transcription.

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.

Gene Expression of Haloferax volcanii on Intermediate and Abundant Sources of Fixed Nitrogen

Haloferax volcanii, a well-developed model archaeon for genomic, transcriptomic, and proteomic analyses, can grow on a defined medium of abundant and intermediate levels of fixed nitrogen. Here we report a global profiling of gene expression of H. volcanii grown on ammonium as an abundant source of fixed nitrogen compared to l-alanine, the latter of which exemplifies an intermediate source of nitrogen that can be obtained from dead cells in natural habitats. By comparing the two growth conditions, 30 genes were found to be differentially expressed, including 16 genes associated with amino acid metabolism and transport. The gene expression profiles contributed to mapping ammonium and l-alanine usage with respect to transporters and metabolic pathways. In addition, conserved DNA motifs were identified in the putative promoter regions and transcription factors were found to be in synteny with the differentially expressed genes, leading us to propose regulons of transcriptionally co-regulated operons. This study provides insight to how H. volcanii responds to and utilizes intermediate vs. abundant sources of fixed nitrogen for growth, with implications for conserved functions in related halophilic archaea.

Global Transcriptional Programs in Archaea Share Features with the Eukaryotic Environmental Stress Response

The environmental stress response (ESR), a global transcriptional program originally identified in yeast, is characterized by a rapid and transient transcriptional response composed of large, oppositely regulated gene clusters. Genes induced during the ESR encode core components of stress tolerance, macromolecular repair, and maintenance of homeostasis. In this review, we investigate the possibility for conservation of the ESR across the eukaryotic and archaeal domains of life. We first re-analyze existing transcriptomics data sets to illustrate that a similar transcriptional response is identifiable in Halobacterium salinarum, an archaeal model organism. To substantiate the archaeal ESR, we calculated gene-by-gene correlations, gene function enrichment, and comparison of temporal dynamics. We note reported examples of variation in the ESR across fungi, then synthesize high-level trends present in expression data of other archaeal species. In particular, we emphasize the need for additional high-throughput time series expression data to further characterize stress-responsive transcriptional programs in the Archaea. Together, this review explores an open question regarding features of global transcriptional stress response programs shared across domains of life.

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.

Characterization of the transcriptome of Haloferax volcanii, grown under four different conditions, with mixed RNA-Seq

Haloferax volcanii is a well-established model species for haloarchaea. Small scale RNomics and bioinformatics predictions were used to identify small non-coding RNAs (sRNAs), and deletion mutants revealed that sRNAs have important regulatory functions. A recent dRNA-Seq study was used to characterize the primary transcriptome. Unexpectedly, it was revealed that, under optimal conditions, H. volcanii contains more non-coding sRNAs than protein-encoding mRNAs. However, the dRNA-Seq approach did not contain any length information. Therefore, a mixed RNA-Seq approach was used to determine transcript length and to identify additional transcripts, which are not present under optimal conditions. In total, 50 million paired end reads of 150 nt length were obtained. 1861 protein-coding RNAs (cdRNAs) were detected, which encoded 3092 proteins. This nearly doubled the coverage of cdRNAs, compared to the previous dRNA-Seq study. About 2/3 of the cdRNAs were monocistronic, and 1/3 covered more than one gene. In addition, 1635 non-coding sRNAs were identified. The highest fraction of non-coding RNAs were cis antisense RNAs (asRNAs). Analysis of the length distribution revealed that sRNAs have a median length of about 150 nt. Based on the RNA-Seq and dRNA-Seq results, genes were chosen to exemplify characteristics of the H. volcanii transcriptome by Northern blot analyses, e.g. 1) the transcript patterns of gene clusters can be straightforward, but also very complex, 2) many transcripts differ in expression level under the four analyzed conditions, 3) some genes are transcribed into RNA isoforms of different length, which can be differentially regulated, 4) transcripts with very long 5'-UTRs and with very long 3'-UTRs exist, and 5) about 30% of all cdRNAs have overlapping 3'-ends, which indicates, together with the asRNAs, that H. volcanii makes ample use of sense-antisense interactions. Taken together, this RNA-Seq study, together with a previous dRNA-Seq study, enabled an unprecedented view on the H. volcanii transcriptome.

Synergistic Impacts of Organic Acids and pH on Growth of Pseudomonas aeruginosa: A Comparison of Parametric and Bayesian Non-parametric Methods to Model Growth

Different weak organic acids have significant potential as topical treatments for wounds infected by opportunistic pathogens that are recalcitrant to standard treatments. These acids have long been used as bacteriostatic compounds in the food industry, and in some cases are already being used in the clinic. The effects of different organic acids vary with pH, concentration, and the specific organic acid used, but no studies to date on any opportunistic pathogens have examined the detailed interactions between these key variables in a controlled and systematic way. We have therefore comprehensively evaluated the effects of several different weak organic acids on growth of the opportunistic pathogen Pseudomonas aeruginosa. We used a semi-automated plate reader to generate growth profiles for two different strains (model laboratory strain PAO1 and clinical isolate PA1054 from a hospital burns unit) in a range of organic acids at different concentrations and pH, with a high level of replication for a total of 162,960 data points. We then compared two different modeling approaches for the interpretation of this time-resolved dataset: parametric logistic regression (with or without a component to include lag phase) vs. non-parametric Gaussian process (GP) regression. Because GP makes no prior assumptions about the nature of the growth, this method proved to be superior in cases where growth did not follow a standard sigmoid functional form, as is common when bacteria grow under stress. Acetic, propionic and butyric acids were all more detrimental to growth than the other acids tested, and although PA1054 grew better than PAO1 under non-stress conditions, this difference largely disappeared as the levels of stress increased. As expected from knowledge of how organic acids behave, their effect was significantly enhanced in combination with low pH, with this interaction being greatest in the case of propionic acid. Our approach lends itself to the characterization of combinatorial interactions between stressors, especially in cases where their impacts on growth render logistic growth models unsuitable.

Stochasticity in transcriptional expression of a negative regulator of Arabidopsis ABA network

Stably heritable spatiotemporal co/over-expression of distinct transcriptional regulators across generations is a desired target as they signal traffic in the cell. Here, the stability and expression pattern of AtHB7 (Arabidopsis homeodomain-leucine zipper class I) cDNA was characterized in 220 random population of transformed tomato clones where no AtHB7 orthologous has been identified in to date. Integration of p35S::AtHB7 casette was tested by the amplification of the stretches (700/425 bp) in the target by NPT II/AtHB7 oligos. Transcriptional expression pattern for the amplicons of the specific transcripts in the leaf tissues of transformants were determined by qRT-PCR. Transgene copy number was negatively correlated with transgene expression level, yet a majority of transformants (78%) carried single-copy of transgene. About 1:3 of the lines containing two-copy inserts showed less transcript expression. Heterologous CaMV 35S promoter drove AtHB7, illuminated no penalty on transgene expression levels, stability or plant phenotype under drought stress. Integration and expression analysis of transcription factors is of great significance for reliable prediction of gene dosing/functions in plant genomes so as to sustain breeding under abiotic stress to guarantee food security.

Heat shock response in archaea

An adequate response to a sudden temperature rise is crucial for cellular fitness and survival. While heat shock response (HSR) is well described in bacteria and eukaryotes, much less information is available for archaea, of which many characterized species are extremophiles thriving in habitats typified by large temperature gradients. Here, we describe known molecular aspects of archaeal heat shock proteins (HSPs) as key components of the protein homeostasis machinery and place this in a phylogenetic perspective with respect to bacterial and eukaryotic HSPs. Particular emphasis is placed on structure–function details of the archaeal thermosome, which is a major element of the HSR and of which subunit composition is altered in response to temperature changes. In contrast with the structural response, it is largely unclear how archaeal cells sense temperature fluctuations and which molecular mechanisms underlie the corresponding regulation. We frame this gap in knowledge by discussing emerging questions related to archaeal HSR and by proposing methodologies to address them. Additionally, as has been shown in bacteria and eukaryotes, HSR is expected to be relevant for the control of physiology and growth in various stress conditions beyond temperature stress. A better understanding of this essential cellular process in archaea will not only provide insights into the evolution of HSR and of its sensing and regulation, but also inspire the development of biotechnological applications, by enabling transfer of archaeal heat shock components to other biological systems and for the engineering of archaea as robust cell factories.

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.