# First evidence of the presence and activity of archaeal C3 group members in an Atlantic intertidal mudflat

Céline Lavergne, Mylène Hugoni, Christine Dupuy, Hélène Agogué

► To cite this version:

Céline Lavergne, Mylène Hugoni, Christine Dupuy, Hélène Agogué. First evidence of the presence and activity of archaeal C3 group members in an Atlantic intertidal mudflat. Scientific Reports, 2018, 8 (1), pp.11790. 10.1038/s41598-018-30222-1 . hal-01858404

HAL Id: hal-01858404

https://hal.science/hal-01858404v1

Submitted on 20 Aug 2018

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_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
SCIENTIFIC REPORTS

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

OPEN

# First evidence of the presence and activity of archaeal C3 group members in an Atlantic intertidal mudflat

Received: 26 January 2018

Accepted: 18 July 2018

Published online: 07 August 2018

Céline Lavergne<fnref n="1" />,<fnref n="2" />, Mylène Hugoni<fnref n="3" />, Christine Dupuy<fnref n="1" /> & Hélène Agogué<fnref n="1" />

The phylogenetic assignment of archaeal communities is constantly evolving, and the recent discovery of new phyla that grouped into superphyla has provided novel insights into archaeal ecology and evolution in ecosystems. In intertidal sediments, archaea are known to be involved in key functional processes such as organic matter turnover, but the ecological relevance of the rarest archaeal groups is poorly investigated, due partly to the lack of cultivated members. The high resolution of microbial diversity provided by high-throughput sequencing technologies now allows the rare biosphere to be described. In this work, we focused on the archaeal C3 group, showing that this phylum is not only present (at the DNA level) independently of sediment depth but also active (at the RNA level) in specific sediment niches depending on vertical physicochemical gradients. Moreover, we highlight the ambiguous phylogenetic affiliation of this group, indicating the need of further research to get new insights into the role of the C3 group.

Recent developments in metagenomic approaches have led to a clear demonstration of the existence of a rare biosphere, defined as low-abundance taxa present in ecosystems (accounting for 0.1 to 0.001% of sequence abundance depending on the study)<fnref n="1" />,<fnref n="2" />. While previous works had focused on eukaryota<fnref n="3" />,⁴ and bacteria⁵,⁶, only a few studies have considered the rare archaeal biosphere⁷,⁸. However, recent work has demonstrated the existence of different fractions among the rare archaeal biosphere. Some of these rare archaea are taxonomically close to abundant members, and some are always active, while others are inactive, which might represent dormant microorganisms that are able to react to seasonal fluctuations and serve as a local seed bank. A third fraction consists of rare archaea that are uncommon in public databases, inactive, and alien to the studied ecosystem, representing a non local seed bank of potential colonizers⁷. These rare taxa could therefore represent a gene reservoir and might play key roles in biogeochemical cycles⁷,⁹.

In intertidal mudflats, sharp vertical gradients in salinity, oxidative state and nutrients¹⁰ shape microbial establishment (at the DNA level)¹¹ as well as microbial activity (at the RNA level)¹². While archaea represent a small fraction of the active prokaryotic community compared with their bacterial counterparts, vertical shifts in major archaeal phyla have been recorded¹². The taxonomic affiliation of several of these rare archaeal taxa that have been newly discovered in many environmental samples is still under active debate¹³,¹⁴. In the present work, we demonstrated that the potentially active archaeal fraction retrieved from coastal sediments was enriched in a classically non-abundant phylum representing an emergent class of *Bathyarchaeota*. These findings raise many questions about the physiology of this group and their adaptation to this peculiar environment.

We were interested in the ecological relevance of the archaeal C3 group, which is still poorly described, as none of its members have yet been cultivated. This group was formerly referred to as subclass MCG-15 of *Bathyarchaeota* (previously known as the *Miscellaneous Crenarchaeotic Group*, MCG)¹⁵ and has been reported to be widespread in anoxic sediments¹⁶, hydrothermal vents¹⁷, and lakes¹⁸.

<footnote n="1">¹Université de La Rochelle - CNRS, UMR 7266, LIENS, 2 rue Olympe de Gouges, 17000, La Rochelle, France.</footnote>

<footnote n="2">²Escuela de Ingeniería Bioquímica, Pontificia Universidad Católica Valparaíso, Avenida Brasil 2085, Valparaíso, Chile.</footnote>

