Multicellular magnetotactic bacteria are genetically heterogeneous consortia with metabolically differentiated cells.

Consortia of multicellular magnetotactic bacteria (MMB) are currently the only known example of bacteria without a unicellular stage in their life cycle. Because of their recalcitrance to cultivation, most previous studies of MMB have been limited to microscopic observations. To study the biology of these unique organisms in more detail, we use multiple culture-independent approaches to analyze the genomics and physiology of MMB consortia at single-cell resolution. We separately sequenced the metagenomes of 22 individual MMB consortia, representing 8 new species, and quantified the genetic diversity within each MMB consortium. This revealed that, counter to conventional views, cells within MMB consortia are not clonal. Single consortia metagenomes were then used to reconstruct the species-specific metabolic potential and infer the physiological capabilities of MMB. To validate genomic predictions, we performed stable isotope probing (SIP) experiments and interrogated MMB consortia using fluorescence in situ hybridization (FISH) combined with nanoscale secondary ion mass spectrometry (NanoSIMS). By coupling FISH with bioorthogonal noncanonical amino acid tagging (BONCAT), we explored their in situ activity as well as variation of protein synthesis within cells. We demonstrate that MMB consortia are mixotrophic sulfate reducers and that they exhibit metabolic differentiation between individual cells, suggesting that MMB consortia are more complex than previously thought. These findings expand our understanding of MMB diversity, ecology, genomics, and physiology, as well as offer insights into the mechanisms underpinning the multicellular nature of their unique lifestyle.

Spatial co-transcriptomics reveals discrete stages of the arbuscular mycorrhizal symbiosis.

The symbiotic interaction of plants with arbuscular mycorrhizal (AM) fungi is ancient and widespread. Plants provide AM fungi with carbon in exchange for nutrients and water, making this interaction a prime target for crop improvement. However, plant-fungal interactions are restricted to a small subset of root cells, precluding the application of most conventional functional genomic techniques to study the molecular bases of these interactions. Here we used single-nucleus and spatial RNA sequencing to explore both Medicago truncatula and Rhizophagus irregularis transcriptomes in AM symbiosis at cellular and spatial resolution. Integrated, spatially registered single-cell maps revealed infected and uninfected plant root cell types. We observed that cortex cells exhibit distinct transcriptome profiles during different stages of colonization by AM fungi, indicating dynamic interplay between both organisms during establishment of the cellular interface enabling successful symbiosis. Our study provides insight into a symbiotic relationship of major agricultural and environmental importance and demonstrates a paradigm combining single-cell and spatial transcriptomics for the analysis of complex organismal interactions.

scMicrobe PTA: Near Complete Genomes from Single Bacterial Cells.

Microbial genomes produced by single-cell amplification are largely incomplete. Here, we show that primary template amplification (PTA), a novel single-cell amplification technique, generated nearly complete genomes from three bacterial isolate species. Furthermore, taxonomically diverse genomes recovered from aquatic and soil microbiomes using PTA had a median completeness of 81%, whereas genomes from standard amplification approaches were usually <30% complete. PTA-derived genomes also included more associated viruses and biosynthetic gene clusters.

Multicellular magnetotactic bacterial consortia are metabolically differentiated and not clonal.

Consortia of multicellular magnetotactic bacteria (MMB) are currently the only known example of bacteria without a unicellular stage in their life cycle. Because of their recalcitrance to cultivation, most previous studies of MMB have been limited to microscopic observations. To study the biology of these unique organisms in more detail, we use multiple culture-independent approaches to analyze the genomics and physiology of MMB consortia at single cell resolution. We separately sequenced the metagenomes of 22 individual MMB consortia, representing eight new species, and quantified the genetic diversity within each MMB consortium. This revealed that, counter to conventional views, cells within MMB consortia are not clonal. Single consortia metagenomes were then used to reconstruct the species-specific metabolic potential and infer the physiological capabilities of MMB. To validate genomic predictions, we performed stable isotope probing (SIP) experiments and interrogated MMB consortia using fluorescence in situ hybridization (FISH) combined with nano-scale secondary ion mass spectrometry (NanoSIMS). By coupling FISH with bioorthogonal non-canonical amino acid tagging (BONCAT) we explored their in situ activity as well as variation of protein synthesis within cells. We demonstrate that MMB consortia are mixotrophic sulfate reducers and that they exhibit metabolic differentiation between individual cells, suggesting that MMB consortia are more complex than previously thought. These findings expand our understanding of MMB diversity, ecology, genomics, and physiology, as well as offer insights into the mechanisms underpinning the multicellular nature of their unique lifestyle.

