Comprehensive de novo structure prediction in a systems-biology context for the archaea Halobacterium sp. NRC-1

Comparative Analysis of rRNA Removal Methods for RNA-Seq Differential Expression in Halophilic Archaea

Despite intense recent research interest in archaea, the scientific community has experienced a bottleneck in the study of genome-scale gene expression experiments by RNA-seq due to the lack of commercial and specifically designed rRNA depletion kits. The high rRNA:mRNA ratio (80–90%: ~10%) in prokaryotes hampers global transcriptomic analysis. Insufficient ribodepletion results in low sequence coverage of mRNA, and therefore, requires a substantially higher number of replicate samples and/or sequencing reads to achieve statistically reliable conclusions regarding the significance of differential gene expression between case and control samples. Here, we show that after the discontinuation of the previous version of RiboZero (Illumina, San Diego, CA, USA) that was useful in partially or completely depleting rRNA from archaea, archaeal transcriptomics studies have experienced a slowdown. To overcome this limitation, here, we analyze the efficiency for four different hybridization-based kits from three different commercial suppliers, each with two sets of sequence-specific probes to remove rRNA from four different species of halophilic archaea. We conclude that the key for transcriptomic success with the currently available tools is the probe-specificity for the rRNA sequence hybridization. With this paper, we provide insights into the archaeal community for selecting certain reagents and strategies over others depending on the archaeal species of interest. These methods yield improved RNA-seq sensitivity and enhanced detection of low abundance transcripts.

Bacterial 'Grounded' Prophages: Hotspots for Genetic Renovation and Innovation

Bacterial genomes are highly plastic allowing the generation of variants through mutations and acquisition of genetic information. The fittest variants are then selected by the econiche thereby allowing the bacterial adaptation and colonization of the habitat. Larger genomes, however, may impose metabolic burden and hence bacterial genomes are optimized by the loss of frivolous genetic information. The activity of temperate bacteriophages has acute consequences on the bacterial population as well as the bacterial genome through lytic and lysogenic cycles. Lysogeny is a selective advantage as the prophage provides immunity to the lysogen against secondary phage attack. Since the non-lysogens are eliminated by the lytic phages, lysogens multiply and colonize the habitat. Nevertheless, all lysogens have an imminent risk of lytic cycle activation and cell lysis. However, a mutation in the attachment sites or in the genes that encode the specific recombinase responsible for prophage excision could result in 'grounding' of the prophage. Since the lysogens with grounded prophage are immune to respective phage infection as well as dodge the induction of lytic cycle, we hypothesize that the selection of these mutant lysogens is favored relative to their normal lysogenic counterparts. These grounded prophages offer several advantages to the bacterial genome evolution through propensity for genetic variations including inversions, deletions, and insertions via horizontal gene transfer. We propose that the grounded prophages expedite bacterial genome evolution by acting as 'genetic buffer zones' thereby increasing the frequency as well as the diversity of variations on which natural selection favors the beneficial variants. The grounded prophages are also hotspots for horizontal gene transfer wherein several ecologically significant genes such as those involved in stress tolerance, antimicrobial resistance, and novel metabolic pathways, are integrated. Moreover, the high frequency of genetic changes within prophages also allows proportionate probability for the de novo genesis of genetic information. Through sequence analyses of well-characterized E. coli prophages we exemplify various roles of grounded prophages in E. coli ecology and evolution. Therefore, the temperate prophages are one of the most significant drivers of bacterial genome evolution and sites of biogenesis of genetic information.

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

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

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

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.

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 role for programmed cell death in the microbial loop

The microbial loop is the conventional model by which nutrients and minerals are recycled in aquatic eco-systems. Biochemical pathways in different organisms become metabolically inter-connected such that nutrients are utilized, processed, released and re-utilized by others. The result is that unrelated individuals end up impacting each others' fitness directly through their metabolic activities. This study focused on the impact of programmed cell death (PCD) on a population's growth as well as its role in the exchange of carbon between two naturally co-occurring halophilic organisms. Flow cytometric, biochemical, ¹⁴C radioisotope tracing assays, and global transcriptomic analyses show that organic algal photosynthate released by Dunalliela salina cells undergoing PCD complements the nutritional needs of other non-PCD D. salina cells. This occurs in vitro in a carbon limited environment and enhances the growth of the population. In addition, a co-occurring heterotroph Halobacterium salinarum re-mineralizes the carbon providing elemental nutrients for the mixoheterotrophic chlorophyte. The significance of this is uncertain and the archaeon can also subsist entirely on the lysate of apoptotic algae. PCD is now well established in unicellular organisms; however its ecological relevance has been difficult to decipher. In this study we found that PCD in D. salina causes the release of organic nutrients such as glycerol, which can be used by others in the population as well as a co-occurring halophilic archaeon. H. salinarum also re-mineralizes the dissolved material promoting algal growth. PCD in D. salina was the mechanism for the flow of dissolved photosynthate between unrelated organisms. Ironically, programmed death plays a central role in an organism's own population growth and in the exchange of nutrients in the microbial loop.

