High Temporal but Low Spatial Heterogeneity of Bacterioplankton in the Chesapeake Bay

Compared to freshwater and the open ocean, less is known about bacterioplankton community structure and spatiotemporal dynamics in estuaries, particularly those with long residence times. The Chesapeake Bay is the largest estuary in the United States, but despite its ecological and economic significance, little is known about its microbial community composition. A rapid screening approach, ITS (internal transcribed spacer)-LH (length heterogeneity)-PCR, was used to screen six rRNA operon (16S rRNA-ITS-23S rRNA) clone libraries constructed from bacterioplankton collected in three distinct regions of the Chesapeake Bay over two seasons. The natural length variation of the 16S-23S rRNA gene ITS region, as well as the presence and location of tRNA-alanine coding regions within the ITS, was determined for 576 clones. Clones representing unique ITS-LH-PCR sizes were sequenced and identified. Dramatic shifts in bacterial composition (changes within subgroups or clades) were observed for the Alphaproteobacteria (Roseobacter clade, SAR11), Cyanobacteria (Synechococcus), and Actinobacteria, suggesting strong seasonal variation within these taxonomic groups. Despite large gradients in salinity and phytoplankton parameters, a remarkably homogeneous bacterioplankton community was observed in the bay in each season. Stronger seasonal, rather than spatial, variation of the bacterioplankton population was also supported by denaturing gradient gel electrophoresis and LH-PCR analyses, indicating that environmental parameters with stronger seasonal, rather than regional, dynamics, such as temperature, might determine bacterioplankton community composition in the Chesapeake Bay.     

Low intraspecific diversity in a polynucleobacter subcluster population numerically dominating bacterioplankton of a freshwater pond

Cultivation-dependent and -independent methods were combined to investigate the microdiversity of a Polynucleobacter subcluster population (Betaproteobacteria) numerically dominating the bacterioplankton of a small, humic freshwater pond. Complete coverage of the population by cultivation allowed the analysis of microdiversity beyond the phylogenetic resolution of ribosomal markers. Fluorescent in situ hybridization with two probes specific for the narrow subcluster C (PnecC bacteria) of the Polynucleobacter cluster revealed that this population contributed up to 60% to the total number of bacterioplankton cells. Microdiversity was investigated for a date at which the highest relative numbers of PnecC were observed. A clone library of fragments of the ribosomal operon (16S rRNA genes, complete 16S-23S internal transcribed spacer 1 [ITS1], partial 23S rRNA genes) amplified with universal bacterial primers was constructed. The library was stepwise screened for fragments from PnecC bacteria and for different ITS genotypes of PnecC bacteria. The isolated PnecC strains were characterized by sequencing of the 16S rRNA genes and the ITS1. Both the clone library and the established culture collection contained only the same three ITS genotypes, and one of them contributed 46% to the entire number of clones. Genomic fingerprinting of the isolates with several methods always resulted in the detection of only one fingerprint per ITS genotype. We conclude that a Polynucleobacter population with an extremely low intraspecific diversity and an uneven structure numerically dominated the bacterioplankton community in the investigated habitat. This low intraspecific diversity is in strong contrast to the high intraspecific diversities found in marine bacterial populations.     

Low Intraspecific Diversity in a Polynucleobacter Subcluster Population Numerically Dominating Bacterioplankton of a Freshwater Pond

Cultivation-dependent and -independent methods were combined to investigate the microdiversity of a Polynucleobacter subcluster population (Betaproteobacteria) numerically dominating the bacterioplankton of a small, humic freshwater pond. Complete coverage of the population by cultivation allowed the analysis of microdiversity beyond the phylogenetic resolution of ribosomal markers. Fluorescent in situ hybridization with two probes specific for the narrow subcluster C (PnecC bacteria) of the Polynucleobacter cluster revealed that this population contributed up to 60% to the total number of bacterioplankton cells. Microdiversity was investigated for a date at which the highest relative numbers of PnecC were observed. A clone library of fragments of the ribosomal operon (16S rRNA genes, complete 16S-23S internal transcribed spacer 1 [ITS1], partial 23S rRNA genes) amplified with universal bacterial primers was constructed. The library was stepwise screened for fragments from PnecC bacteria and for different ITS genotypes of PnecC bacteria. The isolated PnecC strains were characterized by sequencing of the 16S rRNA genes and the ITS1. Both the clone library and the established culture collection contained only the same three ITS genotypes, and one of them contributed 46% to the entire number of clones. Genomic fingerprinting of the isolates with several methods always resulted in the detection of only one fingerprint per ITS genotype. We conclude that a Polynucleobacter population with an extremely low intraspecific diversity and an uneven structure numerically dominated the bacterioplankton community in the investigated habitat. This low intraspecific diversity is in strong contrast to the high intraspecific diversities found in marine bacterial populations.     

