Determination of DNA content of aquatic bacteria by flow cytometry

The distribution of DNA among bacterioplankton and bacterial isolates was determined by flow cytometry of DAPI (4',6'-diamidino-2-phenylindole)-stained organisms. Conditions were optimized to minimize error from nonspecific staining, AT bias, DNA packing, changes in ionic strength, and differences in cell permeability. The sensitivity was sufficient to characterize the small 1- to 2-Mb-genome organisms in freshwater and seawater, as well as low-DNA cells ("dims"). The dims could be formed from laboratory cultivars; their apparent DNA content was 0.1 Mb and similar to that of many particles in seawater. Preservation with formaldehyde stabilized samples until analysis. Further permeabilization with Triton X-100 facilitated the penetration of stain into stain-resistant lithotrophs. The amount of DNA per cell determined by flow cytometry agreed with mean values obtained from spectrophotometric analyses of cultures. Correction for the DNA AT bias of the stain was made for bacterial isolates with known G+C contents. The number of chromosome copies per cell was determined with pure cultures, which allowed growth rate analyses based on cell cycle theory. The chromosome ratio was empirically related to the rate of growth, and the rate of growth was related to nutrient concentration through specific affinity theory to obtain a probe for nutrient kinetics. The chromosome size of a Marinobacter arcticus isolate was determined to be 3.0 Mb by this method. In a typical seawater sample the distribution of bacterial DNA revealed two major populations based on DNA content that were not necessarily similar to populations determined by using other stains or protocols. A mean value of 2.5 fg of DNA cell(-1) was obtained for a typical seawater sample, and 90% of the population contained more than 1.1 fg of DNA cell(-1).     

Ongoing evolution of strand composition in bacterial genomes

We tried to identify the substitutions involved in the establishment of replication strand bias, which has been recognized as an important evolutionary factor in the evolution of bacterial genomes. First, we analyzed the composition asymmetry of 28 complete bacterial genomes and used it to test the possibility that asymmetric deamination of cytosine might be at the origin of the bias. The model showed significant correlation to the data but left unexplained a significant portion of the variance and indicated a systematic underestimation of GC skews in comparison with TA skews. Second, we analyzed the substitutions acting on the genes from five fully sequenced Chlamydia genomes that had not suffered strand switch since speciation. This analysis showed that substitutions were not at equilibrium in Chlamydia trachomatis or in C. muridarum and that strand bias is still an on-going process in these genes. Third, we identified substitutions involved in the adaptation of genes that had switched strands after speciation. These genes adapted quickly to the skewed composition of the new strand, mostly due to C-->T, A-->G, and C-->G asymmetric substitutions. This observation was reinforced by the analysis of genes that switched strands after divergence between Bacillus subtilis and B. halodurans. Finally, we propose a more extended model based on the analysis of the substitution asymmetries of CHLAMYDIA: This model fits well with the data provided by bacterial genomes presenting strong strand bias.     

Bias in Template-to-Product Ratios in Multitemplate PCR

Bias introduced by the simultaneous amplification of specific genes from complex mixtures of templates remains poorly understood. To explore potential causes and the extent of bias in PCR amplification of 16S ribosomal DNAs (rDNAs), genomic DNAs of two closely and one distantly related bacterial species were mixed and amplified with universal, degenerate primers. Quantification and comparison of template and product ratios showed that there was considerable and reproducible overamplification of specific templates. Variability between replicates also contributed to the observed bias but in a comparatively minor way. Based on these initial observations, template dosage and differences in binding energies of permutations of the degenerate, universal primers were tested as two likely causes of this template-specific bias by using 16S rDNA templates modified by site-directed mutagenesis. When mixtures of mutagenized templates containing AT- and GC-rich priming sites were used, templates containing the GC-rich permutation amplified with higher efficiency, indicating that different primer binding energies may to a large extent be responsible for overamplification. In contrast, gene copy number was found to be an unlikely cause of the observed bias. Similarly, amplification from DNA extracted from a natural community to which different amounts of genomic DNA of a single bacterial species were added did not affect relative product ratios. Bias was reduced considerably by using high template concentrations, by performing fewer cycles, and by mixing replicate reaction preparations.     

Effect of natamycin on the enumeration, genetic structure and composition of bacterial community isolated from soils and soybean rhizosphere