The Proteome Folding Project: proteome-scale prediction of structure and function

The incompleteness of proteome structure and function annotation is a critical problem for biologists and, in particular, severely limits interpretation of high-throughput and next-generation experiments. We have developed a proteome annotation pipeline based on structure prediction, where function and structure annotations are generated using an integration of sequence comparison, fold recognition, and grid-computing-enabled de novo structure prediction. We predict protein domain boundaries and three-dimensional (3D) structures for protein domains from 94 genomes (including human, Arabidopsis, rice, mouse, fly, yeast, Escherichia coli, and worm). De novo structure predictions were distributed on a grid of more than 1.5 million CPUs worldwide (World Community Grid). We generated significant numbers of new confident fold annotations (9% of domains that are otherwise unannotated in these genomes). We demonstrate that predicted structures can be combined with annotations from the Gene Ontology database to predict new and more specific molecular functions.

Synergistic use of plant-prokaryote comparative genomics for functional annotations

BACKGROUND

Identifying functions for all gene products in all sequenced organisms is a central challenge of the post-genomic era. However, at least 30-50% of the proteins encoded by any given genome are of unknown or vaguely known function, and a large number are wrongly annotated. Many of these 'unknown' proteins are common to prokaryotes and plants. We set out to predict and experimentally test the functions of such proteins. Our approach to functional prediction integrates comparative genomics based mainly on microbial genomes with functional genomic data from model microorganisms and post-genomic data from plants. This approach bridges the gap between automated homology-based annotations and the classical gene discovery efforts of experimentalists, and is more powerful than purely computational approaches to identifying gene-function associations.

RESULTS

Among Arabidopsis genes, we focused on those (2,325 in total) that (i) are unique or belong to families with no more than three members, (ii) occur in prokaryotes, and (iii) have unknown or poorly known functions. Computer-assisted selection of promising targets for deeper analysis was based on homology-independent characteristics associated in the SEED database with the prokaryotic members of each family. In-depth comparative genomic analysis was performed for 360 top candidate families. From this pool, 78 families were connected to general areas of metabolism and, of these families, specific functional predictions were made for 41. Twenty-one predicted functions have been experimentally tested or are currently under investigation by our group in at least one prokaryotic organism (nine of them have been validated, four invalidated, and eight are in progress). Ten additional predictions have been independently validated by other groups. Discovering the function of very widespread but hitherto enigmatic proteins such as the YrdC or YgfZ families illustrates the power of our approach.

CONCLUSIONS

Our approach correctly predicted functions for 19 uncharacterized protein families from plants and prokaryotes; none of these functions had previously been correctly predicted by computational methods. The resulting annotations could be propagated with confidence to over six thousand homologous proteins encoded in over 900 bacterial, archaeal, and eukaryotic genomes currently available in public databases.

Ab initio modeling led annotation suggests nucleic acid binding function for many DUFs

Expansions of sequence databases driven by new sequencing technology continue apace. These result in a continuous supply of protein sequences and domains that cannot be straightforwardly annotated by simple homology methods. For these, structure-based function prediction may contribute to an improved annotation. Here, short Domains of Unknown Function (DUFs) are ab initio modeled with ROSETTA and screened for likely nucleic acid binding function. Thirty-two DUFs are thereby predicted to have a nucleic acid binding function. In most cases, additional evidence supporting that function could be obtained from structure comparison, domain architectures, distant evolutionary relationships, genome context or protein-protein interaction data. These predictions contribute to the function annotation of thousands of proteins.

Coordination of frontline defense mechanisms under severe oxidative stress

Inference of an environmental and gene regulatory influence network (EGRIN_OS) by integrating transcriptional responses to H₂O₂ and paraquat (PQ) has revealed a multi-tiered oxidative stress (OS)-management program to transcriptionally coordinate three peroxidase/catalase enzymes, two superoxide dismutases, production of rhodopsins, carotenoids and gas vesicles, metal trafficking, and various other aspects of metabolism.

ChIP-chip, microarray, and survival assays have validated important architectural aspects of this network, identified novel defense mechanisms (including two evolutionarily distant peroxidase enxymes), and showed that general transcription factors of the transcription factor B family have an important function in coordinating the OS response (OSR) despite their inability to directly sense ROS.

A comparison of transcriptional responses to sub-lethal doses of H₂O₂ and PQ with predictions of these responses made by an EGRIN model generated earlier from responses to other environmental factors has confirmed that a significant fraction of the OSR is made up of a generalized component that is also observed in response to other stressors.

