Marine planktonic cyanobacteria capable of fixing molecular nitrogen (termed 'diazotrophs') are key in biogeochemical cycling, and the nitrogen fixed is one of the major external sources of nitrogen to the open ocean. Candidatus Atelocyanobacterium thalassa (UCYN-A) is a diazotrophic cyanobacterium known for its widespread geographic distribution in tropical and subtropical oligotrophic oceans, unusually reduced genome and symbiosis with a single-celled prymnesiophyte alga. Recently a novel strain of this organism was also detected in coastal waters sampled from the Scripps Institute of Oceanography pier. We analyzed the metagenome of this UCYN-A2 population by concentrating cells by flow cytometry. Phylogenomic analysis provided strong bootstrap support for the monophyly of UCYN-A (here called UCYN-A1) and UCYN-A2 within the marine Crocosphaera sp. and Cyanothece sp. clade. UCYN-A2 shares 1159 of the 1200 UCYN-A1 protein-coding genes (96.6%) with high synteny, yet the average amino-acid sequence identity between these orthologs is only 86%. UCYN-A2 lacks the same major pathways and proteins that are absent in UCYN-A1, suggesting that both strains can be grouped at the same functional and ecological level. Our results suggest that UCYN-A1 and UCYN-A2 had a common ancestor and diverged after genome reduction. These two variants may reflect adaptation of the host to different niches, which could be coastal and open ocean habitats.

SUPPLEMENTARY INFORMATION is available at Bioinformatics online.

Crocosphaera watsonii, a unicellular nitrogen-fixing cyanobacterium found in oligotrophic oceans, is important in marine carbon and nitrogen cycles. Isolates of C. watsonii can be separated into at least two phenotypes with environmentally important differences, indicating possibly distinct ecological roles and niches. To better understand the evolutionary history and variation in metabolic capabilities among strains and phenotypes, this study compared the genomes of six C. watsonii strains, three from each phenotypic group, which had been isolated over several decades from multiple ocean basins. While a substantial portion of each genome was nearly identical to sequences in the other strains, a few regions were identified as specific to each strain and phenotype, some of which help explain observed phenotypic features. Overall, the small-cell type strains had smaller genomes and a relative loss of genetic capabilities, while the large-cell type strains were characterized by larger genomes, some genetic redundancy, and potentially increased adaptations to iron and phosphorus limitation. As such, strains with shared phenotypes were evolutionarily more closely related than those with the opposite phenotype, regardless of isolation location or date. Unexpectedly, the genome of the type-strain for the species, C. watsonii WH8501, was quite unusual even among strains with a shared phenotype, indicating it may not be an ideal representative of the species. The genome sequences and analyses reported in this study will be important for future investigations of the proposed differences in adaptation of the two phenotypes to nutrient limitation, and to identify phenotype-specific distributions in natural Crocosphaera populations.

The Cytochrome C Oxidase subunit I gene (“COI”) is the de facto standard for animal DNA barcoding. Organism identification based on COI requires an accurate and extensive annotated database of COI sequences. Such a database can also be of value in reconstructing evolutionary history and in diversity studies. Two COI databases are currently available: BOLD and Midori. BOLD’s submissions conform to stringent sequence and metadata requirements; BOLD is specific to COI but makes no attempt to be comprehensive. Midori, derived from GenBank, has more sequences but less stringent standards than BOLD, resulting in higher error rates. To address the need for a comprehensive and accurate COI database, we adapted the ARBitrator algorithm, which classifies based only on sequence properties and has successfully auto-curated bacterial genes mined from GenBank. The adapted algorithm, which we call CO-ARBitrator, built a database of over a million metazoan COI sequences. Sensitivity and specificity are significantly higher than Midori. Specificity is comparable to what BOLD achieves with data quality prerequisites. Results and software are publicly available.