Pronounced summer to winter differences and higher wintertime richness in coastal Antarctic marine bacterioplankton

Marine bacterioplankton studies over the annual cycle in polar systems are limited due to logistic constraints in site access and support. Here, we conducted a comparative study of marine bacterioplankton sampled at several time points over the annual cycle (12 occasions each) at sub-Antarctic Kerguelen Islands (KI) and Antarctic Peninsula (AP) coastal sites in order to establish a better understanding of the extent and nature of variation in diversity and community structure at these different latitudes (49-64S). Molecular methods targeting the 16S rRNA gene (DGGE, CE-SSCP and tag pyrosequencing) suggest a strong seasonal pattern with higher richness in winter and a clear influence of phytoplankton bloom events on bacterioplankton community structure and diversity in both locations. The distribution of sequence tags within Gammaproteobacteria, Alphaproteobacteria and Bacteriodetes differed between the two regions. At both sites, several abundant Rhodobacteraceae, uncultivated Gammaproteobacteria and Bacteriodetes-associated tags displayed intense seasonal variation often with similar trends at both sites. This enhanced understanding of variability in dominant groups of bacterioplankton over the annual cycle contributes to an expanding baseline to understand climate change impacts in the coastal zone of polar oceans and provides a foundation for comparison with open ocean polar systems.     

Disentangling seasonal bacterioplankton population dynamics by high‐frequency sampling

Multiyear comparisons of bacterioplankton succession reveal that environmental conditions drive community shifts with repeatable patterns between years. However, corresponding insight into bacterioplankton dynamics at a temporal resolution relevant for detailed examination of variation and characteristics of specific populations within years is essentially lacking. During 1 year, we collected 46 samples in the Baltic Sea for assessing bacterial community composition by 16S rRNA gene pyrosequencing (nearly twice weekly during productive season). Beta‐diversity analysis showed distinct clustering of samples, attributable to seemingly synchronous temporal transitions among populations (populations defined by 97% 16S rRNA gene sequence identity). A wide spectrum of bacterioplankton dynamics was evident, where divergent temporal patterns resulted both from pronounced differences in relative abundance and presence/absence of populations. Rates of change in relative abundance calculated for individual populations ranged from 0.23 to 1.79 day^?1. Populations that were persistently dominant, transiently abundant or generally rare were found in several major bacterial groups, implying evolution has favoured a similar variety of life strategies within these groups. These findings suggest that high temporal resolution sampling allows constraining the timescales and frequencies at which distinct populations transition between being abundant or rare, thus potentially providing clues about physical, chemical or biological forcing on bacterioplankton community structure.     

Nearly a decade‐long repeatable seasonal diversity patterns of bacterioplankton communities in the eutrophic Lake Donghu (Wuhan,China)

Uncovering which environmental factors govern community diversity patterns and how ecological processes drive community turnover are key questions related to understand the community assembly. However, the ecological mechanisms regulating long‐term variations of bacterioplankton communities in lake ecosystems remain poorly understood. Here we present nearly a decade‐long study of bacterioplankton communities from the eutrophic Lake Donghu (Wuhan, China) using 16S rRNA gene amplicon sequencing with MiSeq platform. We found strong repeatable seasonal diversity patterns in terms of both common (detected in more than 50% samples) and dominant (relative abundance >1%) bacterial taxa turnover. Moreover, community composition tracked the seasonal temperature gradient, indicating that temperature is a key environmental factor controlling observed diversity patterns. Total phosphorus also contributed significantly to the seasonal shifts in bacterioplankton composition. However, any spatial pattern of bacterioplankton communities across the main lake areas within season was overwhelmed by their temporal variabilities. Phylogenetic analysis further indicated that 75%–82% of community turnover was governed by homogeneous selection due to consistent environmental conditions within seasons, suggesting that the microbial communities in Lake Donghu are mainly controlled by niche‐based processes. Therefore, dominant niches available within seasons might be occupied by similar combinations of bacterial taxa with modest dispersal rates throughout different lake areas.     