Natamycin is commonly used to control fungal growth on agar media used for bacterial enumeration or strain isolation. However, there is no conclusive report on the possible effect of this antibiotic on bacterial growth or on the diversity of the recovered soil bacteria. Therefore, the possible effects of natamycin on the numbers of bacteria isolated at 12 degrees C from three different soils and soybean rhizosphere soil were investigated using natamycin concentrations ranging from 0 to 200 mg l(-1). Our results demonstrate that natamycin concentrations, which inhibit the growth of fungi on the media, have a small but significant inhibitory effect on the number of bacterial colony forming units. A natamycin concentration of 50-200 mg l(-1) is required for an efficient control of fungal growth on media in our experimental conditions depending on the soil type. Bacterial community structure was assessed on culturable cells (cells washed from enumeration plates: plate-wash approach) obtained at 12 degrees C from soybean rhizosphere soil by performing Ribosomal Intergenic Spacer Analysis (RISA) fingerprinting. We demonstrate that all natamycin concentrations used alter the structure of the recovered, culturable bacterial community, compared to control without natamycin. Using ARDRA (amplification of the 16S rDNA gene and restriction analysis) genotyping of individual isolates, some differences were observed between the bacterial isolates obtained in the presence or absence of natamycin. Bacterial isolates recovered in the presence of natamycin are more tolerant (maximal growth rate and lag phase) to this compound than those isolated without natamycin, indicating a possible selection of resistant strains. Therefore, high concentration of natamycin cannot be used for isolation of bacterial strains with the aim of studying biodiversity and could bias a selection of strains for practical applications.     

Measuring and mitigating PCR bias in microbiota datasets

PCR amplification plays an integral role in the measurement of mixed microbial communities via high-throughput DNA sequencing of the 16S ribosomal RNA (rRNA) gene. Yet PCR is also known to introduce multiple forms of bias in 16S rRNA studies. Here we present a paired modeling and experimental approach to characterize and mitigate PCR NPM-bias (PCR bias from non-primer-mismatch sources) in microbiota surveys. We use experimental data from mock bacterial communities to validate our approach and human gut microbiota samples to characterize PCR NPM-bias under real-world conditions. Our results suggest that PCR NPM-bias can skew estimates of microbial relative abundances by a factor of 4 or more, but that this bias can be mitigated using log-ratio linear models.     

Biased Gene Conversion, Copy Number, and Apparent Mutation Rate Differences within Chloroplast and Bacterial Genomes

We investigate the possibility that differences between synonymous substitution rates of organelle and bacterial genes differing only in copy number may be due to conversion bias. We find that the rather large observed difference in the synonymous rates between genes in the single copy and inverted-repeat regions of chloroplasts can be accounted for by a very small bias against new mutants. More generally, differences in the within-organelle fixation probability result in different apparent mutation rates as measured by the expected rate of appearance of cells homoplasmic for new mutants. Thus, differences in intracellular population parameters rather than molecular mechanisms can account for some variation in the apparent mutation rates of organelle genes, and possibly in other systems with variable numbers of gene copies. On the other hand, our analysis suggests that conversion bias is not a likely explanation for relatively low mutation rates observed near the replication origin of bacterial chromosomes.     

Noise underlies switching behavior of the bacterial flagellum

We report the switching behavior of the full bacterial flagellum system that includes the filament and the motor in wild-type Escherichia coli cells. In sorting the motor behavior by the clockwise bias, we find that the distributions of the clockwise (CW) and counterclockwise (CCW) intervals are either exponential or nonexponential with long tails. At low bias, CW intervals are exponentially distributed and CCW intervals exhibit long tails. At intermediate CW bias (0.5) both CW and CCW intervals are mainly exponentially distributed. A simple model suggests that these two distinct switching behaviors are governed by the presence of signaling noise within the chemotaxis network. Low noise yields exponentially distributed intervals, whereas large noise yields nonexponential behavior with long tails. These drastically different motor statistics may play a role in optimizing bacterial behavior for a wide range of environmental conditions.     

Extraction and purification of DNA in rhizosphere soil samples for PCR-DGGE analysis of bacterial consortia

