Life history implications of rRNA gene copy number in Escherichia coli

Variability of rRNA Operon Copy Number and Growth Rate Dynamics of Bacillus Isolated from an Extremely Oligotrophic Aquatic Ecosystem

The ribosomal RNA (rrn) operon is a key suite of genes related to the production of protein synthesis machinery and thus to bacterial growth physiology. Experimental evidence has suggested an intrinsic relationship between the number of copies of this operon and environmental resource availability, especially the availability of phosphorus (P), because bacteria that live in oligotrophic ecosystems usually have few rrn operons and a slow growth rate. The Cuatro Ciénegas Basin (CCB) is a complex aquatic ecosystem that contains an unusually high microbial diversity that is able to persist under highly oligotrophic conditions. These environmental conditions impose a variety of strong selective pressures that shape the genome dynamics of their inhabitants. The genus Bacillus is one of the most abundant cultivable bacterial groups in the CCB and usually possesses a relatively large number of rrn operon copies (6–15 copies). The main goal of this study was to analyze the variation in the number of rrn operon copies of Bacillus in the CCB and to assess their growth-related properties as well as their stoichiometric balance (N and P content). We defined 18 phylogenetic groups within the Bacilli clade and documented a range of from six to 14 copies of the rrn operon. The growth dynamic of these Bacilli was heterogeneous and did not show a direct relation to the number of operon copies. Physiologically, our results were not consistent with the Growth Rate Hypothesis, since the copies of the rrn operon were decoupled from growth rate. However, we speculate that the diversity of the growth properties of these Bacilli as well as the low P content of their cells in an ample range of rrn copy number is an adaptive response to oligotrophy of the CCB and could represent an ecological mechanism that allows these taxa to coexist. These findings increase the knowledge of the variability in the number of copies of the rrn operon in the genus Bacillus and give insights about the physiology of this bacterial group under extreme oligotrophic conditions.

Differential Response of Acidobacteria Subgroups to Forest-to-Pasture Conversion and Their Biogeographic Patterns in the Western Brazilian Amazon

Members of the phylum Acidobacteria are among the most abundant soil bacteria on Earth, but little is known about their response to environmental changes. We asked how the relative abundance and biogeographic patterning of this phylum and its subgroups responded to forest-to-pasture conversion in soils of the western Brazilian Amazon. Pyrosequencing of 16S rRNA genes was employed to assess the abundance and composition of the Acidobacteria community across 54 soil samples taken using a spatially nested sampling scheme at the landscape level. Numerically, Acidobacteria represented 20% of the total bacterial community in forest soils and 11% in pasture soils. Overall, 15 different Acidobacteria subgroups of the current 26 subgroups were detected, with Acidobacteria subgroups 1, 3, 5, and 6 accounting together for 87% of the total Acidobacteria community in forest soils and 75% in pasture soils. Concomitant with changes in soil chemistry after forest-to-pasture conversion—particularly an increase in properties linked to soil acidity and nutrient availability—we observed an increase in the relative abundances of Acidobacteria subgroups 4, 10, 17, and 18, and a decrease in the relative abundances of other Acidobacteria subgroups in pasture relative to forest soils. The composition of the total Acidobacteria community as well as the most abundant Acidobacteria subgroups (1, 3, 5, and 6) was significantly more similar in composition across space in pasture soils than in forest soils. These results suggest that preponderant responses of Acidobacteria subgroups, especially subgroups 1, 3, 4, 5, and 6, to forest-to-pasture conversion effects in soils could be used to define management-indicators of agricultural practices in the Amazon Basin. These acidobacterial responses are at least in part through alterations on acidity- and nutrient-related properties of the Amazon soils.

Decreases in average bacterial community rRNA operon copy number during succession

Trait-based studies can help clarify the mechanisms driving patterns of microbial community assembly and coexistence. Here, we use a trait-based approach to explore the importance of rRNA operon copy number in microbial succession, building on prior evidence that organisms with higher copy numbers respond more rapidly to nutrient inputs. We set flasks of heterotrophic media into the environment and examined bacterial community assembly at seven time points. Communities were arrayed along a geographic gradient to introduce stochasticity via dispersal processes and were analyzed using 16 S rRNA gene pyrosequencing, and rRNA operon copy number was modeled using ancestral trait reconstruction. We found that taxonomic composition was similar between communities at the beginning of the experiment and then diverged through time; as well, phylogenetic clustering within communities decreased over time. The average rRNA operon copy number decreased over the experiment, and variance in rRNA operon copy number was lowest both early and late in succession. We then analyzed bacterial community data from other soil and sediment primary and secondary successional sequences from three markedly different ecosystem types. Our results demonstrate that decreases in average copy number are a consistent feature of communities across various drivers of ecological succession. Importantly, our work supports the scaling of the copy number trait over multiple levels of biological organization, ranging from cells to populations and communities, with implications for both microbial ecology and evolution.

Single-cell kinetics of a repressilator when implemented in a single-copy plasmid

Synthetic genetic clocks, such as the Elowitz-Leibler repressilator, will be key regulatory components of future synthetic circuits. We constructed a single-copy repressilator (SCR) by implementing the original repressilator circuit on a single-copy F-plasmid. After verifying its functionality, we studied its behaviour as a function of temperature and compared it with that of the original low-copy-number repressilator (LCR). Namely, we compared the period of oscillations, functionality (the fraction of cells exhibiting oscillations) and robustness to internal fluctuations (the fraction of expected oscillations that would occur). We found that, under optimal temperature conditions, the dynamics of the two systems differs significantly, although qualitatively they respond similarly to temperature changes. Exception to this is in the functionality, in which the SCR is higher at lower temperatures but lower at higher temperatures. Next, by adding IPTG to the medium at low and high concentrations during microscopy sessions, we showed that the functionality of the SCR is more robust to external perturbations, which indicates that the oscillatory behaviour of the LCR can be disrupted by affecting only a few of the copies in a cell. We conclude that the SCR, the first functional, synthetic, single-copy, ring-type genetic clock, is more robust to lower temperatures and to external perturbations than the original LCR. The SCR will be of use in future synthetic circuits, since it complements the array of tasks that the LCR can perform.

Genomic location of the major ribosomal protein gene locus determines Vibrio cholerae global growth and infectivity

The effects on cell physiology of gene order within the bacterial chromosome are poorly understood. In silico approaches have shown that genes involved in transcription and translation processes, in particular ribosomal protein (RP) genes, localize near the replication origin (oriC) in fast-growing bacteria suggesting that such a positional bias is an evolutionarily conserved growth-optimization strategy. Such genomic localization could either provide a higher dosage of these genes during fast growth or facilitate the assembly of ribosomes and transcription foci by keeping physically close the many components of these macromolecular machines. To explore this, we used novel recombineering tools to create a set of Vibrio cholerae strains in which S10-spec-α (S10), a locus bearing half of the ribosomal protein genes, was systematically relocated to alternative genomic positions. We show that the relative distance of S10 to the origin of replication tightly correlated with a reduction of S10 dosage, mRNA abundance and growth rate within these otherwise isogenic strains. Furthermore, this was accompanied by a significant reduction in the host-invasion capacity in Drosophila melanogaster. Both phenotypes were rescued in strains bearing two S10 copies highly distal to oriC, demonstrating that replication-dependent gene dosage reduction is the main mechanism behind these alterations. Hence, S10 positioning connects genome structure to cell physiology in Vibrio cholerae. Our results show experimentally for the first time that genomic positioning of genes involved in the flux of genetic information conditions global growth control and hence bacterial physiology and potentially its evolution.