Physiological potential and evolutionary trajectories of syntrophic sulfate-reducing bacterial partners of anaerobic methanotrophic archaea.

Sulfate-coupled anaerobic oxidation of methane (AOM) is performed by multicellular consortia of anaerobic methanotrophic archaea (ANME) in obligate syntrophic partnership with sulfate-reducing bacteria (SRB). Diverse ANME and SRB clades co-associate but the physiological basis for their adaptation and diversification is not well understood. In this work, we used comparative metagenomics and phylogenetics to investigate the metabolic adaptation among the 4 main syntrophic SRB clades (HotSeep-1, Seep-SRB2, Seep-SRB1a, and Seep-SRB1g) and identified features associated with their syntrophic lifestyle that distinguish them from their non-syntrophic evolutionary neighbors in the phylum Desulfobacterota. We show that the protein complexes involved in direct interspecies electron transfer (DIET) from ANME to the SRB outer membrane are conserved between the syntrophic lineages. In contrast, the proteins involved in electron transfer within the SRB inner membrane differ between clades, indicative of convergent evolution in the adaptation to a syntrophic lifestyle. Our analysis suggests that in most cases, this adaptation likely occurred after the acquisition of the DIET complexes in an ancestral clade and involve horizontal gene transfers within pathways for electron transfer (CbcBA) and biofilm formation (Pel). We also provide evidence for unique adaptations within syntrophic SRB clades, which vary depending on the archaeal partner. Among the most widespread syntrophic SRB, Seep-SRB1a, subclades that specifically partner ANME-2a are missing the cobalamin synthesis pathway, suggestive of nutritional dependency on its partner, while closely related Seep-SRB1a partners of ANME-2c lack nutritional auxotrophies. Our work provides insight into the features associated with DIET-based syntrophy and the adaptation of SRB towards it.

BONCAT-FACS-Seq reveals the active fraction of a biocrust community undergoing a wet-up event.

Determining which microorganisms are active within soil communities remains a major technical endeavor in microbial ecology research. One promising method to accomplish this is coupling bioorthogonal non-canonical amino acid tagging (BONCAT) with fluorescence activated cell sorting (FACS) which sorts cells based on whether or not they are producing new proteins. Combined with shotgun metagenomic sequencing (Seq), we apply this method to profile the diversity and potential functional capabilities of both active and inactive microorganisms in a biocrust community after being resuscitated by a simulated rain event. We find that BONCAT-FACS-Seq is capable of discerning the pools of active and inactive microorganisms, especially within hours of applying the BONCAT probe. The active and inactive components of the biocrust community differed in species richness and composition at both 4 and 21 h after the wetting event. The active fraction of the biocrust community is marked by taxa commonly observed in other biocrust communities, many of which play important roles in species interactions and nutrient transformations. Among these, 11 families within the Firmicutes are enriched in the active fraction, supporting previous reports indicating that the Firmicutes are key early responders to biocrust wetting. We highlight the apparent inactivity of many Actinobacteria and Proteobacteria through 21 h after wetting, and note that members of the Chitinophagaceae, enriched in the active fraction, may play important ecological roles following wetting. Based on the enrichment of COGs in the active fraction, predation by phage and other bacterial members, as well as scavenging and recycling of labile nutrients, appear to be important ecological processes soon after wetting. To our knowledge, this is the first time BONCAT-FACS-Seq has been applied to biocrust samples, and therefore we discuss the potential advantages and shortcomings of coupling metagenomics to BONCAT to intact soil communities such as biocrust. In all, by pairing BONCAT-FACS and metagenomics, we are capable of highlighting the taxa and potential functions that typifies the microbes actively responding to a rain event.

Hybridization breaks species barriers in long-term coevolution of a cyanobacterial population.

Bacterial species often undergo rampant recombination yet maintain cohesive genomic identity. Ecological differences can generate recombination barriers between species and sustain genomic clusters in the short term. But can these forces prevent genomic mixing during long-term coevolution? Cyanobacteria in Yellowstone hot springs comprise several diverse species that have coevolved for hundreds of thousands of years, providing a rare natural experiment. By analyzing more than 300 single-cell genomes, we show that despite each species forming a distinct genomic cluster, much of the diversity within species is the result of hybridization driven by selection, which has mixed their ancestral genotypes. This widespread mixing is contrary to the prevailing view that ecological barriers can maintain cohesive bacterial species and highlights the importance of hybridization as a source of genomic diversity.

