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.

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.

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.

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.

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.

Autoantibody levels in myositis patients correlate with clinical response during B cell depletion with rituximab.

OBJECTIVES: To determine the longitudinal trends in serum levels of four myositis-associated autoantibodies: anti-Jo-1, -transcription intermediary factor 1 γ (TIF1-γ), -signal recognition particle (SRP) and -Mi-2, after B cell depletion with rituximab, and to determine the longitudinal association of these autoantibody levels with disease activity as measured by myositis core-set measures (CSMs). METHODS: Treatment-resistant adult and pediatric myositis subjects (n = 200) received rituximab in the 44-week Rituximab in Myositis Trial. CSMs [muscle enzymes, manual muscle testing (MMT), physician and patient global disease activity, HAQ, and extramuscular disease activity] were evaluated monthly and anti-Jo-1 (n = 28), -TIF1-γ (n = 23), -SRP (n = 25) and -Mi-2 (n = 26) serum levels were measured using validated quantitative ELISAs. Temporal trends and the longitudinal relationship between myositis-associated autoantibodies levels and CSM were estimated using linear mixed models. RESULTS: Following rituximab, anti-Jo-1 levels decreased over time (P < 0.001) and strongly correlated with all CSMs (P < 0.008). Anti-TIF1-γ levels also decreased over time (P < 0.001) and were only associated with HAQ, MMT and physician and patient global disease activity. Anti-SRP levels did not change significantly over time, but were significantly associated with serum muscle enzymes. Anti-Mi-2 levels significantly decreased over time and were associated with muscle enzymes, MMT and the physician global score. CONCLUSION: Anti-Jo-1, anti-TIF1-γ and anti-Mi-2 levels in myositis subjects decreased after B cell depletion and were correlated with changes in disease activity, whereas anti-SRP levels were only associated with longitudinal muscle enzyme levels. The strong association of anti-Jo-1 levels with clinical outcomes suggests that anti-Jo-1 autoantibodies may be a good biomarker for disease activity.

Anti-signal recognition particle autoantibody ELISA validation and clinical associations.

OBJECTIVE: The aim of this study was to develop and validate a quantitative anti-signal recognition particle (SRP) autoantibody serum ELISA in patients with myositis and longitudinal association with myositis disease activity. METHODS: We developed a serum ELISA using recombinant purified full-length human SRP coated on ELISA plates and a secondary antibody that bound human IgG to detect anti-SRP binding. Protein immunoprecipitation was used as the gold standard for the presence of anti-SRP. Serum samples from three groups were analysed: SRP(+) myositis subjects by immunoprecipitation, SRP(-) myositis subjects by immunoprecipitation and non-myositis controls. The ELISA's sensitivity, specificity, positive predictive value and negative predictive value were evaluated. Percentage agreement and test-retest reliability were assessed. Serial samples from seven SRP immunoprecipitation-positive subjects were also tested, along with serum muscle enzymes and manual muscle testing. RESULTS: Using immunoprecipitation, we identified 26 SRP(+) myositis patients and 77 SRP(-) controls (including 38 patients with necrotizing myopathy). Non-myositis control patients included SLE (n = 4) and SSc (n = 7) patients. Anti-SRP positivity by ELISA showed strong agreement (97.1%) with immunoprecipitation (κ = 0.94). The sensitivity, specificity, positive predictive value, and negative predictive value of the anti-SRP ELISA were 88, 100, 100 and 96, respectively. The area under the curve was 0.94, and test-retest reliability was strong (r = 0.91, P < 0.001). Serial samples showed that anti-SRP levels paralleled changes in muscle enzymes and manual muscle testing. CONCLUSION: We developed a quantitative ELISA for detecting serum anti-SRP autoantibodies and validated the assay in myositis. Longitudinal assessment of SRP levels by ELISA may be a useful biomarker for disease activity.

Anti-transcription intermediary factor 1-gamma autoantibody ELISA development and validation.