Increasing evidence indicates that nucleoid spatiotemporal organization is crucial for bacterial physiology since these microorganism lack compartmentalized nucleus. However, it is still unclear how gene order within the chromosome can influence cell physiology. Here, by systematically relocating ribosomal protein genes to different genomic positions in Vibrio cholerae, we revealed drastic differences in growth rate and infectivity of this isogenic strain set. We show that genomic positioning of ribosomal protein genes is crucial for physiology by providing replication-dependent higher dosage. Therefore it might play a key role in genome evolution of bacterial species. This work will contribute to discover genomic rules governing cell physiology which will be essential in the context of the creation of new artificial life forms.

Long-term forest soil warming alters microbial communities in temperate forest soils

Soil microbes are major drivers of soil carbon cycling, yet we lack an understanding of how climate warming will affect microbial communities. Three ongoing field studies at the Harvard Forest Long-term Ecological Research (LTER) site (Petersham, MA) have warmed soils 5°C above ambient temperatures for 5, 8, and 20 years. We used this chronosequence to test the hypothesis that soil microbial communities have changed in response to chronic warming. Bacterial community composition was studied using Illumina sequencing of the 16S ribosomal RNA gene, and bacterial and fungal abundance were assessed using quantitative PCR. Only the 20-year warmed site exhibited significant change in bacterial community structure in the organic soil horizon, with no significant changes in the mineral soil. The dominant taxa, abundant at 0.1% or greater, represented 0.3% of the richness but nearly 50% of the observations (sequences). Individual members of the Actinobacteria, Alphaproteobacteria and Acidobacteria showed strong warming responses, with one Actinomycete decreasing from 4.5 to 1% relative abundance with warming. Ribosomal RNA copy number can obfuscate community profiles, but is also correlated with maximum growth rate or trophic strategy among bacteria. Ribosomal RNA copy number correction did not affect community profiles, but rRNA copy number was significantly decreased in warming plots compared to controls. Increased bacterial evenness, shifting beta diversity, decreased fungal abundance and increased abundance of bacteria with low rRNA operon copy number, including Alphaproteobacteria and Acidobacteria, together suggest that more or alternative niche space is being created over the course of long-term warming.

Engineered ribosomal RNA operon copy-number variants of E. coli reveal the evolutionary trade-offs shaping rRNA operon number

Ribosomal RNA (rrn) operons, characteristically present in several copies in bacterial genomes (7 in E. coli), play a central role in cellular physiology. We investigated the factors determining the optimal number of rrn operons in E. coli by constructing isogenic variants with 5–10 operons. We found that the total RNA and protein content, as well as the size of the cells reflected the number of rrn operons. While growth parameters showed only minor differences, competition experiments revealed a clear pattern: 7–8 copies were optimal under conditions of fluctuating, occasionally rich nutrient influx and lower numbers were favored in stable, nutrient-limited environments. We found that the advantages of quick adjustment to nutrient availability, rapid growth and economic regulation of ribosome number all contribute to the selection of the optimal rrn operon number. Our results suggest that the wt rrn operon number of E. coli reflects the natural, ‘feast and famine’ life-style of the bacterium, however, different copy numbers might be beneficial under different environmental conditions. Understanding the impact of the copy number of rrn operons on the fitness of the cell is an important step towards the creation of functional and robust genomes, the ultimate goal of synthetic biology.

How many initiator tRNA genes does Escherichia coli need?

Multiple copies of a gene require enhanced investment on the part of the cell and, as such, call for an explanation. The observation that Escherichia coli has four copies of initiator tRNA (tRNAi) genes, encoding a special tRNA (tRNA(fMet)) required to start protein synthesis, is puzzling particularly because the cell appears to be unaffected by the removal of one copy. However, the fitness of an organism has both absolute and relative connotations. Thus, we carried out growth competition experiments between E. coli strains that differ in the number of tRNAi genes they contain. This has enabled us to uncover an unexpected link between the number of tRNAi genes and protein synthesis, nutritional status, and fitness. Wild-type strains with the canonical four tRNAi genes are favored in nutrient-rich environments, and those carrying fewer are favored in nutrient-poor environments. Auxotrophs behave as if they have a nutritionally poor internal environment. A heuristic model that links tRNAi gene copy number, genetic stress, and growth rate accounts for the findings. Our observations provide strong evidence that natural selection can work through seemingly minor quantitative variations in gene copy number and thereby impact organismal fitness.

Intragenic homogenization and multiple copies of prey-wrapping silk genes in Argiope garden spiders

Background

Spider silks are spectacular examples of phenotypic diversity arising from adaptive molecular evolution. An individual spider can produce an array of specialized silks, with the majority of constituent silk proteins encoded by members of the spidroin gene family. Spidroins are dominated by tandem repeats flanked by short, non-repetitive N- and C-terminal coding regions. The remarkable mechanical properties of spider silks have been largely attributed to the repeat sequences. However, the molecular evolutionary processes acting on spidroin terminal and repetitive regions remain unclear due to a paucity of complete gene sequences and sampling of genetic variation among individuals. To better understand spider silk evolution, we characterize a complete aciniform spidroin gene from an Argiope orb-weaving spider and survey aciniform gene fragments from congeneric individuals.

Results

We present the complete aciniform spidroin (AcSp1) gene from the silver garden spider Argiope argentata (Aar_AcSp1), and document multiple AcSp1 loci in individual genomes of A. argentata and the congeneric A. trifasciata and A. aurantia. We find that Aar_AcSp1 repeats have >98% pairwise nucleotide identity. By comparing AcSp1 repeat amino acid sequences between Argiope species and with other genera, we identify regions of conservation over vast amounts of evolutionary time. Through a PCR survey of individual A. argentata, A. trifasciata, and A. aurantia genomes, we ascertain that AcSp1 repeats show limited variation between species whereas terminal regions are more divergent. We also find that average dN/dS across codons in the N-terminal, repetitive, and C-terminal encoding regions indicate purifying selection that is strongest in the N-terminal region.

Conclusions

Using the complete A. argentata AcSp1 gene and spidroin genetic variation between individuals, this study clarifies some of the molecular evolutionary processes underlying the spectacular mechanical attributes of aciniform silk. It is likely that intragenic concerted evolution and functional constraints on A. argentata AcSp1 repeats result in extreme repeat homogeneity. The maintenance of multiple AcSp1 encoding loci in Argiope genomes supports the hypothesis that Argiope spiders require rapid and efficient protein production to support their prolific use of aciniform silk for prey-wrapping and web-decorating. In addition, multiple gene copies may represent the early stages of spidroin diversification.

Genome degeneration and adaptation in a nascent stage of symbiosis

