In this study, we aimed to investigate, through high-resolution metagenomics and metatranscriptomics, the composition and the trajectories of microbial communities originating from a natural sample, fed exclusively with methane, over 14 weeks of laboratory incubation. This study builds on our prior data, suggesting that multiple functional guilds feed on methane, likely through guild-to-guild carbon transfer, and potentially through intraguild and intraspecies interactions. We observed that, under two simulated dioxygen partial pressures—low versus high—community trajectories were different, with considerable variability among the replicates. In all microcosms, four major functional guilds were prominently present, representing Methylococcaceae (the true methanotrophs), Methylophilaceae (the nonmethanotrophic methylotrophs), Burkholderiales, and Bacteroidetes. Additional functional guilds were detected in multiple samples, such as members of Opitutae, as well as the predatory species, suggesting additional complexity for methane-oxidizing communities. Metatranscriptomic analysis suggested simultaneous expression of the two alternative types of methanol dehydrogenases in both Methylococcaceae and Methylophilaceae, while high expression of the oxidative/nitrosative stress response genes suggested competition for dioxygen among the community members. The transcriptomic analysis further suggested that Burkholderiales likely feed on acetate that is produced by Methylococcaceae under hypoxic conditions, while Bacteroidetes likely feed on biopolymers produced by both Methylococcaceae and Methylophilaceae.

Changes in the sequence of an organism’s genome, i.e. mutations, are the raw material of evolution ¹ . The frequency and location of mutations can be constrained by specific molecular mechanisms, such as Diversity-generating retroelements (DGRs) ^2–4 . DGRs introduce mutations in specific target genes, and were characterized from several cultivated bacteria and bacteriophages ² . Whilst a larger diversity of DGR loci has been identified in genomic data from environmental samples, i.e. metagenomes, the ecological role of these DGRs and their associated evolutionary drivers remain poorly understood ^5–7 . Here we built and analyzed an extensive dataset of >30,000 metagenome-derived DGRs, and determine that DGRs have a single evolutionary origin and a universal bias towards adenine mutations. We further identified six major lineages of DGRs, each associated with a specific ecological niche defined as a genome type, i.e. whether the DGR is encoded on a viral or cellular genome, a limited set of taxa and environments, and a distinct type of target. Finally, we leverage read mapping and metagenomic time series to demonstrate that DGRs are consistently and broadly active, and responsible for >10% of all amino acid changes in some organisms at a conservative estimate. Overall, these results highlight the strong constraints under which DGRs diversify and expand, and elucidate several distinct roles these elements play in natural communities and in shaping microbial community structure and function in our environment.