Ecological Genomics of Marine Picocyanobacteria

Summary: Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus numerically dominate the picophytoplankton of the world ocean, making a key contribution to global primary production. Prochlorococcus was isolated around 20 years ago and is probably the most abundant photosynthetic organism on Earth. The genus comprises specific ecotypes which are phylogenetically distinct and differ markedly in their photophysiology, allowing growth over a broad range of light and nutrient conditions within the 45°N to 40°S latitudinal belt that they occupy. Synechococcus and Prochlorococcus are closely related, together forming a discrete picophytoplankton clade, but are distinguishable by their possession of dissimilar light-harvesting apparatuses and differences in cell size and elemental composition. Synechococcus strains have a ubiquitous oceanic distribution compared to that of Prochlorococcus strains and are characterized by phylogenetically discrete lineages with a wide range of pigmentation. In this review, we put our current knowledge of marine picocyanobacterial genomics into an environmental context and present previously unpublished genomic information arising from extensive genomic comparisons in order to provide insights into the adaptations of these marine microbes to their environment and how they are reflected at the genomic level.     

Differential Sunlight Sensitivity of Picophytoplankton from Surface Mediterranean Coastal Waters

We tested the sensitivity of coastal picophytoplankton exposed to natural sunlight in short-term experiments. Cell abundance and cell-specific chlorophyll fluorescence were significantly reduced in Prochlorococcus spp. but not in Synechococcus, whereas picoeukaryotes had an intermediate response. These results are the first direct evidence of a differential sensitivity to sunlight of these ubiquitous marine members of unicellular phytoplankton.     

Phylogenetic Relationships of the Glycolytic Enzyme, Glyceraldehyde-3-Phosphate Dehydrogenase, from Parabasalid Flagellates

Over 90% of the open reading frame of gap genes for glycolytic glyceraldehyde-3-phosphate dehydrogenase (GAPDH; EC 1.2.1.12) was obtained with PCR from five species of Parabasala. With gap1 from Trichomonas vaginalis obtained earlier, the data include two sequences each for three species. All sequences were colinear with T. vaginalis gap1 and shared with it as a synapomorphy a 10- to 11-residue insertion not found in any other gap and an S-loop with characteristic features of eubacterial GAPDH. All residues known to be highly conserved in this enzyme were present. The parabasalid sequences formed a robust monophyletic group in phylogenetic reconstructions with distance-based, maximum-parsimony, and maximum-likelihood methods. The two genes of the amphibian commensal, Trichomitus batrachorum, shared a common ancestor with the rest, which separate into two well-supported lineages. T. vaginalis and Tetratrichomonas gallinarum (both representatives of Trichomonadinae) formed one, while Monocercomonas sp. and Tritrichomonas foetus formed the other. These data agreed with and/or were close to published reconstructions based on other macromolecules. They did not support the ancestral position of Monocercomonas sp. proposed on the basis of morphological characteristics but confirmed an early emergence of Trichomitus batrachorum. The sequence pairs obtained from three species indicated either gene duplications subsequent to the divergence of the corresponding lineages or a strong gene conversion later in these lineages. The parabasalid clade was a robust part of the eubacterial radiation of GAPDH and showed no relationships to the clade that contained all other eukaryotic gap genes. The data clearly reveal that the members of this lineage use in their glycolytic pathway a GAPDH species with properties and an evolutionary history that are unique among all eukaryotes studied so far. Received: 28 April 1997 / Accepted: 14 October 1997     

Ancyromonadida: A New Phylogenetic Lineage Among the Protozoa Closely Related to the Common Ancestor of Metazoans, Fungi, and Choanoflagellates (Opisthokonta)