Symbiotic associations between animals and microbes are ubiquitous in nature, with an estimated 15% of all insect species harboring intracellular bacterial symbionts. Most bacterial symbionts share many genomic features including small genomes, nucleotide composition bias, high coding density, and a paucity of mobile DNA, consistent with long-term host association. In this study, we focus on the early stages of genome degeneration in a recently derived insect-bacterial mutualistic intracellular association. We present the complete genome sequence and annotation of Sitophilus oryzae primary endosymbiont (SOPE). We also present the finished genome sequence and annotation of strain HS, a close free-living relative of SOPE and other insect symbionts of the Sodalis-allied clade, whose gene inventory is expected to closely resemble the putative ancestor of this group. Structural, functional, and evolutionary analyses indicate that SOPE has undergone extensive adaptation toward an insect-associated lifestyle in a very short time period. The genome of SOPE is large in size when compared with many ancient bacterial symbionts; however, almost half of the protein-coding genes in SOPE are pseudogenes. There is also evidence for relaxed selection on the remaining intact protein-coding genes. Comparative analyses of the whole-genome sequence of strain HS and SOPE highlight numerous genomic rearrangements, duplications, and deletions facilitated by a recent expansion of insertions sequence elements, some of which appear to have catalyzed adaptive changes. Functional metabolic predictions suggest that SOPE has lost the ability to synthesize several essential amino acids and vitamins. Analyses of the bacterial cell envelope and genes encoding secretion systems suggest that these structures and elements have become simplified in the transition to a mutualistic association.

Intragenomic heterogeneity of 16S rRNA genes causes overestimation of prokaryotic diversity

Ever since Carl Woese introduced the use of 16S rRNA genes for determining the phylogenetic relationships of prokaryotes, this method has been regarded as the "gold standard" in both microbial phylogeny and ecology studies. However, intragenomic heterogeneity within 16S rRNA genes has been reported in many investigations and is believed to bias the estimation of prokaryotic diversity. In the current study, 2,013 completely sequenced genomes of bacteria and archaea were analyzed and intragenomic heterogeneity was found in 952 genomes (585 species), with 87.5% of the divergence detected being below the 1% level. In particular, some extremophiles (thermophiles and halophiles) were found to harbor highly divergent 16S rRNA genes. Overestimation caused by 16S rRNA gene intragenomic heterogeneity was evaluated at different levels using the full-length and partial 16S rRNA genes usually chosen as targets for pyrosequencing. The result indicates that, at the unique level, full-length 16S rRNA genes can produce an overestimation of as much as 123.7%, while at the 3% level, an overestimation of 12.9% for the V6 region may be introduced. Further analysis showed that intragenomic heterogeneity tends to concentrate in specific positions, with the V1 and V6 regions suffering the most intragenomic heterogeneity and the V4 and V5 regions suffering the least intragenomic heterogeneity in bacteria. This is the most up-to-date overview of the diversity of 16S rRNA genes within prokaryotic genomes. It not only provides general guidance on how much overestimation can be introduced when applying 16S rRNA gene-based methods, due to its intragenomic heterogeneity, but also recommends that, for bacteria, this overestimation be minimized using primers targeting the V4 and V5 regions.

High-throughput sequencing: a roadmap toward community ecology

High-throughput sequencing is becoming increasingly important in microbial ecology, yet it is surprisingly under-used to generate or test biogeographic hypotheses. In this contribution, we highlight how adding these methods to the ecologist toolbox will allow the detection of new patterns, and will help our understanding of the structure and dynamics of diversity. Starting with a review of ecological questions that can be addressed, we move on to the technical and analytical issues that will benefit from an increased collaboration between different disciplines.

Intra- and intergenomic variation of ribosomal RNA operons in concurrent Alteromonas macleodii strains

Biodiversity estimates based on ribosomal operon sequence diversity rely on the premise that a sequence is characteristic of a single specific taxon or operational taxonomic unit (OTU). Here, we have studied the sequence diversity of 14 ribosomal RNA operons (rrn) contained in the genomes of two isolates (five operons in each genome) and four metagenomic fosmids, all from the same seawater sample. Complete sequencing of the isolate genomes and the fosmids establish that they represent strains of the same species, Alteromonas macleodii, with average nucleotide identity (ANI) values >97 %. Nonetheless, we observed high levels of intragenomic heterogeneity (i.e., variability between operons of a single genome) affecting multiple regions of the 16S and 23S rRNA genes as well as the internally transcribed spacer 1 (ITS-1) region. Furthermore, the ribosomal operons exhibited intergenomic heterogeneity (i.e., variability between operons located in separate genomes) in each of these regions, compounding the variability. Our data reveal the extensive heterogeneity observed in natural populations of A. macleodii at a single point in time and support the idea that distinct lineages of A. macleodii exist in the deep Mediterranean. These findings highlight the potential of rRNA fingerprinting methods to misrepresent species diversity while simultaneously failing to recognize the ecological significance of individual strains.

Detecting laboratory DNA contamination using polyester-rayon wipes: a method validation study

Due to the high sensitivity of many PCR assays, extraneous target DNA in a laboratory setting can lead to false positive results. To assess the presence of extraneous DNA, many laboratories use gauze wipes to sample laboratory surfaces. The accuracy, precision, limits of detection, linearity, and robustness of a wipe test method and each associated wipe processing step were evaluated using E. coli genomic DNA. The method demonstrated a limit of detection of 108 copies of DNA, which equates to detectable surface concentration of 4.5×10(5) copies of DNA per area sampled. Recovery efficiency or accuracy is 22±10% resulting from a >58% loss of DNA occurring at the wipe wash step. The method is robust, performing consistently despite deliberate modifications of the protocol.

Fundamentals of microbial community resistance and resilience

Microbial communities are at the heart of all ecosystems, and yet microbial community behavior in disturbed environments remains difficult to measure and predict. Understanding the drivers of microbial community stability, including resistance (insensitivity to disturbance) and resilience (the rate of recovery after disturbance) is important for predicting community response to disturbance. Here, we provide an overview of the concepts of stability that are relevant for microbial communities. First, we highlight insights from ecology that are useful for defining and measuring stability. To determine whether general disturbance responses exist for microbial communities, we next examine representative studies from the literature that investigated community responses to press (long-term) and pulse (short-term) disturbances in a variety of habitats. Then we discuss the biological features of individual microorganisms, of microbial populations, and of microbial communities that may govern overall community stability. We conclude with thoughts about the unique insights that systems perspectives – informed by meta-omics data – may provide about microbial community stability.

Dynamic expression of the translational machinery during Bacillus subtilis life cycle at a single cell level

The ability of bacteria to responsively regulate the expression of translation components is crucial for rapid adaptation to fluctuating environments. Utilizing Bacillus subtilis (B. subtilis) as a model organism, we followed the dynamics of the translational machinery at a single cell resolution during growth and differentiation. By comprehensive monitoring the activity of the major rrn promoters and ribosomal protein production, we revealed diverse dynamics between cells grown in rich and poor medium, with the most prominent dissimilarities exhibited during deep stationary phase. Further, the variability pattern of translational activity varied among the cells, being affected by nutrient availability. We have monitored for the first time translational dynamics during the developmental process of sporulation within the two distinct cellular compartments of forespore and mother-cell. Our study uncovers a transient forespore specific increase in expression of translational components. Finally, the contribution of each rrn promoter throughout the bacterium life cycle was found to be relatively constant, implying that differential expression is not the main purpose for the existence of multiple rrn genes. Instead, we propose that coordination of the rrn operons serves as a strategy to rapidly fine tune translational activities in a synchronized fashion to achieve an optimal translation level for a given condition.