Analysis of active regulons within environment and gene regulatory influence network for OS (EGRIN_OS) across diverse environmental conditions has identified the specialized component of oxidative stress response (OSR) that is triggered by sub-lethal OS, but not by other stressors, including sub-inhibitory levels of redox-active metals, extreme changes in oxygen tension, and a sub-lethal dose of γ rays.

Reactive oxygen species (ROS), such as hydrogen peroxide (H₂O₂), superoxide (O₂⁻), and hydroxyl (OH⁻) radicals, are normal by-products of aerobic metabolism. Evolutionarily conserved mechanisms including detoxification enzymes (peroxidase/catalase and superoxide dismutase (SOD)) and free radical scavengers manage this endogenous production of ROS. OS is a condition reached when certain environmental stresses or genetic defects cause the production of ROS to exceed the management capacity. The damage to diverse cellular components including DNA, proteins, lipids, and carbohydrates resulting from OS (Imlay, 2003; Apel and Hirt, 2004; Perrone et al, 2008) is recognized as an important player in many diseases and in the aging process (Finkel, 2005).

We have applied a systems approach to characterize the OSR of an archaeal model organism, Halobacterium salinarum NRC-1. This haloarchaeon grows aerobically at 4.3 M salt concentration in which it routinely faces cycles of desiccation and rehydration, and increased ultraviolet radiation—both of which can increase the production of ROS (Farr and Kogoma, 1991; Oliver et al, 2001). We have reconstructed the physiological adjustments associated with management of excessive OS through the analysis of global transcriptional changes elicited by step exposure to growth sub-inhibitory and sub-lethal levels of H₂O₂ and PQ (a redox-cycling drug that produces O₂⁻; Hassan and Fridovich, 1979) as well as during subsequent recovery from these stresses. We have integrated all of these data into a unified model for OSR to discover conditional functional links between protective mechanisms and normal aspects of metabolism. Subsequent phenotypic analysis of gene deletion strains has verified the conditional detoxification functions of three putative peroxidase/catalase enzymes, two SODs, and the protective function of rhodopsins under increased levels of H₂O₂ and PQ. Similarly, we have also validated ROS scavenging by carotenoids and flotation by gas vesicles as secondary mechanisms that may minimize OS.

Given the ubiquitous nature of OS, it is not entirely surprising that most organisms have evolved similar multiple lines of defense—both passive and active. Although such mechanisms have been extensively characterized using other model organisms, our integrated systems approach has uncovered additional protective mechanisms in H. salinarum (e.g. two evolutionarily distant peroxidase/catalase enzymes) and revealed a structure and hierarchy to the OSR through conditional regulatory associations among various components of the response. We have validated some aspects of the architecture of the regulatory network for managing OS by confirming physical protein–DNA interactions of six transcription factors (TFs) with promoters of genes they were predicted to influence in EGRIN_OS. Furthermore, we have also shown the consequence of deleting two of these TFs on transcript levels of genes they control and survival rate under OS. It is notable that these TFs are not directly associated with sensing ROS, but, rather, they have a general function in coordinating the overall response. This insight would not have been possible without constructing EGRIN_OS through systems integration of diverse datasets.

Although it has been known that OS is a component of diverse environmental stress conditions, we quantitatively show for the first time that much of the transcriptional responses induced by the two treatments could indeed have been predicted using a model constructed from the analysis of transcriptional responses to changes in other environmental factors (UV and γ-radiation, light, oxygen, and six metals). However, using specific examples we also reveal the specific components of the OSR that are triggered only under severe OS. Notably, this model of OSR gives a unified perspective of the interconnections among all of these generalized and OS-specific regulatory mechanisms.

Complexity of cellular response to oxidative stress (OS) stems from its wide-ranging damage to nucleic acids, proteins, carbohydrates, and lipids. We have constructed a systems model of OS response (OSR) for Halobacterium salinarum NRC-1 in an attempt to understand the architecture of its regulatory network that coordinates this complex response. This has revealed a multi-tiered OS-management program to transcriptionally coordinate three peroxidase/catalase enzymes, two superoxide dismutases, production of rhodopsins, carotenoids and gas vesicles, metal trafficking, and various other aspects of metabolism. Through experimental validation of interactions within the OSR regulatory network, we show that despite their inability to directly sense reactive oxygen species, general transcription factors have an important function in coordinating this response. Remarkably, a significant fraction of this OSR was accurately recapitulated by a model that was earlier constructed from cellular responses to diverse environmental perturbations—this constitutes the general stress response component. Notwithstanding this observation, comparison of the two models has identified the coordination of frontline defense and repair systems by regulatory mechanisms that are triggered uniquely by severe OS and not by other environmental stressors, including sub-inhibitory levels of redox-active metals, extreme changes in oxygen tension, and a sub-lethal dose of γ rays.