OBJECTIVES: A quantitative anti-transcription intermediary factor 1-gamma (anti-TIF1-γ) ELISA may improve the detection of cancer-associated myositis (CAM). The aims of this study were the development and validation of a quantitative anti-TIF1-γ autoantibody ELISA in patients with myositis. METHODS: We developed an ELISA using recombinant purified full-length human TIF1-γ. Patient serum was incubated with TIF1-γ-coated ELISA plates, and secondary antibody that bound human IgG was used to detect anti-TIF1-γ binding. Protein immunoprecipitation (IP) was used as the gold standard for the presence of anti-TIF1-γ. Serum samples from myositis patients with positive and negative anti-TIF1-γ by IP, from non-myositis autoimmune patients (SSc, SLE and RA) and from healthy controls were analysed. The ELISA's sensitivity, specificity, positive predictive value (PPV) and negative predictive value (NPV) were evaluated. Agreement between the ELISA and IP results was determined using chi-squared and kappa tests. Test-retest reliability of the ELISA was assessed. RESULTS: We identified 55 myositis patients with and 111 controls without anti-TIF1-γ by IP. Anti-TIF1-γ positivity by ELISA showed strong agreement (93.9%) with IP results (κ = 0.87). The sensitivity, specificity, PPV, NPV and overall accuracy of the anti-TIF1-γ ELISA were 91%, 96%, 93%, 95% and 94%, respectively. The area under the curve (AUC) of a receiver operating characteristic (ROC) curve was 0.938. Test-retest reliability was strong (Pearson r = 0.913, P < 0.001). CONCLUSION: We developed a quantitative ELISA for detecting serum anti-TIF1-γ autoantibodies and validated the assay in myositis and other connective tissue disease patients. The availability of a validated, quantitative ELISA should improve the detection of anti-TIF1-γ autoantibodies and may improve the detection of CAM.

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.

scMicrobe PTA: near complete genomes from single bacterial cells.

Microbial genomes produced by standard single-cell amplification methods are largely incomplete. Here, we show that primary template-directed 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 multiple displacement amplification-based approaches were usually <30% complete. PTA-derived genomes also included more associated viruses and biosynthetic gene clusters.

The salmonella transcriptome in lettuce and cilantro soft rot reveals a niche overlap with the animal host intestine.

Fresh vegetables have been recurrently associated with salmonellosis outbreaks, and Salmonella contamination of retail produce has been correlated positively with the presence of soft rot disease. We observed that population sizes of Salmonella enterica serovar Typhimurium SL1344 increased 56-fold when inoculated alone onto cilantro leaves, versus 2,884-fold when coinoculated with Dickeya dadantii, a prevalent pathogen that macerates plant tissue. A similar trend in S. enterica populations was observed for soft-rotted lettuce leaves. Transcriptome analysis of S. enterica cells that colonized D. dadantii-infected lettuce and cilantro leaves revealed a clear shift toward anaerobic metabolism and catabolism of substrates that are available due to the degradation of plant cells by the pectinolytic pathogen. Twenty-nine percent of the genes that were upregulated in cilantro macerates were also previously observed to have increased expression levels in the chicken intestine. Furthermore, multiple genes induced in soft rot lesions are also involved in the colonization of mouse, pig, and bovine models of host infection. Among those genes, the operons for ethanolamine and propanediol utilization as well as for the synthesis of cobalamin, a cofactor in these pathways, were the most highly upregulated genes in lettuce and cilantro lesions. In S. Typhimurium strain LT2, population sizes of mutants deficient in propanediol utilization or cobalamin synthesis were 10- and 3-fold lower, respectively, than those of the wild-type strain in macerated cilantro (P < 0.0002); in strain SL1344, such mutants behaved similarly to the parental strain. Anaerobic conditions and the utilization of nutrients in macerated plant tissue that are also present in the animal intestine indicate a niche overlap that may explain the high level of adaptation of S. enterica to soft rot lesions, a common postharvest plant disease.

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.

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.

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.

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.

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.

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.

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.