A centimeter-long bacterium with DNA contained in metabolically active, membrane-bound organelles.

Cells of most bacterial species are around 2 micrometers in length, with some of the largest specimens reaching 750 micrometers. Using fluorescence, x-ray, and electron microscopy in conjunction with genome sequencing, we characterized Candidatus (Ca.) Thiomargarita magnifica, a bacterium that has an average cell length greater than 9000 micrometers and is visible to the naked eye. These cells grow orders of magnitude over theoretical limits for bacterial cell size, display unprecedented polyploidy of more than half a million copies of a very large genome, and undergo a dimorphic life cycle with asymmetric segregation of chromosomes into daughter cells. These features, along with compartmentalization of genomic material and ribosomes in translationally active organelles bound by bioenergetic membranes, indicate gain of complexity in the Thiomargarita lineage and challenge traditional concepts of bacterial cells.

Community Structure and Microbial Associations in Sediment-Free Methanotrophic Enrichment Cultures from a Marine Methane Seep.

Syntrophic consortia of anaerobic methanotrophic archaea (ANME) and sulfate-reducing bacteria (SRB) consume large amounts of methane and serve as the foundational microorganisms in marine methane seeps. Despite their importance in the carbon cycle, research on the physiology of ANME-SRB consortia has been hampered by the slow growth and complex physicochemical environment the consortia inhabit. Here, we report successful sediment-free enrichment of ANME-SRB consortia from deep-sea methane seep sediments in the Santa Monica Basin, California. Anoxic Percoll density gradients and size-selective filtration were used to separate ANME-SRB consortia from sediment particles and single cells to accelerate the cultivation process. Over a 3-year period, a subset of the sediment-associated ANME and SRB lineages, predominantly comprised of ANME-2a/2b ("Candidatus Methanocomedenaceae") and their syntrophic bacterial partners, SEEP-SRB1/2, adapted and grew under defined laboratory conditions. Metagenome-assembled genomes from several enrichments revealed that ANME-2a, SEEP-SRB1, and Methanococcoides in different enrichments from the same inoculum represented distinct species, whereas other coenriched microorganisms were closely related at the species level. This suggests that ANME, SRB, and Methanococcoides are more genetically diverse than other members in methane seeps. Flow cytometry sorting and sequencing of cell aggregates revealed that Methanococcoides, Anaerolineales, and SEEP-SRB1 were overrepresented in multiple ANME-2a cell aggregates relative to the bulk metagenomes, suggesting they were physically associated and possibly interacting. Overall, this study represents a successful case of selective cultivation of anaerobic slow-growing microorganisms from sediments based on their physical characteristics, introducing new opportunities for detailed genomic, physiological, biochemical, and ecological analyses. IMPORTANCE Biological anaerobic oxidation of methane (AOM) coupled with sulfate reduction represents a large methane sink in global ocean sediments. Methane consumption is carried out by syntrophic archaeal-bacterial consortia and fuels a unique ecosystem, yet the interactions in these slow-growing syntrophic consortia and with other associated community members remain poorly understood. The significance of this study is the establishment of sediment-free enrichment cultures of anaerobic methanotrophic archaea and sulfate-reducing bacteria performing AOM with sulfate using selective cultivation approaches based on size, density, and metabolism. By reconstructing microbial genomes and analyzing community composition of the enrichment cultures and cell aggregates, we shed light on the diversity of microorganisms physically associated with AOM consortia beyond the core syntrophic partners. These enrichment cultures offer simplified model systems to extend our understanding of the diversity of microbial interactions within marine methane seeps.

Subsurface hydrocarbon degradation strategies in low- and high-sulfate coal seam communities identified with activity-based metagenomics.

Environmentally relevant metagenomes and BONCAT-FACS derived translationally active metagenomes from Powder River Basin coal seams were investigated to elucidate potential genes and functional groups involved in hydrocarbon degradation to methane in coal seams with high- and low-sulfate levels. An advanced subsurface environmental sampler allowed the establishment of coal-associated microbial communities under in situ conditions for metagenomic analyses from environmental and translationally active populations. Metagenomic sequencing demonstrated that biosurfactants, aerobic dioxygenases, and anaerobic phenol degradation pathways were present in active populations across the sampled coal seams. In particular, results suggested the importance of anaerobic degradation pathways under high-sulfate conditions with an emphasis on fumarate addition. Under low-sulfate conditions, a mixture of both aerobic and anaerobic pathways was observed but with a predominance of aerobic dioxygenases. The putative low-molecular-weight biosurfactant, lichysein, appeared to play a more important role compared to rhamnolipids. The methods used in this study-subsurface environmental samplers in combination with metagenomic sequencing of both total and translationally active metagenomes-offer a deeper and environmentally relevant perspective on community genetic potential from coal seams poised at different redox conditions broadening the understanding of degradation strategies for subsurface carbon.

