Global gene expression analysis of the Myxococcus xanthus developmental time course.

To accurately identify the genes and pathways involved in the initiation of the Myxococcus xanthus multicellular developmental program, we have previously reported a method of growing vegetative populations as biofilms within a controllable environment. Using a modified approach to remove up to ~90% rRNAs, we report a comprehensive transcriptional analysis of the M. xanthus developmental cycle while comparing it with the vegetative biofilms grown in rich and poor nutrients. This study identified 1522 differentially regulated genes distributed within eight clusters during development. It also provided a comprehensive overview of genes expressed during a nutrient-stress response, specific development time points, and during development initiation and regulation. We identified several differentially expressed genes involved in key central metabolic pathways suggesting their role in regulating myxobacterial development. Overall, this study will prove an important resource for myxobacterial researchers to delineate the regulatory and functional pathways responsible for development from those of the general nutrient stress response.

Phylogenetically driven sequencing of extremely halophilic archaea reveals strategies for static and dynamic osmo-response.

Organisms across the tree of life use a variety of mechanisms to respond to stress-inducing fluctuations in osmotic conditions. Cellular response mechanisms and phenotypes associated with osmoadaptation also play important roles in bacterial virulence, human health, agricultural production and many other biological systems. To improve understanding of osmoadaptive strategies, we have generated 59 high-quality draft genomes for the haloarchaea (a euryarchaeal clade whose members thrive in hypersaline environments and routinely experience drastic changes in environmental salinity) and analyzed these new genomes in combination with those from 21 previously sequenced haloarchaeal isolates. We propose a generalized model for haloarchaeal management of cytoplasmic osmolarity in response to osmotic shifts, where potassium accumulation and sodium expulsion during osmotic upshock are accomplished via secondary transport using the proton gradient as an energy source, and potassium loss during downshock is via a combination of secondary transport and non-specific ion loss through mechanosensitive channels. We also propose new mechanisms for magnesium and chloride accumulation. We describe the expansion and differentiation of haloarchaeal general transcription factor families, including two novel expansions of the TATA-binding protein family, and discuss their potential for enabling rapid adaptation to environmental fluxes. We challenge a recent high-profile proposal regarding the evolutionary origins of the haloarchaea by showing that inclusion of additional genomes significantly reduces support for a proposed large-scale horizontal gene transfer into the ancestral haloarchaeon from the bacterial domain. The combination of broad (17 genera) and deep (≥5 species in four genera) sampling of a phenotypically unified clade has enabled us to uncover both highly conserved and specialized features of osmoadaptation. Finally, we demonstrate the broad utility of such datasets, for metagenomics, improvements to automated gene annotation and investigations of evolutionary processes.

Identification of an archaeal mercury regulon by chromatin immunoprecipitation.

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

Stable closure of the cytoplasmic half-channel is required for efficient proton transport at physiological membrane potentials in the bacteriorhodopsin catalytic cycle.

The bacteriorhodopsin (BR) Asp96Gly/Phe171Cys/Phe219Leu triple mutant has been shown to translocate protons 66% as efficiently as the wild-type protein. Light-dependent ATP synthesis in haloarchaeal cells expressing the triple mutant is 85% that of the wild-type BR expressing cells. Therefore, the functional activity of BR seems to be largely preserved in the triple mutant despite the observations that its ground-state structure resembles that of the wild-type M state (i.e., the so-called cytoplasmically open state) and that the mutant shows no significant structural changes during its photocycle, in sharp contrast to what occurs in the wild-type protein in which a large structural opening and closing occurs on the cytoplasmic side. To resolve the contradiction between the apparent functional robustness of the triple mutant and the presumed importance of the opening and closing that occurs in the wild-type protein, we conducted additional experiments to compare the behavior of wild-type and mutant proteins under different operational loads. Specifically, we characterized the ability of the two proteins to generate light-driven proton currents against a range of membrane potentials. The wild-type protein showed maximal conductance between -150 and -50 mV, whereas the mutant showed maximal conductance at membrane potentials >+50 mV. Molecular dynamics (MD) simulations of the triple mutant were also conducted to characterize structural changes in the protein and in solvent accessibility that might help to functionally contextualize the current-voltage data. These simulations revealed that the cytoplasmic half-channel of the triple mutant is constitutively open and dynamically exchanges water with the bulk. Collectively, the data and simulations help to explain why this mutant BR does not mediate photosynthetic growth of haloarchaeal cells, and they suggest that the structural closing observed in the wild-type protein likely plays a key role in minimizing substrate back flow in the face of electrochemical driving forces present at physiological membrane potentials.