Molecular and morphological evidence points to the ancyromonad Ancyromonas as a plausible candidate for the closest relative to the common ancestor of metazoans, fungi, and choanoflagellates (the Opisthokonta). Using 18S rDNA sequences from most of the major eukaryotic lineages, maximum-likelihood, minimum-evolution, and maximum-parsimony analyses yielded congruent phylogenies supporting this hypothesis. Combined with ultrastructural similarities between Ancyromonas and opisthokonts, the evidence presented here suggests that Ancyromonas may form an independent lineage, the Ancyromonadida Cavalier-Smith 1997, closer in its relationship to the opisthokonts than is its nearest protist relatives, the Apusomonadida. However, the very low bootstrap support for deep nodes and hypothesis testing indicate that the resolving power of 18S rDNA sequences is limited for examining this aspect of eukaryotic phylogeny. Alternate branching positions for the Ancyromonas lineage cannot be robustly rejected, revealing the importance of ultrastructure when examining the origins of multicellularity. The future use of a multigene approach may additionally be needed to resolve this aspect of eukaryotic phylogeny. Received: 27 March 2000 / Accepted: 12 June 2000     

Accelerated Evolution of Cytochrome b in Simian Primates: Adaptive Evolution in Concert with Other Mitochondrial Proteins?

We have sequenced the cytochrome b gene of Horsfield's tarsier, Tarsius bancanus, to complete a data set of sequences for this gene from representatives of each primate infraorder. These primate cytochrome b sequences were combined with those from representatives of three other mammalian orders (cat, whale, and rat) in an analysis of relative evolutionary rates. The nonsynonymous nucleotide substitution rate of the cytochrome b gene has increased approximately twofold along lineages leading to simian primates compared to that of the tarsier and other primate and nonprimate mammalian species. However, the rate of transversional substitutions at fourfold degenerate sites has remained uniform among all lineages. This increase in the evolutionary rate of cytochrome b is similar in character and magnitude to that described previously for the cytochrome c oxidase subunit II gene. We propose that the evolutionary rate increase observed for cytochrome b and cytochrome c oxidase subunit II may underlie an episode of coadaptive evolution of these two proteins in the mitochondria of simian primates. Received: 15 December 1997 / Accepted: 24 February 1998     

The Origin of Chlorarachniophyte Plastids, as Inferred from Phylogenetic Comparisons of Amino Acid Sequences of EF-Tu

A molecular phylogenetic analysis of elongation factor Tu (EF-Tu) proteins from plastids was performed in an attempt to identify the origin of chlorarachniophyte plastids, which are considered to have evolved from the endosymbiont of a photosynthetic eukaryote. Partial sequences of the genes for plastid EF-Tu proteins (1,080–1,089 bp) were determined for three algae that contain chlorophyll b, namely, Gymnochlora stellata (Chlorarachniophyceae), Bryopsis maxima (Ulvophyceae), and Pyramimonas disomata (Prasinophyceae). The deduced amino acid sequences were used to construct phylogenetic trees of the plastid and bacterial EF-Tu proteins by the maximum likelihood, the maximum parsimony, and the neighbor joining methods. The trees obtained in the present analysis suggest that all plastids that contain chlorophyll b are monophyletic and that the chlorarachniophyte plastids are closely related to those of the Ulvophyceae. The phylogenetic trees also suggest that euglenophyte plastids are closely related to prasinophycean plastids. The results indicate that the chlorarachniophyte plastids evolved from a green algal endosymbiont that was closely related to the Ulvophyceae and that at least two secondary endosymbiotic events have occurred in the lineage of algae with plastids that contain chlorophyll b. Received: 10 March 1997 / Accepted: 28 July 1997     

Modeling Selective Pressures on Phytoplankton in the Global Ocean

