Resolution of Prochlorococcus and Synechococcus Ecotypes by Using 16S-23S Ribosomal DNA Internal Transcribed Spacer Sequences

Cultured isolates of the marine cyanobacteria Prochlorococcus and Synechococcus vary widely in their pigment compositions and growth responses to light and nutrients, yet show greater than 96% identity in their 16S ribosomal DNA (rDNA) sequences. In order to better define the genetic variation that accompanies their physiological diversity, sequences for the 16S-23S rDNA internal transcribed spacer (ITS) region were determined in 32 Prochlorococcus isolates and 25 Synechococcus isolates from around the globe. Each strain examined yielded one ITS sequence that contained two tRNA genes. Dramatic variations in the length and G+C content of the spacer were observed among the strains, particularly among Prochlorococcus strains. Secondary-structure models of the ITS were predicted in order to facilitate alignment of the sequences for phylogenetic analyses. The previously observed division of Prochlorococcus into two ecotypes (called high and low-B/A after their differences in chlorophyll content) were supported, as was the subdivision of the high-B/A ecotype into four genetically distinct clades. ITS-based phylogenies partitioned marine cluster A Synechococcus into six clades, three of which can be associated with a particular phenotype (motility, chromatic adaptation, and lack of phycourobilin). The pattern of sequence divergence within and between clades is suggestive of a mode of evolution driven by adaptive sweeps and implies that each clade represents an ecologically distinct population. Furthermore, many of the clades consist of strains isolated from disparate regions of the world's oceans, implying that they are geographically widely distributed. These results provide further evidence that natural populations of Prochlorococcus and Synechococcus consist of multiple coexisting ecotypes, genetically closely related but physiologically distinct, which may vary in relative abundance with changing environmental conditions.     

High vertical and low horizontal diversity of Prochlorococcus ecotypes in the Mediterranean Sea in summer

Natural populations of the marine cyanobacterium Prochlorococcus exist as two main ecotypes, inhabiting different layers of the ocean's photic zone. These so-called high light- (HL-) and low light (LL-) adapted ecotypes are both physiologically and genetically distinct. HL strains can be separated into two major clades (HLI and HLII), whereas LL strains are more diverse. Here, we used several molecular techniques to study the genetic diversity of natural Prochlorococcus populations during the Prosope cruise in the Mediterranean Sea in the summer of 1999. Using a dot blot hybridization technique, we found that HLI was the dominant HL group and was confined to the upper mixed layer. In contrast, LL ecotypes were only found below the thermocline. Secondly, a restriction fragment length polymorphism analysis of PCR-amplified pcb genes (encoding the major light-harvesting proteins of Prochlorococcus) suggested that there were at least four genetically different ecotypes, occupying distinct but overlapping light niches in the photic zone. At comparable depths, similar banding patterns were observed throughout the sampled area, suggesting a horizontal homogenization of ecotypes. Nevertheless, environmental pcb gene sequences retrieved from different depths at two stations proved all different at the nucleotide level, suggesting a large genetic microdiversity within those ecotypes.     

Phylogenetic diversity of Synechococcus strains isolated from the East China Sea and the East Sea

Phylogenetic relationships among 33 Synechococcus strains isolated from the East China Sea (ECS) and the East Sea (ES) were studied based on 16S rRNA gene sequences and 16S–23S rRNA gene internal transcribed spacer (ITS) sequences. Pigment patterns of the culture strains were also examined. Based on 16S rRNA gene and ITS sequence phylogenies, the Synechococcus isolates were clustered into 10 clades, among which eight were previously identified and two were novel. Half of the culture strains belonged to clade V or VI. All strains that clustered into novel clades exhibited both phycoerythrobilin and phycourobilin. Interestingly, the pigment compositions of isolates belonging to clades V and VI differed from those reported for other oceanic regions. None of the isolates in clade V showed phycourobilin, whereas strains in clade VI exhibited both phycourobilin and phycoerythrobilin, which is in contrast to previous studies. The presence of novel lineages and the different pigment patterns in the ECS and the ES suggests the possibility that some Synechococcus lineages are distributed only in geographically restricted areas and have evolved in these regions. Therefore, further elucidation of the physiological, ecological, and genetic characteristics of the diverse Synechococcus strains is required to understand their spatial and geographical distribution.     

Culture isolation and culture-independent clone libraries reveal new marine Synechococcus ecotypes with distinctive light and N physiologies