Microscale sulfur cycling in the phototrophic pink berry consortia of the Sippewissett Salt Marsh.

Microbial metabolism is the engine that drives global biogeochemical cycles, yet many key transformations are carried out by microbial consortia over short spatiotemporal scales that elude detection by traditional analytical approaches. We investigate syntrophic sulfur cycling in the 'pink berry' consortia of the Sippewissett Salt Marsh through an integrative study at the microbial scale. The pink berries are macroscopic, photosynthetic microbial aggregates composed primarily of two closely associated species: sulfide-oxidizing purple sulfur bacteria (PB-PSB1) and sulfate-reducing bacteria (PB-SRB1). Using metagenomic sequencing and (34) S-enriched sulfate stable isotope probing coupled with nanoSIMS, we demonstrate interspecies transfer of reduced sulfur metabolites from PB-SRB1 to PB-PSB1. The pink berries catalyse net sulfide oxidation and maintain internal sulfide concentrations of 0-500 µm. Sulfide within the berries, captured on silver wires and analysed using secondary ion mass spectrometer, increased in abundance towards the berry interior, while δ(34) S-sulfide decreased from 6‰ to -31‰ from the exterior to interior of the berry. These values correspond to sulfate-sulfide isotopic fractionations (15-53‰) consistent with either sulfate reduction or a mixture of reductive and oxidative metabolisms. Together this combined metagenomic and high-resolution isotopic analysis demonstrates active sulfur cycling at the microscale within well-structured macroscopic consortia consisting of sulfide-oxidizing anoxygenic phototrophs and sulfate-reducing bacteria.

IRF-1 and miRNA126 modulate VCAM-1 expression in response to a high-fat meal.

RATIONALE: A high-fat diet accompanied by hypertriglyceridemia increases an individual's risk for development of atherosclerosis. An early event in this process is monocyte recruitment through binding to vascular cell adhesion molecule 1 (VCAM-1) upregulated on inflamed arterial endothelium. Diets high in polyunsaturated fatty acids (PUFAs) may provide athero-protection by ameliorating this effect. OBJECTIVE: We investigated the acute regulation of VCAM-1 expression in human aortic endothelial cells (HAEC) in response to triglyceride-rich lipoproteins (TGRL) isolated from subjects after consumption of a high-fat meal. METHODS AND RESULTS: Postprandial TGRL isolated from 38 subjects were categorized as proatherogenic or antiatherogenic according to their capacity to alter the inflammatory response of HAEC. Proatherogenic TGRL increased expression of VCAM-1, intercellular adhesion molecule 1 (ICAM-1), and E-selectin by ≈20% compared with stimulation with tumor necrosis factor-α alone, whereas antiatherogenic TGRL decreased VCAM-1 expression by ≈20% while still upregulating ICAM-1. The relative atherogenicity of TGRL positively correlated with particle density of TG, apolipoprotein (Apo)CIII, ApoE, and cholesterol. Ω3-PUFA mimicked the effect of antiatherogenic TGRL by downregulating VCAM-1 expression. TGRL exerted this differential regulation of VCAM-1 by reciprocally modulating expression and activity of the transcription factor interferon regulatory factor 1 (IRF-1) and expression of microRNA 126 (miR-126). Overexpression or silencing of IRF-1 or miR-126 expression recapitulated the proatherogenic or antiatherogenic regulation of VCAM-1. CONCLUSIONS: In response to a high-fat meal, TGRL bias the inflammatory response of endothelium via transcriptional and posttranscriptional editing of VCAM-1. Subjects with an anti-inflammatory response to a meal produced TGRL that was enriched in nonesterified fatty acids, decreased IRF-1 expression, increased miR-126 activity, and diminished monocyte arrest.

Conserved substitution patterns around nucleosome footprints in eukaryotes and Archaea derive from frequent nucleosome repositioning through evolution.

Nucleosomes, the basic repeat units of eukaryotic chromatin, have been suggested to influence the evolution of eukaryotic genomes, both by altering the propensity of DNA to mutate and by selection acting to maintain or exclude nucleosomes in particular locations. Contrary to the popular idea that nucleosomes are unique to eukaryotes, histone proteins have also been discovered in some archaeal genomes. Archaeal nucleosomes, however, are quite unlike their eukaryotic counterparts in many respects, including their assembly into tetramers (rather than octamers) from histone proteins that lack N- and C-terminal tails. Here, we show that despite these fundamental differences the association between nucleosome footprints and sequence evolution is strikingly conserved between humans and the model archaeon Haloferax volcanii. In light of this finding we examine whether selection or mutation can explain concordant substitution patterns in the two kingdoms. Unexpectedly, we find that neither the mutation nor the selection model are sufficient to explain the observed association between nucleosomes and sequence divergence. Instead, we demonstrate that nucleosome-associated substitution patterns are more consistent with a third model where sequence divergence results in frequent repositioning of nucleosomes during evolution. Indeed, we show that nucleosome repositioning is both necessary and largely sufficient to explain the association between current nucleosome positions and biased substitution patterns. This finding highlights the importance of considering the direction of causality between genetic and epigenetic change.

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

