MiST 4.0: a new release of the microbial signal transduction database, now with a metagenomic component.

Signal transduction systems in bacteria and archaea link environmental stimuli to specific adaptive cellular responses. They control gene expression, motility, biofilm formation, development and other processes that are vital to survival. The microbial signal transduction (MiST) database is an online resource that stores tens of thousands of genomes and allows users to explore their signal transduction profiles, analyze genomes in bulk using the database application programming interface (API) and make testable hypotheses about the functions of newly identified signaling systems. However, signal transduction in metagenomes remained completely unexplored. To lay the foundation for research in metagenomic signal transduction, we have prepared a new release of the MiST database, MiST 4.0, which features over 10 000 metagenome-assembled genomes (MAGs), a scaled representation of proteins and detailed BioSample information. In addition, several thousands of new genomes have been processed and stored in the database. A new interface has been developed that allows users to seamlessly switch between genomes and MAGs. MiST 4.0 is freely available at https://mistdb.com; metagenomes and MAGs can also be explored using the API available on the same page.

MiST 3.0: an updated microbial signal transduction database with an emphasis on chemosensory systems.

Bacteria and archaea employ dedicated signal transduction systems that modulate gene expression, second-messenger turnover, quorum sensing, biofilm formation, motility, host-pathogen and beneficial interactions. The updated MiST database provides a comprehensive classification of microbial signal transduction systems. This update is a result of a substantial scaling to accommodate constantly growing microbial genomic data. More than 125 000 genomes, 516 million genes and almost 100 million unique protein sequences are currently stored in the database. For each bacterial and archaeal genome, MiST 3.0 provides a complete signal transduction profile, thus facilitating theoretical and experimental studies on signal transduction and gene regulation. New software infrastructure and distributed pipeline implemented in MiST 3.0 enable regular genome updates based on the NCBI RefSeq database. A novel MiST feature is the integration of unique profile HMMs to link complex chemosensory systems with corresponding chemoreceptors in bacterial and archaeal genomes. The data can be explored online or via RESTful API (freely available at https://mistdb.com).

SeqDepot: streamlined database of biological sequences and precomputed features.

UNLABELLED: Assembling and/or producing integrated knowledge of sequence features continues to be an onerous and redundant task despite a large number of existing resources. We have developed SeqDepot-a novel database that focuses solely on two primary goals: (i) assimilating known primary sequences with predicted feature data and (ii) providing the most simple and straightforward means to procure and readily use this information. Access to >28.5 million sequences and 300 million features is provided through a well-documented and flexible RESTful interface that supports fetching specific data subsets, bulk queries, visualization and searching by MD5 digests or external database identifiers. We have also developed an HTML5/JavaScript web application exemplifying how to interact with SeqDepot and Perl/Python scripts for use with local processing pipelines. AVAILABILITY: Freely available on the web at http://seqdepot.net/. RESTaccess via http://seqdepot.net/api/v1. Database files and scripts maybe downloaded from http://seqdepot.net/download.

Analysis of activator and repressor functions reveals the requirements for transcriptional control by LuxR, the master regulator of quorum sensing in Vibrio harveyi.

UNLABELLED: LuxR-type transcription factors are the master regulators of quorum sensing in vibrios. LuxR proteins are unique members of the TetR superfamily of transcription factors because they activate and repress large regulons of genes. Here, we used chromatin immunoprecipitation and nucleotide sequencing (ChIP-seq) to identify LuxR binding sites in the Vibrio harveyi genome. Bioinformatics analyses showed that the LuxR consensus binding site at repressed promoters is a symmetric palindrome, whereas at activated promoters it is asymmetric and contains only half of the palindrome. Using a genetic screen, we isolated LuxR mutants that separated activation and repression functions at representative promoters. These LuxR mutants exhibit sequence-specific DNA binding defects that restrict activation or repression activity to subsets of target promoters. Altering the LuxR DNA binding site sequence to one more closely resembling the ideal LuxR consensus motif can restore in vivo function to a LuxR mutant. This study provides a mechanistic understanding of how a single protein can recognize a variety of binding sites to differentially regulate gene expression. IMPORTANCE: Bacteria use the cell-cell communication process called quorum sensing to regulate collective behaviors. In vibrios, LuxR-type transcription factors control the quorum-sensing gene expression cascade. LuxR-type proteins are structural homologs of TetR-type transcription factors. LuxR proteins were assumed to function analogously to TetR proteins, which typically bind to a single conserved binding site to repress transcription of one or two genes. We find here that unlike TetR proteins, LuxR acts a global regulator, directly binding upstream of and controlling more than 100 genes. Again unlike TetR, LuxR functions as both an activator and a repressor, and these two activities can be separated by mutagenesis. Finally, the consensus binding motifs driving LuxR-activated and -repressed genes are distinct. This work shows that LuxR, although structurally similar to TetR, has evolved unique features enabling it to differentially control a large regulon of genes in response to quorum-sensing cues.

PAS domain containing chemoreceptor couples dynamic changes in metabolism with chemotaxis.

Chemoreceptors provide sensory specificity and sensitivity that enable motile bacteria to seek optimal positions for growth and metabolism in gradients of various physicochemical cues. Despite the abundance of chemoreceptors, little is known regarding the sensory specificity and the exact contribution of individual chemoreceptors to the lifestyle of bacteria. Azospirillum brasilense are motile bacteria that can fix atmospheric nitrogen under microaerophilic conditions. Here, we characterized a chemoreceptor in this organism, named AerC, which functions as a redox sensor that enables the cells to seek microaerophilic conditions that support optimum nitrogen fixation. AerC is a representative of a widespread class of soluble chemoreceptors that monitor changes in the redox status of the electron transport system via the FAD cofactor associated with its PAS domains. In A. brasilense, AerC clusters at the cell poles. Its cellular localization and contribution to the behavioral response correlate with its expression pattern and with changes in the overall cellular FAD content under nitrogen-fixing conditions. AerC-mediated energy taxis in A. brasilense prevails under conditions of nitrogen fixation, illustrating a strategy by which cells optimize chemosensing to signaling cues that directly affect current metabolic activities and thus revealing a mechanism by which chemotaxis is coordinated with dynamic changes in cell physiology.

The MiST2 database: a comprehensive genomics resource on microbial signal transduction.

The MiST2 database (http://mistdb.com) identifies and catalogs the repertoire of signal transduction proteins in microbial genomes. Signal transduction systems regulate the majority of cellular activities including the metabolism, development, host-recognition, biofilm production, virulence, and antibiotic resistance of human pathogens. Thus, knowledge of the proteins and interactions that comprise these communication networks is an essential component to furthering biomedical discovery. These are identified by searching protein sequences for specific domain profiles that implicate a protein in signal transduction. Compared to the previous version of the database, MiST2 contains a host of new features and improvements including the following: draft genomes; extracytoplasmic function (ECF) sigma factor protein identification; enhanced classification of signaling proteins; novel, high-quality domain models for identifying histidine kinases and response regulators; neighboring two-component genes; gene cart; better search capabilities; enhanced taxonomy browser; advanced genome browser; and a modern, biologist-friendly web interface. MiST2 currently contains 966 complete and 157 draft bacterial and archaeal genomes, which collectively contain more than 245 000 signal transduction proteins. The majority (66%) of these are one-component systems, followed by two-component proteins (26%), chemotaxis (6%), and finally ECF factors (2%).

The third pillar of bacterial signal transduction: classification of the extracytoplasmic function (ECF) sigma factor protein family.

The ability of a bacterial cell to monitor and adaptively respond to its environment is crucial for survival. After one- and two-component systems, extracytoplasmic function (ECF) sigma factors - the largest group of alternative sigma factors - represent the third fundamental mechanism of bacterial signal transduction, with about six such regulators on average per bacterial genome. Together with their cognate anti-sigma factors, they represent a highly modular design that primarily facilitates transmembrane signal transduction. A comprehensive analysis of the ECF sigma factor protein family identified more than 40 distinct major groups of ECF sigma factors. The functional relevance of this classification is supported by the sequence similarity and domain architecture of cognate anti-sigma factors, genomic context conservation, and potential target promoter motifs. Moreover, this phylogenetic analysis revealed unique features indicating novel mechanisms of ECF-mediated signal transduction. This classification, together with the web tool ECFfinder and the information stored in the Microbial Signal Transduction (MiST) database, provides a comprehensive resource for the analysis of ECF sigma factor-dependent gene regulation.

Genomic characterization of methanomicrobiales reveals three classes of methanogens.

BACKGROUND: Methanomicrobiales is the least studied order of methanogens. While these organisms appear to be more closely related to the Methanosarcinales in ribosomal-based phylogenetic analyses, they are metabolically more similar to Class I methanogens. METHODOLOGY/PRINCIPAL FINDINGS: In order to improve our understanding of this lineage, we have completely sequenced the genomes of two members of this order, Methanocorpusculum labreanum Z and Methanoculleus marisnigri JR1, and compared them with the genome of a third, Methanospirillum hungatei JF-1. Similar to Class I methanogens, Methanomicrobiales use a partial reductive citric acid cycle for 2-oxoglutarate biosynthesis, and they have the Eha energy-converting hydrogenase. In common with Methanosarcinales, Methanomicrobiales possess the Ech hydrogenase and at least some of them may couple formylmethanofuran formation and heterodisulfide reduction to transmembrane ion gradients. Uniquely, M. labreanum and M. hungatei contain hydrogenases similar to the Pyrococcus furiosus Mbh hydrogenase, and all three Methanomicrobiales have anti-sigma factor and anti-anti-sigma factor regulatory proteins not found in other methanogens. Phylogenetic analysis based on seven core proteins of methanogenesis and cofactor biosynthesis places the Methanomicrobiales equidistant from Class I methanogens and Methanosarcinales. CONCLUSIONS/SIGNIFICANCE: Our results indicate that Methanomicrobiales, rather than being similar to Class I methanogens or Methanomicrobiales, share some features of both and have some unique properties. We find that there are three distinct classes of methanogens: the Class I methanogens, the Methanomicrobiales (Class II), and the Methanosarcinales (Class III).

The complete genome sequence of Staphylothermus marinus reveals differences in sulfur metabolism among heterotrophic Crenarchaeota.

BACKGROUND: Staphylothermus marinus is an anaerobic, sulfur-reducing peptide fermenter of the archaeal phylum Crenarchaeota. It is the third heterotrophic, obligate sulfur reducing crenarchaeote to be sequenced and provides an opportunity for comparative analysis of the three genomes. RESULTS: The 1.57 Mbp genome of the hyperthermophilic crenarchaeote Staphylothermus marinus has been completely sequenced. The main energy generating pathways likely involve 2-oxoacid:ferredoxin oxidoreductases and ADP-forming acetyl-CoA synthases. S. marinus possesses several enzymes not present in other crenarchaeotes including a sodium ion-translocating decarboxylase likely to be involved in amino acid degradation. S. marinus lacks sulfur-reducing enzymes present in the other two sulfur-reducing crenarchaeotes that have been sequenced -- Thermofilum pendens and Hyperthermus butylicus. Instead it has three operons similar to the mbh and mbx operons of Pyrococcus furiosus, which may play a role in sulfur reduction and/or hydrogen production. The two marine organisms, S. marinus and H. butylicus, possess more sodium-dependent transporters than T. pendens and use symporters for potassium uptake while T. pendens uses an ATP-dependent potassium transporter. T. pendens has adapted to a nutrient-rich environment while H. butylicus is adapted to a nutrient-poor environment, and S. marinus lies between these two extremes. CONCLUSION: The three heterotrophic sulfur-reducing crenarchaeotes have adapted to their habitats, terrestrial vs. marine, via their transporter content, and they have also adapted to environments with differing levels of nutrients. Despite the fact that they all use sulfur as an electron acceptor, they are likely to have different pathways for sulfur reduction.

Genome sequence of Thermofilum pendens reveals an exceptional loss of biosynthetic pathways without genome reduction.

We report the complete genome of Thermofilum pendens, a deeply branching, hyperthermophilic member of the order Thermoproteales in the archaeal kingdom Crenarchaeota. T. pendens is a sulfur-dependent, anaerobic heterotroph isolated from a solfatara in Iceland. It is an extracellular commensal, requiring an extract of Thermoproteus tenax for growth, and the genome sequence reveals that biosynthetic pathways for purines, most amino acids, and most cofactors are absent. In fact, T. pendens has fewer biosynthetic enzymes than obligate intracellular parasites, although it does not display other features that are common among obligate parasites and thus does not appear to be in the process of becoming a parasite. It appears that T. pendens has adapted to life in an environment rich in nutrients. T. pendens was known previously to utilize peptides as an energy source, but the genome revealed a substantial ability to grow on carbohydrates. T. pendens is the first crenarchaeote and only the second archaeon found to have a transporter of the phosphotransferase system. In addition to fermentation, T. pendens may obtain energy from sulfur reduction with hydrogen and formate as electron donors. It may also be capable of sulfur-independent growth on formate with formate hydrogen lyase. Additional novel features are the presence of a monomethylamine:corrinoid methyltransferase, the first time that this enzyme has been found outside the Methanosarcinales, and the presence of a presenilin-related protein. The predicted highly expressed proteins do not include proteins encoded by housekeeping genes and instead include ABC transporters for carbohydrates and peptides and clustered regularly interspaced short palindromic repeat-associated proteins.

MiST: a microbial signal transduction database.

Signal transduction pathways control most cellular activities in living cells ranging from regulation of gene expression to fine-tuning enzymatic activity and controlling motile behavior in response to extracellular and intracellular signals. Because of their extreme sequence variability and extensive domain shuffling, signal transduction proteins are difficult to identify, and their current annotation in most leading databases is often incomplete or erroneous. To overcome this problem, we have developed the microbial signal transduction (MiST) database (http://genomics.ornl.gov/mist), a comprehensive library of the signal transduction proteins from completely sequenced bacterial and archaeal genomes. By searching for domain profiles that implicate a particular protein as participating in signal transduction, we have systematically identified 69 270 two- and one-component proteins in 365 bacterial and archaeal genomes. We have designed a user-friendly website to access and browse the predicted signal transduction proteins within various organisms. Further capabilities include gene/protein sequence retrieval, visualized domain architectures, interactive chromosomal views for exploring gene neighborhood, advanced querying options and cross-species comparison. Newly available, complete genomes are loaded into the database each month. MiST is the only comprehensive and up-to-date electronic catalog of the signaling machinery in microbial genomes.

Burkholderia xenovorans LB400 harbors a multi-replicon, 9.73-Mbp genome shaped for versatility.

Burkholderia xenovorans LB400 (LB400), a well studied, effective polychlorinated biphenyl-degrader, has one of the two largest known bacterial genomes and is the first nonpathogenic Burkholderia isolate sequenced. From an evolutionary perspective, we find significant differences in functional specialization between the three replicons of LB400, as well as a more relaxed selective pressure for genes located on the two smaller vs. the largest replicon. High genomic plasticity, diversity, and specialization within the Burkholderia genus are exemplified by the conservation of only 44% of the genes between LB400 and Burkholderia cepacia complex strain 383. Even among four B. xenovorans strains, genome size varies from 7.4 to 9.73 Mbp. The latter is largely explained by our findings that >20% of the LB400 sequence was recently acquired by means of lateral gene transfer. Although a range of genetic factors associated with in vivo survival and intercellular interactions are present, these genetic factors are likely related to niche breadth rather than determinants of pathogenicity. The presence of at least eleven "central aromatic" and twenty "peripheral aromatic" pathways in LB400, among the highest in any sequenced bacterial genome, supports this hypothesis. Finally, in addition to the experimentally observed redundancy in benzoate degradation and formaldehyde oxidation pathways, the fact that 17.6% of proteins have a better LB400 paralog than an ortholog in a different genome highlights the importance of gene duplication and repeated acquirement, which, coupled with their divergence, raises questions regarding the role of paralogs and potential functional redundancies in large-genome microbes.

Four-helix bundle: a ubiquitous sensory module in prokaryotic signal transduction.

MOTIVATION: Transmembrane chemoreceptors in Escherichia coli utilize ligand-binding domains for detecting various external signals. The structure of this domain in the E.coli aspartate receptor, Tar, is known and its signal transduction mechanism is under investigation. Current domain models for this important sensory module are inaccurate and, therefore, cannot reveal the distribution of this domain within the current genomic landscape. RESULTS: We carried out sensitive and exhaustive PSI-BLAST searches initiated with the sequence corresponding to a known structure of the four-helix, ligand-binding domain of the aspartate chemoreceptor. From the resulting sequences, we built a multiple sequence alignment for this domain family, which confirmed that the current TarH model is erroneous and fails to detect most of the domain homologs. In the process, we developed a technique that visualizes the secondary structure prediction of each protein sequence in order to improve the multiple sequence alignment. We found that the four-helix up-and-down bundle represents a large domain family and includes representatives of all major classes of prokaryotic signal transduction, namely histidine kinases, di-guanylate cyclases and chemotaxis receptors.

Life in hot carbon monoxide: the complete genome sequence of Carboxydothermus hydrogenoformans Z-2901.

We report here the sequencing and analysis of the genome of the thermophilic bacterium Carboxydothermus hydrogenoformans Z-2901. This species is a model for studies of hydrogenogens, which are diverse bacteria and archaea that grow anaerobically utilizing carbon monoxide (CO) as their sole carbon source and water as an electron acceptor, producing carbon dioxide and hydrogen as waste products. Organisms that make use of CO do so through carbon monoxide dehydrogenase complexes. Remarkably, analysis of the genome of C. hydrogenoformans reveals the presence of at least five highly differentiated anaerobic carbon monoxide dehydrogenase complexes, which may in part explain how this species is able to grow so much more rapidly on CO than many other species. Analysis of the genome also has provided many general insights into the metabolism of this organism which should make it easier to use it as a source of biologically produced hydrogen gas. One surprising finding is the presence of many genes previously found only in sporulating species in the Firmicutes Phylum. Although this species is also a Firmicutes, it was not known to sporulate previously. Here we show that it does sporulate and because it is missing many of the genes involved in sporulation in other species, this organism may serve as a "minimal" model for sporulation studies. In addition, using phylogenetic profile analysis, we have identified many uncharacterized gene families found in all known sporulating Firmicutes, but not in any non-sporulating bacteria, including a sigma factor not known to be involved in sporulation previously.

One-component systems dominate signal transduction in prokaryotes.

Two-component systems that link environmental signals to cellular responses are viewed as the primary mode of signal transduction in prokaryotes. By analyzing information encoded by 145 prokaryotic genomes, we found that the majority of signal transduction systems consist of a single protein that contains input and output domains but lacks phosphotransfer domains typical of two-component systems. One-component systems are evolutionarily older, more widely distributed among bacteria and archaea, and display a greater diversity of domains than two-component systems.

The NIT domain: a predicted nitrate-responsive module in bacterial sensory receptors.

A nitrate- and nitrite-sensing (NIT) domain as found in the NasR protein, has been detected in various receptor components of signal transduction pathways in different bacterial lineages. Cellular functions controlled by receptors that contain this novel domain include regulation of gene expression (transcription anti-terminators and histidine kinases), cell motility (chemotaxis receptors) and enzyme activity (diguanylate cyclases and phosphodiesterases).