<footnote n="3">³Université Lyon 1 - UMR CNRS 5557/INRA 1418, Ecologie Microbienne, Villeurbanne, France. Céline Lavergne and Mylène Hugoni contributed equally to this work. Correspondence and requests for materials should be addressed to C.L. (email: celine.lavergne@pucv.cl)</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Relative abundance of archaeal 16S rRNA genes and transcripts in coastal sediments. The left panel shows the relative abundance of archaeal 16S rRNA genes recovered *in situ* on the 5th of July 2012 in an intertidal mudflat of Marennes-Oléron Bay. The right panel shows the relative abundance of archaeal 16S rRNA transcripts recovered in a tidal mesocosm using sediments from the same Marennes-Oléron mudflat. The relative abundance of the 16S rRNA genes or transcripts and the affiliation of the total archaeal OTUs are presented at the class level among the five layers below the sediment surface (bsf).</figcaption></figure>

# Material and Methods

**Sampling and nucleic acid extraction.** Marennes-Oléron Bay is located in the center of the French Atlantic coast. Sampling was performed along ridges at low tide. Triplicate 15-cm-diameter cores were sliced into five layers from 0 to 10 cm below the sediment surface (bsf) then homogenized and stored at -20 °C. DNA was extracted from ~1 g of thawed sediment using the UltraClean Soil kit® (MOBIO, USA) according to the manufacturer's recommendations to maximize yields.

As described in detail in a previous study¹², RNA extracts were obtained from five sediment layers of a tidal mesocosm using the Power Soil™ Total RNA Isolation Kit (MOBIO, USA). cDNAs were generated using SuperScript III and random hexamers (Invitrogen, Life Technologies, USA).

**High-throughput sequencing.** In the current study, three sets of metabarcoding data were used. First, a set of 16S rRNA genes (DNA level) and a set of 16S rRNA transcripts (RNA level), sequenced via 454 pyrosequencing, as partially described in a previous study¹², were used to determine the distribution of archaeal taxa in relation to sediment depth. The 454 pyrosequencing procedure is fully described in previous work¹². Second, to confirm the enrichment of archaeal C3 group sequences in specific niches and obtain longer and more robust sequences for phylogenetic analysis, we re-sequenced the two enriched samples (from 5 to 10 cm below sediment surface) using the Illumina HiSeq platform. Briefly, the Illumina HiSeq approach consisted of amplification of the V3–V5 region of the archaeal 16S rRNA genes using the 519F and 915R primers¹⁸. High-throughput sequencing was achieved after a multiplexing step using HiSeq Rapid Run 300 bp PE technology in an Illumina HiSeq2500 system (GATC Biotech, Konstanz, Germany). Archaeal 16S rRNA paired-end reads were merged with a maximum of 10% mismatches in the overlap region using FLASH¹⁹. Denoising procedures consisted of discarding reads outside of the expected length range (i.e., expected size between 370 and 580 bp) and reads containing ambiguous bases (N). After dereplication, the sequences were clustered into operational taxonomic units (OTUs) using SWARM²⁰ with a local clustering threshold. In the present work, the aggregation distance was equal to 3. Chimeras were then removed using VSEARCH²¹, and low-abundance sequences were filtered at 0.005% (i.e., OTUs with at least 0.005% of all sequences were retained²²), discarding singletons from the datasets. Taxonomic affiliation was evaluated using both RDP Classifier²³ and BLASTN+²⁴ against the 119 SILVA database²⁵. This procedure was automated in the FROGS pipeline²⁶. Raw data are available under Sequence Read Archive (SRA) format within the BioProject PRJNA477428.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Phylogenetic tree constructed with the most abundant C3 group-affiliated OTUs (representing 82 sequences) retrieved in our study. Sequences from the current study are highlighted in green in the tree and the circles indicated the number of OTUs they represent. The tree was built using neighbor-joining calculations in ARB at 423 positions and applying PVP filtering and 1000 bootstraps. Bootstrap values >50% are displayed.</figcaption></figure>