Comparative genomics reveals electron transfer and syntrophic mechanisms differentiating methanotrophic and methanogenic archaea.

The anaerobic oxidation of methane coupled to sulfate reduction is a microbially mediated process requiring a syntrophic partnership between anaerobic methanotrophic (ANME) archaea and sulfate-reducing bacteria (SRB). Based on genome taxonomy, ANME lineages are polyphyletic within the phylum Halobacterota, none of which have been isolated in pure culture. Here, we reconstruct 28 ANME genomes from environmental metagenomes and flow sorted syntrophic consortia. Together with a reanalysis of previously published datasets, these genomes enable a comparative analysis of all marine ANME clades. We review the genomic features that separate ANME from their methanogenic relatives and identify what differentiates ANME clades. Large multiheme cytochromes and bioenergetic complexes predicted to be involved in novel electron bifurcation reactions are well distributed and conserved in the ANME archaea, while significant variations in the anabolic C1 pathways exists between clades. Our analysis raises the possibility that methylotrophic methanogenesis may have evolved from a methanotrophic ancestor.

Dissecting the dominant hot spring microbial populations based on community-wide sampling at single-cell genomic resolution.

With advances in DNA sequencing and miniaturized molecular biology workflows, rapid and affordable sequencing of single-cell genomes has become a reality. Compared to 16S rRNA gene surveys and shotgun metagenomics, large-scale application of single-cell genomics to whole microbial communities provides an integrated snapshot of community composition and function, directly links mobile elements to their hosts, and enables analysis of population heterogeneity of the dominant community members. To that end, we sequenced nearly 500 single-cell genomes from a low diversity hot spring sediment sample from Dewar Creek, British Columbia, and compared this approach to 16S rRNA gene amplicon and shotgun metagenomics applied to the same sample. We found that the broad taxonomic profiles were similar across the three sequencing approaches, though several lineages were missing from the 16S rRNA gene amplicon dataset, likely the result of primer mismatches. At the functional level, we detected a large array of mobile genetic elements present in the single-cell genomes but absent from the corresponding same species metagenome-assembled genomes. Moreover, we performed a single-cell population genomic analysis of the three most abundant community members, revealing differences in population structure based on mutation and recombination profiles. While the average pairwise nucleotide identities were similar across the dominant species-level lineages, we observed differences in the extent of recombination between these dominant populations. Most intriguingly, the creek's Hydrogenobacter sp. population appeared to be so recombinogenic that it more closely resembled a sexual species than a clonally evolving microbe. Together, this work demonstrates that a randomized single-cell approach can be useful for the exploration of previously uncultivated microbes from community composition to population structure.

Activity-based, genome-resolved metagenomics uncovers key populations and pathways involved in subsurface conversions of coal to methane.

Microbial metabolisms and interactions that facilitate subsurface conversions of recalcitrant carbon to methane are poorly understood. We deployed an in situ enrichment device in a subsurface coal seam in the Powder River Basin (PRB), USA, and used BONCAT-FACS-Metagenomics to identify translationally active populations involved in methane generation from a variety of coal-derived aromatic hydrocarbons. From the active fraction, high-quality metagenome-assembled genomes (MAGs) were recovered for the acetoclastic methanogen, Methanothrix paradoxum, and a novel member of the Chlorobi with the potential to generate acetate via the Pta-Ack pathway. Members of the Bacteroides and Geobacter also encoded Pta-Ack and together, all four populations had the putative ability to degrade ethylbenzene, phenylphosphate, phenylethanol, toluene, xylene, and phenol. Metabolic reconstructions, gene analyses, and environmental parameters also indicated that redox fluctuations likely promote facultative energy metabolisms in the coal seam. The active "Chlorobi PRB" MAG encoded enzymes for fermentation, nitrate reduction, and multiple oxygenases with varying binding affinities for oxygen. "M. paradoxum PRB" encoded an extradiol dioxygenase for aerobic phenylacetate degradation, which was also present in previously published Methanothrix genomes. These observations outline underlying processes for bio-methane from subbituminous coal by translationally active populations and demonstrate activity-based metagenomics as a powerful strategy in next generation physiology to understand ecologically relevant microbial populations.