Our view of marine microbes is transforming, as culture-independent methods facilitate rapid characterization of microbial diversity. It is difficult to assimilate this information into our understanding of marine microbe ecology and evolution, because their distributions, traits, and genomes are shaped by forces that are complex and dynamic. Here we incorporate diverse forces—physical, biogeochemical, ecological, and mutational—into a global ocean model to study selective pressures on a simple trait in a widely distributed lineage of picophytoplankton: the nitrogen use abilities of Synechococcus and Prochlorococcus cyanobacteria. Some Prochlorococcus ecotypes have lost the ability to use nitrate, whereas their close relatives, marine Synechococcus, typically retain it. We impose mutations for the loss of nitrogen use abilities in modeled picophytoplankton, and ask: in which parts of the ocean are mutants most disadvantaged by losing the ability to use nitrate, and in which parts are they least disadvantaged? Our model predicts that this selective disadvantage is smallest for picophytoplankton that live in tropical regions where Prochlorococcus are abundant in the real ocean. Conversely, the selective disadvantage of losing the ability to use nitrate is larger for modeled picophytoplankton that live at higher latitudes, where Synechococcus are abundant. In regions where we expect Prochlorococcus and Synechococcus populations to cycle seasonally in the real ocean, we find that model ecotypes with seasonal population dynamics similar to Prochlorococcus are less disadvantaged by losing the ability to use nitrate than model ecotypes with seasonal population dynamics similar to Synechococcus. The model predictions for the selective advantage associated with nitrate use are broadly consistent with the distribution of this ability among marine picocyanobacteria, and at finer scales, can provide insights into interactions between temporally varying ocean processes and selective pressures that may be difficult or impossible to study by other means. More generally, and perhaps more importantly, this study introduces an approach for testing hypotheses about the processes that underlie genetic variation among marine microbes, embedded in the dynamic physical, chemical, and biological forces that generate and shape this diversity.     

Identification of Multiple Gypsy LTR-Retrotransposon Lineages in Vertebrate Genomes

Gypsy LTR-retrotransposons have been identified in the genomes of many organisms, but only a small number of vertebrate examples have been reported to date. Here we show that members of this family are likely to be widespread in many vertebrate classes with the possible exceptions of mammals and birds. Phylogenetic analyses demonstrate that although there are several distinct lineages of vertebrate gypsy LTR-retrotransposons, the majority clusters into one monophyletic clade. Groups of fungal, plant, and insect elements were also observed, suggesting horizontal transfer between phyla may be infrequent. However, in contrast to this, there was little evidence to support sister relationships between elements derived from vertebrate and insect hosts. In fact, the majority of the vertebrate elements appeared to be most closely related to a group of gypsy LTR-retrotransposons present within fungi. This implies either that at least one horizontal transmission between these two phyla has occurred previously or that a gypsy LTR-retrotransposon lineage has been lost from insect taxa. Received: 22 December 1998 / Accepted: 6 April 1999     

Detection of Seven Major Evolutionary Lineages in Cyanobacteria Based on the 16S rRNA Gene Sequence Analysis with New Sequences of Five Marine Synechococcus Strains

Although molecular phylogenetic studies of cyanobacteria on the basis of the 16S rRNA gene sequence have been reported, the topologies were unstable, especially in the inner branchings. Our analysis of 16S rRNA gene phylogeny by the maximum-likelihood and neighbor-joining methods combined with rate homogeneous and heterogeneous models revealed seven major evolutionary lineages of the cyanobacteria, including prochlorophycean organisms. These seven lineages are always stable on any combination of these methods and models, fundamentally corresponding to phylogenetic relationships based on other genes, e.g., psbA, rbcL, rnpB, rpoC, and tufA. Moreover, although known genotypic and phenotypic characters sometimes appear paralleled in independent lineages, many characters are not contradictory within each group. Therefore we propose seven evolutionary groups as a working hypothesis for successive taxonomic reconstruction. New 16S rRNA sequences of five unicellular cyanobacterial strains, PCC 7001, PCC 7003, PCC 73109, PCC 7117, and PCC 7335 of Synechococcus sp., were determined in this study. Although all these strains have been assigned to ``marine clusters B and C,' they were separated into three lineages. This suggests that the organisms classified in the genus Synechococcus evolved diversely and should be reclassified in several independent taxonomic units. Moreover, Synechococcus strains and filamentous cyanobacteria make a monophyletic group supported by a comparatively high statistical confidence value (80 to 100%) in each of the two independent lineages; therefore, these monophylies probably reflect the convergent evolution of a multicellular organization. Received: 3 September 1998 / Accepted: 30 November 1998     