**Phylogenetic analysis.** From the high-throughput sequencing step using the Illumina HiSeq platform, 82 archaeal 16S rRNA gene sequences associated with the C3 group were obtained. All the sequences were analyzed using the ARB software package²⁶ (version July 2014) and the corresponding SILVA SSURef 99 database²⁷, version 128, released on 07.09.2016. All the sequences were automatically aligned using the SINA tool online²⁸. The aligned sequences were added to the tree of the SILVA database. Tree reconstruction was first performed with 241 full 16S rRNA gene sequences using two methods: neighbor joining (NJ, ARB) and maximum likelihood (ML, RAxML²⁹ v8.0.24 using the CIPRES gateway³⁰) to analyze the topology of the tree. Only the 25 most abundant OTUs (Supplementary Information, Table S1) were selected for the final tree. The final tree was computed with 129 sequences (bacteria *Tenericutes* as the outgroup, the Ancient archaeal group, *Lokiarchaeota*, *Thorarchaeota*, *Thaumarchaeota* Marine Group I, *Bathyarchaeota* and the C3 group) using the neighbor-joining method (Jukes Cantor distance, bootstrapping: 1000 replications) and applying the archaeal positional variability by parsimony (PVP) filter to partial sequences (423 positions). For clarity, only selected subsets of the sequences used to build the tree are shown in the figure.

**Results and Discussion**

*Bathyarchaeota* is a recently identified phylum composed of 17 subclasses with a diverse metabolic spectrum, mainly reported as anaerobic³¹. In the last two years, two partial genomes have been obtained for the C3 group (MCG-15 subclass): E09³² and B23³³. These findings agreed with the capacity of the C3 group to fix inorganic CO₂ via the reductive acetyl-CoA (Wood-Ljungdahl, WL) pathway³³. The C3 group is the first archaeal taxa capable of homoacetogenesis³³, as its members use H₂ as an electron donor to form acetate autotrophically³¹ via a more energy-efficient pathway than bacteria. This lower efficiency of bacteria is due to their requirement for ATP to generate formate³³. Thus, these archaea appear to be better adapted to low-energy environments³³. Although previous work indicated that members of *Bathyarchaeota* are also implicated in the methane cycle and could carry out methanogenesis³⁴, only subclass MCG-8 appears to exhibit this peculiarity thus far. The reported partial genomes for the C3 group did not contain a methyl coenzyme M reductase gene (*mcrA* gene) related to methanogenesis. Recent findings suggest that all *Bathyarchaeota* may have been methanogens in the past but that most of them appear to have lost this capacity³⁵. The authors of this study also suggest that subclass MCG-8 may be the first methanogens to lose methanogenesis. The metabolic diversity of *Bathyarchaeota* makes it difficult to appreciate the ecological role of each *Bathyarchaeota* subclass, and more environmental records and information are needed to better understand one of the most important benthic archaeal phyla.

In the current study, we demonstrated the ubiquity of the archaeal C3 group at different vertical sediment depths in an Atlantic intertidal mudflat through an *in situ* study conducted at the DNA level. The abundance of this group was evaluated along a vertical gradient from 0 to 10 cm below sediment surface (bsf), and this abundance was stable throughout the sediment cores, comprising up to 13% of the obtained 16S rRNA archaeal gene sequences (Fig. 1). Despite the stability in terms of 16S rRNA gene sequence numbers at different sediment depths, a mesocosm experiment corresponding to the same sampling sites and depths clearly demonstrated for the first time that this rare group is highly active in mudflat sediments in particular (at the RNA level, Fig. 1). Indeed, the diversity of archaeal 16S rRNA transcripts showed that a considerable proportion of archaeal C3 group sequences varied greatly according to depth and associated vertical gradients (i.e., 13% of total archaeal
16S rRNA transcript sequences were found from 2 to 5 cm bsf, and up to 70% were in the deeper zone, from 5 to 10 cm bsf; Fig. 1). These findings suggest the presence of ecological niches that favor the activity of these archaea, although they remain rare along the vertical gradient. These findings were possible because of the use of archaeal 16S rRNA transcript sequencing, whereas most other related studies have described only microbial diversity at the 16S rRNA gene level or have lacked this specific focus on archaeal communities, thus, potentially overlooking active communities and archaeal diversity.

While the C3 group is considered a class of *Bathyarchaeota*¹⁵, the SILVA database considers C3 as an emergent group that is more closely related to *Thaumarchaeota*. To assess the phylogenetic position of the C3 group transcripts retrieved in the present work, we re-sequenced the two enriched samples previously identified through 454 pyrosequencing (from 5 to 10 cm bsf) using the Illumina HiSeq 2* 300 bp platform to increase sequencing depth and quality. The phylogenetic analysis (Fig. 2) of the 82 sequences of C3 group-affiliated 16S rRNA transcripts retrieved from Marennes-Oléron bay confirmed that the C3 group is closely related to *Bathyarchaeota*.