Marine microbial communities often contain multiple closely related phylogenetic clades, but in many cases, it is still unclear what physiological traits differentiate these putative ecotypes. The numerically abundant marine cyanobacterium Synechococcus can be divided into at least 14 clades. In order to better understand ecotype differentiation in this genus, we assessed the diversity of a Synechococcus community from a well-mixed water column in the Sargasso Sea during March 2002, a time of year when this genus typically reaches its annual peak in abundance. Diversity was estimated from water sampled at three depths (approximately 5, 70, and 170 m) using both culture isolation and construction of cyanobacterial 16S-23S rRNA internal transcribed sequence clone libraries. Clonal isolates were obtained by enrichment with ammonium, nitrite, or nitrate as the sole N source, followed by pour plating. Each method sampled the in situ diversity differently. The combined methods revealed a total of seven Synechococcus phylotypes including two new putative ecotypes, labeled XV and XVI. Although most other isolates grow on nitrate, clade XV exhibited a reduced efficiency in nitrate utilization, and both clade XV and XVI are capable of chromatic adaptation, demonstrating that this trait is more widely distributed among Synechococcus strains than previously known. Thus, as in its sister genus Prochlorococcus, light and nitrogen utilization are important factors in ecotype differentiation in the marine Synechococcus lineage.     

Chromatic adaptation in marine Synechococcus strains

Characterization of two genetically distinct groups of marine Synechococcus sp. strains shows that one, but not the other, increases its phycourobilin/phycoerythrobilin chromophore ratio when growing in blue light. This ability of at least some marine Synechococcus strains to chromatically adapt may help explain their greater abundance in particular ocean environments than cyanobacteria of the genus Prochlorococcus.     

Measurement of Prochlorococcus ecotypes using real-time polymerase chain reaction reveals different abundances of genotypes with similar light physiologies

Prochlorococcus is a marine cyanobacterium which is found at high abundances in world's tropical and subtropical oligotrophic oceans. The genus Prochlorococcus can be divided into two major groups based on light physiology. Both of these groups can be further subdivided into genetically distinct lineages, or ecotypes. Real-time polymerase chain reaction (PCR) assays based on sequence differences in the 16S-23S rDNA internal transcribed spacer or the 23S rDNA were developed to examine the distribution of each ecotype in the field. The real-time PCR assays enabled linear quantification of concentrations ranging from 10 to 4 x 10(5) cells ml(-1). These assays were applied to a stratified water column in the Sargasso Sea. The majority of Prochlorococcus cells above 110 m belonged to the one of the low chlorophyll b/a ratio (high-light adapted) ecotypes, while two types of high chlorophyll b/a ratio (low-light adapted) cells dominated below 110 m. The other three types were found at significantly lower numbers or not detected at all. Differences in the abundance of ecotypes within the major light physiology groupings suggest that other factors, such as nutrient utilization and differential mortality, are driving their relative distributions. Real-time PCR assays will enable further exploration of these factors and temporal and geographic variability in ecotype abundance.     

Ecotype diversity in the marine picoeukaryote Ostreococcus (Chlorophyta, Prasinophyceae)

The importance of the cyanobacteria Prochlorococcus and Synechococcus in marine ecosystems in terms of abundance and primary production can be partially explained by ecotypic differentiation. Despite the dominance of eukaryotes within photosynthetic picoplankton in many areas a similar differentiation has never been evidenced for these organisms. Here we report distinct genetic [rDNA 18S and internal transcribed spacer (ITS) sequencing], karyotypic (pulsed-field gel electrophoresis), phenotypic (pigment composition) and physiological (light-limited growth rates) traits in 12 Ostreococcus strains (Prasinophyceae) isolated from various marine environments and depths, which suggest that the concept of ecotype could also be valid for eukaryotes. Internal transcribed spacer phylogeny grouped together four deep strains isolated between 90 m and 120 m depth from different geographical origins. Three deep strains displayed larger chromosomal bands, different chromosome hybridization patterns, and an additional chlorophyll (chl) c-like pigment. Furthermore, growth rates of deep strains show severe photo-inhibition at high light intensities, while surface strains do not grow at the lowest light intensities. These features strongly suggest distinct adaptation to environmental conditions encountered at surface and the bottom of the oceanic euphotic zone, reminiscent of that described in prokaryotes.     