Lake bacterioplankton dynamics over diurnal timescales

相似文献   

Relationship between bacterioplankton richness, respiration, and production in the Southern North Sea

We investigated the relationship between bacterioplankton production (BP), respiration (BR), and community composition measured by terminal restriction fragment length polymorphism in the southern North Sea over a seasonal cycle. Major changes in bacterioplankton richness were apparent from April to December. While cell-specific BP decreased highly significantly with increasing bacterioplankton richness, cell-specific BR was found to be variable along the richness gradient, suggesting that bacterioplankton respiration is rather independent from shifts in the bacterial community composition. As a consequence, the bacterial growth efficiency [BGE = BP/(BP + BR)] was negatively related to bacterioplankton richness, explaining approximately 43% of the variation in BGE. Our results indicate that despite the observed shifts in the community composition, the main function of the bacterioplankton, the remineralization of dissolved organic carbon to CO(2), is rather stable.     

Relationship between Bacterioplankton Richness, Respiration, and Production in the Southern North Sea

We investigated the relationship between bacterioplankton production (BP), respiration (BR), and community composition measured by terminal restriction fragment length polymorphism in the southern North Sea over a seasonal cycle. Major changes in bacterioplankton richness were apparent from April to December. While cell-specific BP decreased highly significantly with increasing bacterioplankton richness, cell-specific BR was found to be variable along the richness gradient, suggesting that bacterioplankton respiration is rather independent from shifts in the bacterial community composition. As a consequence, the bacterial growth efficiency [BGE = BP/(BP + BR)] was negatively related to bacterioplankton richness, explaining ~43% of the variation in BGE. Our results indicate that despite the observed shifts in the community composition, the main function of the bacterioplankton, the remineralization of dissolved organic carbon to CO₂, is rather stable.     

Contrasting nifD and ribosomal gene relationships among Mesorhizobium from Lotus oroboides in northern Mexico

PCR screens for length variation in a 5' portion of 23S ribosomal RNA and in the 3' end of the 16S rRNA-23S rRNA internal transcribed spacer (ITS) region indicated that nodule bacteria from a Mexican population of Lotus oroboides were diverse on a local scale. Three 23S rRNA length variants and five ITS length variants were detected among the 22 isolates. Sequencing of nearly full-length 16S rRNA genes in three isolates indicated that they fell into the genus Mesorhizobium, but comprised two distinct groups. Two isolates were closely related to M. loti LMG 6125T, while the other isolate clustered with an assemblage of Mesorhizobium taxa that included M. amorphae, M. plurifarium and M. huakuii. However, a phylogenetic tree based on 715 bp of the nitrogenase alpha-subunit (nifD) gene was significantly discordant with the relationships inferred from rRNA sequences. Two isolates that were nearly identical for 16S rRNA had nifD genes that varied at 2% of sites, and one of these nifD sequences was identical to that of another isolate with a strongly divergent 16S rRNA gene. A plasmid screen followed by Southern hybridization indicated that only one of these strains harbored a plasmid-borne nifD gene. These results imply that gene transfer events have altered the distribution of nifD sequences among lineages within this natural population of Mesorhizobium strains.     

Temporal patterns of microbial community structure in the Mid-Atlantic Bight

Although open ocean time-series sites have been areas of microbial research for years, relatively little is known about the population dynamics of bacterioplankton communities in the coastal ocean on kilometer spatial and seasonal temporal scales. To gain a better understanding of microbial community variability, monthly samples of bacterial biomass were collected in 1995-1996 along a 34-km transect near the Long-Term Ecosystem Observatory (LEO-15) off the New Jersey coast. Surface and bottom sampling was performed at seven stations along a transect line with depths ranging from 1 to 35 m (n=178). Microbial populations were fingerprinted using ribosomal 16S rRNA genes and terminal restriction fragment length polymorphism analysis. Results from cluster analysis revealed distinct temporal patterns among the bacterioplankton communities in the Mid-Atlantic Bight rather than grouping by sample location or depth. Principal components analysis models supported the temporal patterns. In addition, partial least squares regression modeling could not discern a significant correlation from traditional oceanographic physical and phytoplankton nutrient parameters on overall bacterial community variability patterns at LEO-15. These results suggest factors not traditionally measured during oceanographic studies are structuring coastal microbial communities.     