Deciphering the structure of gene regulatory networks across the tree of life remains one of the major challenges in postgenomic biology. We present a novel ChIP-seq workflow for the archaea using the model organism Halobacterium salinarum sp. NRC-1 and demonstrate its application for mapping the genome-wide binding sites of natively expressed transcription factors. This end-to-end pipeline is the first protocol for ChIP-seq in archaea, with methods and tools for each stage from gene tagging to data analysis and biological discovery. Genome-wide binding sites for transcription factors with many binding sites (TfbD) are identified with sensitivity, while retaining specificity in the identification the smaller regulons (bacteriorhodopsin-activator protein). Chromosomal tagging of target proteins with a compact epitope facilitates a standardized and cost-effective workflow that is compatible with high-throughput immunoprecipitation of natively expressed transcription factors. The Pique package, an open-source bioinformatics method, is presented for identification of binding events. Relative to ChIP-Chip and qPCR, this workflow offers a robust catalog of protein-DNA binding events with improved spatial resolution and significantly decreased cost. While this study focuses on the application of ChIP-seq in H. salinarum sp. NRC-1, our workflow can also be adapted for use in other archaea and bacteria with basic genetic tools.

JContextExplorer: a tree-based approach to facilitate cross-species genomic context comparison.

BACKGROUND: Cross-species comparisons of gene neighborhoods (also called genomic contexts) in microbes may provide insight into determining functionally related or co-regulated sets of genes, suggest annotations of previously un-annotated genes, and help to identify horizontal gene transfer events across microbial species. Existing tools to investigate genomic contexts, however, lack features for dynamically comparing and exploring genomic regions from multiple species. As DNA sequencing technologies improve and the number of whole sequenced microbial genomes increases, a user-friendly genome context comparison platform designed for use by a broad range of users promises to satisfy a growing need in the biological community. RESULTS: Here we present JContextExplorer: a tool that organizes genomic contexts into branching diagrams. We implement several alternative context-comparison and tree rendering algorithms, and allow for easy transitioning between different clustering algorithms. To facilitate genomic context analysis, our tool implements GUI features, such as text search filtering, point-and-click interrogation of individual contexts, and genomic visualization via a multi-genome browser. We demonstrate a use case of our tool by attempting to resolve annotation ambiguities between two highly homologous yet functionally distinct genes in a set of 22 alpha and gamma proteobacteria. CONCLUSIONS: JContextExplorer should enable a broad range of users to analyze and explore genomic contexts. The program has been tested on Windows, Mac, and Linux operating systems, and is implemented both as an executable JAR file and java WebStart. Program executables, source code, and documentation is available at http://www.bme.ucdavis.edu/facciotti/resources_data/software/.

Role of squalene in the organization of monolayers derived from lipid extracts of Halobacterium salinarum.

We have studied interfacial compressibility and lateral organization in monolayer configurations of total (squalene containing) and polar (squalene-devoid) lipid extracts of Halobacterium salinarum NRC-1, an extremely halophilic archaeon. Pressure-area isotherms derived from Langmuir experiments reveal that packing characteristics and elastic compressibility are strongly influenced by the presence of squalene in the total lipid extract. In conjunction with control experiments using mixtures of DPhPC and squalene, our results establish that the presence of squalene significantly extends elastic area compressibility of total lipid extracts, suggesting it has a role in facilitating tighter packing of archaeal lipid mixtures. Moreover, we find that squalene also influences spatial organization in archaeal membranes. Epifluorescence and atomic force microscopy characterization of Langmuir monolayers transferred onto solid hydrophilic substrates reveal an unusual domain morphology. Individual domains of microscopic dimensions (as well as their extended networks) exhibiting a peculiar bowl-like topography are evident in atomic force microscopy images. The tall rims outlining individual domains indicate that squalene accumulates at the domain periphery in a manner similar to the accumulation of cholesterol at domain boundaries in their mixtures with phospholipids. Taken together, the results presented here support the notion that squalene plays a role in modulating molecular packing and lateral organization (i.e., domain formation) in the membranes of archaea analogous to that of cholesterol in eukaryotic membranes.

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