Basin-scale distribution patterns of picocyanobacterial lineages in the Atlantic Ocean

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus are major contributors to oceanic primary production. The genera are genetically diverse, comprising several known ecotypes or lineages. However, little is known of the distribution of these lineages over large geographic areas. Here, we analysed the relative abundance of Prochlorococcus ecotypes and Synechococcus lineages at the ocean basin scale along an Atlantic Meridional Transect (AMT) using dot blot hybridization and fluorescence in situ hybridization (FISH) techniques. The transect covered several contrasting oceanic provinces (gyres, upwelling, temperate regions) as well as environmentally 'equivalent' regions in the northern and southern hemisphere (northern and southern gyres and temperate regions). Flow cytometric data revealed a discrete separation in abundance of major picocyanobacterial genera. Prochlorococcus reached highest abundance in oligotrophic regions, while more mesotrophic waters were dominated by Synechococcus. Individual genetic lineages of both Prochlorococcus and Synechococcus showed highly similar distributions in corresponding regions in the northern and southern hemisphere. In addition, Prochlorococcus showed a distinctive depth distribution, with HLI and HLII ecotypes near the surface and co-occurring LL ecotypes further down in the water column. Conversely, Synechococcus generally revealed no obvious depth preference, but did show highly specific distribution at the horizontal scale, with clades I and IV particularly dominating temperate, mesotrophic waters in both the northern and southern hemispheres. The data clearly reveal that specific picocyanobacterial lineages proliferate in similar oceanic provinces separated by large spatial scales. Furthermore, comparison with an earlier AMT dataset suggests that basin scale distribution patterns for Prochlorococcus ecotypes are remarkably reproducible from year to year.     

CONSPECIFICITY AND INDO-PACIFIC DISTRIBUTION OF SYMBIODINIUM GENOTYPES (DINOPHYCEAE) FROM GIANT CLAMS

We previously reported the occurrence of genetically‐diverse symbiotic dinoflagellates (zooxanthellae) within and between 7 giant clam species (Tridacnidae) from the Philippines based on the algal isolates' allozyme and random amplified polymorphic DNA (RAPD) patterns. We also reported that these isolates all belong to clade A of the Symbiodinium phylogeny with identical 18S rDNA sequences. Here we extend the genetic characterization of Symbiodinium isolates from giant clams and propose that they are conspecific. We used the combined DNA sequences of the internal transcribed spacer (ITS)1, 5.8S rDNA, and ITS2 regions (rDNA‐ITS region) because the ITS1 and ITS2 regions evolve faster than 18S rDNA and have been shown to be useful in distinguishing strains of other dinoflagellates. DGGE of the most variable segment of the rDNA‐ITS region, ITS1, from clonal representatives of clades A, B, and C showed minimal intragenomic variation. The rDNA‐ITS region shows similar phylogenetic relationships between Symbiodinium isolates from symbiotic bivalves and some cnidarians as does 18S rDNA, and that there are not many different clade A species or strains among cultured zooxanthellae (CZ) from giant clams. The CZ from giant clams had virtually identical sequences, with only a single nucleotide difference in the ITS2 region separating two groups of isolates. These data suggest that there is one CZ species and perhaps two CZ strains, each CZ strain containing individuals that have diverse allozyme and RAPD genotypes. The CZ isolated from giant clams from different areas in the Philippines (21 isolates, 7 clam species), the Australian Great Barrier Reef (1 isolate, 1 clam species), Palau (8 isolates, 7 clam species), and Okinawa, Japan (1 isolate, 1 clam species) shared the same rDNA‐ITS sequences. Furthermore, analysis of fresh isolates from giant clams collected from these geographical areas shows that these bivalves also host indistinguishable clade C symbionts. These data demonstrate that conspecific Symbiodinium genotypes, particularly clade A symbionts, are distributed in giant clams throughout the Indo‐Pacific.     

Effect of temperature and light on growth of and photosynthesis by Synechococcus isolates typical of those predominating in the octopus spring microbial mat community of Yellowstone National Park

