Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2005, p. 4822-4832, Vol. 71, No. 8
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.8.4822-4832.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Jeffrey J. Walker,2,3,
Stephen J. Mojzsis,1,2 and
Norman R. Pace2,3*
Department of Geological Sciences, University of Colorado, Boulder, Colorado 80309-0399,1 Department of Molecular, Cellular, and Developmental Biology, University of Colorado, Boulder, Colorado 80309-0347,3 Center for Astrobiology, University of Colorado, Boulder, Colorado 80309-03922
Received 8 December 2004/ Accepted 7 March 2005
| ABSTRACT |
|---|
|
|
|---|
92%. Clone libraries were
90% bacterial and
10% archaeal, and eucaryotic rRNA genes were not detected in the libraries. The most abundant sequences were representative of novel proteobacteria (
28%), planctomycetes (
17%), and actinobacteria (
14%). Sequences representative of cyanobacteria, long considered to dominate these communities, comprised <5% of clones. Approximately 10% of the sequences were most closely related to those of
-proteobacterial anoxygenic phototrophs. These results provide a framework for understanding the kinds of organisms that build contemporary stromatolites, their ecology, and their relevance to stromatolites preserved in the geological record. | INTRODUCTION |
|---|
|
|
|---|
The stromatolites of Hamelin Pool, shown in Fig. 1, are examples of living stromatolites actively forming as a result of microbial activity. Stromatolites rise up to 0.5 m from the floor of this shallow hypersaline pool and display a variety of morphologies, including columnar, club, spheroidal, domal, nodular, or irregular shape (31). The stromatolites grow by accretion and precipitation of material to the outer surface layer, where the most conspicuous microbiology also occurs, commonly as a thinly stratified microbial mat (43, 49). Stromatolite morphology is often attributed to different kinds of microbial mat communities. Mat communities have been classified mostly by conspicuous cyanobacterial morphotypes and whether mats have a smooth or irregular appearance. Stromatolite communities also are structured vertically. Studies of Bahamian stromatolites, for instance, find that lithified interiors beneath surface mat communities are modified by endolithic microbes (44).
|
The objective of the present study was to determine the molecular phylogenetic makeup of microbial communities that build Hamelin Pool stromatolites and perhaps to shed light on their relevance to Precambrian stromatolites. The composition of such microbial communities cannot be determined by methods that require cultivation because most organisms (>99%) from natural environments are not cultured by standard methods (40). Consequently, we determined the composition and diversity of communities by sequence analysis of small-subunit (SSU) rRNA genes amplified by PCR with universal primers (515F and 1391R) from genomic DNA extracted directly from stromatolite samples. The compositions of communities were compared by using phylogenetic statistics to assess the similarity of sequence collections in phylogenetic trees. With this information, we tested the following hypotheses: (i) the stromatolite communities are comprised primarily of cyanobacteria, (ii) endolithic communities of stromatolites are distinct from surface communities, and (iii) different stromatolite morphotypes are associated with different microbial communities. These analyses collectively provide new perspective on the structure and composition of the Hamelin Pool stromatolite communities.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Sample description.
Samples of two stromatolites were collected from the Hamelin Pool in the area of Hamelin Station in Shark Bay during the late afternoon of July 9, 2003. Samples of the upper
3 cm of stromatolites were collected from the intertidal zone, preserved in 70% ethanol in sterile 50-ml plastic tubes, and stored at 4°C.
Sample HPDOM (for Hamelin Pool Domal) had a smooth, domal morphology (Fig. 1) with a thin (
1-mm) indurated surface colored green, orange, and black. The shape of the HPDOM stromatolite was most similar to the "smooth surface mat" described by Reid et al. (43). It was collected when the stromatolite was directly exposed to air at low-to-mid tide. A chip of HPDOM was used to make a petrographic thin section for observation by optical microscopy. Figures 2A and B show representative images of the crust and interior of HPDOM, respectively. The beige-white porous interior consisted of submillimeter-sized carbonate sand, foraminifera, bivalve shells, shell fragments, and microcrystalline carbonate (micrite). Aragonite needles often coated or partially filled the pore spaces of HPDOM, and some detrital quartz grains and dolomite crystals (grown in situ) were observed. The microfabric of HPDOM most closely resembled the "unlaminated calcarenite" described by Reid et al. (43). HPDOM was subsampled by manual separation into a surface sample ca. 1 to 2 mm thick (HPDOM-S) and an interior sample taken from a vertical interval ca. 10 to 15 mm below the surface (HPDOM-I).
|
5 mm of HPIRR was analyzed in bulk.
Laser scanning confocal microscopy (LSCM).
We imaged fluorescence associated with oxygenic phototrophs with a Leica TCS SP2 AOBS laser scanning confocal microscope. We imaged fluorescent emissions at 680 to 690 nm, which is due to chlorophyll molecules in reaction centers of photosystem II (PS II) of oxygenic phototrophs in vivo (6). Wavelengths that excite PS II fluorescence depend on associated accessory antenna pigments, which vary in different groups of oxygenic phototrophs. Thus, excitation profiles of PS II fluorescence are diagnostic of different kinds of oxygenic phototrophs (6). To determine excitation profiles, we measured the relative fluorescence intensity at 680 to 690 nm when samples were excited by wavelengths of 458, 476, 488, 496, 514, 543, and 633 nm. We also detected fluorescence of the antenna pigment phycoerythrin directly at 570 nm with an excitation wavelength of 543 nm (6).
Genomic DNA extraction.
We extracted duplicate community genomic DNA samples from each stromatolite sample (HPDOM-S, HPSOM-I, and HPIRR) and two negative controls with a bead-beating extraction protocol (16). A sample of each stromatolite (
5 g) was crushed and homogenized in a flame-sterilized steel ore-crusher mortar and pestle (Fisher Scientific), and 0.25 to 0.35 g of crushed stromatolite was used in each DNA extraction. To prevent contamination, procedures were performed in a UV sterilized AirClean 600 PCR workstation (AirClean Systems). Community DNA was resuspended in 40 µl of TE (10 mM Tris, 1 mM EDTA). DNA concentrations were measured with the PicoGreen dsDNA quantitation kit (Molecular Probes) and a fluorometer. Average yields of extracted DNA were
5 µg/g of crushed rock.
PCR amplification of rRNA genes.
SSU rRNA genes were amplified by PCR with community DNA as a template and the universal primers 515F (5'-GTGCCAGCMGCCGCGGTAA-3') and 1391R (5'-GACGGGCGGTGWGTRCA-3') (30). Each 50-µl PCR contained
50 ng of template DNA, 1.25 U of HotMaster Taq DNA Polymerase (Eppendorf), 1x HotMaster buffer (2.5 mM Mg2+), 1 mg of bovine serum albumin/ml, 50 µM concentrations of each deoxynucleoside triphosphate, and 0.2 µM concentrations of each primer. PCR conditions were 94°C for 2 min, followed by 30 cycles of 94°C for 15 s, 52°C for 15 s, and 65°C for 60 s, with a final 20- min extension at 65°C. PCR controls with no template or extraction control as a template were negative.
Clone library construction.
We constructed one clone library, 96 randomly selected rRNA clones, for each community DNA. PCR-amplified products from four independent PCRs were pooled to reduce the chances of PCR artifacts (28), purified by agarose gel electrophoresis with the Montage DNA gel extraction kit (Millipore), and cloned into the pCR4.0 vector by using the TOPO TA cloning kit and TOP10 chemically competent cells (Invitrogen). For each library, 96 colonies were picked, grown for 16 h, and processed into a quick plasmid preparation (1:4 dilution of overnight culture in Tris-EDTA, incubated at 85°C for 10 min and centrifuged).
Restriction fragment length polymorphism analysis and sequencing.
rRNA gene clones were amplified by PCR from plasmid preparations with the vector primers T3 and T7. Amplified rRNA clones were initially screened by restriction fragment length polymorphism analysis to assess the diversity of clones prior to sequencing (25). The results showed that the libraries were sufficiently diverse to sequence all clones. PCR-amplified rRNA clones were prepared for DNA sequencing with the ExoSAP-it PCR cleanup kit (USB) and sequenced by using DYEnamic ET terminator cycle sequencing kit reaction mixtures on a MegaBACE 1000 DNA sequencer (Amersham Biosciences) with T3 and T7 vector primers. A total of 531 clones were analyzed in the present study.
Phylogenetic analysis and chimera detection.
Raw sequence data were analyzed and assembled into contiguous sequences with the PHRED and PHRAP software packages (17, 18). Sequences were compared to a current database of rRNA gene sequences from GenBank with BLAST (4). These steps were performed with the aid of XplorSeq, an in-house software package developed to automate and manage sequence analysis (D. N. Frank, unpublished software). Complete rRNA gene sequences were exported from XplorSeq and aligned to a database of >10,000 known SSU rRNA sequences (25) with the Arb software package (33). We used the method described by Hugenholtz et al. to identify chimeras (25). One chimera was identified and excluded from further analysis. Phylogenetic analysis was performed with the Arb and PAUP* 4.0 software packages (48) as previously described (23).
Richness and coverage estimators.
To estimate richness and coverage, sequences were collected into operational taxonomic units (OTUs) based on sequence identity at different fixed levels from 95 to 100% in 1% increments with a computer script written in the Python programming language that analyzes distance matrices exported from Arb. OTU richness and coverage estimators were calculated with the software program EstimateS version 6b1 (13) with default settings and 100 random sample repetitions.
Comparative phylogenetic methods.
We compared collections of community rRNA gene sequences with the phylogenetic test (p-test), which calculates differences between sequence sets in phylogenetic trees (34). This method calculates the probability that sequences in two sets compared are derived from the same sequence population with randomly permuted trees. The p-test tests whether there are more distinct lineages in a tree than expected by chance. P values were calculated as the ratio of 1,000 randomly permuted trees with fewer parsimony changes than the original tree. We considered P values of
0.05 as significant to reject the hypothesis that two sequence sets were derived from the same sequence population. The p-tests were determined with the Phylonode software module (C. Lozupone and R. Knight, unpublished software). We also used Phylonode to automatically identify internal nodes with significant differences. Each internal node was selected and tested to determine whether removal of that node significantly changed the result from the original tree. We adjusted P values for multiple tests with the Bonferroni correction.
Nucleotide sequence accession numbers.
Sequences of 478 rRNA gene clones determined are available under GenBank accession numbers AY851765 to AY852242. Clone names designate community (HPDOMI, HPDOMS, and HPIRR), replicate library number (1 or 2), and the clone number (e.g., A01).
| RESULTS |
|---|
|
|
|---|
Autosimilarity and sample representativeness.
Microbial communities are highly complex, and we do not sample all unique sequences. Consequently, we tested whether the clone libraries were representative samples of rRNA genes amplified and cloned from community DNA by assessment of the autosimilarity (similarity among replicate samples randomly drawn from a population (11) of duplicate clone libraries from each stromatolite with the p-test (see Materials and Methods)). The results showed high autosimilarities with no significant differences in phylogenetic representation between duplicate libraries from the same stromatolite sample (P > 0.995). This indicates clone libraries were representative samples of at least the abundant phylogenetic diversity of sequences cloned from each community. The similarity of duplicate libraries is illustrated by the phylogenetic trees in Fig. 3, which show a relatively even distribution of sequences from duplicate libraries for each community. Based on these results, we combined duplicate libraries from each community and performed all additional analyses on combined sequence sets unless stated otherwise.
|
|
Comparison to known microbial diversity.
The relatedness of Hamelin Pool sequences to known rRNA gene sequences was determined by comparison to the >200,000 rRNA gene sequences currently in GenBank with BLAST to find the closest match for each Hamelin Pool sequence. Figure 4 shows the distribution of relatedness between Hamelin Pool sequences and their closest match in GenBank. Hamelin Pool sequences were related to known sequences by an average of only 92% sequence identity, which represents microbial diversity substantially different from what is currently known. Sequence relatedness ranged from 85 to 100% identity and varied by phylogenetic group. Table 2 shows the average relatedness of Hamelin Pool sequences to known sequences for each phylogenetic division represented. Only 3% of Hamelin Pool sequences were related to known sequences at
97% sequence identity, which suggests 97% of sequences represent new species (47).
|
|
Hamelin Pool stromatolite communities were comprised of rRNA genes representative of ca. 90% Bacteria and 10% Archaea. We detected no eucaryotic rRNA genes by PCR with universal primers in the present study. The stromatolite communities were complex, comprised of organisms with rRNA genes representative of 19 of the 52 identified divisions of Bacteria (42), as well as representatives of euryarchaeota and crenarchaeota. The most abundant rRNA gene sequences in the combined data set were representative of proteobacteria (
29%), planctomycetes (
17%), and actinobacteria (
15%). In contrast, cyanobacterial sequences comprised <5% of the stromatolite communities. No sequences indicative of novel phylogenetic divisions were detected among the >500 clones analyzed.
Surface compared to endolithic communities.
In order to test the hypothesis that endolithic communities of stromatolites are distinct from surface communities, we compared rRNA gene sequences from HPDOM-I to those of HPDOM-S. The p-test showed a significant difference between the surface and interior communities (P < 0.001). Both communities were comprised of diverse sequences representative of Bacteria and Archaea. Sequences from the interior community represented 19 bacterial divisions. Sequences from the surface community represented 12 of the same bacterial divisions and one additional division, TM6.
With the p-test, we tested (i) whether there was a significant difference between the communities for each lineage in the tree and (ii) which lineages reduced the significant difference between the communities when removed from the tree. Lineages with significant differences are shown in Table 2. Five of twelve bacterial divisions with sequences from both communities showed significant differences (P < 0.05), as did the distribution of euryarchaeota and
-,
-, and
-proteobacteria. Except for the firmicutes, lineages with significant differences also contained most sequences. Removal of lineages did not reduce the significant difference between the communities, which suggests the difference between the communities is distributed throughout the phylogenetic tree and not concentrated in any particular lineage.
Figure 5A shows the relative distribution of clones from the surface and interior communities for the most abundant phylogenetic lineages. There were no identical sequences and few closely related sequences shared between the surface and interior communities. Only ca. 20% of sequences from the two communities were related at 97% or greater sequence identity, which suggests about 80% of each community differed from the others at a minimum of the species level. The most obvious differences in sequence composition between layers were in relative abundances of
-proteobacterial, planctomycete, and cyanobacterial sequences, as shown in Fig. 5A.
|
98% identical. The HPDOM cyanobacterial sequences were most closely related to the filamentous Microcoleus spp., but by only
94% sequence identity, which indicates the most abundant cyanobacteria in HPDOM are relatively novel compared to known species. Sequences representative of
-proteobacteria also were more abundant in the surface community, and all were most closely related to cultured organisms known for anoxygenic photosynthesis, as shown in Fig. 6B. Planctomycete sequences were enriched in the interior community, and there were significant phylogenetic differences from exterior community, as shown in Table 2. Planctomycete sequences from both communities were substantially different from known sequences, with an average of
10% rRNA sequence divergence.
|
The communities were similar in their relatively low abundances of cyanobacterial sequences and high abundances of
-proteobacterial sequences, which were most closely related to known anoxygenic phototrophs (Fig. 6B). HPIRR had two distinct groups of cyanobacterial sequences, one most closely related to Microcoleus spp. sequences and the other to sequences of nonfilamentous Pleurocapsa spp. (Fig. 6A).
There were more differences than similarities in the composition of HPDOM and HPIRR at the bacterial phylogenetic division level, as shown in Table 2 and Fig. 5B. Only about 10% of sequences from the two communities were related at
97% sequence identity, and only about 15% of sequences were related at
95% sequence identity. The most obvious differences in composition of the two communities were in sequences representative of actinobacteria and euryarchaeota, which were more abundant in HPIRR, and planctomycetes, which were more abundant in HPDOM.
Microscopy.
In order to determine the distribution of oxygenic phototrophs in Hamelin Pool stromatolite communities, we imaged fluorescence indicative of pigments associated with oxygenic photosynthesis in thin sections by LSCM. We also analyzed spectral signatures of fluorescent cells by LSCM; these signatures are diagnostic of photosynthetic pigments associated with different groups of oxygenic phototrophs (see Materials and Methods). Fluorescence was induced only at 633 nm and not other wavelengths (458, 476, 488, 496, 514, and 543 nm), which is indicative of cyanobacteria in vivo and contraindicative of cells that contain chloroplasts, which should be excited at 458 nm (6). Figure 7 shows the spatial distribution and spectral signature of fluorescent cells in HPDOM. All observable fluorescence was located in a layer
0.2 mm thick at the surface of the samples. Individual cells with diameters of 2 to 4 µm were observed at the surface. The underlying band of diffuse fluorescence was due to scattered fluorescence emitted by cells embedded in the mineral matrix. The excitation observed in cells at 633 nm and corresponding emission at 685 nm are diagnostic of phycocyanin. A subset of fluorescent cells shown in Fig. 7G to I also exhibited fluorescence diagnostic of phycoerythrin (excitation, 543 nm; emission, 590 nm) (6).
|
| DISCUSSION |
|---|
|
|
|---|
We acknowledge, however, that there are potential experimental artifacts, which can introduce biases in how well clone libraries represent the actual rRNA gene community in a sample. These include variable PCR amplification that results from primer selectivity and differential extraction of genomic DNA from environmental samples (28, 51). Different species may also have variable numbers of rRNA gene copies (2). Therefore, the relative abundance of rRNA genes surveyed from a sample with PCR-based methods may not correspond directly to the relative abundance of organisms in the original sample.
Nonetheless, studies show a significant correspondence between abundances measured with clone libraries and assessed with other methods, such as fluorescent in situ hybridization (27), rRNA hybridization (35), and clone libraries made by reverse transcription of rRNA (38). General correspondence is also observed between libraries made with different PCR primer combinations (20, 46). Collectively, the results show that, when used carefully, these methods accurately identify the most abundant organisms present in the original sample. Although we cannot rule out systematic biases in the present study, the high autosimilarity of our duplicate libraries indicates there was little random bias in the methods.
Sample representativeness.
The statistical treatment of our results suggests that the clone libraries were representative of the population of sequences amplified and cloned from the stromatolite communities. Sample collection is a complex and essential consideration in general ecology (11), which attracts attention in molecular microbial ecology as efforts to understand microbial diversity increase and the cost of DNA sequencing decreases (7, 19, 26). Natural microbial communities tend to be highly complex, and the in situ diversity is seldom entirely sampled with molecular methods. Consequently, for comparison and estimation of community compositions it is essential to assess how well samples represent communities (11, 34). We assessed sample representativeness with two methods that test different measures of autosimilarity. This approach has been proposed as a more objective and efficient means to assess sample representativeness than other methods currently used, such as estimation of rarefaction or coverage (11, 29).
Although results from comparative phylogenetic methods and rarefaction analysis suggested our samples were representative, we propose that measurement of the autosimilarity of sequence data by phylogenetic methods is more robust and objective than measurement by OTU-based methods because the latter approach imposes a taxonomic species concept on phylogenetic data. This can lead to questionable assumptions about the taxonomic relationships of organisms based on rRNA sequence variation. For example, evidence suggests the relationship between taxonomic units and SSU rRNA gene sequence varies by phylogenetic division (19), but the details have not been articulated. Stackebrandt and Goebel (47) studied this relationship in Bacteria and concluded that strains with SSU rRNA genes that are <97% identical are most likely different species, but strains with SSU rRNA genes that are >97% identical could be the same species or different species.
Phylogenetic composition.
The microbial communities of the Hamelin Pool stromatolites are comprised of
90% Bacteria and
10% Archaea, similar to results from other studies of different kinds of microbial mats (45). Much of the novel diversity we found was associated with bacterial phylogenetic divisions that are considered highly sampled, such as the
-proteobacteria, where we found sequences <90% identical to known sequences. This shows there is undiscovered novel diversity even in lineages represented by many sequences. We also obtained sequences characteristic of five candidate divisions of Bacteria, divisions with no cultured representatives that previously had been encountered in other environmental surveys (24, 42). These sequences constitute new and independent evidence that these candidate divisions in fact represent unique bacterial lineages.
We found sequences representative of a relatively low abundance of archaea, but both euryarchaeotes and crenarchaeotes were represented. Little is known about the metabolic diversity of archaea, and it is therefore not possible to assess their potential roles in the stromatolite communities. Euryarchaeal clones were most closely related to halophilic archaea, which is consistent with the hypersaline environment of the Hamelin Pool. Crenarchaeal rRNA sequences were not specifically associated with those of any cultured organisms. Absence of eucaryotic sequences suggests that if they were present in these stromatolites, they were relatively scarce, in contrast to previous observations (5). Our results corroborate the findings of Reid et al. (44), who also found that diatoms and other eucaryotes are scarce in Bahamian stromatolites.
Primary productivity, biomass, and community structure.
The emergent properties of stromatolite communities are governed by how community members contribute to the biomass, energy flow, and mineralization of the community. Our results provide a snapshot of the relative proportions of phylogenetic types in the stromatolite communities. The standing crop biomass of a community is a function of the relative abundance and size of its members. Community members are sustained by a flow of energy, which primary producers ultimately generate from an external source of energy, such as light or chemical energy. Both the relative size and the turnover rate of primary producers influences the number of individuals needed to sustain a community. A distinction between the overall biomass of the stromatolite communities and that of its primary producers is potentially important to understand the processes that form stromatolites.
The low abundance of cyanobacterial rRNA genes in the Hamelin Pool communities suggests that other kinds of organisms comprise much of the biomass. Surveys of hypersaline microbial-mat communities (45) and Bahamian stromatolite communities (L. K. Baumgartner, N. R. Pace, and P. T. Visscher, unpublished data) also find low relative proportions of cyanobacterial rRNA genes in universal libraries. Biogeochemical studies of the Bahamian stromatolites show most oxygenic photosynthesis occurs in a submillimeter layer near the surface (50). This suggests that the distribution of oxygenic phototrophs we observed in Hamelin Pool communities is similar to that of Bahamian communities.
If the main source of primary productivity in Hamelin Pool communities is oxygenic photosynthesis by cyanobacteria, then these stromatolites are supported by a relatively small fraction of the community (
5% of rRNA genes). Such a scenario seems possible considering cyanobacteria can have very high turnover rates (39), and the cyanobacterial cells we observed microscopically had on average about twice the diameter of other cells. This observation implies cyanobacterial cells in the community have cell volumes about eight times larger than other organisms. In addition, there is no obvious external source of chemical energy in Hamelin Pool sufficient to drive new net primary production by metabolisms other than oxygenic photosynthesis.
Nonetheless, we find evidence that potential anoxygenic phototrophs are numerically abundant in the stromatolite community and thus may represent a considerable fraction of the standing biomass. Such organisms undoubtedly contribute to gross primary productivity in the community by photosynthesis, but anoxygenic phototrophs rely on a supply of electron donors such as sulfide, hydrogen, and reduced organics. Energy that typically produces such electron donors in comparable microbial systems originates from oxygenic phototrophs, cyanobacteria, and is transformed by a complex web of metabolisms (14). For example, sulfate reduction occurs in Bahamian stromatolites in an anoxic layer several millimeters below the photosynthetic zone (50). Such a sulfidic layer could support anoxygenic photosynthesis in stromatolite communities. Most other Hamelin Pool sequences were only distantly related to those of cultivated organisms, and so the metabolisms of the organisms they represent are not predictable.
Geological implications.
The significant differences we found between the HPDOM and HPIRR communities show that stromatolite morphology may be attributable to different microbial communities, although more comparisons are required to confirm this hypothesis. Ancient stromatolites have many morphologies (reviewed in reference 22), which may indicate that different microbial communities formed the basis for each unique structure. Variations in stromatolite morphology could result, for instance, from the nature of microbial extracellular materials, which might trap sediments differently or attract different community members.
Although cyanobacteria do not dominate Hamelin Pool stromatolite communities numerically, they likely play the main role in primary production. Our results suggest that novel and diverse microorganisms form much of the structure and fabric of the communities and, therefore, must contribute substantially to stromatolite formation. It is possible that microbial communities fueled by metabolisms other than oxygenic photosynthesis formed some ancient stromatolites. For example, stromatolites in Yellowstone hotsprings are built by communities in which the primary producers are considered anoxygenic photosynthetic chloroflexi fueled by electron donors derived from hydrothermal fluids (15). Therefore, we conclude that ancient stromatolites cannot be taken alone as evidence for oxygenic photosynthesis.
| ACKNOWLEDGMENTS |
|---|
This study was supported by NSF grant DEB-0085490 to N.R.P. and the NASA Cooperative Agreement with the University of Colorado Center for Astrobiology to N.R.P. and S.J.M. D.P. thanks the Fonds de Recherche sur la Nature et les Technologies du Québec and the NAI for graduate support.
| FOOTNOTES |
|---|
D.P. and J.J.W. contributed equally to this study. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| J. Bacteriol. | Microbiol. Mol. Biol. Rev. | Eukaryot. Cell | All ASM Journals |
|---|