BACKGROUND: Discovery of functionally significant short, statistically overrepresented subsequence patterns (motifs) in a set of sequences is a challenging problem in bioinformatics. Oftentimes, not all sequences in the set contain a motif. These non-motif-containing sequences complicate the algorithmic discovery of motifs. Filtering the non-motif-containing sequences from the larger set of sequences while simultaneously determining the identity of the motif is, therefore, desirable and a non-trivial problem in motif discovery research. RESULTS: We describe MotifCatcher, a framework that extends the sensitivity of existing motif-finding tools by employing random sampling to effectively remove non-motif-containing sequences from the motif search. We developed two implementations of our algorithm; each built around a commonly used motif-finding tool, and applied our algorithm to three diverse chromatin immunoprecipitation (ChIP) data sets. In each case, the motif finder with the MotifCatcher extension demonstrated improved sensitivity over the motif finder alone. Our approach organizes candidate functionally significant discovered motifs into a tree, which allowed us to make additional insights. In all cases, we were able to support our findings with experimental work from the literature. CONCLUSIONS: Our framework demonstrates that additional processing at the sequence entry level can significantly improve the performance of existing motif-finding tools. For each biological data set tested, we were able to propose novel biological hypotheses supported by experimental work from the literature. Specifically, in Escherichia coli, we suggested binding site motifs for 6 non-traditional LexA protein binding sites; in Saccharomyces cerevisiae, we hypothesize 2 disparate mechanisms for novel binding sites of the Cse4p protein; and in Halobacterium sp. NRC-1, we discoverd subtle differences in a general transcription factor (GTF) binding site motif across several data sets. We suggest that small differences in our discovered motif could confer specificity for one or more homologous GTF proteins. We offer a free implementation of the MotifCatcher software package at http://www.bme.ucdavis.edu/facciotti/resources_data/software/.

Mauve assembly metrics.

SUMMARY: High-throughput DNA sequencing technologies have spurred the development of numerous novel methods for genome assembly. With few exceptions, these algorithms are heuristic and require one or more parameters to be manually set by the user. One approach to parameter tuning involves assembling data from an organism with an available high-quality reference genome, and measuring assembly accuracy using some metrics. We developed a system to measure assembly quality under several scoring metrics, and to compare assembly quality across a variety of assemblers, sequence data types, and parameter choices. When used in conjunction with training data such as a high-quality reference genome and sequence reads from the same organism, our program can be used to manually identify an optimal sequencing and assembly strategy for de novo sequencing of related organisms. AVAILABILITY: GPL source code and a usage tutorial is at http://ngopt.googlecode.com CONTACT: aarondarling@ucdavis.edu SUPPLEMENTARY INFORMATION: Supplementary data is available at Bioinformatics online.

Regulatory multidimensionality of gas vesicle biogenesis in Halobacterium salinarum NRC-1.

It is becoming clear that the regulation of gas vesicle biogenesis in Halobacterium salinarum NRC-1 is multifaceted and appears to integrate environmental and metabolic cues at both the transcriptional and posttranscriptional levels. The mechanistic details underlying this process, however, remain unclear. In this manuscript, we quantify the contribution of light scattering made by both intracellular and released gas vesicles isolated from Halobacterium salinarum NRC-1, demonstrating that each form can lead to distinct features in growth curves determined by optical density measured at 600 nm (OD(600)). In the course of the study, we also demonstrate the sensitivity of gas vesicle accumulation in Halobacterium salinarum NRC-1 on small differences in growth conditions and reevaluate published works in the context of our results to present a hypothesis regarding the roles of the general transcription factor tbpD and the TCA cycle enzyme aconitase on the regulation of gas vesicle biogenesis.

Prevalence of transcription promoters within archaeal operons and coding sequences.

Despite the knowledge of complex prokaryotic-transcription mechanisms, generalized rules, such as the simplified organization of genes into operons with well-defined promoters and terminators, have had a significant role in systems analysis of regulatory logic in both bacteria and archaea. Here, we have investigated the prevalence of alternate regulatory mechanisms through genome-wide characterization of transcript structures of approximately 64% of all genes, including putative non-coding RNAs in Halobacterium salinarum NRC-1. Our integrative analysis of transcriptome dynamics and protein-DNA interaction data sets showed widespread environment-dependent modulation of operon architectures, transcription initiation and termination inside coding sequences, and extensive overlap in 3' ends of transcripts for many convergently transcribed genes. A significant fraction of these alternate transcriptional events correlate to binding locations of 11 transcription factors and regulators (TFs) inside operons and annotated genes-events usually considered spurious or non-functional. Using experimental validation, we illustrate the prevalence of overlapping genomic signals in archaeal transcription, casting doubt on the general perception of rigid boundaries between coding sequences and regulatory elements.

