Genome signatures in metagenomic datasets

Genome signatures are used to identify and cluster sequences de novo from an acid biofilm microbial community metagenomic dataset, revealing information about the low-abundance community members.

BackgroundMetal sulfide mineral dissolution during bioleaching and acid mine drainage (AMD) formation creates an environment that is inhospitable to most life. Despite dominance by a small number of bacteria, AMD microbial biofilm communities contain a notable variety of coexisting and closely related Euryarchaea, most of which have defied cultivation efforts. For this reason, we used metagenomics to analyze variation in gene content that may contribute to niche differentiation among co-occurring AMD archaea. Our analyses targeted members of the Thermoplasmatales and related archaea. These results greatly expand genomic information available for this archaeal order.ResultsWe reconstructed near-complete genomes for uncultivated, relatively low abundance organisms A-, E-, and Gplasma, members of Thermoplasmatales order, and for a novel organism, Iplasma. Genomic analyses of these organisms, as well as Ferroplasma type I and II, reveal that all are facultative aerobic heterotrophs with the ability to use many of the same carbon substrates, including methanol. Most of the genomes share genes for toxic metal resistance and surface-layer production. Only Aplasma and Eplasma have a full suite of flagellar genes whereas all but the Ferroplasma spp. have genes for pili production. Cryogenic-electron microscopy (cryo-EM) and tomography (cryo-ET) strengthen these metagenomics-based ultrastructural predictions. Notably, only Aplasma, Gplasma and the Ferroplasma spp. have predicted iron oxidation genes and Eplasma and Iplasma lack most genes for cobalamin, valine, (iso)leucine and histidine synthesis.ConclusionThe Thermoplasmatales AMD archaea share a large number of metabolic capabilities. All of the uncultivated organisms studied here (A-, E-, G-, and Iplasma) are metabolically very similar to characterized Ferroplasma spp., differentiating themselves mainly in their genetic capabilities for biosynthesis, motility, and possibly iron oxidation. These results indicate that subtle, but important genomic differences, coupled with unknown differences in gene expression, distinguish these organisms enough to allow for co-existence. Overall this study reveals shared features of organisms from the Thermoplasmatales lineage and provides new insights into the functioning of AMD communities.

The assimilation of organic nutrients by autotrophs, a form of mixotrophy, has been demonstrated in the globally abundant marine picocyanobacterial genera Prochlorococcus and Synechococcus. However, the range of compounds used and the distribution of organic compound uptake genes within picocyanobacteria are unknown. Here we analyze genomic and metagenomic data from around the world to determine the extent and distribution of mixotrophy in these phototrophs. Analysis of 49 Prochlorococcus and 18 Synechococcus isolate genomes reveals that all have the transporters necessary to take up amino acids, peptides and sugars. However, the number and type of transporters and associated catabolic genes differ between different phylogenetic groups, with lowlight IV Prochlorococcus, and 5.1B, 5.2 and 5.3 Synechococcus strains having the largest number. Metagenomic data from 68 stations from the Tara Oceans expedition indicate that the genetic potential for mixotrophy in picocyanobacteria is globally distributed and differs between clades. Phylogenetic analyses indicate gradual organic nutrient transporter gene loss from the low-light IV to the high-light II Prochlorococcus. The phylogenetic differences in genetic capacity for mixotrophy, combined with the ubiquity of picocyanobacterial organic compound uptake genes suggests that mixotrophy has a more central role in picocyanobacterial ecology than was previously thought.