Molecular Phylogenetic Inference of the Woolly Mammoth Mammuthus primigenius, Based on Complete Sequences of Mitochondrial Cytochrome b and 12S Ribosomal RNA Genes

Complete sequences of cytochrome b (1,137 bases) and 12S ribosomal RNA (961 bases) genes in mitochondrial DNA were successfully determined from the woolly mammoth (Mammuthus primigenius), African elephant (Loxodonta africana), and Asian elephant (Elephas maximus). From these sequence data, phylogenetic relationships among three genera were examined. Molecular phylogenetic trees reconstructed by the neighbor-joining and the maximum parsimony methods provided an identical topology both for cytochrome b and 12S rRNA genes. These results support the ``Mammuthus-Loxodonta' clade, which is contrary to some previous morphological reports that Mammuthus is more closely related to Elephas than to Loxodonta. Received: 8 April 1997 / Accepted: 23 July 1997     

A Phylogenetic Assessment of the Eukaryotic Light-Harvesting Antenna Proteins, with Implications for Plastid Evolution

The light-harvesting complexes (LHCs) are a superfamily of chlorophyll-binding proteins present in all photosynthetic eukaryotes. The Lhc genes are nuclear-encoded, yet the pigment–protein complexes are localized to the thylakoid membrane and provide a marker to follow the evolutionary paths of plastids with different pigmentation. The LHCs are divided into the chlorophyll a/b-binding proteins of the green algae, euglenoids, and higher plants and the chlorophyll a/c-binding proteins of various algal taxa. This work examines the phylogenetic position of the LHCs from three additional taxa: the rhodophytes, the cryptophytes, and the chlorarachniophytes. Phylogenetic analysis of the LHC sequences provides strong statistical support for the clustering of the rhodophyte and cryptomonad LHC sequences within the chlorophyll a/c-binding protein lineage, which includes the fucoxanthin–chlorophyll proteins (FCP) of the heterokonts and the intrinsic peridinin–chlorophyll proteins (iPCP) of the dinoflagellates. These associations suggest that plastids from the heterokonts, haptophytes, cryptomonads, and the dinoflagellate, Amphidinium, evolved from a red algal-like ancestor. The Chlorarachnion LHC is part of the chlorophyll a/b-binding protein assemblage, consistent with pigmentation, providing further evidence that its plastid evolved from a green algal secondary endosymbiosis. The Chlorarachnion LHC sequences cluster with the green algal LHCs that are predominantly associated with photosystem II (LHCII). This suggests that the green algal endosymbiont that evolved into the Chlorarachnion plastid was acquired following the emergence of distinct LHCI and LHCII complexes. Received: 25 February 1998 / Accepted: 13 May 1998     

Phylogeny of Nematoda and Cephalorhyncha Derived from 18S rDNA

Phylogenetic relationships of nematodes, nematomorphs, kinorhynchs, priapulids, and some other major groups of invertebrates were studied by 18S rRNA gene sequencing. Kinorhynchs and priapulids form the monophyletic Cephalorhyncha clade that is the closest to the coelomate animals. When phylogenetic trees were generated by different methods, the position of nematomorphs appeared to be unstable. Inclusion of Enoplus brevis, a representative of a slowly evolving nematode lineage, in the set of analyzed species refutes the tree patterns, previously derived from molecular data, where the nematodes appear as a basal bilateral lineage. The nematodes seem to be closer to the coelomate animals than was speculated earlier. According to the results obtained, nematodes, nematomorphs, tardigrades, arthropods, and cephalorhynchs are a paraphyletic association of closely related taxa. Received: 1 December 1997 / Accepted: 9 April 1998     

