# Diversity of Methanogens in Animals' Gut

Cheick Oumar Guindo, Bernard Davoust, Michel Drancourt, Ghiles Grine

► To cite this version:

Cheick Oumar Guindo, Bernard Davoust, Michel Drancourt, Ghiles Grine. Diversity of Methanogens in Animals' Gut. Microorganisms, 2021, 9 (1), pp.13. 10.3390/microorganisms9010013. hal-03216617

HAL Id: hal-03216617

https://amu.hal.science/hal-03216617

Submitted on 26 Jun 2023

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 Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Diversity of Methanogens in Animals' Gut

Cheick Oumar Guindo <fnref n="1" />,<fnref n="2" /> ID, Bernard Davoust <fnref n="2" /> ID, Michel Drancourt <fnref n="1" />,<fnref n="2" /> and Ghiles Grine <fnref n="2" />,<fnref n="3" />,*

**Abstract:** Methanogens are members of anaerobe microbiota of the digestive tract of mammals, including humans. However, the sources, modes of acquisition, and dynamics of digestive tract methanogens remain poorly investigated. In this study, we aimed to expand the spectrum of animals that could be sources of methanogens for humans by exploring methanogen carriage in animals. We used real-time PCR, PCR-sequencing, and multispacer sequence typing to investigate the presence of methanogens in 407 fecal specimens collected from nine different mammalian species investigated here. While all the negative controls remained negative, we obtained by PCR-sequencing seven different species of methanogens, of which three (*Methanobrevibacter smithii*, *Methanobrevibacter millerae* and *Methanomassiliicoccus luminyensis*) are known to be part of the methanogens present in the human digestive tract. *M. smithii* was found in 24 cases, including 12/24 (50%) in pigs, 6/24 (25%) in dogs, 4/24 (16.66%) in cats, and 1/24 (4.16%) in both sheep and horses. Genotyping these 24 *M. smithii* revealed five different genotypes, all known in humans. Our results are fairly representative of the methanogen community present in the digestive tract of certain animals domesticated by humans, and other future studies must be done to try to cultivate methanogens here detected by molecular biology to better understand the dynamics of methanogens in animals and also the likely acquisition of methanogens in humans through direct contact with these animals or through consumption of the meat and/or milk of certain animals, in particular cows.

**Keywords:** mammals' digestive tract; dynamics of methanogens; sources of methanogens and zoonotic methanogens

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

check for
updates

**Citation:** Guindo, C.O.; Davoust, B.; Drancourt, M.; Grine, G. Diversity of Methanogens in Animals' Gut. *Microorganisms* **2021**, *9*, 13. https://dx.doi.org/10.3390/microorganisms9010013

Received: 8 December 2020

Accepted: 20 December 2020

Published: 23 December 2020

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Methanogens are archaea characterized by their unique capability in producing methane from byproducts of bacterial anaerobe fermentations, being members of anaerobe microbiota of the digestive tract microbiota of several mammals [1]. Methanogens as strict anaerobes are classified to be limited to anoxic habitats. However, some studies have shown that some methanogens are able to produce methane in soils rich in oxygen [2] and even in human microbiota [3]. Methanogens were classified into three biochemical groups based on the substrates used for hydrogen production: hydrogenotrophic, aceticlastic, and methylotrrophic [2,4,5]. The group most described in human microbiota is hydrogenotrophic methanogens, which oxidize $H_2$, formate or a few simple alcohols, and reduce $CO_2$ to $CH_4$.

Accordingly, methanogens gained interest in the clinical microbiology over the past years after they were detected by PCR-based methods and cultured from the gut microbiota [6,7]; their translocation in milk and urines has been further observed [8]. Moreover, methanogens have been associated with dysbiosis such as in the case of vaginosis [9], urinary tract infections [10], and anaerobe abscesses of the brain [11,12], the muscle [13], the oral cavity in the case of periodontitis, and periimplantitis [14,15] in the case of refractory sinusitis [16]. Recently, we observed blood-borne methanogens associated with endocarditis [17]. In all these situations, anaerobe bacteria were associated in the methanogen-disease process, and this observation was probably reflecting methanogen specificities,

<footnote n="1">¹ IHU Méditerranée Infection, 13005 Marseille, France; cheicko86@gmail.com (C.O.G.); michel.drancourt@univ-amu.fr (M.D.)</footnote>

<footnote n="2">² IRD, MEPHI, Aix-Marseille Université, 13005 Marseille, France; bernard.davoust@gmail.com</footnote>

<footnote n="3">³ Faculty of Odontology, Aix-Marseille Université, 13005 Marseille, France</footnote>

<footnote>* Correspondence: grineghiles@univ-amu.fr; Tel.: +33-(0)4-13-73-24-01; Fax: +33-(0)-13-73-24-02</footnote>
including the absolute oxygen intolerance and the necessity of bacterial fermentative products to produce methane [4,18].

Currently, 16 different methanogens have been cultured from digestive-tract microbiota of animals [19–22], and PCR-based methods of detecting species-specific sequences traced an additional 4 species [21] (Table 1).

<table><caption>**Table 1.** Methanogens found in digestive tract microbiota of animals.</caption><thead><tr><th>Methanogens Species Obtained by Culture</th><th>Additional Species Detected by PCR-Based Methods But Not by Culture</th></tr></thead><tbody><tr><td>Methanosarcina sp.</td><td></td></tr><tr><td>Methanobacterium formicicum</td><td></td></tr><tr><td>Methanomicrobium mobile</td><td></td></tr><tr><td>Methanosarcina barkeri</td><td></td></tr><tr><td>Methanobacterium bryantii</td><td></td></tr><tr><td>Methanobrevibacter, ruminantium</td><td></td></tr><tr><td>Methanobrevibacter millerae</td><td>Methanobrevibacter smithii</td></tr><tr><td>Methanobrevibacter olleyae</td><td>Methanimicrococcus spp.</td></tr><tr><td>Methanoculleus olentangyi</td><td>Methanosphaera spp.</td></tr><tr><td>Methanobrevibacter woesei</td><td>Methanobacterium spp.</td></tr><tr><td>Methanobrevibacter gottschalkii</td><td></td></tr><tr><td>Methanobrevibacter thaueri</td><td></td></tr><tr><td>Methanobrevibacter wolinii</td><td></td></tr><tr><td>Methanobrevibacter cuticularis</td><td></td></tr><tr><td>Methanobrevibacter curvatus</td><td></td></tr><tr><td>Methanobrevibacter filiformi</td><td></td></tr></tbody></table>