The Exploration of Novel Regulatory Relationships Drives Haloarchaeal Operon-Like Structural Dynamics over Short Evolutionary Distances.

Operons are a dominant feature of bacterial and archaeal genome organization. Numerous investigations have related aspects of operon structure to operon function, making operons exemplars for studies aimed at deciphering Nature's design principles for genomic organization at a local scale. We consider this understanding to be both fundamentally important and ultimately useful in the de novo design of increasingly complex synthetic circuits. Here we analyze the evolution of the genomic context of operon-like structures in a set of 76 sequenced and annotated species of halophilic archaea. The phylogenetic depth and breadth of this dataset allows insight into changes in operon-like structures over shorter evolutionary time scales than have been studied in previous cross-species analysis of operon evolution. Our analysis, implemented in the updated software package JContextExplorer finds that operon-like context as measured by changes in structure frequently differs from a sequence divergence model of whole-species phylogeny and that changes seem to be dominated by the exploration of novel regulatory relationships.

Combining Microbial Culturing With Mathematical Modeling in an Introductory Course-Based Undergraduate Research Experience.

Quantitative techniques are a critical part of contemporary biology research, but students interested in biology enter college with widely varying quantitative skills and attitudes toward mathematics. Course-based undergraduate research experiences (CUREs) may be an early way to build student competency and positive attitudes. Here we describe the design, implementation, and assessment of an introductory quantitative CURE focused on halophilic microbes. In this CURE, students culture and isolate halophilic microbes from environmental and food samples, perform growth assays, then use mathematical modeling to quantify the growth rate of strains in different salinities. To assess how the course may impact students' future academic plans and attitudes toward the use of math in biology, we used pre- and post-quarter surveys. Students who completed the course showed more positive attitudes toward science learning and an increased interest in pursuing additional quantitative biology experiences. We argue that the classroom application of microbiology methods, combined with mathematical modeling using student-generated data, provides a degree of student ownership, collaboration, iteration, and discovery that makes quantitative learning both relevant and exciting to students.

Draft Genome Sequences of 16 Halophilic Prokaryotes Isolated from Diverse Environments.

Halophile-specific enzymes have wide-ranging industrial and commercial applications. Despite their importance, there is a paucity of available halophile whole-genome sequences. Here, we report the draft genome sequences of 16 diverse salt-tolerant strains of bacteria and archaea isolated from a variety of high-salt environments.

Model-based deconvolution of genome-wide DNA binding.

MOTIVATION: Chromatin immunoprecipitation followed by hybridization to a genomic tiling microarray (ChIP-chip) is a routinely used protocol for localizing the genomic targets of DNA-binding proteins. The resolution to which binding sites in this assay can be identified is commonly considered to be limited by two factors: (1) the resolution at which the genomic targets are tiled in the microarray and (2) the large and variable lengths of the immunoprecipitated DNA fragments. RESULTS: We have developed a generative model of binding sites in ChIP-chip data and an approach, MeDiChI, for efficiently and robustly learning that model from diverse data sets. We have evaluated MeDiChI's performance using simulated data, as well as on several diverse ChIP-chip data sets collected on widely different tiling array platforms for two different organisms (Saccharomyces cerevisiae and Halobacterium salinarium NRC-1). We find that MeDiChI accurately predicts binding locations to a resolution greater than that of the probe spacing, even for overlapping peaks, and can increase the effective resolution of tiling array data by a factor of 5x or better. Moreover, the method's performance on simulated data provides insights into effectively optimizing the experimental design for increased binding site localization accuracy and efficacy. AVAILABILITY: MeDiChI is available as an open-source R package, including all data, from http://baliga.systemsbiology.net/medichi.

Training interdisciplinary scientists for systems biology.

The rise of interdisciplinary research programs in recent years has spawned numerous questions regarding the best way to organize interdisciplinary programs and how to best train new interdisciplinary scientists in ways that will catalyze novel discoveries in biology. Systems biology, a relatively new branch of science, can be considered in many respects a poster child for modern interdisciplinary science. It not only requires that people from different traditional disciplines work closely together but it also requires the development of unique training environments to educate the next generation of systems biologists. The unique scientific and training challenges associated with the development of systems biology are certainly faced across the spectrum of interdisciplinary endeavors. Therefore, it is useful for scientists interested in building interdisciplinary research programs to consider the merits of successful systems biology initiatives. Institute for Systems Biology is one such example, and in the following, several key aspects of Institute for Systems Biology that make it both a unique and successful interdisciplinary science and training center are discussed.