Two transcription factors are necessary for iron homeostasis in a salt-dwelling archaeon

Because iron toxicity and deficiency are equally life threatening, maintaining intracellular iron levels within a narrow optimal range is critical for nearly all known organisms. However, regulatory mechanisms that establish homeostasis are not well understood in organisms that dwell in environments at the extremes of pH, temperature, and salinity. Under conditions of limited iron, the extremophile Halobacterium salinarum, a salt-loving archaeon, mounts a specific response to scavenge iron for growth. We have identified and characterized the role of two transcription factors (TFs), Idr1 and Idr2, in regulating this important response. An integrated systems analysis of TF knockout gene expression profiles and genome-wide binding locations in the presence and absence of iron has revealed that these TFs operate collaboratively to maintain iron homeostasis. In the presence of iron, Idr1 and Idr2 bind near each other at 24 loci in the genome, where they are both required to repress some genes. By contrast, Idr1 and Idr2 are both necessary to activate other genes in a putative a feed forward loop. Even at loci bound independently, the two TFs target different genes with similar functions in iron homeostasis. We discuss conserved and unique features of the Idr1-Idr2 system in the context of similar systems in organisms from other domains of life.

Innate immune detection of the type III secretion apparatus through the NLRC4 inflammasome

The mammalian innate immune system uses Toll-like receptors (TLRs) and Nod-LRRs (NLRs) to detect microbial components during infection. Often these molecules work in concert; for example, the TLRs can stimulate the production of the proforms of the cytokines IL-1beta and IL-18, whereas certain NLRs trigger their subsequent proteolytic processing via caspase 1. Gram-negative bacteria use type III secretion systems (T3SS) to deliver virulence factors to the cytosol of host cells, where they modulate cell physiology to favor the pathogen. We show here that NLRC4/Ipaf detects the basal body rod component of the T3SS apparatus (rod protein) from S. typhimurium (PrgJ), Burkholderia pseudomallei (BsaK), Escherichia coli (EprJ and EscI), Shigella flexneri (MxiI), and Pseudomonas aeruginosa (PscI). These rod proteins share a sequence motif that is essential for detection by NLRC4; a similar motif is found in flagellin that is also detected by NLRC4. S. typhimurium has two T3SS: Salmonella pathogenicity island-1 (SPI1), which encodes the rod protein PrgJ, and SPI2, which encodes the rod protein SsaI. Although PrgJ is detected by NLRC4, SsaI is not, and this evasion is required for virulence in mice. The detection of a conserved component of the T3SS apparatus enables innate immune responses to virulent bacteria through a single pathway, a strategy that is divergent from that used by plants in which multiple NB-LRR proteins are used to detect T3SS effectors or their effects on cells. Furthermore, the specific detection of the virulence machinery permits the discrimination between pathogenic and nonpathogenic bacteria.

A single transcription factor regulates evolutionarily diverse but functionally linked metabolic pathways in response to nutrient availability

During evolution, enzyme-coding genes are acquired and/or replaced through lateral gene transfer and compiled into metabolic pathways. Gene regulatory networks evolve to fine tune biochemical fluxes through such metabolic pathways, enabling organisms to acclimate to nutrient fluctuations in a competitive environment. Here, we demonstrate that a single TrmB family transcription factor in Halobacterium salinarum NRC-1 globally coordinates functionally linked enzymes of diverse phylogeny in response to changes in carbon source availability. Specifically, during nutritional limitation, TrmB binds a cis-regulatory element to activate or repress 113 promoters of genes encoding enzymes in diverse metabolic pathways. By this mechanism, TrmB coordinates the expression of glycolysis, TCA cycle, and amino-acid biosynthesis pathways with the biosynthesis of their cognate cofactors (e.g. purine and thiamine). Notably, the TrmB-regulated metabolic network includes enzyme-coding genes that are uniquely archaeal as well as those that are conserved across all three domains of life. Simultaneous analysis of metabolic and gene regulatory network architectures suggests an ongoing process of co-evolution in which TrmB integrates the expression of metabolic enzyme-coding genes of diverse origins.

Fold prediction of VP24 protein of Ebola and Marburg viruses using de novo fragment assembly