The sources, modes of acquisition, and dynamics of digestive-tract methanogens remain poorly investigated. We previously reported that one-day newborns exhibited culturable *Methanobrevibacter smithii* (M. smithii) in the gastric fluid [23], suggesting a perinatal source of acquisition. Accordingly, we reported that mother milk did contain culturable M. smithii and culturable *Methanobrevibacter oralis* (M. oralis) [8]. Yet, it is unclear whether these one-day methanogens do persist along with the digestive tract of the newborns or whether this is just one of several waves of acquisition of methanogens along the first months of life [23–27]. Therefore, the search for methanogens sources other than mother milk is of interest.

Certain mammals (cow, sheep, donkey, horse, cat, pig, rabbit, rat, rhinoceros, baboon, monkey, and hippopotamus); birds (goose, turkey, and chicken) and insects (termites) are acknowledged to harbor digestive tract methanogens, and *M. smithii* in particular has already been detected from bovine and also from Wistar rats [21,22,28–43]. In this study, we aimed to expand the spectrum of animals that could be sources of methanogens for humans, by exploring methanogen carriage in animals.

## 2. Materials and Methods

### 2.1. Feces Samples

After the obtention of verbal consent from animals' owners, feces samples were collected from nine different animal species, namely cat, dog, horse, sheep, rabbit, cow, pig, goat, and donkey from animals living in metropolitan France, more precisely in the Marseille metropolitan area, Southeastern France (Table 2). Dogs and cats were fed industrial dry-kibble feed; horses were fed hay + straw + pellets; sheep and goats were fed pasture (grass) and dry supplementary feed; rabbits were fed dehydrated alfalfa + hay + pellets (other vegetables, cereals, mineral salts, and vitamins); cows were fed hay + straw + pasture (grass) and whole plant maize silage; pigs were fed straw + dry pelleted feed (formula consisting mainly of maize, wheat, oats, peas, soybeans, cereals, oilseeds, and minerals), and donkeys were fed grass and hay. Feces samples were stored at +4 °C for five weeks before being processed for DNA extraction as reported below.
<table><caption>**Table 2.** Details of 407 feces samples here investigated for the presence of methanogens.</caption><thead><tr><td>Origin of Samples</td><td>Species</td><td>Collected Number per Sample</td><td>Collection Sites</td></tr></thead><tbody><tr><td>Cat</td><td>Felis silvestris catus</td><td>105</td><td>Marseille</td></tr><tr><td>Dog</td><td>Canis lupus</td><td>52</td><td>Marseille</td></tr><tr><td>Horse</td><td>Equus caballus</td><td>89</td><td>Marseille and Carnoux</td></tr><tr><td>Sheep</td><td>Ovis aries</td><td>29</td><td>Bourganeuf</td></tr><tr><td>Rabbit</td><td>Oryctolagus cuniculus</td><td>2</td><td>Allauch</td></tr><tr><td>Cow</td><td>Bos taurus</td><td>57</td><td>Bourganeuf and Allauch</td></tr><tr><td>Pig</td><td>Sus scrofa domesticus</td><td>64</td><td>Avignon</td></tr><tr><td>Goat</td><td>Capra aegagrus hircus</td><td>5</td><td>Allauch</td></tr><tr><td>Donkey</td><td>Equus asinus</td><td>4</td><td>Allauch</td></tr></tbody></table>

## 2.2. DNA Extraction and PCR Assays