Variable copy number, intra-genomic heterogeneities and lateral transfers of the 16S rRNA gene in Pseudomonas

Even though the 16S rRNA gene is the most commonly used taxonomic marker in microbial ecology, its poor resolution is still not fully understood at the intra-genus level. In this work, the number of rRNA gene operons, intra-genomic heterogeneities and lateral transfers were investigated at a fine-scale resolution, throughout the Pseudomonas genus. In addition to nineteen sequenced Pseudomonas strains, we determined the 16S rRNA copy number in four other Pseudomonas strains by Southern hybridization and Pulsed-Field Gel Electrophoresis, and studied the intra-genomic heterogeneities by Denaturing Gradient Gel Electrophoresis and sequencing. Although the variable copy number (from four to seven) seems to be correlated with the evolutionary distance, some close strains in the P. fluorescens lineage showed a different number of 16S rRNA genes, whereas all the strains in the P. aeruginosa lineage displayed the same number of genes (four copies). Further study of the intra-genomic heterogeneities revealed that most of the Pseudomonas strains (15 out of 19 strains) had at least two different 16S rRNA alleles. A great difference (5 or 19 nucleotides, essentially grouped near the V1 hypervariable region) was observed only in two sequenced strains. In one of our strains studied (MFY30 strain), we found a difference of 12 nucleotides (grouped in the V3 hypervariable region) between copies of the 16S rRNA gene. Finally, occurrence of partial lateral transfers of the 16S rRNA gene was further investigated in 1803 full-length sequences of Pseudomonas available in the databases. Remarkably, we found that the two most variable regions (the V1 and V3 hypervariable regions) had probably been laterally transferred from another evolutionary distant Pseudomonas strain for at least 48.3 and 41.6% of the 16S rRNA sequences, respectively. In conclusion, we strongly recommend removing these regions of the 16S rRNA gene during the intra-genus diversity studies.

BactQuant: an enhanced broad-coverage bacterial quantitative real-time PCR assay

BACKGROUND

Bacterial load quantification is a critical component of bacterial community analysis, but a culture-independent method capable of detecting and quantifying diverse bacteria is needed. Based on our analysis of a diverse collection of 16 S rRNA gene sequences, we designed a broad-coverage quantitative real-time PCR (qPCR) assay--BactQuant--for quantifying 16 S rRNA gene copy number and estimating bacterial load. We further utilized in silico evaluation to complement laboratory-based qPCR characterization to validate BactQuant.

METHODS

The aligned core set of 4,938 16 S rRNA gene sequences in the Greengenes database were analyzed for assay design. Cloned plasmid standards were generated and quantified using a qPCR-based approach. Coverage analysis was performed computationally using >670,000 sequences and further evaluated following the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) guidelines.

RESULTS

A bacterial TaqMan® qPCR assay targeting a 466 bp region in V3-V4 was designed. Coverage analysis showed that 91% of the phyla, 96% of the genera, and >80% of the 89,537 species analyzed contained at least one perfect sequence match to the BactQuant assay. Of the 106 bacterial species evaluated, amplification efficiencies ranged from 81 to 120%, with r2-value of >0.99, including species with sequence mismatches. Inter- and intra-run coefficient of variance was <3% and <16% for Ct and copy number, respectively.

CONCLUSIONS

The BactQuant assay offers significantly broader coverage than a previously reported universal bacterial quantification assay BactQuant in vitro performance was better than the in silico predictions.

Large variations in bacterial ribosomal RNA genes

Ribosomal RNA (rRNA) genes, essential to all forms of life, have been viewed as highly conserved and evolutionarily stable, partly because very little is known about their natural variations. Here, we explored large-scale variations of rRNA genes through bioinformatic analyses of available complete bacterial genomic sequences with an emphasis on formation mechanisms and biological significance. Interestingly, we found bacterial genomes in which no 16S rRNA genes harbor the conserved core of the anti–Shine-Dalgarno sequence (5′-CCTCC-3′). This loss was accompanied by elimination of Shine-Dalgarno–like sequences upstream of their protein-coding genes. Those genomes belong to 1 or 2 of the following categories: primary symbionts, hemotropic Mycoplasma, and Flavobacteria. We also found many rearranged rRNA genes and reconstructed their history. Conjecturing the underlying mechanisms, such as inversion, partial duplication, transposon insertion, deletion, and substitution, we were able to infer their biological significance, such as co-orientation of rRNA transcription and chromosomal replication, lateral transfer of rRNA gene segments, and spread of rRNA genes with an apparent structural defect through gene conversion. These results open the way to understanding dynamic evolutionary changes of rRNA genes and the translational machinery.

The variable subdomain of Escherichia coli SecA functions to regulate SecA ATPase activity and ADP release

Bacterial SecA proteins can be categorized by the presence or absence of a variable subdomain (VAR) located within nucleotide-binding domain II of the SecA DEAD motor. Here we show that VAR is dispensable for SecA function, since the VAR deletion mutant secAΔ519-547 displayed a wild-type rate of cellular growth and protein export. Loss or gain of VAR is extremely rare in the history of bacterial evolution, indicating that it appears to contribute to secA function within the relevant species in their natural environments. VAR removal also results in additional secA phenotypes: azide resistance (Azi(r)) and suppression of signal sequence defects (PrlD). The SecAΔ(519-547) protein was found to be modestly hyperactive for SecA ATPase activities and displayed an accelerated rate of ADP release, consistent with the biochemical basis of azide resistance. Based on our findings, we discuss models whereby VAR allosterically regulates SecA DEAD motor function at SecYEG.

Phylogenetic clustering of soil microbial communities by 16S rRNA but not 16S rRNA genes

We evaluated phylogenetic clustering of bacterial and archaeal communities from redox-dynamic subtropical forest soils that were defined by 16S rRNA and rRNA gene sequences. We observed significant clustering for the RNA-based communities but not the DNA-based communities, as well as increasing clustering over time of the highly active taxa detected by only rRNA.

Bacterial cells carrying synthetic dual-function operon survived starvation

A synthetic dual-function operon with a bistable structure was designed and successfully integrated into the bacterial genome. Bistability was generated by the mutual inhibitory structure comprised of the promoters P(tet) and P(lac) and the repressors LacI and TetR. Dual function essential for cell growth was introduced by replacing the genes (i.e., hisC and leuB) encoding proteins involved in the biosynthesis of histidine and leucine from their native chromosomal locations to the synthetic operon. Both colony formation and population dynamics of the cells carrying this operon showed that the cells survived starvation and the newly formed population transited between the two stable states, representing the induced hisC and leuB levels, in accordance with the nutritional status. The results strongly suggested that the synthetic design of proto-operons sensitive to external perturbations is practical and functional in native cells.

Combined phylogenetic and genomic approaches for the high-throughput study of microbial habitat adaptation