Virus particle 24 (VP24) is the smallest protein of the Ebola and Marburg virus genomes. Recent experiments show that Ebola VP24 blocks binding of tyrosine-phosphorylated STAT-1 homodimer (PY-STAT1) to the NPI-1 subfamily of importin alpha, thereby preventing nuclear accumulation of this interferon-promoting transcription factor which, in turn, reduces the innate immune response of the host target. Lacking an experimental structure for VP24, we applied de novo protein structure prediction using the fragment assembly-based Rosetta method to classify its fold topology and better understand its biological function. Filtering and ranking of models were performed with the DFIRE all-atom statistical potential and the CHARMM22 force field with a generalized Born solvent model. From 40,000 Rosetta-generated structures and selective comparisons with the SCOP database, a structural match to two of our top 10-ranking models was the Armadillo repeat fold topology. Specific members of this fold family include importin alpha, importin beta, and exportin. We propose that, unlike the nuclear import of host cargo, VP24 lacks a classical nuclear localization signal (NLS) and targets importin alpha in a similar manner to the observed heterodimeric complex with exportin, thereby interfering with the auto-inhibitory NLS on importin alpha and blocking peripheral docking sites for PY-STAT1 assembly.

Halobacterium salinarum NRC-1 PeptideAtlas: toward strategies for targeted proteomics and improved proteome coverage

The relatively small numbers of proteins and fewer possible post-translational modifications in microbes provide a unique opportunity to comprehensively characterize their dynamic proteomes. We have constructed a PeptideAtlas (PA) covering 62.7% of the predicted proteome of the extremely halophilic archaeon Halobacterium salinarum NRC-1 by compiling approximately 636 000 tandem mass spectra from 497 mass spectrometry runs in 88 experiments. Analysis of the PA with respect to biophysical properties of constituent peptides, functional properties of parent proteins of detected peptides, and performance of different mass spectrometry approaches has highlighted plausible strategies for improving proteome coverage and selecting signature peptides for targeted proteomics. Notably, discovery of a significant correlation between absolute abundances of mRNAs and proteins has helped identify low abundance of proteins as the major limitation in peptide detection. Furthermore, we have discovered that iTRAQ labeling for quantitative proteomic analysis introduces a significant bias in peptide detection by mass spectrometry. Therefore, despite identifying at least one proteotypic peptide for almost all proteins in the PA, a context-dependent selection of proteotypic peptides appears to be the most effective approach for targeted proteomics.

A predictive model for transcriptional control of physiology in a free living cell

The environment significantly influences the dynamic expression and assembly of all components encoded in the genome of an organism into functional biological networks. We have constructed a model for this process in Halobacterium salinarum NRC-1 through the data-driven discovery of regulatory and functional interrelationships among approximately 80% of its genes and key abiotic factors in its hypersaline environment. Using relative changes in 72 transcription factors and 9 environmental factors (EFs) this model accurately predicts dynamic transcriptional responses of all these genes in 147 newly collected experiments representing completely novel genetic backgrounds and environments-suggesting a remarkable degree of network completeness. Using this model we have constructed and tested hypotheses critical to this organism's interaction with its changing hypersaline environment. This study supports the claim that the high degree of connectivity within biological and EF networks will enable the construction of similar models for any organism from relatively modest numbers of experiments.

A personal journey of discovery: developing technology and changing biology

This autobiographical article describes my experiences in developing chemically based, biological technologies for deciphering biological information: DNA, RNA, proteins, interactions, and networks. The instruments developed include protein and DNA sequencers and synthesizers, as well as ink-jet technology for synthesizing DNA chips. Diverse new strategies for doing biology also arose from novel applications of these instruments. The functioning of these instruments can be integrated to generate powerful new approaches to cloning and characterizing genes from a small amount of protein sequence or to using gene sequences to synthesize peptide fragments so as to characterize various properties of the proteins. I also discuss the five paradigm changes in which I have participated: the development and integration of biological instrumentation; the human genome project; cross-disciplinary biology; systems biology; and predictive, personalized, preventive, and participatory (P4) medicine. Finally, I discuss the origins, the philosophy, some accomplishments, and the future trajectories of the Institute for Systems Biology.

Superfamily assignments for the yeast proteome through integration of structure prediction with the gene ontology

Saccharomyces cerevisiae is one of the best-studied model organisms, yet the three-dimensional structure and molecular function of many yeast proteins remain unknown. Yeast proteins were parsed into 14,934 domains, and those lacking sequence similarity to proteins of known structure were folded using the Rosetta de novo structure prediction method on the World Community Grid. This structural data was integrated with process, component, and function annotations from the Saccharomyces Genome Database to assign yeast protein domains to SCOP superfamilies using a simple Bayesian approach. We have predicted the structure of 3,338 putative domains and assigned SCOP superfamily annotations to 581 of them. We have also assigned structural annotations to 7,094 predicted domains based on fold recognition and homology modeling methods. The domain predictions and structural information are available in an online database at http://rd.plos.org/10.1371_journal.pbio.0050076_01.

The anatomy of microbial cell state transitions in response to oxygen