Host population diversity as a driver of viral infection cycle in wild populations of green sulfur bacteria with long standing virus-host interactions.

Temperate phages are viruses of bacteria that can establish two types of infection: a lysogenic infection in which the virus replicates with the host cell without producing virions, and a lytic infection where the host cell is eventually destroyed, and new virions are released. While both lytic and lysogenic infections are routinely observed in the environment, the ecological and evolutionary processes regulating these viral dynamics are still not well understood, especially for uncultivated virus-host pairs. Here, we characterized the long-term dynamics of uncultivated viruses infecting green sulfur bacteria (GSB) in a model freshwater lake (Trout Bog Lake, TBL). As no GSB virus has been formally described yet, we first used two complementary approaches to identify new GSB viruses from TBL; one in vitro based on flow cytometry cell sorting, the other in silico based on CRISPR spacer sequences. We then took advantage of existing TBL metagenomes covering the 2005-2018 period to examine the interactions between GSB and their viruses across years and seasons. From our data, GSB populations in TBL were constantly associated with at least 2-8 viruses each, including both lytic and temperate phages. The dominant GSB population in particular was consistently associated with two prophages with a nearly 100% infection rate for >10 years. We illustrate with a theoretical model that such an interaction can be stable given a low, but persistent, level of prophage induction in low-diversity host populations. Overall, our data suggest that lytic and lysogenic viruses can readily co-infect the same host population, and that host strain-level diversity might be an important factor controlling virus-host dynamics including lytic/lysogeny switch.

A pipeline for targeted metagenomics of environmental bacteria.

BACKGROUND: Metagenomics and single cell genomics provide a window into the genetic repertoire of yet uncultivated microorganisms, but both methods are usually taxonomically untargeted. The combination of fluorescence in situ hybridization (FISH) and fluorescence activated cell sorting (FACS) has the potential to enrich taxonomically well-defined clades for genomic analyses. METHODS: Cells hybridized with a taxon-specific FISH probe are enriched based on their fluorescence signal via flow cytometric cell sorting. A recently developed FISH procedure, the hybridization chain reaction (HCR)-FISH, provides the high signal intensities required for flow cytometric sorting while maintaining the integrity of the cellular DNA for subsequent genome sequencing. Sorted cells are subjected to shotgun sequencing, resulting in targeted metagenomes of low diversity. RESULTS: Pure cultures of different taxonomic groups were used to (1) adapt and optimize the HCR-FISH protocol and (2) assess the effects of various cell fixation methods on both the signal intensity for cell sorting and the quality of subsequent genome amplification and sequencing. Best results were obtained for ethanol-fixed cells in terms of both HCR-FISH signal intensity and genome assembly quality. Our newly developed pipeline was successfully applied to a marine plankton sample from the North Sea yielding good quality metagenome assembled genomes from a yet uncultivated flavobacterial clade. CONCLUSIONS: With the developed pipeline, targeted metagenomes at various taxonomic levels can be efficiently retrieved from environmental samples. The resulting metagenome assembled genomes allow for the description of yet uncharacterized microbial clades. Video abstract.

Probing the active fraction of soil microbiomes using BONCAT-FACS.

The ability to link soil microbial diversity to soil processes requires technologies that differentiate active microbes from extracellular DNA and dormant cells. Here, we use BONCAT (bioorthogonal non-canonical amino acid tagging) to measure translationally active cells in soils. We compare the active population of two soil depths from Oak Ridge (Tennessee, USA) and find that a maximum of 25-70% of the extractable cells are active. Analysis of 16S rRNA sequences from BONCAT-positive cells recovered by fluorescence-activated cell sorting (FACS) reveals that the phylogenetic composition of the active fraction is distinct from the total population of extractable cells. Some members of the community are found to be active at both depths independently of their abundance rank, suggesting that the incubation conditions favor the activity of similar organisms. We conclude that BONCAT-FACS is effective for interrogating the active fraction of soil microbiomes in situ and provides a new approach for uncovering the links between soil processes and specific microbial groups.

Optimizing de novo genome assembly from PCR-amplified metagenomes.