Application of DNA fingerprinting methods enables the detection of diverse members of soil bacterial consortia, even including those bacteria not yet cultivated. However, extraction and purification of DNA from soil samples without bias is difficult. We compared five different DNA isolation methods and three purification methods for rhizosphere soil samples. Purified DNA extracts were amplified in PCR using universal bacterial primers and the PCR products were analysed with denaturing gradient gel electrophoresis (DGGE) for the visualisation of DNA bands representing dominant bacterial species. Both the isolation and purification methods affected the apparent bacterial community structure of the samples.     

Determination of DNA Content of Aquatic Bacteria by Flow Cytometry

The distribution of DNA among bacterioplankton and bacterial isolates was determined by flow cytometry of DAPI (4′,6′-diamidino-2-phenylindole)-stained organisms. Conditions were optimized to minimize error from nonspecific staining, AT bias, DNA packing, changes in ionic strength, and differences in cell permeability. The sensitivity was sufficient to characterize the small 1- to 2-Mb-genome organisms in freshwater and seawater, as well as low-DNA cells (“dims”). The dims could be formed from laboratory cultivars; their apparent DNA content was 0.1 Mb and similar to that of many particles in seawater. Preservation with formaldehyde stabilized samples until analysis. Further permeabilization with Triton X-100 facilitated the penetration of stain into stain-resistant lithotrophs. The amount of DNA per cell determined by flow cytometry agreed with mean values obtained from spectrophotometric analyses of cultures. Correction for the DNA AT bias of the stain was made for bacterial isolates with known G+C contents. The number of chromosome copies per cell was determined with pure cultures, which allowed growth rate analyses based on cell cycle theory. The chromosome ratio was empirically related to the rate of growth, and the rate of growth was related to nutrient concentration through specific affinity theory to obtain a probe for nutrient kinetics. The chromosome size of a Marinobacter arcticus isolate was determined to be 3.0 Mb by this method. In a typical seawater sample the distribution of bacterial DNA revealed two major populations based on DNA content that were not necessarily similar to populations determined by using other stains or protocols. A mean value of 2.5 fg of DNA cell⁻¹ was obtained for a typical seawater sample, and 90% of the population contained more than 1.1 fg of DNA cell⁻¹.     

Deletional bias and the evolution of bacterial genomes

Although bacteria increase their DNA content through horizontal transfer and gene duplication, their genomes remain small and, in particular, lack nonfunctional sequences. This pattern is most readily explained by a pervasive bias towards higher numbers of deletions than insertions. When selection is not strong enough to maintain them, genes are lost in large deletions or inactivated and subsequently eroded. Gene inactivation and loss are particularly apparent in obligate parasites and symbionts, in which dramatic reductions in genome size can result not from selection to lose DNA, but from decreased selection to maintain gene functionality. Here we discuss the evidence showing that deletional bias is a major force that shapes bacterial genomes.     

Populational Heterogeneity vs. Temporal Fluctuation in Escherichia coli Flagellar Motor Switching

The dynamic switching of the bacterial flagellar motor regulates cell motility in bacterial chemotaxis. It has been reported under physiological conditions that the switching bias of the flagellar motor undergoes large temporal fluctuations, which reflects noise propagating in the chemotactic signaling network. On the other hand, nongenetic heterogeneity is also observed in flagellar motor switching, as a large group of switching motors show different switching bias and frequency under the same physiological condition. In this work, we present simultaneous measurement of groups of Escherichia coli flagellar motor switching and compare them to long time recording of single switching motors. Consistent with previous studies, we observed temporal fluctuations in switching bias in long time recording experiments. However, the variability in switching bias at the populational level showed much higher volatility than its temporal fluctuation. These results suggested stable individuality in E. coli motor switching. We speculate that uneven expression of key regulatory proteins with amplification by the ultrasensitive response of the motor can account for the observed populational heterogeneity and temporal fluctuations.     

Populational Heterogeneity vs. Temporal Fluctuation in Escherichia coli Flagellar Motor Switching