This phylogenetic analysis separated the C3 OTUs obtained in the current study into two different groups. The first contained environmental sequences from the subseafloor, shallow sediment and extreme environments such as the Napoli mud volcano, while the second group might correspond to a coastal group. Interestingly, 16 OTUs were placed in the tree near benthic environmental sequences originating in the Aber-Benoît tidal basin, which is located in a neighboring region of the current study site³⁶.

The relationship of the archaeal C3 group with both *Thaumarchaeota* and *Bathyarchaeota* as a part of the TACK superphylum is clear, but whether the C3 group is a class of *Bathyarchaeota* cannot be determined using only the V3–V4 region of the 16S rRNA gene marker. A metagenomic and metatranscriptomic survey in C3 group-enriched ecosystems would help to reveal the metabolic functions of this mysterious, non-abundant archaeal group. Thus, we highlight the necessity of (1) exploring potentially active biodiversity at the RNA level to detect key ecological players among nonabundant taxa; and (2) studying further the C3 group, as it could be a key taxon for better understanding biogeochemical processes in coastal sediments.

## References

1. Pedrós-Alió, C. The rare bacterial biosphere. *Ann. Rev. Mar. Sci.* 4, 449–466 (2012).

2. Lynch, M. D. J. & Neufeld, J. D. Ecology and exploration of the rare biosphere. *Nat. Rev. Microbiol.* 13, 217–229 (2015).

3. Bachy, C. & Worden, A. Z. Microbial Ecology: Finding Structure in the Rare Biosphere. *Curr. Biol.* 24, R315–R317 (2014).

4. Logares, R. et al. Patterns of rare and abundant marine microbial eukaryotes. *Curr. Biol. CB* 24, 813–821 (2014).

5. Campbell, B. J., Yu, L., Heidelberg, J. F. & Kirchman, D. L. Activity of abundant and rare bacteria in a coastal ocean. *Proc. Natl. Acad. Sci. USA* 108, 12776–12781 (2011).

6. Sogin, M. L. et al. Microbial diversity in the deep sea and the underexplored "rare biosphere". *Proc. Natl. Acad. Sci.* 103, 12115–12120 (2006).

7. Hugoni, M. et al. Structure of the rare archaeal biosphere and seasonal dynamics of active ecotypes in surface coastal waters. *Proc. Natl. Acad. Sci.* 110, 6004–6009 (2013).

8. Galand, P. E., Casamayor, E. O., Kirchman, D. L. & Lovejoy, C. Ecology of the rare microbial biosphere of the Arctic Ocean. *Proc. Natl. Acad. Sci. USA* 106, 22427–22432 (2009).

9. Wang, Y. et al. Quantifying the Importance of the Rare Biosphere for Microbial Community Response to Organic Pollutants in a Freshwater Ecosystem. *Appl. Environ. Microbiol.* 83 (2017).

10. Herlemann, D. P. et al. Transitions in bacterial communities along the 2000 km salinity gradient of the Baltic Sea. *ISME J.* 5, 1571–1579 (2011).

11. Böer, S. I. et al. Time- and sediment depth-related variations in bacterial diversity and community structure in subtidal sands. *ISME J.* 3, 780–791 (2009).

12. Lavergne, C. et al. Diel Rhythm Does Not Shape the Vertical Distribution of Bacterial and Archaeal 16S rRNA Transcript Diversity in Intertidal Sediments: a Mesocosm Study. *Microb. Ecol.*, https://doi.org/10.1007/s00248-017-1048-1 (2017).

13. Castelle, C. J. et al. Genomic expansion of domain archaea highlights roles for organisms from new phyla in anaerobic carbon cycling. *Curr. Biol. CB* 25, 690–701 (2015).

14. Zaremba-Niedzwiedzka, K. et al. Asgard archaea illuminate the origin of eukaryotic cellular complexity. *Nature* 541, 353–358 (2017).

15. Meng, J. et al. Genetic and functional properties of uncultivated MCG archaea assessed by metagenome and gene expression analyses. *ISME J* 8, 650–659 (2014).

16. Kubo, K. et al. Archaea of the Miscellaneous Crenarchaeotal Group are abundant, diverse and widespread in marine sediments. *ISME J* 6, 1949–1965 (2012).