BACKGROUND: Metagenomics has transformed our understanding of microbial diversity across ecosystems, with recent advances enabling de novo assembly of genomes from metagenomes. These metagenome-assembled genomes are critical to provide ecological, evolutionary, and metabolic context for all the microbes and viruses yet to be cultivated. Metagenomes can now be generated from nanogram to subnanogram amounts of DNA. However, these libraries require several rounds of PCR amplification before sequencing, and recent data suggest these typically yield smaller and more fragmented assemblies than regular metagenomes. METHODS: Here we evaluate de novo assembly methods of 169 PCR-amplified metagenomes, including 25 for which an unamplified counterpart is available, to optimize specific assembly approaches for PCR-amplified libraries. We first evaluated coverage bias by mapping reads from PCR-amplified metagenomes onto reference contigs obtained from unamplified metagenomes of the same samples. Then, we compared different assembly pipelines in terms of assembly size (number of bp in contigs ≥ 10 kb) and error rates to evaluate which are the best suited for PCR-amplified metagenomes. RESULTS: Read mapping analyses revealed that the depth of coverage within individual genomes is significantly more uneven in PCR-amplified datasets versus unamplified metagenomes, with regions of high depth of coverage enriched in short inserts. This enrichment scales with the number of PCR cycles performed, and is presumably due to preferential amplification of short inserts. Standard assembly pipelines are confounded by this type of coverage unevenness, so we evaluated other assembly options to mitigate these issues. We found that a pipeline combining read deduplication and an assembly algorithm originally designed to recover genomes from libraries generated after whole genome amplification (single-cell SPAdes) frequently improved assembly of contigs ≥10 kb by 10 to 100-fold for low input metagenomes. CONCLUSIONS: PCR-amplified metagenomes have enabled scientists to explore communities traditionally challenging to describe, including some with extremely low biomass or from which DNA is particularly difficult to extract. Here we show that a modified assembly pipeline can lead to an improved de novo genome assembly from PCR-amplified datasets, and enables a better genome recovery from low input metagenomes.

Hidden diversity of soil giant viruses.

Known giant virus diversity is currently skewed towards viruses isolated from aquatic environments and cultivated in the laboratory. Here, we employ cultivation-independent metagenomics and mini-metagenomics on soils from the Harvard Forest, leading to the discovery of 16 novel giant viruses, chiefly recovered by mini-metagenomics. The candidate viruses greatly expand phylogenetic diversity of known giant viruses and either represented novel lineages or are affiliated with klosneuviruses, Cafeteria roenbergensis virus or tupanviruses. One assembled genome with a size of 2.4 Mb represents the largest currently known viral genome in the Mimiviridae, and others encode up to 80% orphan genes. In addition, we find more than 240 major capsid proteins encoded on unbinned metagenome fragments, further indicating that giant viruses are underexplored in soil ecosystems. The fact that most of these novel viruses evaded detection in bulk metagenomes suggests that mini-metagenomics could be a valuable approach to unearth viral giants.

Visualizing in situ translational activity for identifying and sorting slow-growing archaeal-bacterial consortia.

To understand the biogeochemical roles of microorganisms in the environment, it is important to determine when and under which conditions they are metabolically active. Bioorthogonal noncanonical amino acid tagging (BONCAT) can reveal active cells by tracking the incorporation of synthetic amino acids into newly synthesized proteins. The phylogenetic identity of translationally active cells can be determined by combining BONCAT with rRNA-targeted fluorescence in situ hybridization (BONCAT-FISH). In theory, BONCAT-labeled cells could be isolated with fluorescence-activated cell sorting (BONCAT-FACS) for subsequent genetic analyses. Here, in the first application, to our knowledge, of BONCAT-FISH and BONCAT-FACS within an environmental context, we probe the translational activity of microbial consortia catalyzing the anaerobic oxidation of methane (AOM), a dominant sink of methane in the ocean. These consortia, which typically are composed of anaerobic methane-oxidizing archaea (ANME) and sulfate-reducing bacteria, have been difficult to study due to their slow in situ growth rates, and fundamental questions remain about their ecology and diversity of interactions occurring between ANME and associated partners. Our activity-correlated analyses of >16,400 microbial aggregates provide the first evidence, to our knowledge, that AOM consortia affiliated with all five major ANME clades are concurrently active under controlled conditions. Surprisingly, sorting of individual BONCAT-labeled consortia followed by whole-genome amplification and 16S rRNA gene sequencing revealed previously unrecognized interactions of ANME with members of the poorly understood phylum Verrucomicrobia This finding, together with our observation that ANME-associated Verrucomicrobia are found in a variety of geographically distinct methane seep environments, suggests a broader range of symbiotic relationships within AOM consortia than previously thought.