The dynamic switching of the bacterial flagellar motor regulates cell motility in bacterial chemotaxis. It has been reported under physiological conditions that the switching bias of the flagellar motor undergoes large temporal fluctuations, which reflects noise propagating in the chemotactic signaling network. On the other hand, nongenetic heterogeneity is also observed in flagellar motor switching, as a large group of switching motors show different switching bias and frequency under the same physiological condition. In this work, we present simultaneous measurement of groups of Escherichia coli flagellar motor switching and compare them to long time recording of single switching motors. Consistent with previous studies, we observed temporal fluctuations in switching bias in long time recording experiments. However, the variability in switching bias at the populational level showed much higher volatility than its temporal fluctuation. These results suggested stable individuality in E. coli motor switching. We speculate that uneven expression of key regulatory proteins with amplification by the ultrasensitive response of the motor can account for the observed populational heterogeneity and temporal fluctuations.     

Cloning and sequencing of an Agrobacterium tumefaciens beta-glucosidase gene involved in modifying a vir-inducing plant signal molecule.

Induction of Agrobacterium tumefaciens virulence genes by plant phenolic compounds is essential for successful T-DNA transfer to a host plant. In Douglas fir needles, the major virulence region inducer is the glycoside coniferin (J. W. Morris and R. O. Morris, Proc. Natl. Acad. Sci. USA 87:3612-3618, 1990). Agrobacterium strains with high beta-glucosidase activity respond to coniferin and infect Douglas fir seedlings, whereas most strains with low beta-glucosidase activity fail to respond to coniferin and are avirulent on this host. We have cloned two beta-glucosidase genes from A. tumefaciens B3/73 and sequenced one of them, cbg1. It appears to be part of a polycistronic unit and shows a high bias for GC-rich codons. When expressed in Escherichia coli, Cbg1 beta-glucosidase hydrolyzes coniferin but not cellobiose. The 88-kDa predicted product of cbg1 is highly similar to one other bacterial beta-glucosidase and several fungal beta-glucosidases. There is little homology between Cbg1 and other bacterial beta-glucosidases, including an Agrobacterium cellobiase.     

Effects of organic antagonists of Ca(2+), Na(+), and K(+) on chemotaxis and motility of escherichia coli

Various Ca(2+) antagonists used in animal research, many of them known to be Ca(2+) channel blockers, inhibited Escherichia coli chemotaxis (measured as entry of cells into a capillary containing attractant). The most effective of these, acting in the nanomolar range, was omega-conotoxin GVIA. The next most effective were gallopamil and verapamil. At concentrations around 100-fold higher than that needed for inhibition of chemotaxis, each of these antagonists inhibited motility (measured as entry of cells into a capillary lacking attractant). Various other Ca(2+) antagonists were less effective, though chemotaxis was almost always more sensitive to inhibition than was motility. Cells treated with each of these Ca(2+) antagonists swam with a running bias, i.e., tumbling was inhibited. Similarly, some Na(+) antagonists used in animal research inhibited bacterial chemotaxis. E. coli chemotaxis was inhibited by saxitoxin at concentrations above 10(-7) M, while more than 10(-4) M was needed to inhibit motility. Cells treated with saxitoxin swam with a tumbling bias. In the case of other Na(+) antagonists in animals, aconitine inhibited bacterial chemotaxis 10 times more effectively than it inhibited motility, and two others inhibited chemotaxis and motility at about the same concentration. In the case of K(+) antagonists used in animal research, 4-aminopyridine blocked E. coli chemotaxis between 10(-3) M and, totally, 10(-2) M, while motility was not affected at 10(-2) M; on the other hand, tetraethylammonium chloride failed to inhibit either chemotaxis or motility at 10(-2) M.     

Survey of bromodeoxyuridine uptake among environmental bacteria and variation in uptake rates in a taxonomically diverse set of bacterial isolates

Incorporation of 5-Bromo-2′-Deoxyuridine (BrdU) into DNA can be used to target replicating bacteria in the environment, but differential uptake capacity is a potential bias. Among 23 bacterial isolates commonly found in soils, most took up BrdU, but at up to 10-fold different cell-specific rates. Combined with results from an in silico analysis of 1000 BrdU-labeled 16S rRNA gene sequences, our results demonstrate a BrdU uptake bias with no apparent relationship between taxa affiliation and ability to incorporate BrdU.     

A large-scale analysis of codon usage bias in 4868 bacterial genomes shows association of codon adaptation index with GC content,protein functional domains and bacterial phenotypes