High-throughput sequencing technologies provide new opportunities to address longstanding questions about habitat adaptation in microbial organisms. How have microbes managed to adapt to such a wide range of environments, and what genomic features allow for such adaptation? We review recent large-scale studies of habitat adaptation, with emphasis on those that utilize phylogenetic techniques. On the basis of current trends, we summarize methodological challenges faced by investigators, and the tools, techniques and analytical approaches available to overcome them. Phylogenetic approaches and detailed information about each environmental sample will be crucial as the ability to collect genome sequences continues to expand.

Quantitative metagenomic analyses based on average genome size normalization

Over the past quarter-century, microbiologists have used DNA sequence information to aid in the characterization of microbial communities. During the last decade, this has expanded from single genes to microbial community genomics, or metagenomics, in which the gene content of an environment can provide not just a census of the community members but direct information on metabolic capabilities and potential interactions among community members. Here we introduce a method for the quantitative characterization and comparison of microbial communities based on the normalization of metagenomic data by estimating average genome sizes. This normalization can relieve comparative biases introduced by differences in community structure, number of sequencing reads, and sequencing read lengths between different metagenomes. We demonstrate the utility of this approach by comparing metagenomes from two different marine sources using both conventional small-subunit (SSU) rRNA gene analyses and our quantitative method to calculate the proportion of genomes in each sample that are capable of a particular metabolic trait. With both environments, to determine what proportion of each community they make up and how differences in environment affect their abundances, we characterize three different types of autotrophic organisms: aerobic, photosynthetic carbon fixers (the Cyanobacteria); anaerobic, photosynthetic carbon fixers (the Chlorobi); and anaerobic, nonphotosynthetic carbon fixers (the Desulfobacteraceae). These analyses demonstrate how genome proportionality compares to SSU rRNA gene relative abundance and how factors such as average genome size and SSU rRNA gene copy number affect sampling probability and therefore both types of community analysis.

Optimizing the analysis of human intestinal microbiota with phylogenetic microarray

Phylogenetic microarrays present an attractive strategy to high-throughput interrogation of complex microbial communities. In this work, we present several approaches to optimize the analysis of intestinal microbiota with the recently developed Microbiota Array. First, we determined how 16S rDNA-specific PCR amplification influenced bacterial detection and the consistency of measured abundance values. Bacterial detection improved with an increase in the number of PCR amplification cycles, but 25 cycles were sufficient to achieve the maximum possible detection. A PCR-caused deviation in the measured abundance values was also observed. We also developed two mathematical algorithms that aimed to account for a predicted cross-hybridization of 16S rDNA fragments among different species, and to adjust the measured hybridization signal based on the number of 16S rRNA gene copies per species genome. The 16S rRNA gene copy adjustment indicated that the presence of members of the class Clostridia might be overestimated in some 16S rDNA-based studies. Finally, we show that the examination of total community RNA with phylogenetic microarray can provide estimates of the relative metabolic activity of individual community members. Complementary profiling of genomic DNA and total RNA isolated from the same sample presents an opportunity to assess population structure and activity in the same microbial community.

Shifts in the host range of a promiscuous plasmid through parallel evolution of its replication initiation protein

The ability of bacterial plasmids to adapt to novel hosts and thereby shift their host range is key to their long-term persistence in bacterial communities. Promiscuous plasmids of the incompatibility group P (IncP)-1 can colonize a wide range of hosts, but it is not known if and how they can contract, shift or further expand their host range. To understand the evolutionary mechanisms of host range shifts of IncP-1 plasmids, an IncP-1β mini-replicon was experimentally evolved in four hosts in which it was initially unstable. After 1000 generations in serial batch cultures under antibiotic selection for plasmid maintenance (kanamycin resistance), the stability of the mini-plasmid dramatically improved in all coevolved hosts. However, only plasmids evolved in Shewanella oneidensis showed improved stability in the ancestor, indicating that adaptive mutations had occurred in the plasmid itself. Complete genome sequence analysis of nine independently evolved plasmids showed seven unique plasmid genotypes that had various kinds of single mutations at one locus, namely, the N-terminal region of the replication initiation protein TrfA. Such parallel evolution indicates that this region was under strong selection. In five of the seven evolved plasmids, these trfA mutations resulted in a significantly higher plasmid copy number. Evolved plasmids were found to be stable in four other naive hosts, but could no longer replicate in Pseudomonas aeruginosa. This study shows that plasmids can specialize to a novel host through trade-offs between improved stability in the new host and the ability to replicate in a previously permissive host.

Refined method for the genomic integration of complex synthetic circuits

Genetic reconstruction of regulatory gene circuits is currently applied in systematic dynamics and structure-function studies of intact cellular networks in systems biology. We present a modified procedure for the integration of complex genetic circuits into the Escherichia coli genome, to provide an efficient synthetic approach for stochastic study and the artificial engineering of genetic networks. Linear artificial sequences of various lengths were easily integrated into the bacterial genome at one time. Comparison of the cellular concentrations of proteins encoded by genes carried on plasmids or the genome indicated that genome recombination could minimize the copy number noise in the genetic circuit, allowing precise design and interpretation of the cellular network. The refined recombination procedure allowed efficient construction of a single copy of a complex genetic circuit in cells, and the resultant reduced fluctuation in copy number led to accurate phenotypic behaviour of the genome-integrated synthetic switch corresponding to the design principle. The availability of long-fragment insertions makes the reconstruction of complex networks easy on the genome, and provides a powerful tool for precise engineering in synthetic and systems biology.

The systemic imprint of growth and its uses in ecological (meta)genomics

Microbial minimal generation times range from a few minutes to several weeks. They are evolutionarily determined by variables such as environment stability, nutrient availability, and community diversity. Selection for fast growth adaptively imprints genomes, resulting in gene amplification, adapted chromosomal organization, and biased codon usage. We found that these growth-related traits in 214 species of bacteria and archaea are highly correlated, suggesting they all result from growth optimization. While modeling their association with maximal growth rates in view of synthetic biology applications, we observed that codon usage biases are better correlates of growth rates than any other trait, including rRNA copy number. Systematic deviations to our model reveal two distinct evolutionary processes. First, genome organization shows more evolutionary inertia than growth rates. This results in over-representation of growth-related traits in fast degrading genomes. Second, selection for these traits depends on optimal growth temperature: for similar generation times purifying selection is stronger in psychrophiles, intermediate in mesophiles, and lower in thermophiles. Using this information, we created a predictor of maximal growth rate adapted to small genome fragments. We applied it to three metagenomic environmental samples to show that a transiently rich environment, as the human gut, selects for fast-growers, that a toxic environment, as the acid mine biofilm, selects for low growth rates, whereas a diverse environment, like the soil, shows all ranges of growth rates. We also demonstrate that microbial colonizers of babies gut grow faster than stabilized human adults gut communities. In conclusion, we show that one can predict maximal growth rates from sequence data alone, and we propose that such information can be used to facilitate the manipulation of generation times. Our predictor allows inferring growth rates in the vast majority of uncultivable prokaryotes and paves the way to the understanding of community dynamics from metagenomic data.