Adjustment of physiology in response to changes in oxygen availability is critical for the survival of all organisms. However, the chronology of events and the regulatory processes that determine how and when changes in environmental oxygen tension result in an appropriate cellular response is not well understood at a systems level. Therefore, transcriptome, proteome, ATP, and growth changes were analyzed in a halophilic archaeon to generate a temporal model that describes the cellular events that drive the transition between the organism's two opposing cell states of anoxic quiescence and aerobic growth. According to this model, upon oxygen influx, an initial burst of protein synthesis precedes ATP and transcription induction, rapidly driving the cell out of anoxic quiescence, culminating in the resumption of growth. This model also suggests that quiescent cells appear to remain actively poised for energy production from a variety of different sources. Dynamic temporal analysis of relationships between transcription and translation of key genes suggests several important mechanisms for cellular sustenance under anoxia as well as specific instances of post-transcriptional regulation.

Discovering functional linkages and uncharacterized cellular pathways using phylogenetic profile comparisons: a comprehensive assessment

BACKGROUND

A widely-used approach for discovering functional and physical interactions among proteins involves phylogenetic profile comparisons (PPCs). Here, proteins with similar profiles are inferred to be functionally related under the assumption that proteins involved in the same metabolic pathway or cellular system are likely to have been co-inherited during evolution.

RESULTS

Our experimentation with E. coli and yeast proteins with 16 different carefully composed reference sets of genomes revealed that the phyletic patterns of proteins in prokaryotes alone could be adequate enough to make reasonably accurate functional linkage predictions. A slight improvement in performance is observed on adding few eukaryotes into the reference set, but a noticeable drop-off in performance is observed with increased number of eukaryotes. Inclusion of most parasitic, pathogenic or vertebrate genomes and multiple strains of the same species into the reference set do not necessarily contribute to an improved sensitivity or accuracy. Interestingly, we also found that evolutionary histories of individual pathways have a significant affect on the performance of the PPC approach with respect to a particular reference set. For example, to accurately predict functional links in carbohydrate or lipid metabolism, a reference set solely composed of prokaryotic (or bacterial) genomes performed among the best compared to one composed of genomes from all three super-kingdoms; this is in contrast to predicting functional links in translation for which a reference set composed of prokaryotic (or bacterial) genomes performed the worst. We also demonstrate that the widely used random null model to quantify the statistical significance of profile similarity is incomplete, which could result in an increased number of false-positives.

CONCLUSION

Contrary to previous proposals, it is not merely the number of genomes but a careful selection of informative genomes in the reference set that influences the prediction accuracy of the PPC approach. We note that the predictive power of the PPC approach, especially in eukaryotes, is heavily influenced by the primary endosymbiosis and subsequent bacterial contributions. The over-representation of parasitic unicellular eukaryotes and vertebrates additionally make eukaryotes less useful in the reference sets. Reference sets composed of highly non-redundant set of genomes from all three super-kingdoms fare better with pathways showing considerable vertical inheritance and strong conservation (e.g. translation apparatus), while reference sets solely composed of prokaryotic genomes fare better for more variable pathways like carbohydrate metabolism. Differential performance of the PPC approach on various pathways, and a weak positive correlation between functional and profile similarities suggest that caution should be exercised while interpreting functional linkages inferred from genome-wide large-scale profile comparisons using a single reference set.

The Rab6 effector Bicaudal D1 associates with Chlamydia trachomatis inclusions in a biovar-specific manner

Chlamydia species are obligate intracellular bacteria that replicate within a membrane-bound vacuole, the inclusion, which is trafficked to the peri-Golgi region by processes that are dependent on early chlamydial gene expression. Although neither the host nor the chlamydial proteins that regulate the intracellular trafficking have been clearly defined, several enhanced green fluorescent protein (EGFP)-tagged Rab GTPases, including Rab6, are recruited to Chlamydia trachomatis inclusions. To further characterize the association of Rab6 with C. trachomatis inclusions, we examined the intracellular localization of guanine nucleotide-binding mutants of Rab6 and demonstrated that only active GTP-bound and not inactive GDP-bound EGFP-Rab6 mutants were recruited to the inclusion, suggesting that EGFP-Rab6 interacts with the inclusion via a host Rab6 effector or a chlamydial protein that mimics a Rab6 effector. Using EGFP-tagged fusion proteins, we also demonstrated that the Rab6 effector Bicaudal D1 (BICD1) localized to C. trachomatis inclusions in a biovar-specific manner. In addition, we demonstrated that EGFP-Rab6 and its effector EGFP-BICD1 are recruited to the inclusion in a microtubule- and Golgi apparatus-independent but chlamydial gene expression-dependent mechanism. Finally, in contrast to the Rab6-dependent Golgi apparatus localization of endogenous BICD1, EGFP-BICD1 was recruited to the inclusion by a Rab6-independent mechanism. Collectively, these data demonstrate that neither Rab6 nor BICD1 is trafficked to the inclusion via a Golgi apparatus-localized intermediate, suggesting that each protein is trafficked to the C. trachomatis serovar L2 inclusion by a unique, but as-yet-undefined, mechanism.