Phylogenetic Depth of the Bacterial Genera Aquifex and Thermotoga Inferred from Analysis of Ribosomal Protein, Elongation Factor, and RNA Polymerase Subunit Sequences

The phylogenetic placement of the Aquifex and Thermotoga lineages has been inferred from (i) the concatenated ribosomal proteins S10, L3, L4, L23, L2, S19, L22, and S3 encoded in the S10 operon (833 aa positions); (ii) the joint sequences of the elongation factors Tu(1α) and G(2) coded by the str operon tuf and fus genes (733 aa positions); and (iii) the joint RNA polymerase β- and β′-type subunits encoded in the rpoBC operon (1130 aa positions). Phylogenies of r-protein and EF sequences support with moderate (r-proteins) to high statistical confidence (EFs) the placement of the two hyperthermophiles at the base of the bacterial clade in agreement with phylogenies of rRNA sequences. In the more robust EF-based phylogenies, the branching of Aquifex and Thermotoga below the successive bacterial lineages is given at bootstrap proportions of 82% (maximum likelihood; ML) and 85% (maximum parsimony; MP), in contrast to the trees inferred from the separate EF-Tu(1α) and EF-G(2) data sets, which lack both resolution and statistical robustness. In the EF analysis MP outperforms ML in discriminating (at the 0.05 level) trees having A. pyrophilus and T. maritima as the most basal lineages from competing alternatives that have (i) mesophiles, or the Thermus genus, as the deepest bacterial radiation and (ii) a monophyletic A. pyrophilus–T. maritima cluster situated at the base of the bacterial clade. RNAP-based phylogenies are equivocal with respect to the Aquifex and Thermotoga placements. The two hyperthermophiles fall basal to all other bacterial phyla when potential artifacts contributed by the compositionally biased and fast-evolving Mycoplasma genitalium and Mycoplasma pneumoniae sequences are eschewed. However, the branching order of the phyla is tenuously supported in ML trees inferred by the exhaustive search method and is unresolved in ML trees inferred by the quartet puzzling algorithm. A rooting of the RNA polymerase-subunit tree at the mycoplasma level seen in both the MP trees and the ML trees reconstructed with suboptimal amino acid substitution models is not supported by the EF-based phylogenies which robustly affiliate mycoplasmas with low-G+C gram-positives and, most probably, reflects a ``long branch attraction' artifact. Received: 22 September 1999 / Accepted: 11 January 2000     

Inhibition by ultraviolet and photosynthetically available radiation lowers model estimates of depth‐integrated picophytoplankton photosynthesis: global predictions for Prochlorococcus and Synechococcus

Phytoplankton photosynthesis is often inhibited by ultraviolet (UV) and intense photosynthetically available radiation (PAR), but the effects on ocean productivity have received little consideration aside from polar areas subject to periodic enhanced UV‐B due to depletion of stratospheric ozone. A more comprehensive assessment is important for understanding the contribution of phytoplankton production to the global carbon budget, present and future. Here, we consider responses in the temperate and tropical mid‐ocean regions typically dominated by picophytoplankton including the prokaryotic lineages, Prochlorococcus and Synechococcus. Spectral models of photosynthetic response for each lineage were constructed using model strains cultured at different growth irradiances and temperatures. In the model, inhibition becomes more severe once exposure exceeds a threshold (E_max) related to repair capacity. Model parameters are presented for Prochlorococcus adding to those previously presented for Synechococcus. The models were applied to estimate midday, water column photosynthesis based on an atmospheric model of spectral radiation, satellite‐derived spectral water transparency and temperature. Based on a global survey of inhibitory exposure severity, a full‐latitude section of the mid‐Pacific and near‐equatorial region of the east Pacific were identified as representative regions for prediction of responses over the entire water column. Comparing predictions integrated over the water column including versus excluding inhibition, production was 7–28% lower due to inhibition depending on strain and site conditions. Inhibition was consistently greater for Prochlorococcus compared to two strains of Synechococcus. Considering only the surface mixed layer, production was inhibited 7–73%. On average, including inhibition lowered estimates of midday productivity around 20% for the modeled region of the Pacific with UV accounting for two‐thirds of the reduction. In contrast, most other productivity models either ignore inhibition or only include PAR inhibition. Incorporation of E_max model responses into an existing spectral model of depth‐integrated, daily production will enable efficient global predictions of picophytoplankton productivity including inhibition.     