DNA extraction was performed by mixing 0.2 g of each feces sample with 500 µL of G2 buffer (QIAGEN, Hilden, Germany) in an Eppendorf tube (Fisher Scientific, Illkirch, France). Then, 0.3 g of acid-washed beads ≤106 µm (Sigma-Aldrich, Saint-Quentin Fallavier, France) was added in each tube and shaken in a FastPrep BIO 101 device (MP Biomedicals, Illkirch, France) for 45 s for mechanical lysis before 10-min incubation at 100 °C. A 180 µL volume of the mixture was then incubated with 20 µL of proteinase K (QIAGEN) at 56 °C overnight before a second mechanical lysis was performed. Total DNA was finally extracted with the EZ1 Advanced XL extraction kit (QIAGEN) and 200 µL eluted volume. Sterile phosphate-buffered saline (PBS) was used as a negative control in each DNA extraction run. Extracted DNA was incorporated into real-time PCR performed using Metha_16S_2_MBF: 5'-CGAACCGGATTAGATACCCG-3' and Metha_16S_2_MBR: 5'- CCCGCCAATTCCTTTAAGTT-3' primers (Eurogentec, Angers, France) and a FAM_Metha_16S_2_MBP 6FAM- CCTGGGAAGTACGGTCGCAAG probe targeting the 16S DNA gene of methanogens, designed in our laboratory (Eurogentec). PCR amplification was done in a 20 µL volume including 15 µL of mix and 5 µL of extracted DNA. Five µL of ultrapure water (Fisher Scientific) was used instead of DNA in the negative controls. The amplification reaction was performed in a CFX96 thermocycler (BioRad, Marnes-la-Coquette, France) incorporating a protocol with a cycle of 50 °C for 2-min, followed by 39 cycles of 95 °C for 5-min, 95 °C for 5 s and finally 60 °C for 30 s. The PCR-sequencing was done in a 20 µL volume, including 15 µL of mix and 5 µL of extracted DNA. Five µL of ultrapure water (Fisher Scientific, Illkirch, France) was used instead of DNA in the negative controls. The amplification reaction was performed in a CFX96 thermocycler (BioRad, Marnes-la-Coquette, France) incorporating a protocol with a cycle of 50 °C for 2-min, followed by 39 cycles of 95 °C for 5-min, 95 °C for 5 s and finally 60 °C for 30 s. Amplification of the archaeal 16S rRNA gene (primers used: SDArch0333aS15, 5-TCCAGGCCCTACGGG-3 and SDArch0958aA19, 5-YCCGGCGTTGAMTCCAATT-3) was performed as previously described [8,9,33,34]. Sequencing reactions (Sangers' method) were carried out using the BigDye Terminator, version 1.1, cycle sequencing kit DNA according to the manufacturer's instructions (Applied Biosystems, Foster City, USA). Nucleotide sequences were assembled using Chromas Pro software, version 1.7 (Technelysium Pty Ltd., Tewantin, Australia) and compared with sequences available in the GenBank database using the online NCBI BLAST program (http://blast.ncbi.nlm.nih.gov.gate1.inist.fr/Blast.cgi). We considered the sequences as belonging to the same species if the percentage of identity was >98.7%; as different species if between 95-98.7%, and different genera if this threshold was < 95% with respect to the first hit obtained by BLAST [44].
## 2.3. Multispacer Sequence Typing

We carried out a multispacer sequence typing (MST) technique on all fecal specimens positive by PCR-sequencing as previously described in our laboratory [23,45]. PCRs were realized in a 2720 Thermal Cycler (Applied Biosystems, Foster City, California, USA) and followed all the steps described for standard PCR used for the molecular analysis of fecal specimens. Negative controls consisting of PCR mixture without DNA template were included in each PCR run. All PCR products were sequenced in both directions using the same primers as used for PCRs in a 2720 Thermal Cycler (Applied Biosystems) with an initial 1-min denaturation step at 96 °C, followed by 25 cycles denaturation for 10 s each at 96 °C, a 20 s annealing step at 50 °C, and a 4-min extension step at 60 °C. Sequencing products were purified using the MultiScreen 96-well plates Millipore (Merck, Molsheim, France), containing 5% of Sephadex G-50 (Sigma-Aldrich), and sequences were analyzed on an ABI PRISM 31309 Genetic Analyzer (Applied Biosystems, Foster City, California, USA) and edited using the ChromasPro software (version 1.42; Technelysium Pty Ltd., Tewantin, Australia). For each intergenic spacer, a spacer type (ST) was defined as a sequence exhibiting unique genetic polymorphism (SNPs and indels). MST genotypes were defined as a unique combination of the four spacer sequences [23,45].

## 2.4. Phylogenetic Analyses

Sequences were edited using ChromasPro software (ChromasPro 1.7, Technelysium Pty Ltd., Tewantin, Australia). Molecular phylogenetic and evolutionary analyses were conducted in MEGA7 as previously described [46].

## 2.5. Statistical Analyses

We used R software for data analysis (https://www.r-project.org/). The Chi 2 test was used to compare the prevalence between the different animal species with a threshold $\alpha = 0.05$.

# 3. Results

In this study, a total of 407 fecal specimens collected from nine different mammalian species were investigated by real-time PCR and PCR-sequencing for the presence of methanogens using primers targeting the broad-range archaeal 16S rRNA gene.

Firstly, incorporating the 16S rRNA archaeal gene PCR primers newly designed in our laboratory into real-time PCR, we detected the presence of methanogen DNA in all animals here investigated and none of the negative controls. We found that 100.0% of cat feces specimens were positive with Ct values of 33.51 ± 1.28; 78.8% of dog feces specimens were positive with Ct values of 27.71 ± 0.94; 84.4% of horse feces specimens were positive with Ct values of 25 ± 2.95; 96.6% of sheep feces specimens were positive with Ct values of 27.19 ± 3.11; 100% of rabbit feces specimens were positive with Ct values of 27.1 ± 1.36; 100% of cow feces specimens were positive with Ct values of 24.11 ± 1.94; 100% of pig feces specimens were positive with Ct values of 22.15 ± 2.75; 80% of goat feces specimens were positive with Ct values of 19.18 ± 2.46; and 100% of donkey feces specimens were positive with Ct values of 18.82 ± 1.44 (Table 3).

Secondly, sequencing the standard PCR products was used for the precise identification of methanogens at the genus and species levels in each sample. In cats, 50/105 successfully sequenced samples yielded 20 *Methanocorpusculum aggregans* (*M. aggregans*), 13 *Methanocorpusculum labreanum* (*M. labreanum*), 09 *Methanobrevibacter millerae* (*M. millerae*), 04 *M. smithii*, 02 *Methanobrevibacter thaueri* (*M. thaueri*), and 02 *Methanobrevibacter olleyae* (*M. olleyae*). In dogs, 30/52 successfully sequenced samples yielded 13 *M. labreanum*, 06 *M. smithii*, 05 *M. aggregans*, 03 *M. thaueri*, 02 *M. millerae*, and 01 *M. olleyae*. In horses, 24/89 successfully sequenced samples yielded 11 *M. aggregans*, 10 *M. olleyae*, 01 *M. smithii*, 01 *M. millerae*, and 01 *M. labreanum*. In sheep, 28/29 successfully sequenced samples yielded 23 *M. labreanum*, 03 *M. millerae*, 01 *M. smithii*, and 01 *M. aggregans*. In rabbits, 2/2 successfully sequenced samples yielded 02 *M. thaueri*. In cows, 44/57 successfully sequenced samples
yielded 22 *M. aggregans*, 11 *M. millerae*, 06 *M. labreanum*, and 05 *M. thaueri*. In pigs, 25/64 successfully sequenced samples yielded 12 *M. smithii*, 09 *M. millerae*, 03 *Methanomassiliicoccus luminiensis* (*M. luminiensis*), and 01 *M. olleyae*. In goats, 4/5 successfully sequenced samples yielded 03 *M. labreanum* and 01 *M. aggregans*. Finally, in donkeys, 4/4 successfully sequenced samples yielded 04 *M. aggregans* (Table 4).

<table><caption>**Table 3.** Comparison of prevalence based on real-time PCR between animal species.</caption><thead><tr><th>Animal Species</th><th>Number of Samples Analyzed</th><th>Number of Positive Samples by RT-PCR</th><th>Prevalence [IC 95%]</th><th>p-Value</th></tr></thead><tbody><tr><td>Cat</td><td>105</td><td>105</td><td>100.0 [96.5–100.0]</td><td rowspan="8">9.9 × 10<sup>-8</sup></td></tr><tr><td>Dog</td><td>52</td><td>41</td><td>78.8 [65.3–88.9]</td></tr><tr><td>Horse</td><td>89</td><td>75</td><td>84.4 [75.0–91.1]</td></tr><tr><td>Sheep</td><td>29</td><td>28</td><td>96.6 [82.2–99.9]</td></tr><tr><td>Rabbit</td><td>2</td><td>2</td><td>100.0 [15.8–100.0]</td></tr><tr><td>Cow</td><td>57</td><td>57</td><td>100.0 [93.7–100.0]</td></tr><tr><td>Pig</td><td>64</td><td>64</td><td>100.0 [94.4–100.0]</td></tr><tr><td>Goat</td><td>5</td><td>4</td><td>80.0 [28.4–99.5]</td></tr><tr><td>Donkey</td><td>4</td><td>4</td><td>100.0 [39.8–100.0]</td></tr></tbody></table>

<table><caption>**Table 4.** Comparison of prevalence based on PCR-sequencing between animal species.</caption><thead><tr><th>Animal Species</th><th>Number of Samples Analyzed</th><th>Number of Positive Samples by PCR-Sequencing</th><th>Prevalence [IC 95%]</th><th>p-Value</th></tr></thead><tbody><tr><td>Cat</td><td>105</td><td>50</td><td>47.6 [37.8–57.6]</td><td rowspan="8">1.4 × 10<sup>-12</sup></td></tr><tr><td>Dog</td><td>52</td><td>30</td><td>57.7 [43.2–71.3]</td></tr><tr><td>Horse</td><td>89</td><td>24</td><td>27.0 [18.1–37.4]</td></tr><tr><td>Sheep</td><td>29</td><td>28</td><td>96.6 [82.2–99.9]</td></tr><tr><td>Rabbit</td><td>2</td><td>2</td><td>100.0 [15.8–100.0]</td></tr><tr><td>Cow</td><td>57</td><td>44</td><td>77.2 [64.2–87.3]</td></tr><tr><td>Pig</td><td>64</td><td>25</td><td>39.1 [27.1–52.1]</td></tr><tr><td>Goat</td><td>5</td><td>4</td><td>80.0 [28.4–99.5]</td></tr><tr><td>Donkey</td><td>4</td><td>4</td><td>100.0 [39.8–100.0]</td></tr></tbody></table>

We obtained a total of seven different species of methanogens in our study as illustrated with Venn diagrams (Figure 1).). The Venn diagram shows which species of methanogens are found in common in humans and in animal samples analyzed in this study, and which species of methanogens are found exclusively in animals and exclusively in humans. Indeed, three methanogens species (*M. smithii*, *M. millerae*, and *M. luminyensis*) are known to be part of the methanogens present in the human digestive tract. The remaining four (*M. thaueri*, *M. olleyae*, *M. labreanum* and *M. aggregans*) are not known to date in humans. However, we did not find in our study the other 10 species of methanogens present in the human digestive tract, including *Methanobrevibacter arboriphilicus*, *M. oralis*, *Methanosphaera stadtmanae* (*M. stadtmanae*), *Candidatus Methanomethylophilus alvus* (*Ca. Methanomethylophilus alvus*), *Candidatus Methanomassiliicoccus intestinalis* (*Ca. Methanomassiliicoccus intestinalis*), *Methanoculleus chikugoensis* (*M. chikugoensis*), *Methanobacterium congolense* (*M. congolense*), *Methanoculleus bourgensis* (*M. bourgensis*), *Candidatus Nitrososphaera evergladensis* (*Ca. Nitrososphaera evergladensis*), and *Methanosarcinia mazei* (*M. mazei*).
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Venn diagram between the methanogens found in our study and those known from the human digestive tract. This Venn diagram shows which species of methanogens are found in common in humans and in animal samples analyzed in this study and which species of methanogens are found exclusively in animals and exclusively in humans.</figcaption></figure>

Among the 211 sequences obtained, 153 (72.51%) of them have an identity percentage greater than 99%, 43 (20.37%) have an identity percentage lower than 98.7%, and 15 (7.10%) have an identity percentage lower than 95% (Table 5). The phylogenetic trees of sequences obtained with a percentage identity lower than 98.7% and sequences with a percentage identity lower than 95% indicated new species and new genera, respectively (Figure 2 and Supplementary Figures). We obtained 24 *M. smithii* by PCR-sequencing including 12/24 (50%) in pigs, 6/24 (25%) in dogs, 4/24 (16.66%) in cats, and 1/24 (4.16%) in both sheep and horses. Genotyping the 24 *M. smithii* revealed five different genotypes. Genotype 1 was found in 8/24 (33.33%); genotype 2 in 10/24 (41.66%); genotype 3 in 4/24 (16.66%); and genotypes 4 and 5 in 1/24 (4.16%) each (Table 6).

<table><caption>**Table 5.** Percentage of identity among the sequences obtained.</caption><thead><tr><th>Animal Species</th><th>Percentage > 99%</th><th>Percentage < 98.7%</th><th>Percentage < 95%</th></tr></thead><tbody><tr><td>Cat</td><td>50</td><td>0</td><td>0</td></tr><tr><td>Dog</td><td>30</td><td>0</td><td>0</td></tr><tr><td>Horse</td><td>7</td><td>15</td><td>2</td></tr><tr><td>Sheep</td><td>18</td><td>7</td><td>3</td></tr><tr><td>Rabbit</td><td>2</td><td>0</td><td>0</td></tr><tr><td>Cow</td><td>29</td><td>12</td><td>3</td></tr><tr><td>Pig</td><td>16</td><td>4</td><td>5</td></tr><tr><td>Goat</td><td>0</td><td>2</td><td>2</td></tr><tr><td>Donkey</td><td>1</td><td>3</td><td>0</td></tr></tbody></table>
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Molecular phylogenetic analysis, based on 16S rRNA partial gene, showed the position of *Methanomassiliicoccus-like* sequences detected in feces of pigs. The evolutionary history was inferred using the Neighbor-Joining method. The optimal tree with the sum of branch length = 0.82722721 is shown. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1.000 replicates) are shown next to the branches. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Maximum Composite Likelihood method and are in the units of the number of base substitutions per site. The analysis involved 14 nucleotide sequences. All positions containing gaps and missing data were eliminated. There were a total of 415 positions in the final dataset. Evolutionary analyses were conducted in MEGA7. Bootstrap values ≥ 95% are indicated at nodes. In red: out of group. Species highlighted: methanogens species detected in this study.</figcaption></figure>

<table><caption>**Table 6.** Summary of the results of multispacer sequence typing.</caption><thead><tr><th rowspan="2">Samples</th><th rowspan="2">Origin</th><th rowspan="2">Collection Sites</th><th colspan="5">Genotypes</th></tr><tr><th>Spacer 1</th><th>Spacer 2</th><th>Spacer 3</th><th>Spacer 4</th><th>Spacer Type *</th></tr></thead><tbody><tr><td>1</td><td>Sheep</td><td>Bouganeuf</td><td></td><td></td><td>×</td><td></td><td>1</td></tr><tr><td>2</td><td>Horse</td><td>Marseille</td><td></td><td></td><td>×</td><td></td><td>1</td></tr><tr><td>3</td><td>Pig</td><td>Avignon</td><td>×</td><td>×</td><td>×</td><td>×</td><td>2</td></tr><tr><td>4</td><td>Pig</td><td>Avignon</td><td>×</td><td>×</td><td>×</td><td>×</td><td>2</td></tr><tr><td>5</td><td>Pig</td><td>Avignon</td><td>×</td><td>×</td><td>×</td><td>×</td><td>2</td></tr><tr><td>6</td><td>Pig</td><td>Avignon</td><td>×</td><td>×</td><td>×</td><td>×</td><td>2</td></tr><tr><td>7</td><td>Pig</td><td>Avignon</td><td>×</td><td>×</td><td>×</td><td></td><td>3</td></tr><tr><td>8</td><td>Pig</td><td>Avignon</td><td>×</td><td>×</td><td></td><td></td><td>4</td></tr><tr><td>9</td><td>Pig</td><td>Avignon</td><td></td><td>×</td><td>×</td><td></td><td>5</td></tr><tr><td>10</td><td>Pig</td><td>Avignon</td><td></td><td></td><td>×</td><td></td><td>1</td></tr><tr><td>11</td><td>Pig</td><td>Avignon</td><td>×</td><td>×</td><td>×</td><td>×</td><td>2</td></tr><tr><td>12</td><td>Pig</td><td>Avignon</td><td>×</td><td>×</td><td>×</td><td></td><td>3</td></tr><tr><td>13</td><td>Pig</td><td>Avignon</td><td>×</td><td>×</td><td>×</td><td>×</td><td>2</td></tr><tr><td>14</td><td>Pig</td><td>Avignon</td><td>×</td><td>×</td><td>×</td><td>×</td><td>2</td></tr><tr><td>15</td><td>Dog</td><td>Marseille</td><td>×</td><td>×</td><td>×</td><td>×</td><td>2</td></tr><tr><td>16</td><td>Dog</td><td>Marseille</td><td>×</td><td>×</td><td>×</td><td>×</td><td>2</td></tr><tr><td>17</td><td>Dog</td><td>Marseille</td><td>×</td><td>×</td><td>×</td><td></td><td>3</td></tr><tr><td>18</td><td>Dog</td><td>Marseille</td><td>×</td><td>×</td><td>×</td><td></td><td>3</td></tr><tr><td>19</td><td>Dog</td><td>Marseille</td><td>×</td><td>×</td><td>×</td><td>×</td><td>2</td></tr><tr><td>20</td><td>Dog</td><td>Marseille</td><td></td><td></td><td>×</td><td></td><td>1</td></tr><tr><td>21</td><td>Cat</td><td>Marseille</td><td></td><td></td><td>×</td><td></td><td>1</td></tr><tr><td>22</td><td>Cat</td><td>Marseille</td><td></td><td></td><td>×</td><td></td><td>1</td></tr><tr><td>23</td><td>Cat</td><td>Marseille</td><td></td><td></td><td>×</td><td></td><td>1</td></tr><tr><td>24</td><td>Cat</td><td>Marseille</td><td></td><td></td><td>×</td><td></td><td>1</td></tr></tbody></table>

Genotyping these *M. smithii* strains revealed the presence of five different genotypes. * Spacer type was determined according to reference [23,45].
## 4. Discussion

It is known and published that methanogens colonize the gastrointestinal tract of certain mammals, particularly herbivorous ones [35]. Most methanogens identified in mammals belong to the phylum Euryarchaeota, with a high percentage of the species *M. smithii* [36], a species being the most prevalent one in humans [47]. Our report is the largest one showing the presence of methanogens in nine mammals in the same study. Our results confirmed the published data on the presence of methanogens in the digestive tract of cats, dogs, horses, cows, sheep, rabbits, goats, pigs, and donkeys [21, 22, 28–34, 36–38]. In addition, all methanogens found in this study belong to the phylum Euryarchaeota, which is in accordance with the results obtained in studies conducted on the human digestive tract [6,48]. Our results give an insight on the concentration of methanogens present in the intestinal microbiota of each animal species analyzed and on the prevalence of methanogens in domestic animals by humans.

The results of the analysis of the 16S RNA sequences obtained from our samples show that there is a real diversity of methanogenic archaea genera (*Methanosphaera*, *Methanocorpusculum*, *Methanocalculus*, *Methanoculleus*, *Methanogenium*, *Methanoplanus*, *Methanolacinia*, *Methanobacterium*, *Methanomicrobium*, *Methanomassiliicoccus* and *Methanobrevibacter*) in the digestive tract of animals (cats, dogs, horses, sheep, cows, rabbits, goats, pigs, and donkeys) as in humans [6,48]. All sequences with a percentage lower than 98.7% have been deposited in the GenBank database (accession no MT587812 to MT587864) and EBI database (accession no MT793590; MT819603; MT822292; MT822293; and MT822482).

*Methanomassiliicoccus luminiyensis* was known to colonize the human digestive tract, and it has never been detected in animals’ digestive tracts [46]. For the first time, this study demonstrated the presence of the species *M. luminiyensis* in pigs and not in the other animals investigated here. These results could be explained by the fact that the pig is an omnivore, which means that its diet is close to that of humans compared to other animals. In addition, 50% of *M. smithii* in our study was found in pigs, indicating that *M. smithii* was the most prevalent methanogen in the digestive tract of pigs, consistent with work carried out in humans where the high prevalence of *M. smithii* in the digestive tract has been demonstrated [45,49].

These results are representative of the methanogen community present in the digestive tract of certain animals domesticated by humans, and other future studies must be done to try to cultivate methanogens detected here by molecular biology to better understand the dynamics of methanogens in animals. The possible ways of methanogens’ acquisition in humans could be contact with animals and/or through consumption of milk/dairy products of certain animals, in particular cows, since a recent study demonstrated an association between the acquisition of *M. smithii* in children and the consumption of dairy products [50].

**Supplementary Materials:** The following are available online at https://www.mdpi.com/2076-2607/9/1/13/s1.

**Author Contributions:** C.O.G.: Samples collection, study design, manipulations, data information, results interpretation, draft writing. B.D.: samples collection, data information. M.D.: Supervision, study design, results interpretation, draft writing. G.G.: Grant acquisition, supervision, study design, samples collection, results interpretation, draft writing. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by a grant of GDR Archaea, Centre National de la Recherche Scientifique, Paris, France.

**Acknowledgments:** We thank GDR Archaea, Centre National de la Recherche Scientifique, Paris, France for its financial support.

**Conflicts of Interest:** All the authors declare that there is no conflict of interest.
## References

1. Michelland, R.J.; Monteils, V.; Combes, S.; Cauquil, L.; Gidenne, T.; Fortun-Lamothe, L. Comparison of the archaeal community in the fermentative compartment and faeces of the cow and the rabbit. *Anaerobe* **2010**, *16*, 396–401. [CrossRef]

2. Lyu, Z.; Lu, Y. Metabolic shift at the class level sheds light on adaptation of methanogens to oxidative environments. *ISME J.* **2018**, *12*, 411–423. [CrossRef] [PubMed]

3. Grine, G.; Terrer, E.; Boualam, M.A.; Aboudharam, G.; Chaudet, H.; Ruimy, R.; Drancourt, M. Tobacco-smoking-related prevalence of methanogens in the oral fluid microbiota. *Sci. Rep.* **2018**, *8*, 9197. [CrossRef] [PubMed]

4. Sogodogo, E.; Drancourt, M.; Grine, G. Methanogens as emerging pathogens in anaerobic abscesses. *Eur. J. Clin. Microbiol. Infect. Dis.* **2019**, *38*, 811–818. [CrossRef] [PubMed]

5. Ruaud, A.; Esquivel-Elizondo, S.; de la Cuesta-Zuluaga, J.; Waters, J.L.; Angenent, L.T.; Youngblut, N.D.; Ley, R.E. Syntrophy via interspecies H₂ transfer between *Christensenella* and *Methanobrevibacter* underlies their global cooccurrence in the human gut. *mBio* **2020**, *11*, e03235-19. [CrossRef]

6. Nkamga, V.D.; Henrissat, B.; Drancourt, M. Archaea: Essential inhabitants of the human digestive microbiota. *Hum. Microbiome J.* **2017**, *3*, 1–8. [CrossRef]

7. Guindo, C.O.; Drancourt, M.; Grine, G. Digestive tract methanodrome: Physiological roles of human microbiota-associated methanogens. *Microb. Pathog.* **2020**, *149*, 104425. [CrossRef]

8. Togo, A.H.; Grine, G.; Khelaifia, S.; des Robert, C.; Brevaut, V.; Caputo, A.; Baptiste, E.; Bonnet, M.; Levasseur, A.; Drancourt, M.; et al. Culture of methanogenic archaea from human colostrum and milk. *Sci. Rep.* **2019**, *9*, 18653. [CrossRef]

9. Grine, G.; Drouet, H.; Fenollar, F.; Bretelle, F.; Raoult, D.; Drancourt, M. Detection of Methanobrevibacter smithii in vaginal samples collected from women diagnosed with bacterial vaginosis. *Eur. J. Clin. Microbiol. Infect. Dis.* **2019**, *38*, 1643–1649. [CrossRef]

10. Grine, G.; Lotte, R.; Chirio, D.; Chevalier, A.; Raoult, D.; Drancourt, M.; Ruimy, R. Co-culture of Methanobrevibacter smithii with enterobacteria during urinary infection. *EBioMedicine* **2019**, *43*, 333–337. [CrossRef]

11. Drancourt, M.; Nkamga, V.D.; Lakhe, N.A.; Régis, J.-M.; Dufour, H.; Fournier, P.-E.; Bechah, Y.; Michael Scheld, W.; Raoult, D. Evidence of archaeal methanogens in brain abscess. *Clin. Infect. Dis.* **2017**, *65*, 1–5. [CrossRef] [PubMed]

12. Nkamga, V.D.; Lotte, R.; Chirio, D.; Lonjon, M.; Roger, P.-M.; Drancourt, M.; Ruimy, R. Methanobrevibacter oralis detected along with Aggregatibacter actinomycetemcomitans in a series of community-acquired brain abscesses. *Clin. Microbiol. Infect.* **2018**, *24*, 207–208. [CrossRef] [PubMed]

13. Nkamga, V.D.; Lotte, R.; Roger, P.-M.; Drancourt, M.; Ruimy, R. Methanobrevibacter smithii and Bacteroides thetaiotaomicron cultivated from a chronic paravertebral muscle abscess. *Clin. Microbiol. Infect.* **2016**, *22*, 1008–1009. [CrossRef]

14. Faveri, M.; Gonçalves, L.F.H.; Feres, M.; Figueiredo, L.C.; Gouveia, L.A.; Shibli, J.A.; Mayer, M.P.A. Prevalence and microbiological diversity of Archaea in peri-implantitis subjects by 16S ribosomal RNA clonal analysis: Archaea in peri-implantitis subjects. *J. Periodontal Res.* **2011**, *46*, 338–344. [CrossRef]

15. Nguyen-Hieu, T.; Khelaifia, S.; Aboudharam, G.; Drancourt, M. Methanogenic archaea in subgingival sites: A review. *APMIS* **2013**, *121*, 467–477. [CrossRef]

16. Sogodogo, E.; Fellag, M.; Loukil, A.; Nkamga, V.D.; Michel, J.; Dessi, P.; Fournier, P.-E.; Drancourt, M. Nine cases of methanogenic archaea in refractory sinusitis, an emerging clinical entity. *Front. Public Health* **2019**, *7*, 38. [CrossRef]

17. Drancourt, M.; Djemai, K.; Gouriet, F.; Grine, G.; Loukil, A.; Bedotto, M.; Levasseur, A.; Lepidi, H.; Bou-Khalil, J.; Khelaifia, S.; et al. Methanobrevibacter smithii archaemia in febrile patients with bacteremia, including those with endocarditis. *Clin. Infect. Dis.* **2020**, ciaa998. [CrossRef]

18. Guindo, C.O.; Terrer, E.; Chabrière, E.; Aboudharam, G.; Drancourt, M.; Grine, G. Culture of salivary methanogens assisted by chemically produced hydrogen. *Anaerobe* **2020**, *61*, 102128. [CrossRef]

19. Miller, T.L. Description of *Methanobrevibacter gottschalkii* sp. nov., *Methanobrevibacter thaueri* sp. nov., *Methanobrevibacter woesei* sp. nov. and *Methanobrevibacter wolinii* sp. nov. *Int. J. Syst. Evol. Microbiol.* **2002**, *52*, 819–822.

20. Pol, A.; Demeyer, D.I. Fermentation of methanol in the sheep rumen. *Appl. Environ. Microbiol.* **1988**, *54*, 832–834. [CrossRef]

21. Janssen, P.H.; Kirs, M. Structure of the archaeal community of the rumen. *Appl. Environ. Microbiol.* **2008**, *74*, 3619–3625. [CrossRef] [PubMed]

22. Jarvis, G.N.; Strömpl, C.; Burgess, D.M.; Skillman, L.C.; Moore, E.R.B.; Joblin, K.N. Isolation and identification of ruminal methanogens from grazing cattle. *Curr. Microbiol.* **2000**, *40*, 327–332. [CrossRef] [PubMed]

23. Grine, G.; Boualam, M.A.; Drancourt, M. Methanobrevibacter smithii, a methanogen consistently colonising the newborn stomach. *Eur. J. Clin. Microbiol. Infect. Dis.* **2017**, *36*, 2449–2455. [CrossRef] [PubMed]

24. Mihajlovski, A.; Doré, J.; Levenez, F.; Alric, M.; Brugère, J.-F. Molecular evaluation of the human gut methanogenic archaeal microbiota reveals an age-associated increase of the diversity. *Environ. Microbiol. Rep.* **2010**, *2*, 272–280. [CrossRef]

25. Danielsson, R.; Schnürer, A.; Arthurson, V.; Bertilsson, J. Methanogenic population and CH₄ production in swedish dairy cows fed different levels of forage. *Appl. Environ. Microbiol.* **2012**, *78*, 6172–6179. [CrossRef]

26. Odamaki, T.; Kato, K.; Sugahara, H.; Hashikura, N.; Takahashi, S.; Xiao, J.; Abe, F.; Osawa, R. Age-related changes in gut microbiota composition from newborn to centenarian: A cross-sectional study. *BMC Microbiol.* **2016**, *16*, 90. [CrossRef]

27. Palmer, C.; Bik, E.M.; DiGiulio, D.B.; Relman, D.A.; Brown, P.O. Development of the human infant intestinal microbiota. *PLoS Biol.* **2007**, *5*, e177. [CrossRef]
28. Miller, T.L.; Wolin, M.J. Methanogens in human and animal intestinal Tracts. *Syst. Appl. Microbiol.* **1986**, *7*, 223-229. [CrossRef]

29. Miller, T.L.; Wolin, M.J.; Zhao, H.X.; Bryant, M.P. Characteristics of methanogens isolated from bovine rumen. *Appl. Environ. Microbiol.* **1986**, *51*, 201-202. [CrossRef]

30. Skillman, L.C.; Evans, P.N.; Strömpl, C.; Joblin, K.N. 16S rDNA directed PCR primers and detection of methanogens in the bovine rumen. *Lett. Appl. Microbiol.* **2006**, *42*, 222-228. [CrossRef]

31. Zhang, Y.; Zamudio Cañas, E.M.; Zhu, Z.; Linville, J.L.; Chen, S.; He, Q. Robustness of archaeal populations in anaerobic co-digestion of dairy and poultry wastes. *Bioresour. Technol.* **2011**, *102*, 779-785. [CrossRef] [PubMed]

32. Luo, Y.; Su, Y.; Wright, A.-D.G.; Zhang, L.; Smidt, H.; Zhu, W. Lean breed landrace pigs harbor fecal methanogens at higher diversity and density than obese breed erhualian pigs. *Archaea* **2012**, *2012*, 1-9. [CrossRef] [PubMed]

33. Hernandez-Sanabria, E.; Goonewardene, L.A.; Wang, Z.; Zhou, M.; Moore, S.S.; Guan, L.L. Influence of sire breed on the interplay among rumen microbial populations inhabiting the rumen liquid of the progeny in beef cattle. *PLoS ONE* **2013**, *8*, e58461. [CrossRef] [PubMed]

34. Rea, S.; Bowman, J.P.; Popovski, S.; Pimm, C.; Wright, A.-D.G. *Methanobrevibacter millerae* sp. nov. and *Methanobrevibacter olleyae* sp. nov., methanogens from the ovine and bovine rumen that can utilize formate for growth. *Int. J. Syst. Evol. Microbiol.* **2007**, *57*, 450-456. [CrossRef] [PubMed]

35. Murru, F.; Fliegerova, K.; Mura, E.; Mrázek, J.; Kopečný, J.; Moniello, G. A comparison of methanogens of different regions of the equine hindgut. *Anaerobe* **2018**, *54*, 104-110. [CrossRef]

36. Wright, A.-D.G.; Williams, A.J.; Winder, B.; Christophersen, C.T.; Rodgers, S.L.; Smith, K.D. Molecular diversity of rumen methanogens from sheep in western Australia. *Appl. Environ. Microbiol.* **2004**, *70*, 1263-1270. [CrossRef]

37. Lee, J.-H.; Kumar, S.; Lee, G.-H.; Chang, D.-H.; Rhee, M.-S.; Yoon, M.-H.; Kim, B.-C. *Methanobrevibacter boviskoreani* sp. nov., isolated from the rumen of Korean native cattle. *Int. J. Syst. Evol. Microbiol.* **2013**, *63*, 4196-4201. [CrossRef]

38. Tun, H.M.; Brar, M.S.; Khin, N.; Jun, L.; Hui, R.K.-H.; Dowd, S.E.; Leung, F.C.-C. Gene-centric metagenomics analysis of feline intestinal microbiome using 454 junior pyrosequencing. *J. Microbiol. Methods* **2012**, *88*, 369-376. [CrossRef]

39. Beauchemin, K.A.; Ungerfeld, E.M.; Eckard, R.J.; Wang, M. Review: Fifty years of research on rumen methanogenesis: Lessons learned and future challenges for mitigation. *Animal* **2020**, *14*, s2-s16. [CrossRef]

40. Qin, H.; Lang, H.; Yang, H. Characterization of the methanogen community in a household anaerobic digester fed with swine manure in China. *Appl. Microbiol. Biotechnol.* **2013**, *97*, 8163-8171. [CrossRef]

41. Saengkerdsub, S.; Ricke, S.C. Ecology and characteristics of methanogenic archaea in animals and humans. *Crit. Rev. Microbiol.* **2014**, *40*, 97-116. [CrossRef] [PubMed]

42. Whitford, M.F.; Teather, R.M.; Forster, R.J. Phylogenetic analysis of methanogens from the bovine rumen. *BMC Microbiol.* **2001**, *1*, 5. [CrossRef] [PubMed]

43. Mukumba, P.; Makaka, G.; Mamphweli, S. Anaerobic digestion of donkey dung for biogas production. *S. Afr. J. Sci.* **2016**, *112*, 1-4. [CrossRef]

44. Stackebrandt, E.; Ebers, J. Taxonomic parameters revisited: Tarnished gold standards. *Microbiol. Today* **2006**, *33*, 152-155.

45. Nkamga, V.D.; Huynh, H.T.T.; Aboudharam, G.; Ruimy, R.; Drancourt, M. Diversity of human-associated methanobrevibacter smithii isolates revealed by multispacer sequence typing. *Curr. Microbiol.* **2015**, *70*, 810-815. [CrossRef]

46. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. *Mol. Biol. Evol.* **2016**, *33*, 1870-1874. [CrossRef]

47. Dridi, B.; Henry, M.; El Khéchine, A.; Raoult, D.; Drancourt, M. High prevalence of *Methanobrevibacter smithii* and *Methanosphaera stadtmanae* detected in the human gut using an improved DNA detection protocol. *PLoS ONE* **2009**, *4*, e7063. [CrossRef]

48. Dridi, B.; Raoult, D.; Drancourt, M. Archaea as emerging organisms in complex human microbiomes. *Anaerobe* **2011**, *17*, 56-63. [CrossRef]

49. Sogodogo, E.; Doumbo, O.; Aboudharam, G.; Kouriba, B.; Diawara, O.; Koita, H.; Togora, S.; Drancourt, M. First characterization of methanogens in oral cavity in Malian patients with oral cavity pathologies. *BMC Oral Health* **2019**, *19*, 232. [CrossRef]

50. Van de Pol, J.A.A.; van Best, N.; Mbakwa, C.A.; Thijs, C.; Savelkoul, P.H.; Arts, I.C.W.; Hornef, M.W.; Mommers, M.; Penders, J. Gut colonization by methanogenic archaea is associated with organic dairy consumption in children. *Front. Microbiol.* **2017**, *8*, 35. [CrossRef]