Seasonal dynamics of bacterioplankton community structure in a eutrophic lake as determined by 5S rRNA analysis.

Community structure of bacterioplankton was studied during the major growth season for phytoplankton (April to October) in the epilimnion of a temperate eutrophic lake (Lake Plusssee, northern Germany) by using comparative 5S rRNA analysis. Estimates of the relative abundances of single taxonomic groups were made on the basis of the amounts of single 5S rRNA bands obtained after high-resolution electrophoresis of RNA directly from the bacterioplankton. Full-sequence analysis of single environmental 5S rRNAs enabled the identification of single taxonomic groups of bacteria. Comparison of partial 5S rRNA sequences allowed the detection of changes of single taxa over time. Overall, the whole bacterioplankton community showed two to eight abundant (>4% of the total 5S rRNA) taxa. A distinctive seasonal succession was observed in the taxonomic structure of this pelagic community. A rather-stable community structure, with seven to eight different taxonomic units, was observed beginning in April during the spring phytoplankton bloom. A strong reduction in this diversity occurred at the beginning of the clear-water phase (early May), when only two to four abundant taxa were observed, with one taxon dominating (up to 72% of the total 5S rRNA). The community structure during summer stagnation (June and July) was characterized by frequent changes of different dominating taxa. During late summer, a dinoflagellate bloom (Ceratium hirudinella) occurred, with Comamonas acidovorans (beta-subclass of the class Proteobacteria) becoming the dominant bacterial species (average abundance of 43% of the total 5S rRNA). Finally, the seasonal dynamics of the community structure of bacterioplankton were compared with the abundances of other major groups of the aquatic food web, such as phyto- and zooplankton, revealing that strong grazing pressure by zooplankton can reduce microbial diversity substantially in pelagic environments.     

Paenibacillus larvae 16S-23S rDNA intergenic transcribed spacer (ITS) regions: DNA fingerprinting and characterization

Paenibacillus larvae is the causative agent of American foulbrood in honey bee (Apis mellifera) larvae. PCR amplification of the 16S-23S ribosomal DNA (rDNA) intergenic transcribed spacer (ITS) regions, and agarose gel electrophoresis of the amplified DNA, was performed using genomic DNA collected from 134 P. larvae strains isolated in Connecticut, six Northern Regional Research Laboratory stock strains, four strains isolated in Argentina, and one strain isolated in Chile. Following electrophoresis of amplified DNA, all isolates exhibited a common migratory profile (i.e., ITS-PCR fingerprint pattern) of six DNA bands. This profile represented a unique ITS-PCR DNA fingerprint that was useful as a fast, simple, and accurate procedure for identification of P. larvae. Digestion of ITS-PCR amplified DNA, using mung bean nuclease prior to electrophoresis, characterized only three of the six electrophoresis bands as homoduplex DNA and indicating three true ITS regions. These three ITS regions, DNA migratory band sizes of 915, 1010, and 1474 bp, signify a minimum of three types of rrn operons within P. larvae. DNA sequence analysis of ITS region DNA, using P. larvae NRRL B-3553, identified the 3' terminal nucleotides of the 16S rRNA gene, 5' terminal nucleotides of the 23S rRNA gene, and the complete DNA sequences of the 5S rRNA, tRNA(ala), and tRNA(ile) genes. Gene organization within the three rrn operon types was 16S-23S, 16S-tRNA(ala)-23S, and l6S-5S-tRNA(ile)-tRNA(ala)-23S and these operons were named rrnA, rrnF, and rrnG, respectively. The 23S rRNA gene was shown by I-CeuI digestion and pulsed-field gel electrophoresis of genomic DNA to be present as seven copies. This was suggestive of seven rrn operon copies within the P. larvae genome. Investigation of the 16S-23S rDNA regions of this bacterium has aided the development of a diagnostic procedure and has helped genomic mapping investigations via characterization of the ITS regions.     

Comparison of Free-Living and Particle-Associated Bacterial Communities in the Chesapeake Bay by Stable Low-Molecular-Weight RNA Analysis