Picophytoplankton biomass, community structure and productivity in the Great Astrolabe Lagoon, Fiji

 Phytoplankton biomass, community structure and productivity of the Great Astrolabe lagoon and surrounding ocean were studied using measurements of chlorophyll concentration and carbon uptake. The contribution of picophytoplankton to biomass, productivity and community structure was estimated by size fractionation, ¹⁴C-incubation and flow cytometry analysis. Picoplankton red fluorescence was demonstrated to be a proxy for chlorophyll <3 μm. Consequently, the percentage contribution to chl a<3 μm from each picoplankton group could be calculated using regression estimated values of ψ _i (fg chl a per unit of red fluorescence). In the lagoon, average chlorophyll concentration was 0.8 mg m^-3 with 45% of phytoplankton <3 μm. Primary production reached 1.3 g C m^-2 day^-1 with 53% due to phytoplankton <3 μm. Synechococcus was the most abundant group at all stations, followed by Prochlorococcus and picoeukaryotes. At all stations, Prochlorococcus represented less than 4% of the chl a <3 μm, Synechococcus between 85 and 95%, and Picoeukaryotes between 5 and 10%. In the upper 40 m of surrounding oceanic waters, phytoplankton biomass was dominated by the >3 μm size fraction. In deeper water, the <1 μm size fraction dominated. Prochlorococcus was the most abundant picoplankton group and their contributions to the chlorophyll a<3 μm were close to that of the picoeukaryotes (50% each). Accepted: 27 May 1999     

Comparative Genomics of Mitochondrial DNA in Members of the Drosophila melanogaster Subgroup

In this study, a comparative genomics approach is employed to investigate the forces that shape evolutionary change in the mitochondrial DNA (mtDNA) of members of the Drosophila melanogaster subgroup. This approach facilitates differentiation of the patterns of variation resulting from processes acting at a higher level from those acting on a single gene. The mitochondrial genomes of three isofemale lines of D. simulans (siI, -II, and -III), two of D. melanogaster (Oregon R and a line from Zimbabwe), and D. mauritiana (maI and -II), and one of D. sechellia were sequenced and compared with that derived from D. yakuba. Data presented here indicate that at least three broad mechanisms shape the evolutionary dynamics of mtDNA in these taxa. The first set of mechanisms is intrinsic to the molecule. Dominant processes may be interpreted as selection for an increased rate of replication of the mtDNA molecule, biases in DNA repair, and differences in the pattern of nucleotide substitution among strands. In the genes encoded on the major strand (62% of the coding DNA) changes to or from C predominate, whereas on the minor changes to or from G predominate. The second set of mechanisms affects distinct lineages. There are evolutionary rate differences among lineages, possibly owing to population demographic changes or changes in mutational biases. This is supported by the heterogeneity found in synonymous, nonsynonymous, and silent substitutions. The third set of mechanisms differentially affects distinct genes. A maximum-likelihood sliding-window analysis detected four disjunct regions that have a significantly different nucleotide substitution process from that derived from the complete sequence. These data show the potential for comparative genomics to tease apart subtle forces that shape the evolution of DNA. Received: 30 July 1999 / Accepted: 16 March 2000     

Ocean time-series reveals recurring seasonal patterns of virioplankton dynamics in the northwestern Sargasso Sea