Previous molecular analysis of the Octopus Spring cyanobacterial mat revealed numerous genetically distinct 16S rRNA sequences from predominant Synechococcus populations distantly related to the readily cultivated unicellular cyanobacterium Synechococcus lividus. Patterns in genotype distribution relative to temperature and light conditions suggested that the organisms contributing these 16S rRNA sequences may fill distinct ecological niches. To test this hypothesis, Synechococcus isolates were cultivated using a dilution and filtration approach and then shown to be genetically relevant to natural mat populations by comparisons of similarities of 16S rRNA genes and 16S-23S internal transcribed spacer (ITS) regions. Most isolates were identical or nearly identical at both loci to predominant mat genotypes; others showed 1- to 2-nucleotide differences at the 16S rRNA locus and even greater difference in ITS sequences. Isolates with predominant mat genotypes had distinct temperature ranges and optima for growth that were consistent with their distributions in the mat. Isolates with genotypes not previously detected or known to be predominant in the mat exhibited temperature ranges and optima that were not representative of predominant mat populations and also grew more slowly. Temperature effects on photosynthesis did not reflect temperature relations for growth. However, the isolate with the highest temperature optimum and upper limit was capable of performing photosynthesis at a higher temperature than other isolates. Growth rate and photosynthetic responses provided evidence for light acclimation but evidence of, at best, only subtle light adaptation.     

Clade-specific 16S ribosomal DNA oligonucleotides reveal the predominance of a single marine Synechococcus clade throughout a stratified water column in the Red Sea

Phylogenetic relationships among members of the marine Synechococcus genus were determined following sequencing of the 16S ribosomal DNA (rDNA) from 31 novel cultured isolates from the Red Sea and several other oceanic environments. This revealed a large genetic diversity within the marine Synechococcus cluster consistent with earlier work but also identified three novel clades not previously recognized. Phylogenetic analyses showed one clade, containing halotolerant isolates lacking phycoerythrin (PE) and including strains capable, or not, of utilizing nitrate as the sole N source, which clustered within the MC-A (Synechococcus subcluster 5.1) lineage. Two copies of the 16S rRNA gene are present in marine Synechococcus genomes, and cloning and sequencing of these copies from Synechococcus sp. strain WH 7803 and genomic information from Synechococcus sp. strain WH 8102 reveal these to be identical. Based on the 16S rDNA sequence information, clade-specific oligonucleotides for the marine Synechococcus genus were designed and their specificity was optimized. Using dot blot hybridization technology, these probes were used to determine the in situ community structure of marine Synechococcus populations in the Red Sea at the time of a Synechococcus maximum during April 1999. A predominance of genotypes representative of a single clade was found, and these genotypes were common among strains isolated into culture. Conversely, strains lacking PE, which were also relatively easily isolated into culture, represented only a minor component of the Synechococcus population. Genotypes corresponding to well-studied laboratory strains also appeared to be poorly represented in this stratified water column in the Red Sea.     

Light variability illuminates niche-partitioning among marine Picocyanobacteria

Prochlorococcus and Synechococcus picocyanobacteria are dominant contributors to marine primary production over large areas of the ocean. Phytoplankton cells are entrained in the water column and are thus often exposed to rapid changes in irradiance within the upper mixed layer of the ocean. An upward fluctuation in irradiance can result in photosystem II photoinactivation exceeding counteracting repair rates through protein turnover, thereby leading to net photoinhibition of primary productivity, and potentially cell death. Here we show that the effective cross-section for photosystem II photoinactivation is conserved across the picocyanobacteria, but that their photosystem II repair capacity and protein-specific photosystem II light capture are negatively correlated and vary widely across the strains. The differences in repair rate correspond to the light and nutrient conditions that characterize the site of origin of the Prochlorococcus and Synechococcus isolates, and determine the upward fluctuation in irradiance they can tolerate, indicating that photoinhibition due to transient high-light exposure influences their distribution in the ocean.     

Phylogenetic characterization of a new thermoacidophilic bacterium isolated from hot springs in Japan