17. Wemheuer, B., Taube, R., Akyol, P., Wemheuer, F. & Daniel, R. Microbial diversity and biochemical potential encoded by thermal spring metagenomes derived from the Kamchatka peninsula. *Archaea* 2013, 13 (2013).

18. Hugoni, M. et al. Temporal dynamics of active Archaea in oxygen-depleted zones of two deep lakes. *Environ. Microbiol. Rep.* 7, 321–9 (2015).

19. Magoč, T. & Salzberg, S. L. FLASH: fast length adjustment of short reads to improve genome assemblies. *Bioinformatics* 27, 2957–2963 (2011).

20. Escudié, F. et al. FROGS: Find, Rapidly, OTUs with Galaxy Solution. *Bioinformatics*, https://doi.org/10.1093/bioinformatics/btx791 (2017).

21. Rognes, T., Flouri, T., Nichols, B., Quince, C. & Mahé, F. VSEARCH: a versatile open source tool for metagenomics. *PeerJ* 4, e2584 (2016).

22. Bokulich, N. A. et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. *Nat. Methods* 10, 57–59 (2013).

23. Wang, Q., Garrity, G. M., Tiedje, J. M. & Cole, J. R. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. *Appl. Environ. Microbiol.* 73, 5261–5267 (2007).

24. Camacho, C. et al. BLAST+: architecture and applications. *BMC Bioinformatics* 10, 421 (2009).

25. Pruesse, E. et al. SILVA: A comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. *Nucleic Acids Res.* 35, 7188–7196 (2007).

26. Ludwig, W. et al. ARB: a software environment for sequence data. *Nucleic Acids Res.* 32, 1363–1371 (2004).

27. Yilmaz, P. et al. The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. *Nucleic Acids Res.* 42, D643–D648 (2014).

28. Pruesse, E., Peplies, J. & Glöckner, F. O. SINA: Accurate high-throughput multiple sequence alignment of ribosomal RNA genes. *Bioinformatics* 28, 1823–1829 (2012).
29. Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. *Bioinformatics* **30**, 1312–1313 (2014).

30. Miller, M. A., Pfeiffer, W. & Schwartz, T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In *Gateway Computing Environments Workshop (GCE)*, 2010, 1–8, https://doi.org/10.1109/gce.2010.5676129 (2010).

31. Adam, P. S., Borrel, G., Brochier-Armanet, C. & Gribaldo, S. The growing tree of Archaea: new perspectives on their diversity, evolution and ecology. *ISME J.* https://doi.org/10.1038/ismej.2017.122 (2017).

32. Lloyd, K. G. et al. Predominant archaea in marine sediments degrade detrital proteins. *Nature* **496**, 215–218 (2013).

33. He, Y. et al. Genomic and enzymatic evidence for acetogenesis among multiple lineages of the archaeal phylum Bathyarchaeota widespread in marine sediments. *Nat. Microbiol.* **1**, 16035 (2016).

34. Evans, P. N. et al. Methane metabolism in the archaeal phylum Bathyarchaeota revealed by genome-centric metagenomics. *Science* **350**, 434–438 (2015).

35. Borrel, G., Adam, P. S. & Gribaldo, S. Methanogenesis and the Wood–Ljungdahl Pathway: An Ancient, Versatile, and Fragile Association. *Genome Biol. Evol.* **8**, 1706–1711 (2016).

36. Stauffert, M., Duran, R., Gassie, C. & Cravo-Laureau, C. Response of archaeal communities to oil spill in bioturbated mudflat sediments. *Microb. Ecol.* **67**, 108–119 (2014).

## Acknowledgements

This research was supported by a PhD grant from the Charente Maritime department to C.L., by the program CPER of Poitou-Charentes (2006–2013), and by the french national program EC2CO/CNRS-INSU (“CAPABIOC”, 2012–2014) to H.A. We acknowledge the Molecular Core Facilities at LIENS laboratory and are grateful to Jean-Christophe Auguet for his expertise. Authors also thank the editor members and the anonymous reviewers for their highly positive comments and suggestions.

## Author Contributions

The project was designed by C.D. and H.A.; C.D., H.A. and C.L. collected the samples, H.A. and C.L. extracted the nucleic acids and prepared the libraries, C.L. and M.H. analysed the data, C.L. generated the figures and M.H., C.L. and H.A. wrote the manuscript. All authors reviewed and approved the manuscript.

## Additional Information

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-018-30222-1.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018