Systems biology in drug safety and metabolism: integration of microarray, real-time PCR and enzyme approaches

The last decade has seen a rapid expansion in the field of functional genomics, due mainly to the global gene expression profiling capabilities provided by techniques, such as microarray analysis. Application of these technologies in fields as diverse as plant research, to public health and environmental sciences, forensic science and drug research, shows the versatility of these tools and the promise they hold for revolutionizing research in the life sciences. In drug discovery, attempts have been made to use functional genomics in target identification and validation, lead selection and optimization, and in preclinical studies to predict clinical outcome. These studies have provided a plethora of data and undoubtedly expanded our understanding of genetic alterations in diseased and non-diseased states, but the benefits that these technologies hold have not yet been fully realized. This review discusses how a comprehensive approach to gene regulation studies, a 'systems biology' approach, is being applied in a drug development setting to address mechanism-based questions and issues raised by regulatory authorities.

A systems view of haloarchaeal strategies to withstand stress from transition metals

Given that transition metals are essential cofactors in central biological processes, misallocation of the wrong metal ion to a metalloprotein can have resounding and often detrimental effects on diverse aspects of cellular physiology. Therefore, in an attempt to characterize unique and shared responses to chemically similar metals, we have reconstructed physiological behaviors of Halobacterium NRC-1, an archaeal halophile, in sublethal levels of Mn(II), Fe(II), Co(II), Ni(II), Cu(II), and Zn(II). Over 20% of all genes responded transiently within minutes of exposure to Fe(II), perhaps reflecting immediate large-scale physiological adjustments to maintain homeostasis. At steady state, each transition metal induced growth arrest, attempts to minimize oxidative stress, toxic ion scavenging, increased protein turnover and DNA repair, and modulation of active ion transport. While several of these constitute generalized stress responses, up-regulation of active efflux of Co(II), Ni(II), Cu(II), and Zn(II), down-regulation of Mn(II) uptake and up-regulation of Fe(II) chelation, confer resistance to the respective metals. We have synthesized all of these discoveries into a unified systems-level model to provide an integrated perspective of responses to six transition metals with emphasis on experimentally verified regulatory mechanisms. Finally, through comparisons across global transcriptional responses to different metals, we provide insights into putative in vivo metal selectivity of metalloregulatory proteins and demonstrate that a systems approach can help rapidly unravel novel metabolic potential and regulatory programs of poorly studied organisms.

The Inferelator: an algorithm for learning parsimonious regulatory networks from systems-biology data sets de novo

We present a method (the Inferelator) for deriving genome-wide transcriptional regulatory interactions, and apply the method to predict a large portion of the regulatory network of the archaeon Halobacterium NRC-1. The Inferelator uses regression and variable selection to identify transcriptional influences on genes based on the integration of genome annotation and expression data. The learned network successfully predicted Halobacterium's global expression under novel perturbations with predictive power similar to that seen over training data. Several specific regulatory predictions were experimentally tested and verified.

Post-genomics of the model haloarchaeon Halobacterium sp. NRC-1

Halobacteriumsp. NRC-1 is an extremely halophilic archaeon that is easily cultured and genetically tractable. Since its genome sequence was completed in 2000, a combination of genetic, transcriptomic, proteomic, and bioinformatic approaches have provided insights into both its extremophilic lifestyle as well as fundamental cellular processes common to all life forms. Here, we review post-genomic research on this archaeon, including investigations of DNA replication and repair systems, phototrophic, anaerobic, and other physiological capabilities, acidity of the proteome for function at high salinity, and role of lateral gene transfer in its evolution.

From genomes to function: haloarchaea as model organisms

Haloarchaea are adapted to high-salt environments and accumulate equally high salt concentrations in the cytoplasm. The genomes of representatives of six haloarchaeal genera have been fully or partially sequenced, allowing the analysis of haloarchaeal propertiesin silico. Transcriptome and proteome analyses have been established forHalobacterium salinarumandHaloferax volcanii. Genetic systems are available including methods that allow the fast in-frame deletion or modification of chromosomal genes. The high-efficiency transformation system ofHf. volcaniiallows the isolation of genes essential for a biological process by complementation of loss-of-function mutants. For the analysis of haloarchaeal biology many molecular genetic, biochemical, structural and cell biological methods have been adapted to application at high salt concentrations. Recently it has become clear that several different mechanisms allow the adaptation of proteins to the high salt concentration of the cytoplasm. Taken together, the wealth of techniques available make haloarchaea excellent archaeal model species.