Three strains of thermophilic-acidophilic bacteria isolated previously from different hot springs in Japan were characterized by molecular genetic methods. The strategy taken involved PCR amplification, sequencing and restriction pattern analysis of 16S rDNA, 16S-23S rDNA spacer polymorphism analysis and genomic DNA-DNA hybridization. A phylogenetic analysis based on 16S rDNA sequences showed that the new thermoacidophilic isolates formed a genetically coherent group at the species level and fell into a major cluster together with members of the genera Alicyclobacillus and Sulfobacillus with A. acidocaldarius and A. acidoterrestris as their closest relatives. The levels of binary sequence similarity between the isolates and the two Alicyclobacillus species were 97.6 to 97.9%, values considered low enough to warrant placement of the isolates in a distinct species of the genus Alicyclobacillus. The 16S rDNA restriction pattern analysis, but not 16S-23S rDNA spacer polymorphism analysis, was useful for differentiating the isolates from the established Alicyclobacillus species. DNA-DNA hybridization assays demonstrated a distinct phylogenetic position of our isolates as a genospecies within the genus Alicyclobacillus. On the basis of these results, the thermoacidophilic isolates should be classified into a new species of Alicyclobacillus. The results of this study suggest that this new genospecies of Alicyclobacillus is widely distributed in hot springs in Japan.     

Molecular genetic variation within and among isolates of QPX (Thraustochytridae), a parasite of the hard clam Mercenaria mercenaria

The thraustochytrid known as QPX (Quahog Parasite Unknown) has sporadically caused disease in the hard clam Mercenaria mercenaria along the east coast of North America since the 1960s. We hypothesized that genetically distinct QPX strains might be responsible for outbreaks of QPX disease in different areas and tested this hypothesis by comparing several QPX isolates recovered from the recent outbreak in Raritan Bay, New York with QPX strains isolated from 2 outbreaks in Massachusetts, USA. There was no variation in small subunit rDNA (SSU rDNA), 5.8S rDNA, or 4 mitochondrial gene sequences. In contrast, both of the ribosomal ribonucleic acid (rRNA) operon intergenic spacers, internal transcribed spacers 1 and 2 (ITS1 and ITS2), revealed substantial sequence variation. However, strain-specific sequences were not detected because the ITS sequence variation within QPX isolates was comparable to the variation between isolates. ITS1 sequences recovered from an infected clam by amplification with a QPX ITS2-specific primer were identical to those recovered from the QPX isolates.     

Influence of molecular resolution on sequence-based discovery of ecological diversity among Synechococcus populations in an alkaline siliceous hot spring microbial mat

Previous research has shown that sequences of 16S rRNA genes and 16S-23S rRNA internal transcribed spacer regions may not have enough genetic resolution to define all ecologically distinct Synechococcus populations (ecotypes) inhabiting alkaline, siliceous hot spring microbial mats. To achieve higher molecular resolution, we studied sequence variation in three protein-encoding loci sampled by PCR from 60°C and 65°C sites in the Mushroom Spring mat (Yellowstone National Park, WY). Sequences were analyzed using the ecotype simulation (ES) and AdaptML algorithms to identify putative ecotypes. Between 4 and 14 times more putative ecotypes were predicted from variation in protein-encoding locus sequences than from variation in 16S rRNA and 16S-23S rRNA internal transcribed spacer sequences. The number of putative ecotypes predicted depended on the number of sequences sampled and the molecular resolution of the locus. Chao estimates of diversity indicated that few rare ecotypes were missed. Many ecotypes hypothesized by sequence analyses were different in their habitat specificities, suggesting different adaptations to temperature or other parameters that vary along the flow channel.     

Cryptic diversity of the symbiotic cyanobacterium Synechococcus spongiarum among sponge hosts

Cyanobacteria are common members of sponge-associated bacterial communities and are particularly abundant symbionts of coral reef sponges. The unicellular cyanobacterium Synechococcus spongiarum is the most prevalent photosynthetic symbiont in marine sponges and inhabits taxonomically diverse hosts from tropical and temperate reefs worldwide. Despite the global distribution of S. spongiarum , molecular analyses report low levels of genetic divergence among 16S ribosomal RNA (rRNA) gene sequences from diverse sponge hosts, resulting either from the widespread dispersal ability of these symbionts or the low phylogenetic resolution of a conserved molecular marker. Partial 16S rRNA and entire 16S–23S rRNA internal transcribed spacer (ITS) genes were sequenced from cyanobacteria inhabiting 32 sponges (representing 18 species, six families and four orders) from six geographical regions. ITS phylogenies revealed 12 distinct clades of S. spongiarum that displayed 9% mean sequence divergence among clades and less than 1% sequence divergence within clades. Symbiont clades ranged in specificity from generalists to specialists, with most (10 of 12) clades detected in one or several closely related hosts. Although multiple symbiont clades inhabited some host sponges, symbiont communities appear to be structured by both geography and host phylogeny. In contrast, 16S rRNA sequences were highly conserved, exhibiting less than 1% sequence divergence among symbiont clades. ITS gene sequences displayed much higher variability than 16S rRNA sequences, highlighting the utility of ITS sequences in determining the genetic diversity and host specificity of S. spongiarum populations among reef sponges. The genetic diversity of S. spongiarum revealed by ITS sequences may be correlated with different physiological capabilities and environmental preferences that may generate variable host–symbiont interactions.     

