Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2005, p. 6870-6877, Vol. 71, No. 11
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.11.6870-6877.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Department of Microbiology, University of Massachusetts, Amherst, Massachusetts 01003
Received 19 April 2005/ Accepted 29 June 2005
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
An alternative approach might be to correlate levels of mRNA transcripts for key genes involved in Fe(III) reduction with rates of Fe(III) reduction and then measure their transcript levels in the environment (13). This is a possibility because even though there is a wide phylogenetic diversity of dissimilatory Fe(III)-reducing microorganisms with different strategies for Fe(III) reduction (31), molecular studies have demonstrated that Geobacteraceae are the predominant microorganisms in many sedimentary environments in which Fe(III) reduction is an important process (4, 21, 39-41), as well as on the surfaces of electrodes harvesting electricity from aquatic sediments (9, 20, 43). Thus, if it were possible to identify a gene that was highly conserved within the Geobacteraceae family with expression patterns that correlated with rates of Fe(III) reduction, it should theoretically be possible to estimate Fe(III) reduction rates by monitoring the expression of this gene.
Previous studies have suggested that it is possible to monitor the gene expression of Geobacteraceae in subsurface environments (23). For example, Geobacteraceae nifD transcript levels could be measured in Fe(III)-reducing subsurface sediments (23). These nitrogen fixation studies were facilitated by the fact that the genes involved in the N2 fixation pathway are well understood for a diversity of microorganisms. However, Fe(III) respiration is not as well understood, making it difficult to identify the appropriate gene(s) to monitor in the environment. Initial studies of Fe(III) reduction rates focused on Geobacter sulfurreducens, a heavily studied member of the Geobacteraceae. This organism was selected because it is closely related to the Geobacter species that are most prevalent in many Fe(III)-reducing environments (31) and because both this organism's complete genome sequence (32) and a genetic system (15) are available.
Analysis of chemostat cultures of G. sulfurreducens demonstrated that there was a direct correlation between levels of transcripts for fumarate reductase and the rate of fumarate respiration (13). More importantly, levels of transcripts for omcB, a gene which encodes an outer membrane c-type cytochrome required for Fe(III) reduction (27), were directly correlated with rates of Fe(III) reduction (13). However, the level of omcB transcripts was also significantly influenced by whether metabolism was limited by electron donor or Fe(III). These results suggest that environmental conditions could also have an important influence on the expression of omcB (13). Furthermore, although omcB is an essential gene for Fe(III) reduction in G. sulfurreducens, this gene has not been detected in the genome of the closely related species Geobacter metallireducens (J. Butler, D. R. Lovley, et al., unpublished data) or in genomic DNAs extracted from Geobacteraceae-dominated environments. Thus, monitoring levels ofomcB transcripts is unlikely to be a useful strategy for evaluating rates of Fe(III) reduction in microbial communities comprised of multiple Geobacter species.
An alternative to monitoring the expression of genes directly involved in Fe(III) respiration would be to evaluate rates of central metabolism, which provides the electrons that are required for dissimilatory Fe(III) reduction. Acetate is expected to be the primary electron donor for Fe(III) reduction in a variety of sedimentary environments (30) as well as during in situ bioremediation of uranium (5). Citrate synthase is a key gene in the incorporation of acetate into the tricarboxylic acid (TCA) cycle (19, 47), and the genome of G. sulfurreducens was found to contain a novel putative citrate synthase that is more closely related to the citrate synthases of eukaryotic microorganisms than to those of prokaryotes (32). The gene product was purified and found to possess the only citrate synthase activity in G. sulfurreducens (11). All members of the Geobacteraceae examined contained a similar citrate synthase gene (11), designated gltA. We report that the levels of transcripts for gltA correspond with the rates of acetate metabolism and Fe(III) reduction in chemostat cultures of G. sulfurreducens and of current production when G. sulfurreducens is growing with an electrode serving as the sole electron acceptor. Furthermore, Geobacteraceae-specific gltA genes can be detected in a variety of environments in which Geobacteraceae are important members of the microbial community, and levels of gltA transcripts in the environment are related to the availability of acetate to support Fe(III) reduction.
| MATERIALS AND METHODS |
|---|
|
|
|---|
G. sulfurreducens grown under chemostat conditions.
G. sulfurreducens was grown in electron donor-limited chemostats, with acetate (5 mM) provided as the electron donor and Fe(III) citrate (55 mM) provided as the electron acceptor, at 30°C as previously described (18). G. sulfurreducens was grown in triplicate chemostats at dilution rates of 0.03 h1 and 0.07 h1. Analyses of Fe(II), acetate, and protein were performed as previously described (18).
Growth of G. sulfurreducens on an electrode.
G. sulfurreducens cells that had been grown with acetate (10 mM) as the electron donor and Fe(III) citrate (55 mM) as the electron acceptor were harvested by centrifugation at 3,150 x g at 4°C and resuspended in anoxic medium lacking an electron donor or acceptor. This cell suspension served as an inoculum for the anaerobic anodic chamber (200 ml of medium) of a two-chambered electrode system constructed as previously described (9, 10). The electron acceptor provided for growth in the anode chamber consisted of a graphite electrode poised with a potentiostat (AMEL Instruments, Milan, Italy) at a constant potential of +0.52 V (in reference to a standard H2 electrode), and acetate (10 mM) was provided as the electron donor.
A Power Lab 4SP unit connected to a Power Macintosh computer collected current and voltage measurements directly from potentiostat outputs every 10 seconds. The data were logged with Chart 4.0 software (ADInstruments, Mountain View, CA) as levels of current (mA) production over time.
Geobacteraceae-dominated environment 1: U(VI)-contaminated aquifer.
A small-scale in situ bioremediation experiment was conducted on the grounds of a former uranium ore-processing facility in Rifle, Colorado (5), during the months of August and September 2004. This site, designated the Old Rifle site, is part of the Uranium Mill Tailings Remedial Action program of the U.S. Department of Energy. Using previously described procedures (5), acetate (100 mM) and a bromide tracer (10 mM) were injected into the subsurface to provide <10 mM acetate in the groundwater over the course of 28 days via an injection gallery composed of five injection wells which was installed up-gradient from a series of four monitoring wells (each 1.83 m down-gradient of the injection gallery/previous monitoring well). The test plot was adjacent to the 2002 Rifle experiment plot and was one-fourth the size and similar in design (5).
Geobacteraceae-dominated environment 2: site exposed to high concentrations of CMA.
The second site consisted of a highway runoff recharge pool located adjacent to State Route 25 in Plymouth, Massachusetts. The pool was constructed to collect runoff generated by State Route 25, which opened in August 1987 (14). The Massachusetts Department of Environmental Protection enacted restrictions on this area requiring the use of nonchloride deicing agents along a 1,900-m section of highway impacting nearby cranberry bogs. Since the road opened, the primary road deicing agent used on this stretch of highway has been calcium magnesium acetate (CMA). The unconfined aquifer underlying the study site is part of the Wareham Outwash Plain consisting of fine to coarse-grained sand. The concentration of acetate in the groundwater varies between 0 and 5 mM (36).
Geobacteraceae-dominated environment 3: Fe(III) reduction zone of petroleum-contaminated aquifer.
Aquifer sediments at the USGS Groundwater Toxics Site in Bemidji, Minnesota, have been contaminated with crude oil for 18 years as a result of a break in an oil pipeline (7, 16, 24). Fe(III) reduction is an important terminal electron-accepting process in portions of the aquifer (3, 29, 40). Sediments were collected from the Fe(III) reduction zone of the contaminant plume in 2004 by split-spoon sample cores and transported immediately to the laboratory as previously described (23).
Geobacteraceae-dominated environment 4: fuel cell deployed in freshwater pond.
Sediment fuel cells were deployed at the University of Massachusetts Field Station in Nantucket, Massachusetts, as previously described (38, 43). A graphite electrode with a surface area of 189.58 in2 was submerged in the sediment and connected via a 560-
resistor to a cathode of the same size in the overlying aerobic water. Two control electrodes were not connected to a cathode and did not harvest electricity. Current was produced by the sediment fuel cells for 16.3 weeks, and the voltage was logged hourly by a Campbell scientific logger with cellular phone access for remote monitoring. Electrodes were deployed in May and harvested in October, and current generated by the freshwater sediment fuel cell ranged from 0.05 to 1.45 mA (25).
Extraction of mRNA from pure cultures.
Solutions for RNA extraction were made with diethyl pyrocarbonate-treated water (Ambion). As previously described (23), chemostat cultures (200 ml) at a steady state were transferred to prechilled 50-ml conical tubes and centrifuged at 4,000 rpm for 15 min at 4°C. The supernatant was discarded, and pellets were flash frozen in an ethanol-dry ice bath and stored at 80°C. Prior to RNA extraction, the pellets were resuspended in 1.5 ml TPE buffer (100 mM Tris-HCl, 100 mM KH2PO4, 10 mM EDTA; pH 8.0) and divided into aliquots in eight separate 2-ml screw-cap tubes.
Cells were harvested from current-harvesting electrodes at six different rates of current, i.e., 0.5 mA, 0.75 mA, 1.0 mA, 1.2 mA, 1.4 mA, and 2.0 mA. Once the desired current was obtained, electrodes were removed from the anodic chamber, rinsed with RNA Protect (QIAGEN), and vigorously scraped with a sterile razor blade into 100 ml of RNA Protect, producing a graphite slurry. Graphite suspensions were then transferred to prechilled 50-ml conical tubes and centrifuged at 4,000 rpm for 15 min at 4°C. The supernatant was discarded, and pellets were flash frozen in an ethanol-dry ice bath and stored at 80°C. Prior to RNA extraction, each pellet was resuspended in 3 ml TPE buffer and divided into aliquots in eight separate 2-ml screw-cap tubes. RNAs were extracted from chemostat and electrode pellets as previously described (23).
Extraction of mRNA from environmental samples.
Groundwater samples were collected from the same well (well M16) at the U(VI)-contaminated site in Rifle, Colorado, every other day over the course of 27 days in the months of August and September 2004. At the site exposed to CMA in Plymouth, Massachusetts, groundwater samples were collected from well AY in June and October 2004. Groundwater from both sites was concentrated by impact filtration on 293-mm-diameter Supor membrane disk filters (Pall Life Sciences). The filters were placed into Whirl-Pack bags, flash frozen in a dry ice-ethanol bath, and shipped back to the laboratory, where they were stored at 80°C. RNA could only be extracted from one-half of a filter at a time. The frozen filter halves were first crushed into a fine powder, dispensed into eight different 2-ml screw-cap tubes, and suspended in 800 µl of TPE buffer. RNA was extracted from this filter suspension by a previously described acetone precipitation protocol (23).
In order to extract RNA from the surfaces of current-harvesting anodes, the electrodes were first removed from the sediment, placed into Whirl-Pack bags, and flash frozen in liquid nitrogen. Each frozen graphite electrode was then fragmented into smaller pieces and allocated into 15-ml conical tubes. Four milliliters of TE-sucrose buffer (10 mM Tris-HCl, 1 mM EDTA, 6.7% sucrose, pH 8.0), 60 µl 10% sodium dodecyl sulfate, 20 µl lysozyme (50 mg/ml), and 6 µl proteinase K (20 mg/ml) were added to each tube. Samples were then incubated at 37°C for 10 min and centrifuged at 4,000 rpm for 20 min. The supernatants were allocated into 2-ml screw-cap tubes, and a hot acidic phenol-chloroform extraction followed by isopropanol precipitation at 30°C was done as previously described (23).
RNA extracted from all pure culture and environmental RNA samples was further purified with an RNeasy RNA cleanup kit (QIAGEN) and treated with DNA-free DNase (Ambion) according to the manufacturer's instructions. All samples had A260/A280 ratios of 1.8 to 2.0, indicating that they were of high purity (6). In order to confirm the fact that PCR products generated from cDNA templates did not result from the amplification of contaminating DNA, all PCR analyses included negative controls with RNA that had not been subjected to reverse transcription-PCR (RT-PCR).
Testing and design of primers.
Degenerate primers targeting Geobacteraceae gltA genes were designed from nucleotide sequences from the G. sulfurreducens (32), G. metallireducens, Desulfuromonas palmitatis, Pelobacter carbinolicus, and Pelobacter propionicus genomes. Preliminary sequence data from the G. metallireducens, D. palmitatis, P. carbinolicus, and P. propionicus genomes were obtained from the DOE Joint Genome Institute website (www.jgi.doe.gov). These nucleotide sequences were initially aligned in Clustal X (44) and imported into the Genetic Computer Group sequence editor (Wisconsin Package, version 1.0; Madison, WI). This alignment was then examined, and conserved regions were targeted for primer design.
The following primer set was used to amplify a 750-bp fragment of the gltA gene from genomic DNA extracted from pure cultures of P. massiliensis and G. psychrophilus, Plymouth groundwater, Bemidji sediment, and the current-harvesting anode: CS100F (5' TCARTSTATTGGCGGTGC 3') and CS850R (5' CGGAGTCTTSAIGCAGAACTC 3'). The gltA gene was amplified from DNA extracted from Rifle groundwater with CS18F (5' CTCGCGACATCC GGAGTCT 3') and CS821R (5' TGTCCGGCGTTCAGGGTAT 3').
Primers were designed for TaqMan PCR analysis according to the manufacturer's instructions (amplicon size, 100 to 250 bp), and representative products from each of these primer sets were verified by sequencing. For pure-culture quantitative RT-PCR analyses of G. sulfurreducens, the following TaqMan primers that targeted G. sulfurreducens gltA, rpoD, recA, and proC transcripts were designed from the G. sulfurreducens genome sequence (32): CS418F (5' TGC TCTCCGTCGGTATCCTT 3') and CS552R (5' AGATGAATGCGGCGATGAC 3') amplified the citrate synthase gene, and the constitutively expressed housekeeping genes recA, proC, and rpoD were amplified with recA660F (5' GTGAAGGTGGTCAAGAACAAGGT 3') and recA737R (5' GGAAAT GCCCTCACCGTAGTAA 3'), proC2F (5' CATGCTGAAGGGAAGCACTCT 3') and proC77R (5' GGCCAGCAGCCCTTTGAT 3'), and rpoD1132F (5' TCATGAAGGCGGTGGACAA 3') and rpoD1210R (5' GCCTGTCGA ATCCACCAAGT 3'), respectively.
Before TaqMan primers could be designed to quantify environmental RNA transcripts, it was necessary to construct cDNA libraries using primer sets that targeted a larger region of the gene (400 to 800 bp). The gltA gene was amplified from cDNA generated from Plymouth groundwater, Bemidji sediment, and the current-harvesting anode with CS100F/CS850R, and the primer set CS18F/CS821R was used to amplify cDNA from Rifle groundwater. The proC gene was also amplified from cDNA generated from RNAs extracted from Rifle groundwater with proC75f (5' ATWGGIGGIGGIAATATGGC 3') and proC471r (5' TCCCCACCAGGTCGAACA 3'). The dominant sequences within each cDNA library constructed from these sites were then used to design specific Geobacteraceae primers suitable for quantitative RT-PCR analyses.
TaqMan primers specific for Geobacteraceae gltA genes in Rifle and Plymouth groundwater were amplified with CS375F (5' AACAAGATGACCGCCTGG G 3') and CS598R (5' TCGTGGTCGGAGTGGAGAAT 3'). Primers specific for Geobacteraceae gltA mRNA on the freshwater electrode included nanfwcs352F (5' ACCACCTACGAGGATGCCTG 3') and nanfwcs452R (5' TCGGGAGCGATGATGTTTC 3'), and proC cDNA from Rifle groundwater was amplified with proC100F (5' AGATACGGCATACAGGTTAC 3') and proC250R (5' ACGCCCGCCATGATGGAGATG 3').
Quantification of gene expression by TaqMan PCR.
A DuraScript enhanced avian RT single-strand synthesis kit (Sigma) was used to generate cDNA from gltA, recA, proC, and rpoD transcripts as previously described (23). Once the appropriate cDNA fragments were generated by RT-PCR, TaqMan PCR amplification and detection were performed with the GeneAmp 5700 sequence detection system (PE Biosystems, Foster City, CA). Optimal TaqMan PCR conditions were determined by using the manufacturer's guidelines. Each PCR mixture consisted of a total volume of 50 µl and contained 3 to 9 µl of the appropriate primers (stock concentration, 5 µM), 5 µl 10x SYBR green PCR buffer (PE Biosystems, Foster City, CA), 6 µl MgCl2 solution (25 mM; PE Biosystems), 4 µl deoxynucleoside triphosphate mix (2.5 mM; PE Biosystems), 0.5 U AmpErase uracil-N-glycosylase (UNG; PE Biosystems), and 0.25 U AmpliTaq Gold (PE Biosystems). The thermal cycling parameters consisted of a UNG activation step at 50°C for 2 min, a denaturation step at 95°C for 10 min, and 40 cycles of 95°C for 15 s and 63°C for 1 min. Standard curves were constructed as previously described (23) and covered a range of
8 orders of magnitude. To verify amplification and correct amplicon sizes, aliquots from real-time PCR were examined in an ethidium bromide-stained 1% agarose gel.
Amplification of gltA and 16S rRNA RT-PCR products and construction of cDNA libraries.
The optimal amplification conditions for CS100F/850R and CS18F/821R were determined in a gradient thermal cycler (MJ Research Inc., Waltham, MA). The amplification parameters included an initial denaturation step of 95°C for 5 min; 20 cycles of 95°C for 1 min, 55°C (0.5°C per cycle) for 1.5 min, and 72°C for 1.5 min; another 20 cycles of 95°C for 1 min, 50°C for 1.5 min, and 72°C for 1.5 min; and a final extension at 72°C for 10 min.
In an attempt to reduce PCR bias and to provide the most reliable evaluation of microbial community dynamics, more than one bacterial primer set targeting different regions of 16S rRNA was used. The primer sets included 8F (17) with 519R (26) and 338F (2) with 907R (26). Once the appropriate 16S cDNA gene fragments were generated by RT-PCR, the following PCR parameters were used: an initial denaturation step at 95°C for 5 min, followed by 35 cycles of 95°C (45 s), 50°C (1 min), and 72°C (1 min) and a final extension step at 72°C for 7 min. To ensure sterility, the PCR mixtures were exposed to UV radiation for 8 min prior to the addition of cDNA template and Taq polymerase.
16S rRNA and gltA PCR products were purified with a gel extraction kit (QIAGEN), and clone libraries were constructed with a TOPO TA cloning kit, version M (Invitrogen, Carlsbad, CA), according to the manufacturer's instructions. One hundred plasmid inserts from each cDNA clone library were then sequenced with the M13F primer at the University of Massachusetts Sequencing Facility.
Phylogenetic analysis.
16S rRNA and gltA gene sequences were compared to the GenBank nucleotide and protein databases using the BLASTn and BLASTx algorithms (1). Nucleotide and amino acid sequences for each gene were initially aligned in Clustal X (44) and imported into the Genetic Computer Group sequence editor (Wisconsin Package, version 1.0; Madison, WI). These alignments were then imported into Clustal W (45) and Mview (12), where similarity and identity matrices were generated.
Aligned sequences were imported into PAUP 4.0b10 (42), where phylogenetic trees were inferred. Distances and the branching order were determined and compared using maximum parsimony and distance-based algorithms (HKY85 and Jukes-Cantor). Bootstrap values were obtained from 100 replicates.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
In order to evaluate whether transcript levels for gltA might be affected by changes in the rate of metabolism, G. sulfurreducens was grown in chemostats at dilution rates near the lowest (0.03 h1) and highest (0.07 h1) rates at which G. sulfurreducens can be maintained in such systems (18). Acetate served as the electron donor and was the growth-limiting factor in order to mimic conditions that might typically be found in subsurface environments. Fe(III) was the electron acceptor. As expected, the steady-state cell yields (0.166 ± 0.029 mg protein/ml [mean ± standard error; n = 3] for 0.03 h1 and 0.171 ± 0.014 mg protein/ml for 0.07 h1) and Fe(II) concentrations (47.03 ± 8.48 mM for 0.03 h1 and 45.92 ± 3.46 mM for 0.07 h1) were comparable, and steady-state acetate concentrations were below the detection limit of 5 µM.
Quantitative real-time PCR studies demonstrated that the number of gltA transcripts was two times higher in the culture that was grown at 0.07 h1 than in the 0.03 h1 culture (Fig. 1A). In contrast, the levels of transcripts for the constitutively expressed housekeeping genes proC, recA, and rpoD remained relatively constant at both growth rates. These results suggest that levels of transcripts for gltA are related to the rates of growth and metabolism of G. sulfurreducens, and presumably other Geobacteraceae species.
|
Phylogenetic analysis of 16S and gltA cDNA clone libraries.
The phylogeny of the metabolically active members of the microbial community, as determined from analyses of 16S rRNA transcripts, as well as the presence and phylogeny of gltA transcripts, was evaluated in four environments in which Geobacteraceae were expected to be important components of the microbial community. The addition of acetate to a U(VI)-contaminated groundwater site in Rifle, Colorado, resulted in a dramatic increase in the number of Geobacteraceae 16S rRNA transcripts expressed in the groundwater, i.e., only 6.5% of the 16S rRNA bacterial transcripts were most similar to Geobacteraceae sequences on day 0, compared to 85% of the transcripts on day 14 (Fig. 2). It took approximately 7 days for acetate to reach well M16; the acetate concentration was 3.88 mM on day 7. Exposure to acetate resulted in a significant increase in the proportion of Geobacteraceae in the bacterial community; on day 8, Geobacteraceae accounted for 33% of the community. As a control, RNA was extracted from groundwater collected from a well that had not been exposed to acetate at the beginning (day 0) and end (day 27) of the experiment. The increase in Geobacteraceae sequences over time that was associated with the acetate-amended wells was not observed in these nonamended control samples, i.e., Geobacteraceae accounted for only 12.3% and 18.9% of the bacterial population on days 0 and 27, respectively.
|
|
|
gltA expression patterns in the environment.
In order to evaluate gltA expression patterns in the environment, gltA primers that specifically targeted the Geobacteraceae citrate synthase genes from each site were designed for TaqMan RT-PCR analysis. There was a remarkable correspondence between acetate levels in the groundwater and levels of transcripts for gltA during the in situ uranium bioremediation field study at the Old Rifle site. As acetate concentrations rose, gltA transcript levels increased (Fig. 4). A significant increase in the number of gltA transcripts was not observed until day 8, when acetate reached the well. When acetate levels fell due to a rain event that diluted the acetate with recharge water and prevented the injection of acetate for several days, the level of gltA transcripts decreased, and then gltA transcript levels increased concurrent with a renewed increase in acetate over time. In contrast to the gltA gene, the number of transcripts for the constitutively expressed housekeeping gene proC remained relatively constant (Fig. 4).
|
Implications.
The results of this study suggest that the ability to estimate levels of expression of the citrate synthase gene from Geobacteraceae (gltA) may be a useful tool for assessing the activity of these organisms in sedimentary environments. The eukaryotic nature of the citrate synthase gene suggests that few other prokaryotes expected to be found in anoxic subsurface environments will have similar citrate synthase sequences, making it easier to design primers to specifically monitor the metabolism of Geobacteraceae. Pure culture studies suggest that levels of gltA transcripts are directly related to rates of metabolism and electron transfer to Fe(III) and electrodes. Furthermore, in situ levels of gltA transcripts in Geobacteraceae-dominated environments can be monitored, and the finding that transcript levels fluctuate in accordance with changes in acetate availability in subsurface environments suggests that gltA expression patterns are related to in situ metabolic activity. These results demonstrate that monitoring gene expression in subsurface environments can provide an indication not only of the metabolic state of Geobacteraceae (23), but also of their rates of metabolism. Such information is invaluable for making decisions about the concentrations of various nutrients and electron donors that should be added to best stimulate in situ bioremediation. This information will also aid in the development of in silico models that can predict the metabolic responses of microorganisms to various bioremediation or energy-harvesting strategies. Such studies have the potential to change the design of bioremediation and energy harvesting from an empirical practice to more of a science (28).
| ACKNOWLEDGMENTS |
|---|
This research was supported by the Office of Science (BER), U.S. Department of Energy, with funds from the Natural and Accelerated Bioremediation Research (NABIR) program (grant DE-FG02-97ER62475), and by the Genomes to Life Program (cooperative agreement no. DE-FC02-02ER63446).
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| J. Bacteriol. | Microbiol. Mol. Biol. Rev. | Eukaryot. Cell | All ASM Journals |
|---|