Microbial minimal generation times vary from a few minutes to several weeks. The reasons for this disparity have been thought to lie on different life-history strategies: fast-growing microbes grow extremely fast in rich media, but are less capable of dealing with stress and/or poor nutrient conditions. Prokaryotes have evolved a set of genomic traits to grow fast, including biased codon usage and transient or permanent gene multiplication for dosage effects. Here, we studied the relative role of these traits and show they can be used to predict minimal generation times from the genomic data of the vast majority of microbes that cannot be cultivated. We show that this inference can also be made with incomplete genomes and thus be applied to metagenomic data to test hypotheses about the biomass productivity of biotopes and the evolution of microbiota in the human gut after birth. Our results also allow a better understanding of the co-evolution between growth rates and genomic traits and how they can be manipulated in synthetic biology. Growth rates have been a key variable in microbial physiology studies in the last century, and we show how intimately they are linked with genome organization and prokaryotic ecology.

Nonoptimal microbial response to antibiotics underlies suppressive drug interactions

Suppressive drug interactions, in which one antibiotic can actually help bacterial cells to grow faster in the presence of another, occur between protein and DNA synthesis inhibitors. Here, we show that this suppression results from nonoptimal regulation of ribosomal genes in the presence of DNA stress. Using GFP-tagged transcription reporters in Escherichia coli, we find that ribosomal genes are not directly regulated by DNA stress, leading to an imbalance between cellular DNA and protein content. To test whether ribosomal gene expression under DNA stress is nonoptimal for growth rate, we sequentially deleted up to six of the seven ribosomal RNA operons. These synthetic manipulations of ribosomal gene expression correct the protein-DNA imbalance, lead to improved survival and growth, and completely remove the suppressive drug interaction. A simple mathematical model explains the nonoptimal regulation in different nutrient environments. These results reveal the genetic mechanism underlying an important class of suppressive drug interactions.

The biology of environmental mycobacteria

Although the environmental mycobacteria are slow growing relative to other microorganisms in water and soil which would suggest that they are poor competitors, compensating factors permit survival, growth and persistence in natural and human-engineered environments. Factors such as the hydrophobic, lipid-rich impermeable envelope, biofilm formation, acid resistance, anaerobic survival and metabolism of recalcitrant carbon compounds permit survival and growth of the environmental mycobacteria in a wide range of natural and human-engineered habitats. High numbers of environmental mycobacteria are found in coastal swamps and estuaries and boreal, peat-rich forest soils and waters. The hydrophobic surface results in concentration of the environmental mycobacteria at interfaces (air-water and surface-water) and in aerosolized droplets ejected from water. The survival and growth in protozoa and amoebae permit environmental mycobacteria to persist in habitats subject to predation and likely led to survival and growth in phagocytic cells of animals. Finally, slow growth allows time for mycobacterial cells to adapt to changing conditions before loss of viability.

Visualization of ribosomal RNA operon copy number distribution

Background

Results of microbial ecology studies using 16S rRNA sequence information can be deceiving due to differences in rRNA operon copy number and genome size of the detected organisms. It therefore will be useful for investigators to have a better understanding of how these two parameters differ in various organism types. In this study, the number of ribosomal operons and genome size were separately mapped onto a Bacterial phylogenetic tree.

Results

A representative Bacterial tree was constructed using 31 marker genes found in 578 bacterial genome sequences. Organism names are displayed on the trees using graduations of color such that similar colors indicate similar numbers of operons or genome size. The resulting images provide an intuitive understanding of how copy number and genome size vary in different Bacterial phyla.

Conclusion

Once the phylogenetic position of a novel organism is known the number of rRNA operons, and to a lesser extent the genome size, can be estimated by examination of the colored maps. Further detail can then be obtained for members of relevant taxa from the rrnDB database.

Active transcription of rRNA operons condenses the nucleoid in Escherichia coli: examining the effect of transcription on nucleoid structure in the absence of transertion

In Escherichia coli the genome must be compacted approximately 1,000-fold to be contained in a cellular structure termed the nucleoid. It is proposed that the structure of the nucleoid is determined by a balance of multiple compaction forces and one major expansion force. The latter is mediated by transertion, a coupling of transcription, translation, and translocation of nascent membrane proteins and/or exported proteins. In supporting this notion, it has been shown consistently that inhibition of transertion by the translation inhibitor chloramphenicol results in nucleoid condensation due to the compaction forces that remain active in the cell. Our previous study showed that during optimal growth, RNA polymerase is concentrated into transcription foci or "factories," analogous to the eukaryotic nucleolus, indicating that transcription and RNA polymerase distribution affect the nucleoid structure. However, the interpretation of the role of transcription in the structure of the nucleoid is complicated by the fact that transcription is implicated in both compacting forces and the expansion force. In this work, we used a new approach to further examine the effect of transcription, specifically from rRNA operons, on the structure of the nucleoid, when the major expansion force was eliminated. Our results showed that transcription is necessary for the chloramphenicol-induced nucleoid compaction. Further, an active transcription from multiple rRNA operons in chromosome is critical for the compaction of nucleoid induced by inhibition of translation. All together, our data demonstrated that transcription of rRNA operons is a key mechanism affecting genome compaction and nucleoid structure.

Genome-scale reconstruction of Escherichia coli's transcriptional and translational machinery: a knowledge base, its mathematical formulation, and its functional characterization

Metabolic network reconstructions represent valuable scaffolds for ‘-omics’ data integration and are used to computationally interrogate network properties. However, they do not explicitly account for the synthesis of macromolecules (i.e., proteins and RNA). Here, we present the first genome-scale, fine-grained reconstruction of Escherichia coli's transcriptional and translational machinery, which produces 423 functional gene products in a sequence-specific manner and accounts for all necessary chemical transformations. Legacy data from over 500 publications and three databases were reviewed, and many pathways were considered, including stable RNA maturation and modification, protein complex formation, and iron–sulfur cluster biogenesis. This reconstruction represents the most comprehensive knowledge base for these important cellular functions in E. coli and is unique in its scope. Furthermore, it was converted into a mathematical model and used to: (1) quantitatively integrate gene expression data as reaction constraints and (2) compute functional network states, which were compared to reported experimental data. For example, the model predicted accurately the ribosome production, without any parameterization. Also, in silico rRNA operon deletion suggested that a high RNA polymerase density on the remaining rRNA operons is needed to reproduce the reported experimental ribosome numbers. Moreover, functional protein modules were determined, and many were found to contain gene products from multiple subsystems, highlighting the functional interaction of these proteins. This genome-scale reconstruction of E. coli's transcriptional and translational machinery presents a milestone in systems biology because it will enable quantitative integration of ‘-omics’ datasets and thus the study of the mechanistic principles underlying the genotype–phenotype relationship.

Systems biology aims to understand the interactions of cellular components in a systemic manner. Mathematical modeling is critical to the integration and analysis of these components on a conceptual as well as mechanistic level. To date, detailed genome-scale reconstructions of metabolism have become available for a growing number of organisms. Although metabolism has an important role in cells, other cellular functions need to be considered as well, such as signaling, regulation, and macromolecular synthesis. For instance, the cellular machinery required for RNA and protein synthesis consists of a complex set of proteins. Here, we show that one can collect all of the necessary information for a prokaryotic organism to create a gene-specific, fine-grained representation of the macromolecular synthesis machinery. E. coli was chosen as a model organism because of the wealth of available information. The explicit representation of transcription and translation in terms of a mass-balanced network enables a detailed, quantitative accounting of the protein synthesis capabilities of E. coli in silico. Hence, this study demonstrates the feasibility of constructing very large networks and also represents a critical step toward building cellular models of growth that can account for gene-specific protein production in a stoichiometric fashion on the genome scale.