There are an estimated 10³⁰ virioplankton in the world oceans, the majority of which are phages (viruses that infect bacteria). Marine phages encompass enormous genetic diversity, affect biogeochemical cycling of elements, and partially control aspects of prokaryotic production and diversity. Despite their importance, there is a paucity of data describing virioplankton distributions over time and depth in oceanic systems. A decade of high-resolution time-series data collected from the upper 300 m in the northwestern Sargasso Sea revealed recurring temporal and vertical patterns of virioplankton abundance in unprecedented detail. An annual virioplankton maximum developed between 60 and 100 m during periods of summer stratification and eroded during winter convective mixing. The timing and vertical positioning of this seasonal pattern was related to variability in water column stability and the dynamics of specific picophytoplankton and heterotrophic bacterioplankton lineages. Between 60 and 100 m, virioplankton abundance was negatively correlated to the dominant heterotrophic bacterioplankton lineage SAR11, as well as the less abundant picophytoplankton, Synechococcus. In contrast, virioplankton abundance was positively correlated to the dominant picophytoplankton lineage Prochlorococcus, and the less abundant alpha-proteobacteria, Rhodobacteraceae. Seasonally, virioplankton abundances were highly synchronous with Prochlorococcus distributions and the virioplankton to Prochlorococcus ratio remained remarkably constant during periods of water column stratification. The data suggest that a significant fraction of viruses in the mid-euphotic zone of the subtropical gyres may be cyanophages and patterns in their abundance are largely determined by Prochlorococcus dynamics in response to water column stability. This high-resolution, decadal survey of virioplankton abundance provides insight into the possible controls of virioplankton dynamics in the open ocean.     

Co-occurring Synechococcus ecotypes occupy four major oceanic regimes defined by temperature,macronutrients and iron

Marine picocyanobacteria, comprised of the genera Synechococcus and Prochlorococcus, are the most abundant and widespread primary producers in the ocean. More than 20 genetically distinct clades of marine Synechococcus have been identified, but their physiology and biogeography are not as thoroughly characterized as those of Prochlorococcus. Using clade-specific qPCR primers, we measured the abundance of 10 Synechococcus clades at 92 locations in surface waters of the Atlantic and Pacific Oceans. We found that Synechococcus partition the ocean into four distinct regimes distinguished by temperature, macronutrients and iron availability. Clades I and IV were prevalent in colder, mesotrophic waters; clades II, III and X dominated in the warm, oligotrophic open ocean; clades CRD1 and CRD2 were restricted to sites with low iron availability; and clades XV and XVI were only found in transitional waters at the edges of the other biomes. Overall, clade II was the most ubiquitous clade investigated and was the dominant clade in the largest biome, the oligotrophic open ocean. Co-occurring clades that occupy the same regime belong to distinct evolutionary lineages within Synechococcus, indicating that multiple ecotypes have evolved independently to occupy similar niches and represent examples of parallel evolution. We speculate that parallel evolution of ecotypes may be a common feature of diverse marine microbial communities that contributes to functional redundancy and the potential for resiliency.     

Evolution of the mitochondrial DNA control region in the mbuna (Cichlidae) species flock of lake Malawi, East Africa

Considerable controversy has surrounded the application of mitochondrial DNA data to reconstruction of evolutionary relationships among the endemic cichlids of Lake Malawi. Central to this debate has been the issue of whether lineage sorting is complete, and thus whether these data actually reflect species phylogeny, or simply gene genealogy. Review of all mtDNA control region sequences available for members of one monophyletic subset of this species flock, the Malawi rockfishes, or mbuna, strongly indicates that lineage sorting is incomplete: Character-based analyses of these sequences reconstruct gene, not species, interrelationships. Analysis of the pattern of nucleotide substitutions differentiating these mtDNA alleles suggests that pyrimidine residues undergo transition substitutions more often than do purines. Estimation of the magnitude of derived sequence differentiation in light of the reconstructed gene genealogy suggests that the mbuna may be of considerably more recent vintage than previous molecular characterizations have indicated. Received: 6 April 1996 / Accepted: 3 March 1997