Free-living and particle-associated bacterial communities in the Chesapeake Bay estuary were analyzed and compared by using acridine orange direct counts and low-molecular-weight (LMW) RNA analysis. Samples were taken from top and bottom waters at upper- and mid-bay sites in December 1992. Free-living bacteria dominated the bacterial numbers at all sampling sites, although particle-associated bacteria increased in areas with greater particle loads. LMW RNAs (5S rRNA and tRNA) obtained directly from free-living, particle-associated, and total bacterioplankton communities were analyzed by high-resolution electrophoresis. There were distinct differences in the migration distances between LMW RNAs of free-living and particle-associated communities taken from the same site, indicating that the two communities differ in composition. In addition, LMW RNA profiles differed minimally with depth for all of the communities examined, presumably because of vertical mixing. 5S rRNAs of free-living communities from the upper- and mid-bay regions differed considerably. Particle-associated RNAs, on the other hand, were very similar, suggesting consistent environmental conditions on particles that select for similar community members. Lastly, several isolated bacteria had 5S rRNAs that were not detected in their respective extracted community 5S rRNAs, indicating that these isolated organisms were not representative of dominant members.     

Heteroduplex structures in 16S-23S rRNA intergenic transcribed spacer PCR products reveal ribosomal interoperonic polymorphisms within single Frankia strains

AIMS: Detection of polymorphisms in intergenic transcribed spacer (ITS) 16S-23S rRNA within single Frankia strains. METHODS AND RESULTS: Polymorphisms in the 16S-23S rRNA ITS were investigated in single-colony subcultures of seven Frankia isolates. Multiple ITS-polymerase chain reaction (PCR) bands were detected solely in isolates BMG5.5 and BMG5.11. The slow-migrating bands in the ITS-PCR agarose gel electrophoresis profiles of the isolates were revealed to be heteroduplexes on the basis of their migration shift in different electrophoretic matrices, southern hybridization and the single-strand DNA mung bean endonuclease digestion. Laser-scanned capillary electrophoresis detected two ITS-PCR fragments differing in length by three and six nucleotide insertions/deletions in strains BMG5.5 and BMG5.11, respectively. Sequence analysis of the cloned ITS showed that in strain BMG5.5 the two ITS differed by the presence of three to four copies of the 3-bp tandem repeat 5'-TGG-3'. In strain BMG5.11, the two ITS differed by the presence of two to three copies of the 6-bp tandem repeat 5'-CTTGGG-3'. CONCLUSIONS: We demonstrate the occurrence of ITS 16S-23S rRNa polymorphisms within single Frankia strains. SIGNIFICANCE AND IMPACT OF THE STUDY: We reported the occurrence of ITS 16S-23S rRNA polymorphisms within single Frankia strains from Elaeagnus host group recognized as the more flexible strains within Frankia genus. Furthermore, we underscored the applied interest of strains BMG5.11 and BMG5.5 in future ecological studies using ITS 16S-23S rRNA as molecular marker.     

Respiratory succession and community succession of bacterioplankton in seasonally anoxic estuarine waters

Anoxia occurs in bottom waters of stratified estuaries when respiratory consumption of oxygen, primarily by bacteria, outpaces atmospheric and photosynthetic reoxygenation. Once water becomes anoxic, bacterioplankton must change their metabolism to some form of anaerobic respiration. Analysis of redox chemistry in water samples spanning the oxycline of Chesapeake Bay during the summer of 2004 suggested that there was a succession of respiratory metabolism following the loss of oxygen. Bacterial community doubling time, calculated from bacterial abundance (direct counts) and production (anaerobic leucine incorporation), ranged from 0.36 to 0.75 day and was always much shorter than estimates of the time that the bottom water was anoxic (18 to 44 days), indicating that there was adequate time for bacterial community composition to shift in response to changing redox conditions. However, community composition (as determined by PCR-denaturing gradient gel electrophoresis analysis of 16S rRNA genes) in anoxic waters was very similar to that in surface waters in June when nitrate respiration was apparent in the water column and only partially shifted away from the composition of the surface community after nitrate was depleted. Anoxic water communities did not change dramatically until August, when sulfate respiration appeared to dominate. Surface water populations that remained dominant in anoxic waters were Synechococcus sp., Gammaproteobacteria in the SAR86 clade, and Alphaproteobacteria relatives of Pelagibacter ubique, including a putative estuarine-specific Pelagibacter cluster. Populations that developed in anoxic water were most similar (<92% similarity) to uncultivated Firmicutes, uncultivated Bacteroidetes, Gammaproteobacteria in the genus Thioalcalovibrio, and the uncultivated SAR406 cluster. These results indicate that typical estuarine bacterioplankton switch to anaerobic metabolism under anoxic conditions but are ultimately replaced by different organisms under sulfidic conditions.     