Evolution of multiple phosphodiesterase isoforms in stickleback involved in cAMP signal transduction pathway

Background

Duplicate genes are considered to have evolved through the partitioning of ancestral functions among duplicates (subfunctionalization) and/or the acquisition of novel functions from a beneficial mutation (neofunctionalization). Additionally, an increase in gene dosage resulting from duplication may also confer an advantageous effect, as has been suggested for histone, tRNA, and rRNA genes. Currently, there is little understanding of the effect of increased gene dosage on subcellular networks like signal transduction pathways. Addressing this issue may provide further insights into the evolution by gene duplication.

Results

We analyzed the evolution of multiple stickleback phosphodiesterase (PDE, EC: 3.1.4.17) 1C genes involved in the cyclic nucleotide signaling pathway. Stickleback has 8–9 copies of this gene, whereas only one or two loci exist in other model vertebrates. Our phylogenetic and synteny analyses suggested that the multiple PDE1C genes in stickleback were generated by repeated duplications of >100-kbp chromosome segments. Sequence evolution analysis did not provide strong evidence for neofunctionalization in the coding sequences of stickleback PDE1C isoforms. On the other hand, gene expression analysis suggested that the derived isoforms acquired expression in new organs, implying their neofunctionalization in terms of expression patterns. In addition, at least seven isoforms of the stickleback PDE1C were co-expressed with olfactory-type G-proteins in the nose, suggesting that PDE1C dosage is increased in the stickleback olfactory transduction (OT) pathway. In silico simulations of OT implied that the increased PDE1C dosage extends the longevity of the depolarization signals of the olfactory receptor neuron.

Conclusion

The predicted effect of the increase in PDE1C products on the OT pathway may play an important role in stickleback behavior and ecology. However, this possibility should be empirically examined. Our analyses imply that an increase in gene product sometimes has a significant, yet unexpected, effect on the functions of subcellular networks.

Three genomes from the phylum Acidobacteria provide insight into the lifestyles of these microorganisms in soils

The complete genomes of three strains from the phylum Acidobacteria were compared. Phylogenetic analysis placed them as a unique phylum. They share genomic traits with members of the Proteobacteria, the Cyanobacteria, and the Fungi. The three strains appear to be versatile heterotrophs. Genomic and culture traits indicate the use of carbon sources that span simple sugars to more complex substrates such as hemicellulose, cellulose, and chitin. The genomes encode low-specificity major facilitator superfamily transporters and high-affinity ABC transporters for sugars, suggesting that they are best suited to low-nutrient conditions. They appear capable of nitrate and nitrite reduction but not N(2) fixation or denitrification. The genomes contained numerous genes that encode siderophore receptors, but no evidence of siderophore production was found, suggesting that they may obtain iron via interaction with other microorganisms. The presence of cellulose synthesis genes and a large class of novel high-molecular-weight excreted proteins suggests potential traits for desiccation resistance, biofilm formation, and/or contribution to soil structure. Polyketide synthase and macrolide glycosylation genes suggest the production of novel antimicrobial compounds. Genes that encode a variety of novel proteins were also identified. The abundance of acidobacteria in soils worldwide and the breadth of potential carbon use by the sequenced strains suggest significant and previously unrecognized contributions to the terrestrial carbon cycle. Combining our genomic evidence with available culture traits, we postulate that cells of these isolates are long-lived, divide slowly, exhibit slow metabolic rates under low-nutrient conditions, and are well equipped to tolerate fluctuations in soil hydration.

Phenotypic and genotypic analyses of clinical Fusobacterium nucleatum and Fusobacterium periodonticum isolates from the human gut

Fusobacterium nucleatum is a Gram-negative anaerobic rod that is part of the normal human microflora, and has also been associated with various infections. Bacterial strains belonging to the species are typically heterogeneous in both phenotype and genotype, which can hinder their identification in a clinical setting. The majority of F. nucleatum isolates originate from oral sites, however the species is also a resident of the human gastrointestinal tract. The aim of this study was to compare F. nucleatum isolates from human intestinal biopsy samples to try and determine whether isolates from this site are divergent from oral isolates. We used a variety of phenotypic and genotypic markers to compare 21 F. nucleatum and Fusobacterium periodonticum isolates from the GI tract to oral isolates and recognized type strains in order to study heterogeneity within this set. 16S rDNA and rpoB gene sequence analysis allowed us to build phylogenetic trees that consistently placed isolates into distinct clusters. 16S rDNA copy number analyses using Denaturing Gradient Gel Electrophoresis (DGGE) demonstrated potential for use as a method to examine clonality amongst species. Phenotypic analyses gave variable results that were generally unhelpful in distinguishing between phylogenetic clusters. Our results suggest that a) F. periodonticum isolates are not restricted to the oral niche; b) phenotypic classification is not sufficient to subspeciate isolates; c) heterogeneity within the species is extensive but constrained; and d) F. nucleatum isolates from the gut tend to identify with the animalis subspecies.

Small-scale copy number variation and large-scale changes in gene expression

The expression dynamics of interacting genes depends, in part, on the structure of regulatory networks. Genetic regulatory networks include an overrepresentation of subgraphs commonly known as network motifs. In this article, we demonstrate that gene copy number is an omnipresent parameter that can dramatically modify the dynamical function of network motifs. We consider positive feedback, bistable feedback, and toggle switch motifs and show that variation in gene copy number, on the order of a single or few copies, can lead to multiple orders of magnitude change in gene expression and, in some cases, switches in deterministic control. Further, small changes in gene copy number for a 3-gene motif with successive inhibition (the "repressilator") can lead to a qualitative switch in system behavior among oscillatory and equilibrium dynamics. In all cases, the qualitative change in expression is due to the nonlinear nature of transcriptional feedback in which duplicated motifs interact via common pools of transcription factors. We are able to implicitly determine the critical values of copy number which lead to qualitative shifts in system behavior. In some cases, we are able to solve for the sufficient condition for the existence of a bifurcation in terms of kinetic rates of transcription, translation, binding, and degradation. We discuss the relevance of our findings to ongoing efforts to link copy number variation with cell fate determination by viruses, dynamics of synthetic gene circuits, and constraints on evolutionary adaptation.

A coarse-grained biophysical model of E. coli and its application to perturbation of the rRNA operon copy number