Multiple synonymous codons code for the same amino acid, resulting in the degeneracy of the genetic code and in the preferred used of some codons called codon bias usage (CBU). We performed a large-scale analysis of codon usage bias analysing the distribution of the codon adaptation index (CAI) and the codon relative adaptiveness index (RA) in 4868 bacterial genomes. We found that CAI values differ significantly between protein functional domains and part of the protein outside domains and show how CAI, GC content and preferred usage of polymerase III alpha subunits are related. Additionally, we give evidence of the association between CAI and bacterial phenotypes.     

Regulation of Switching Frequency and Bias of the Bacterial Flagellar Motor by CheY and Fumarate

The effect of CheY and fumarate on switching frequency and rotational bias of the bacterial flagellar motor was analyzed by computer-aided tracking of tethered Escherichia coli. Plots of cells overexpressing CheY in a gutted background showed a bell-shaped correlation curve of switching frequency and bias centering at about 50% clockwise rotation. Gutted cells (i.e., with cheA to cheZ deleted) with a low CheY level but a high cytoplasmic fumarate concentration displayed the same correlation of switching frequency and bias as cells overexpressing CheY at the wild-type fumarate level. Hence, a high fumarate level can phenotypically mimic CheY overexpression by simultaneously changing the switching frequency and the bias. A linear correlation of cytoplasmic fumarate concentration and clockwise rotation bias was found and predicts exclusively counterclockwise rotation without switching when fumarate is absent. This suggests that (i) fumarate is essential for clockwise rotation in vivo and (ii) any metabolically induced fluctuation of its cytoplasmic concentration will result in a transient change in bias and switching probability. A high fumarate level resulted in a dose-response curve linking bias and cytoplasmic CheY concentration that was offset but with a slope similar to that for a low fumarate level. It is concluded that fumarate and CheY act additively presumably at different reaction steps in the conformational transition of the switch complex from counterclockwise to clockwise motor rotation.     

On Ribosome Load,Codon Bias and Protein Abundance

Different codons encoding the same amino acid are not used equally in protein-coding sequences. In bacteria, there is a bias towards codons with high translation rates. This bias is most pronounced in highly expressed proteins, but a recent study of synthetic GFP-coding sequences did not find a correlation between codon usage and GFP expression, suggesting that such correlation in natural sequences is not a simple property of translational mechanisms. Here, we investigate the effect of evolutionary forces on codon usage. The relation between codon bias and protein abundance is quantitatively analyzed based on the hypothesis that codon bias evolved to ensure the efficient usage of ribosomes, a precious commodity for fast growing cells. An explicit fitness landscape is formulated based on bacterial growth laws to relate protein abundance and ribosomal load. The model leads to a quantitative relation between codon bias and protein abundance, which accounts for a substantial part of the observed bias for E. coli. Moreover, by providing an evolutionary link, the ribosome load model resolves the apparent conflict between the observed relation of protein abundance and codon bias in natural sequences and the lack of such dependence in a synthetic gfp library. Finally, we show that the relation between codon usage and protein abundance can be used to predict protein abundance from genomic sequence data alone without adjustable parameters.     

Long-range correlations in prokaryotic chromosomes

One of the fascinating properties of the DNA sequences of prokaryotic and eukaryotic chromosomes is that they possess long-range order. Computational methods like spectral analysis, mutual information and DNA random walks have been used to probe long-range order via-long range correlations. This work attempts to show the advantage of using the Information Theoretic measure of mutual information for this purpose. A number Mu is found which indicates the existence of long-range order. Mu is the ratio between the value of mutual information function between two nucleotides of a DNA sequence separated by a large distance of 100 kilobases to the value expected from a randomized sequence of the same DNA. It is found that in spite of the constant shuffling of nucleotides due to insertion, deletion, inversion and recombination that occur during evolution, the chromosomal structure of prokaryotes is not always mosaic. While all archaeal chromosomes show mosaic structure and lack long-range order, a sizable fraction of the bacterial chromosomes do possess long-range order. A statistical multivariate analysis has been done to find which of the physical variables like genome size or GC% affects the organization of the chromosome or correlates with the long-range order. The existence of long-range order in bacterial chromosomes could be directly correlated to the degree of gene strand bias shown by it. Firmicutes which have low GC content also have pronounced strand bias and show long-range correlations. It is observed that the occurrence of long-range order in bacteria is independent of genome size, but depends on its GC content and gene strand bias.