An eco-informatics tool for microbial community studies: supervised classification of Amplicon Length Heterogeneity (ALH) profiles of 16S rRNA

Support vector machines (SVM) and K-nearest neighbors (KNN) are two computational machine learning tools that perform supervised classification. This paper presents a novel application of such supervised analytical tools for microbial community profiling and to distinguish patterning among ecosystems. Amplicon length heterogeneity (ALH) profiles from several hypervariable regions of 16S rRNA gene of eubacterial communities from Idaho agricultural soil samples and from Chesapeake Bay marsh sediments were separately analyzed. The profiles from all available hypervariable regions were concatenated to obtain a combined profile, which was then provided to the SVM and KNN classifiers. Each profile was labeled with information about the location or time of its sampling. We hypothesized that after a learning phase using feature vectors from labeled ALH profiles, both these classifiers would have the capacity to predict the labels of previously unseen samples. The resulting classifiers were able to predict the labels of the Idaho soil samples with high accuracy. The classifiers were less accurate for the classification of the Chesapeake Bay sediments suggesting greater similarity within the Bay's microbial community patterns in the sampled sites. The profiles obtained from the V1+V2 region were more informative than that obtained from any other single region. However, combining them with profiles from the V1 region (with or without the profiles from the V3 region) resulted in the most accurate classification of the samples. The addition of profiles from the V 9 region appeared to confound the classifiers. Our results show that SVM and KNN classifiers can be effectively applied to distinguish between eubacterial community patterns from different ecosystems based only on their ALH profiles.     

A comparison of phytoplankton assemblages in the Chesapeake and Delaware estuaries (USA), with emphasis on diatoms

The Delaware Bay is characterized as having greater nutrient and turbidity levels than the Chesapeake Bay. In reference to these differences, a one year study was conducted to identify any similarities and differences in the phytoplankton populations in these estuaries. The results indicated patterns of similarity in the diatom composition, with the total phytoplankton assemblage forming two site groups along a salinity gradient in each bay. These site groups were associated with stations located in the tidal fresh-oligohaline and meso-polyhaline regions of both estuaries. The seasonal concentrations of diatoms and total phytoplankton in both of these regions were higher in the Chesapeake Bay.Subtle differences between the two estuaries include a more diversified and abundant assemblage of neritic phytoplankters (including dinoflagellates) are present in the lower Chesapeake Bay. In contrast, a diatom dominated community is more characteristic of Delaware Bay. It is suggested the entry of neritic species into lower regions of the estuaries was enhanced by the reduced amount of rainfall and flow rates that occurred during the study period. The greater success of neritic species in the Chesapeake Bay is attributed to the lower turbidity of that estuary compared to Delaware Bay.     

Is the 16S-23S rRNA internal transcribed spacer region a good tool for use in molecular systematics and population genetics? A case study in cyanobacteria

We amplified, TA-cloned, and sequenced the 16S-23S internal transcribed spacer (ITS) regions from single isolates of several cyanobacterial species, Calothrix parietina, Scytonema hyalinum, Coelodesmium wrangelii, Tolypothrix distorta, and a putative new genus (isolates SRS6 and SRS70), to investigate the potential of this DNA sequence for phylogenetic and population genetic studies. All isolates carried ITS regions containing the sequences coding for two tRNA molecules (tRNA and tRNA). We retrieved additional sequences without tRNA features from both C. parietina and S. hyalinum. Furthermore, in S. hyalinum, we found two of these non-tRNA-encoding regions to be identical in length but different in sequence. This is the first report of ITS regions from a single cyanobacterial isolate not only different in configuration, but also, within one configuration, different in sequence. The potential of the ITS region as a tool for studying molecular systematics and population genetics is significant, but the presence of multiple nonidentical rRNA operons poses problems. Multiple nonidentical rRNA operons may impact both studies that depend on comparisons of phylogenetically homologous sequences and those that employ restriction enzyme digests of PCR products. We review current knowledge of the numbers and kinds of 16S-23S ITS regions present across bacterial groups and plastids, and we discuss broad patterns congruent with higher-level systematics of prokaryotes.