# Long-Term Incubation of Lake Water Enables Genomic Sampling of Consortia Involving Planctomycetes and Candidate Phyla Radiation Bacteria

Alexander L Jaffe, Maxime Fuster, Marie C Schoelmerich, Lin-Xing Chen, Jonathan Colombet, Hermine Billard, Télesphore Sime-Ngando, Jillian F Banfield

<figure><img src="image_1.png" /></figure>

## To cite this version:

Alexander L Jaffe, Maxime Fuster, Marie C Schoelmerich, Lin-Xing Chen, Jonathan Colombet, et al.. Long-Term Incubation of Lake Water Enables Genomic Sampling of Consortia Involving Planctomycetes and Candidate Phyla Radiation Bacteria. mSystems, 2022, 7 (2), 10.1128/msystems.00223-22. [https://doi.org/10.1128/msystems.00223-22](https://doi.org/10.1128/msystems.00223-22)

HAL Id: hal-03865093

https://hal.science/hal-03865093

Submitted on 22 Nov 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

# Long-Term Incubation of Lake Water Enables Genomic Sampling of Consortia Involving *Planctomycetes* and Candidate Phyla Radiation Bacteria

Alexander L. Jaffe,ª Maxime Fuster,ᵇ Marie C. Schoelmerich,ᶜ Lin-Xing Chen,ᵈ Jonathan Colombet,ᵇ Hermine Billard,ᵇ
Télesphore Sime-Ngando,ᵇ Jillian F. Banfieldᶜ,d,e,f

Alexander L. Jaffe and Maxime Fuster contributed equally to this work. Author order was determined based on seniority.

**ABSTRACT** Microbial communities in lakes can profoundly impact biogeochemical processes through their individual activities and collective interactions. However, the complexity of these communities poses challenges, particularly for studying rare organisms such as Candidate Phyla Radiation bacteria (CPR) and enigmatic entities such as aster-like nanoparticles (ALNs). Here, a reactor was inoculated with water from Lake Fargette, France, and maintained under dark conditions at 4°C for 31 months and enriched for ALNs, diverse *Planctomycetes*, and CPR bacteria. We reconstructed draft genomes and predicted metabolic traits for 12 diverse *Planctomycetes* and 9 CPR bacteria, some of which are likely representatives of undescribed families or genera. One CPR genome representing the little-studied lineage *"Candidatus Peribacter"* was curated to completion (1.239 Mbp) and unexpectedly encodes the full gluconeogenesis pathway. Metatranscriptomic data indicate that some planctomycetes and CPR bacteria were active under the culture conditions, accounting for ~30% and ~1% of RNA reads mapping to the genome set, respectively. We also reconstructed genomes and obtained transmission electron microscope images for numerous viruses, including one with a >300-kbp genome and several predicted to infect *Planctomycetes*. Together, our analyses suggest that freshwater *Planctomycetes* are central players in a subsystem that includes ALNs, symbiotic CPR bacteria, and viruses.

**IMPORTANCE** Laboratory incubations of natural microbial communities can aid in the study of member organisms and their networks of interaction. This is particularly important for understudied lineages for which key elements of basic biology are still emerging. Using genomics and microscopy, we found that members of the bacterial lineage *Planctomycetes* may be central players in a subset of a freshwater lake microbiome that includes other bacteria, archaea, viruses, and mysterious entities, called aster-like nanoparticles (ALNs), whose origin is unknown. Our results help constrain the possible origins of ALNs and provide insight into possible interactions within a complex lake ecosystem.

**KEYWORDS** CPR bacteria, aster-like nanoparticles, lake microbiome, metagenomics, viruses

Freshwater lakes host diverse microbial communities that likely control ecosystem biogeochemistry (1, 2). Here, we established a laboratory culture based on an inoculum from Lake Fargette, France, a site chosen as part of a parallel study of enigmatic aster-like nanoparticles (ALNs) (3, 4). ALNs are enigmatic, organic, femtoplankton entities that exhibit bloom-like behavior in various freshwater and coastal environments (3, 4). In this experiment,

**Editor** Jack A. Gilbert, University of California
San Diego

**Copyright** © 2022 Jaffe et al. This is an open-
access article distributed under the terms of
the Creative Commons Attribution 4.0
International license.

Address correspondence to Jillian F. Banfield,
jbanfield@berkeley.edu.

The authors declare a conflict of interest. J.F.B.
is a founder of Metagenomi. The other authors
declare no competing interests.

Received 11 March 2022

Accepted 11 March 2022

<footnote>ªDepartment of Plant and Microbial Biology, University of California, Berkeley, California, USA</footnote>

<footnote>ᵇLaboratoire Microorganismes: Génome et Environnement (LMGE), UMR CNRS 6023, Université Clermont-Auvergne, Clermont-Ferrand, France</footnote>

<footnote>ᶜInnovative Genomics Institute, University of California, Berkeley, California, USA</footnote>

<footnote>ᵈDepartment of Earth and Planetary Science, University of California, Berkeley, California, USA</footnote>

<footnote>ᵉDepartment of Environmental Science, Policy, and Management, University of California, Berkeley, California, USA</footnote>

<footnote>fChan Zuckerberg Biohub, San Francisco, California, USA</footnote>
microscopy showed that the proportion of ALNs increased substantially during incubation, from 26% to 36% of imaged objects (see Table S1 at https://zenodo.org/record/5362898#.YiuiRBPML5Y). To seek clues to the origins of ALNs and the organisms that they might associate with, and to better understand the lake ecosystem overall, we studied the culture using a combination of metagenomics, metatranscriptomics, and microscopy. In doing so, we recovered draft genomes for abundant and transcriptionally active *Planctomycetes* as well as CPR bacteria, phages, and eukaryotic viruses. Overall, we provide clues to associations and potential interactions among microbial groups in a lake ecosystem.

**Community composition and genome reconstruction.** Concentrate from the 0.2-µm to 25-µm size fraction of the highly eutrophic Lake Targette, France, was incubated at 4°C in the dark with filtered and sterilized lake water (<20 kDa). The incubation was performed under dark conditions to favor heterotrophic organisms associated with ALNs instead of phototrophic eukaryotes. After 31 months, DNA and RNA were extracted for metagenomic and transcriptomic analyses. DNA reads were assembled and the scaffolds profiled to identify ribosomal protein S3 (rpS3). Profiling using the predicted rpS3 protein sequences revealed that the enrichment was bacterially dominated, although several members of the *Thaumarchaeota* were present (Fig. 1a; see also Table S2 at https://zenodo.org/record/5362898#.YiuiRBPML5Y). The most abundant organisms overall were *Planctomycetes* (~27% of overall rpS3 coverage), including the most abundant singular organism (Fig. 1 and Table S2). CPR bacteria were 5 of the top 25 most abundant organism groups (~9% overall rpS3 coverage; Fig. 1 and Table S2). Compared to baseline abundances in Lake Targette, these results indicated enrichment of these groups, particularly CPR bacteria, which were barely detectable in the lake (<1% relative abundance; Fig. 1b and Table S3).

Where possible, scaffolds were assigned to genome bins that ranged in quality from draft to nearly complete. The 48 genomes captured most of the phylogenetic diversity (76% of the most abundant rpS3 genes; Fig. 1a and Tables S2 and S4 at https://zenodo.org/record/5362898#.YiuiRBPML5Y); 12 genomes represent phylogenetically diverse *Planctomycetes*, including several from the *Planctomycetes* and *Phycisphaerae* classes (see Fig. S1 and Table S4). Metabolic reconstructions suggested that the *Planctomycetes* are primarily heterotrophs with the potential to oxidize nitrite or reduce nitrate in three cases (Fig. 1e). We also recovered 9 genomes of CPR bacteria, 8 of which were classified as "Candidatus Parcubacteria" (Fig. S1). The "Candidatus Parcubacteria" genomes encoded minimal metabolic capacities, consistent with symbiotic lifestyles (5). However, several had a phylogenetically distinct *nirK* gene that may play a role in denitrification or energy conservation (6) (Fig. 1e). Read mapping from a metatranscriptome collected contemporaneously with the metagenome suggested that some *Planctomycetes* and, to a lesser extent, CPR bacteria were actively transcribing under the culture conditions, with *Planctomycetes* and CPR accounting for about 27% and 1%, respectively, of RNA reads stringently mapping to the nonredundant set of genomes (Fig. 1d and Table S4). However, only a small proportion of total RNA reads (201,809 reads, or ~0.2% of the total) mapped stringently to the genome set and passed filtering for non-mRNAs (see the supplemental material at https://zenodo.org/record/5362898#.YiuiRBPML5Y).

The ninth draft CPR genome was for a member of the undersampled CPR lineage "Candidatus Peribacteria." We manually curated the original bin of 6 fragments into a single fragment that was circularized by a small, unbinned contig, and all scaffolding errors and local misassemblies were fixed (see Fig. S2 at https://zenodo.org/record/5362898#.YiuiRBPML5Y). The newly reported, fully curated 1,239,242-bp genome shares ~89% similarity in its 16S rRNA gene and is largely syntenous with a closely related "Candidatus Peribacteria" genome from Rifle, Colorado (7), supporting the accuracy of both assemblies (see Fig. S3). Unlike the Rifle genomes, this peribacter likely cannot synthesize purines de novo. Other notable differences include the presence of vacuolar-type H⁺/Na⁺-transporting ATPase complex and the lack of genes for biosynthesis of mevalonate. Based on biosynthetic deficiencies, we conclude that this bacterium probably relies on other organisms for many building block but to a lesser extent than the "Candidatus Parcubacteria." Supporting this idea, we observed that the peribacter genome encodes all gluconeogenesis enzymes, including
<figure><img src="image_4.png" /><figcaption>**FIG 1** Long-term incubation enriched for members of the *Planctomycetes* and CPR bacteria. (a) Rank abundance curve based on ribosomal protein S3 (rps3) coverage for organisms recovered at the end of incubation. Asterisks indicate marker genes that were binned into genomes. (b) Relative abundance of CPR bacteria and *Planctomycetes* during monthly sampling of Lake Fargette based on 16S rRNA amplicon sequencing. Each point represents the relative abundance of CPR bacteria or *Planctomycetes* in a given month. (c) Overall community composition at the end of incubation, based on cumulative coverage of rpS3. Panel d displays the fraction of RNA reads from the end of incubation that could be mapped to genomes. (e) Sequence characteristics, metabolic predictions, and CRISPR-Cas loci information for genomes affiliated with the *Planctomycetes* and CPR bacteria. Cells with color fill indicate the fraction of key genes for each pathway (as defined by KEGGDecoder) that are present. X indicates genomes with CRISPR loci and if those loci contained spacers targeting at least one curated phage genome from the sample.</figcaption></figure>

fructose-1,6-bisphosphatase I, which was not found in the "Candidatus Parcubacteria" genomes (Fig. 1e).

We reconstructed 12 phage sequences and 5 phage-like sequences, including 3 circularized genomes exceeding 100 kbp and 2 incomplete phage/phage-like fragments of >300 kbp (see Table S5 at https://zenodo.org/record/5362898#.YiuiRBPML5Y). Phylogenetic analyses of encoded terminase and capsid proteins suggested that the phages likely fall within the *Caudovirales*, which are known to include numerous tailed phages with large capsids (8). Additionally, we used phage gene content and analyses of bacterial CRISPR loci to infer that phages infect *Planctomycetes*, *Proteobacteria*, and *Bacteroidetes* (Fig. 2 and Table S5).
A.

<figure><img src="image_5.png" /></figure>

B. viral id

<figure><img src="image_6.png" /><figcaption>**FIG 2** TEM imaging and viral genomics in the enrichment culture. (A) Featured in this image is a cell inferred to be a *Planctomycetes* (pl) with a characteristic stalk and holdfast (h). Attached to the cell are four phage particles (v) (two different sizes; thus, likely different phages). Also visible is one large tailed jumbo phage (ph) that is 145 nm in diameter with a 214-nm tail as well as several aster-like nanoparticles (ALNs) (a). (B) Genome sizes and predicted hosts for phages and eukaryotic viruses.</figcaption></figure>

We also reconstructed large fragments of 17 eukaryotic viruses. Based on the phylo-
genetic placement of the major capsid protein and homologs of the poxvirus late transcription
factor VLTF3, some viruses belonged to *Iridoviridae* (Betairidoviridae), extended *Mimiviridae*,
*Phycodnaviridae*, and Pitho-like viruses (see Fig. S4 and S5 and Table S5 at https://zenodo.org/
record/5362898#.YiuiRBPML5Y). Other viruses could only be classified at the superclade level,
including those within a potentially novel *Phycodnaviridae*, *Asfarviridae*, *Megavirales* (PAM)
clade, or did not contain either marker protein (see Fig. S4 and S5). Interestingly, we detected
very little transcription of genes from phage or eukaryotic viruses (mean coverage, ≪1×), sug-
gesting that they were not actively replicating in the incubated community (see Table S5). The
lack of eukaryotes in the enrichment suggests that some of these particles have derived from
the inoculum and persisted for over 2 years.

**Imaging of incubation community.** We imaged a diversity of cellular, cellular-like, and noncellular particles, as well as many ALNs, using transmission electron microscopy (TEM) (Fig. 2 and Fig. S6 at https://zenodo.org/record/5362898#.YiuiRBPML5Y). Based on the presence of the extracellular holdfast (9, 10), we infer that the cell imaged in Fig. 2 is likely a member of the *Planctomycetes*, with at least two attached distinct types of virus-like particles (VLPs). This finding is consistent with the high abundance of *Planctomycetes* in the enrichment as well as multiple phages predicted to infect them (~18 to 100× coverage). We also imaged numerous tailed phages (see Fig. S6), including one with a 145-nm-diameter capsid and 214-nm tail (Fig. 2, ph). This large capsid size is consistent with those of jumbo phages with genomes in the 300-kbp range (8), two of which were reconstructed here (see Table S5). Despite these observations, TEM image counts suggested that the overall proportion of tailed phages decreased from 26% to less than 1% of imaged objects over the course of the incubation (see Table S1).

DISCUSSION

*Planctomycetes* are globally distributed across freshwater ecosystems, where they are thought to play important roles in nitrogen and carbon cycling (11, 12). Our analyses expand genomic sampling for these organisms, and coenrichment suggests that they interact
with CPR bacteria. Although we cannot establish a direct association between *Planctomycetes* and CPR bacteria from the current data, reported lifestyles for other CPR suggest that they are host cell-attached, at least at sometimes (13). Furthermore, the enrichment contains a similar level of diversity of both *Planctomycetes* and CPR, raising the possibility of species-specific associations (see Fig. S1 at https://zenodo.org/record/5362898#.YiuiRBPML5Y). As episymbionts, CPR would almost certainly influence the physiology of *Planctomycetes* cells, either via mutualistic interactions or parasitism, and, thus, the biogeochemical cycles that *Planctomycetes* mediate. The approach used in the current study could guide future coevolution of CPR bacteria of the "Candidatus Parcubacteria" and "Candidatus Peribacteria" lineages to test this hypothesis.

Like CPR bacteria, phages clearly impact their hosts and, thus, ecosystem structure and performance. It is intriguing that phages were maintained over the long incubation period, as they are often lost from laboratory cultures (14). A subset of the phages clearly infect *Planctomycetes*, the diversity of which may have enabled their sustained replication.

One motivation of our study was to seek clues to the origin of ALNs, and it is also intriguing to find them coenriched with *Planctomycetes* and CPR bacteria. Their defined morphologies suggest that they develop under genetic control; however, the genetic system responsible for producing ALNs remains unknown. We found no evidence for nucleic acid sequences from completely unknown organisms or mobile genetic elements or were not able to extract them with the methods used here. Thus, one possible interpretation is that ALNs originate from cooccurring bacteria and that *Planctomycetes* and CPR should be among candidates for further study.

Materials and methods can be found in the supplemental material at https://zenodo.org/record/5362898#.YiuiRBPML5Y.

**Data availability.** Read data and draft genomes from this study are available through NCBI at PRJNA757735. Genome accession information for the 48 bacterial and archaeal genomes is also listed in Table S3 at https://zenodo.org/record/5362898#.YiuiRBPML5Y. Custom codes for the described analyses are also available on GitHub (github.com/alexanderjaffe/alm-enrichment). All supplementary figures, tables, extended data files, and genomes (including phage/viral genomes) are also available through Zenodo (https://doi.org/10.5281/zenodo.5362897).

## ACKNOWLEDGMENTS

We thank Alex Crits-Christoph, Adair Borges, Rohan Sachdeva, and Yue Clare Lou for informatics support and helpful discussions. We also thank Plateforme CYSTEM-UCA PARTNER (Clermont-Ferrand, France) for their technical support and expertise, the Innovative Genomics Institute, and the Vincent J. Coates Genomics Sequencing Laboratory at UC Berkeley.

M.F., H.B., J.C., and T.S.-N. performed sample collection, incubation experiments, and microscopy. A.L.J., L.-X.C., M.C.S., and J.F.B. performed metagenomics and sequence analysis, and J.F.B carried out the manual genome curation. All authors contributed to project design and manuscript writing.

J.F.B. is a founder of Metagenomi. The other authors declare no competing interests.

Funding was provided by the Innovative Genomics Institute at UC Berkeley, the Chan Zuckerberg Biohub, and a Moore Foundation grant 71785. M.F. was supported by a Ph.D. fellowship from the CPER 2015-2020 SYMBIOSE challenge program (French Ministry of Research, UCA, CNRS, INRA, Auvergne-Rhône-Alpes Region, FEDER). M.C.S. was supported by a DFG fellowship. This study is a contribution to the "C NO LIMIT" project funded by the Interdisciplinary Mission of the French National Center of Scientific Research (CNRS) Program X-life, 2018 edition. This research was also financed by the French government IDEX-ISITE initiative 16-IDEX-0001 (CAP 20-25).

## REFERENCES

1. Grossart H-P, Massana R, McMahon KD, Walsh DA. 2020. Linking metagenomics to aquatic microbial ecology and biogeochemical cycles. Limnol Oceanogr 65:e11382. https://doi.org/10.1002/Ino.11382.

2. Tran PQ, McIntyre PB, Kraemer BM, Vadeboncoeur Y, Kimirei IA, Tamatamah R, et al. 2019. Depth-discrete eco-genomics of Lake Tanganyika reveals roles of diverse microbes, including candidate phyla, in tropical freshwater nutrient cycling. bioRxiv https://doi.org/10.1101/834861.

3. Colombet J, Billard H, Viguès B, Balor S, Boulé C, Geay L, Benzerara K, Menguy N, Ilango G, Fuster M, Enault F, Bardot C, Gautier V, Pradeep Ram AS, Sime-Ngando T. 2019. Discovery of high abundances of Aster-like Downloaded from https://journals.asm.org/journal/msystems on 01 April 2022 by 193.54.49.188.
nanoparticles in pelagic environments. *Front Microbiol* 10:2376. https://doi.org/10.3389/fmicb.2019.02376.

4. Fuster M, Billard H, Mandart M, Steiger J, Sime-Ngando T, Colombet J. 2020. Trophic conditions influence widespread distribution of Aster-like nanoparticles within aquatic environments. *Microb Ecol* 80:741-745. https://doi.org/10.1007/s00248-020-01541-6.

5. Castelle CJ, Brown CT, Anantharaman K, Probst AJ, Huang RH, Banfield JF. 2018. Biosynthetic capacity, metabolic variety and unusual biology in the CPR and DPANN radiations. *Nat Rev Microbiol* 16:629-645. https://doi.org/10.1038/s41579-018-0076-2.

6. Danczak RE, Johnston MD, Kenah C, Slattery M, Wrighton KC, Wilkins MJ. 2017. Members of the Candidate Phyla Radiation are functionally differentiated by carbon- and nitrogen-cycling capabilities. *Microbiome* 5:112. https://doi.org/10.1186/s40168-017-0331-1.

7. Anantharaman K, Brown CT, Burstein D, Castelle CJ, Probst AJ, Thomas BC, Williams KH, Banfield JF. 2016. Analysis of five complete genome sequences for members of the class *Peribacteria* in the recently recognized *Peregrinibacteria* bacterial phylum. *PeerJ* 4:e1607. https://doi.org/10.7717/peerj.1607.

8. Yuan Y, Gao M. 2017. Jumbo bacteriophages: an overview. *Front Microbiol* 8:403.

9. Boedeker C, Schüler M, Reintjes G, Jeske O, van Teeseling MCF, Jogler M, Rast P, Borchert D, Devos DP, Kucklick M, Schaffer M, Kolter R, van Niftrik L, Engelmann S, Amann R, Rohde M, Engelhardt H, Jogler C. 2017. Determining the bacterial cell biology of Planctomycetes. *Nat Commun* 8:14853. https://doi.org/10.1038/ncomms14853.

10. Jogler C, Glöckner FO, Kolter R. 2011. Characterization of Planctomyces limnophilus and development of genetic tools for its manipulation establish it as a model species for the phylum Planctomycetes. *Appl Environ Microbiol* 77:5826-5829. https://doi.org/10.1128/AEM.05132-11.

11. Wiegand S, Jogler M, Jogler C. 2018. On the maverick Planctomycetes. *FEMS Microbiol Rev* 42:739-760. https://doi.org/10.1093/femsre/fuy029.

12. Andrei A, Salcher MM, Mehrshad M, Rychtecký P, Znachor P, Ghai R. 2019. Niche-directed evolution modulates genome architecture in freshwater Planctomycetes. *ISME J* 13:1056-1071. https://doi.org/10.1038/s41396-018-0332-5.

13. He X, McLean JS, Edlund A, Yooseph S, Hall AP, Liu S-Y, Dorrestein PC, Esquenazi E, Hunter RC, Cheng G, Nelson KE, Lux R, Shi W. 2015. Cultivation of a human-associated TM7 phylotype reveals a reduced genome and epi-biotic parasitic lifestyle. *Proc Natl Acad Sci U S A* 112:244-249. https://doi.org/10.1073/pnas.1419038112.

14. Paez-Espino D, Sharon I, Morovic W, Stahl B, Thomas BC, Barrangou R, Banfield JF. 2015. CRISPR immunity drives rapid phage genome evolution in *Streptococcus thermophilus*. *mBio* 6:e00262-15. https://doi.org/10.1128/mBio.00262-15.