Multi-locus sequence analysis, taxonomic resolution and biogeography of marine Synechococcus

Conserved markers such as the 16S rRNA gene do not provide sufficient molecular resolution to identify spatially structured populations of marine Synechococcus, or 'ecotypes' adapted to distinct ecological niches. Multi-locus sequence analysis targeting seven 'core' genes was employed to taxonomically resolve Synechococcus isolates and correlate previous phylogenetic analyses encompassing a range of markers. Despite the recognized importance of lateral gene transfer in shaping the genomes of marine cyanobacteria, multi-locus sequence analysis of more than 120 isolates reflects a clonal population structure of major lineages and subgroups. A single core genome locus, petB, encoding the cytochrome b(6) subunit of the cytochrome b(6) f complex, was selected to expand our understanding of the diversity and ecology of marine Synechococcus populations. Environmental petB sequences cloned from contrasting sites highlight numerous genetically and ecologically distinct clusters, some of which represent novel, environmentally abundant clades without cultured representatives. With a view to scaling ecological analyses, the short sequence, taxonomic resolution and accurate automated alignment of petB is ideally suited to high-throughput and high-resolution sequencing projects to explore links between the ecology, evolution and biology of marine Synechococcus.     

Glucosylglycerate: a secondary compatible solute common to marine cyanobacteria from nitrogen-poor environments

The synthesis and accumulation of compatible solutes represent an essential part of the salt acclimation strategy of microorganisms. Glucosylglycerol is considered to be the typical compatible solute among marine cyanobacteria. However, genes that encode enzymes for the synthesis of glucosylglycerol were not detected in the genome sequences of marine picoplanktonic Prochlorococcus strains. Instead, we noticed the presence of genes that putatively encode for glucosylglycerate (GGA) synthesis among Prochlorococcus and most other closely related marine picocyanobacteria. Recombinant proteins from Prochlorococcus marinus SS120 and Synechococcus sp. PCC 7002 exhibited glucosyl-phosphoglycerate synthase (GpgS) activity, and GpgS is a key enzyme of GGA synthesis. GGA accumulation was found to be salt- as well as nitrogen-regulated in the coastal strain Synechococcus sp. PCC 7002. Moreover, GGA was also detected in all picoplanktonic Prochlorococcus and Synechococcus strains harbouring gpgS genes, especially under N-limiting conditions. These results suggest that marine picocyanobacteria acquired the capacity to synthesize the negatively charged compound GGA during their evolution. Our results establish GGA as the fifth most widespread compatible solute among cyanobacteria. Additionally, GGA appears to replace glutamate as an anion to counter monovalent cations in marine picocyanobacteria from N-poor environments.     

Taxonomic resolution, ecotypes and the biogeography of Prochlorococcus

In order to expand our understanding of the diversity and biogeography of Prochlorococcus ribotypes, we PCR-amplified, cloned and sequenced the 16S/23S rRNA ITS region from sites in the Atlantic and Pacific oceans. Ninety-three per cent of the ITS sequences could be assigned to existing Prochlorococcus clades, although many novel subclades were detected. We assigned the sequences to operational taxonomic units using a graduated scale of sequence identity from 80% to 99.5% and correlated Prochlorococcus diversity with respect to environmental variables and dispersal time between the sites. Dispersal time was estimated using a global ocean circulation model. The significance of specific environmental variables was dependent on the degree of sequence identity used to define a taxon: light correlates with broad-scale diversity (90% cut-off), temperature with intermediate scale (95%) whereas no correlation with phosphate was observed. Community structure was correlated with dispersal time between sample sites only when taxa were defined using the finest sequence similarity cut-off. Surprisingly, the concentration of nitrate, which cannot be used as N source by the Prochlorococcus strains in culture, explains some variation in community structure for some definitions of taxa. This study suggests that the spatial distribution of Prochlorococcus ecotypes is shaped by a hierarchy of environmental factors as well dispersal limitation.