We propose a biophysical model of Escherichia coli that predicts growth rate and an effective cellular composition from an effective, coarse-grained representation of its genome. We assume that E. coli is in a state of balanced exponential steady-state growth, growing in a temporally and spatially constant environment, rich in resources. We apply this model to a series of past measurements, where the growth rate and rRNA-to-protein ratio have been measured for seven E. coli strains with an rRNA operon copy number ranging from one to seven (the wild-type copy number). These experiments show that growth rate markedly decreases for strains with fewer than six copies. Using the model, we were able to reproduce these measurements. We show that the model that best fits these data suggests that the volume fraction of macromolecules inside E. coli is not fixed when the rRNA operon copy number is varied. Moreover, the model predicts that increasing the copy number beyond seven results in a cytoplasm densely packed with ribosomes and proteins. Assuming that under such overcrowded conditions prolonged diffusion times tend to weaken binding affinities, the model predicts that growth rate will not increase substantially beyond the wild-type growth rate, as indicated by other experiments. Our model therefore suggests that changing the rRNA operon copy number of wild-type E. coli cells growing in a constant rich environment does not substantially increase their growth rate. Other observations regarding strains with an altered rRNA operon copy number, such as nucleoid compaction and the rRNA operon feedback response, appear to be qualitatively consistent with this model. In addition, we discuss possible design principles suggested by the model and propose further experiments to test its validity.

A bacterium like E. coli can be thought of as a self-replicating factory, where inventory synthesis, degradation, and management is concerted according to a well-defined set of rules encoded in the organism's genome. Since the organism's survival depends on this set of rules, these rules were most likely optimized by evolution. Therefore, by writing down these rules, what could one learn about Escherichia coli? We examined E. coli growing in the simplest imaginable environment, one constant in space and time and rich in resources, and attempted to identify the rules that relate the genome to the cell composition and self-replication time. With more than 4,400 genes, a full-scale model would be prohibitively complicated, and therefore we “coarse-grained” E. coli by lumping together genes and proteins of similar function. We used this model to examine measurements of strains with reduced copy number of ribosomal-RNA genes, and also to show that increasing this copy number overcrowds the cell with ribosomes and proteins. As a result, there appears to be an optimum copy number with respect to the wild-type genome, in agreement with observation. We hope that this model will improve and further challenge our understanding of bacterial physiology, also in more complicated environments.

Genome of the epsilonproteobacterial chemolithoautotroph Sulfurimonas denitrificans

Sulfur-oxidizing epsilonproteobacteria are common in a variety of sulfidogenic environments. These autotrophic and mixotrophic sulfur-oxidizing bacteria are believed to contribute substantially to the oxidative portion of the global sulfur cycle. In order to better understand the ecology and roles of sulfur-oxidizing epsilonproteobacteria, in particular those of the widespread genus Sulfurimonas, in biogeochemical cycles, the genome of Sulfurimonas denitrificans DSM1251 was sequenced. This genome has many features, including a larger size (2.2 Mbp), that suggest a greater degree of metabolic versatility or responsiveness to the environment than seen for most of the other sequenced epsilonproteobacteria. A branched electron transport chain is apparent, with genes encoding complexes for the oxidation of hydrogen, reduced sulfur compounds, and formate and the reduction of nitrate and oxygen. Genes are present for a complete, autotrophic reductive citric acid cycle. Many genes are present that could facilitate growth in the spatially and temporally heterogeneous sediment habitat from where Sulfurimonas denitrificans was originally isolated. Many resistance-nodulation-development family transporter genes (10 total) are present; of these, several are predicted to encode heavy metal efflux transporters. An elaborate arsenal of sensory and regulatory protein-encoding genes is in place, as are genes necessary to prevent and respond to oxidative stress.

Identification of Aeromonas veronii genes required for colonization of the medicinal leech, Hirudo verbana

Most digestive tracts contain a complex consortium of beneficial microorganisms, making it challenging to tease apart the molecular interactions between symbiont and host. The digestive tract of Hirudo verbana, the medicinal leech, is an ideal model system because it harbors a simple microbial community in the crop, comprising the genetically amenable Aeromonas veronii and a Rikenella-like bacterium. Signature-tagged mutagenesis (STM) was used to identify genes required for digestive tract colonization. Of 3,850 transposon (Tn) mutants screened, 46 were identified as colonization mutants. Previously we determined that the complement system of the ingested blood remained active inside the crop and prevented serum-sensitive mutants from colonizing. The identification of 26 serum-sensitive mutants indicated a successful screen. The remaining 20 serum-resistant mutants are described in this study and revealed new insights into symbiont-host interactions. An in vivo competition assay compared the colonization levels of the mutants to that of a wild-type competitor. Attenuated colonization mutants were grouped into five classes: surface modification, regulatory, nutritional, host interaction, and unknown function. One STM mutant, JG736, with a Tn insertion in lpp, encoding Braun's lipoprotein, was characterized in detail. This mutant had a >25,000-fold colonization defect relative to colonization by the wild-type strain at 72 h and, in vitro, an increased sensitivity to sodium dodecyl sulfate, suggesting the presence of an additional antimicrobial property in the crop. The classes of genes identified in this study are consistent with findings from previous STM studies involving pathogenic bacteria, suggesting parallel molecular requirements for beneficial and pathogenic host colonization.

Isolation and characterization of soil bacteria that define Terriglobus gen. nov., in the phylum Acidobacteria

Bacteria in the phylum Acidobacteria are widely distributed and abundant in soils, but their ecological roles are poorly understood, owing in part to a paucity of cultured representatives. In a molecular survey of acidobacterial diversity at the Michigan State University Kellogg Biological Station Long-Term Ecological Research site, 27% of acidobacterial 16S rRNA gene clones in a never-tilled, successional plant community belonged to subdivision 1, whose relative abundance varied inversely with soil pH. Strains of subdivision 1 were isolated from these never-tilled soils using low-nutrient medium incubated for 3 to 4 weeks under elevated levels of carbon dioxide, which resulted in a slightly acidified medium that matched the pH optima of the strains (between 5 and 6). Colonies were approximately 1 mm in diameter and either white or pink, the latter due to a carotenoid(s) that was synthesized preferentially under 20% instead of 2% oxygen. Strains were gram-negative, aerobic, chemo-organotrophic, nonmotile rods that produced an extracellular matrix. All strains contained either one or two copies of the 16S rRNA encoding gene, which along with a relatively slow doubling time (10 to 15 h at ca. 23 degrees C) is suggestive of an oligotrophic lifestyle. Six of the strains are sufficiently similar to one another, but distinct from previously named Acidobacteria, to warrant creation of a new genus, Terriglobus, with Terriglobus roseus defined as the type species. The physiological and nutritional characteristics of Terriglobus are consistent with its potential widespread distribution in soil.

Performance of the translational apparatus varies with the ecological strategies of bacteria

Protein synthesis is the predominant activity of growing bacteria; the protein synthesis system accounts for more than one-half the cell's dry mass and consumes most of the cell's energy during rapid growth. Translation has been studied extensively using model organisms, and the translational apparatus is qualitatively similar in terms of structure and function across all known forms of life. However, little is known about variation between organisms in translational performance. Using measurements of macromolecular content in a phylogenetically diverse collection of bacteria with contrasting ecological strategies, we found that the translational power (the rate of protein synthesis normalized to the mass of the protein synthesis system) is three- to fourfold higher among bacteria that respond rapidly to nutrient availability than among bacteria that respond slowly. An analysis of codon use in completely sequenced bacterial genomes confirmed that the selective forces acting on translation vary with the ecological strategy. We propose that differences in translational power result from ecologically based variation among microbes in the relative importance of two competing benefits: reducing the biomass invested in the protein synthesis system and reducing the energetic expense of protein synthesis.