From replication to cultivation: hot news from Haloarchaea

Haloarchaea have developed into model organisms that are utilized to study many biological processes. Examples are the mechanisms of chromosome maintenance, gene expression and its regulation, protein export and degradation, and motility and sensing. In addition to the analysis of model species like Halobacterium salinarum and Haloferax volcanii, natural communities have been characterized. Halophilic Archaea were found in low-salt environments and are thus more widespread than previously thought.

Predicted roles for hypothetical proteins in the low-temperature expressed proteome of the Antarctic archaeon Methanococcoides burtonii

Using liquid chromatography-mass spectrometry, 528 proteins were identified that are expressed during growth at 4 degrees C in the cold adapted archaeon, Methanococcoides burtonii. Of those, 135 were annotated previously as unique or conserved hypothetical proteins. We have performed a comprehensive, integrated analysis of the latter proteins using threading, InterProScan, predicted subcellular localization and visualization of conserved gene context across multiple prokaryotic genomes. Functional information was obtained for 55 proteins, providing new insight into the physiology of M. burtonii. Many of the proteins were predicted to be involved in DNA/RNA binding or modification and cell signaling, suggesting a complex, uncharacterized regulatory network controlling cellular processes during growth at low-temperature. Novel enzymatic functions were predicted for several proteins, including a putative candidate gene for the posttranslational modification of the key methanogenesis enzyme coenzyme M methyl reductase. A bacterial-like CRISPR locus was identified as a strong candidate for archaeal-bacterial lateral gene transfer. Gene context analysis proved a valuable augmentation to the other predictive methods in several cases, by revealing conserved gene associations and annotations in other microbial genomes. Our results underscore the importance of addressing the "hypothetical protein problem" for a complete understanding of cell physiology.

The Haloferax volcanii FtsY homolog is critical for haloarchaeal growth but does not require the A domain

The targeting of many Sec substrates to the membrane-associated translocation pore requires the cytoplasmic signal recognition particle (SRP). In Eukarya and Bacteria it has been shown that membrane docking of the SRP-substrate complex occurs via the universally conserved SRP receptor (Sralpha/beta and FtsY, respectively). While much has been learned about the archaeal SRP in recent years, few studies have examined archaeal Sralpha/FtsY homologs. In the present study the FtsY homolog of Haloferax volcanii was characterized in its native host. Disruption of the sole chromosomal copy of ftsY in H. volcanii was possible only under conditions where either the full-length haloarchaeal FtsY or an amino-terminally truncated version of this protein lacking the A domain, was expressed in trans. Subcellular fractionation analysis of H. volcanii ftsY deletion strains expressing either one of the complementing proteins revealed that in addition to a cytoplasmic pool, both proteins cofractionate with the haloarchaeal cytoplasmic membrane. Moreover, membrane localization of the universally conserved SRP subunit SRP54, the key binding partner of FtsY, was detected in both H. volcanii strains. These analyses suggest that the H. volcanii FtsY homolog plays a crucial role but does not require its A domain for haloarchaeal growth.

Genome sequence of Haloarcula marismortui: a halophilic archaeon from the Dead Sea

We report the complete sequence of the 4,274,642-bp genome of Haloarcula marismortui, a halophilic archaeal isolate from the Dead Sea. The genome is organized into nine circular replicons of varying G+C compositions ranging from 54% to 62%. Comparison of the genome architectures of Halobacterium sp. NRC-1 and H. marismortui suggests a common ancestor for the two organisms and a genome of significantly reduced size in the former. Both of these halophilic archaea use the same strategy of high surface negative charge of folded proteins as means to circumvent the salting-out phenomenon in a hypersaline cytoplasm. A multitiered annotation approach, including primary sequence similarities, protein family signatures, structure prediction, and a protein function association network, has assigned putative functions for at least 58% of the 4242 predicted proteins, a far larger number than is usually achieved in most newly sequenced microorganisms. Among these assigned functions were genes encoding six opsins, 19 MCP and/or HAMP domain signal transducers, and an unusually large number of environmental response regulators-nearly five times as many as those encoded in Halobacterium sp. NRC-1--suggesting H. marismortui is significantly more physiologically capable of exploiting diverse environments. In comparing the physiologies of the two halophilic archaea, in addition to the expected extensive similarity, we discovered several differences in their metabolic strategies and physiological responses such as distinct pathways for arginine breakdown in each halophile. Finally, as expected from the larger genome, H. marismortui encodes many more functions and seems to have fewer nutritional requirements for survival than does Halobacterium sp. NRC-1.