Genomic convergence between Akkermansia muciniphila in different mammalian hosts
BMC Microbiology volume 21, Article number: 298 (2021)
Akkermansia muciniphila is a member of the human gut microbiota where it resides in the mucus layer and uses mucin as the sole carbon, nitrogen and energy source. A. muciniphila is the only representative of the Verrucomicrobia phylum in the human gut. However, A. muciniphila 16S rRNA gene sequences have also been found in the intestines of many vertebrates.
We detected A. muciniphila-like bacteria in the intestines of animals belonging to 15 out of 16 mammalian orders. In addition, other species belonging to the Verrucomicrobia phylum were detected in fecal samples. We isolated 10 new A. muciniphila strains from the feces of chimpanzee, siamang, mouse, pig, reindeer, horse and elephant. The physiology and genome of these strains were highly similar in comparison to the type strain A. muciniphila MucT. Overall, the genomes of the new strains showed high average nucleotide identity (93.9 to 99.7%). In these genomes, we detected considerable conservation of at least 75 of the 78 mucin degradation genes that were previously detected in the genome of the type strain MucT.
The low genomic divergence observed in the new strains may indicate that A. muciniphila favors mucosal colonization independent of the differences in hosts. In addition, the conserved mucus degradation capability points towards a similar beneficial role of the new strains in regulating host metabolic health.
The gastrointestinal (GI) tract of vertebrates is colonized with a dense and diverse microbiota . Several factors affect the gut microbiota composition of vertebrates including diet, host phylogeny and gut morphology . The microbiota has had a large influence on animal evolution, and can be seen as an obligate and beneficial symbiont . The main phyla representing the gut microbiota in mammals are Firmicutes, Bacteroidetes, Proteobacteria, Actinobacteria, and Verrucomicrobia . The gut microbiota produces short chain fatty acids from degradation of otherwise indigestible components providing the host with the ability to digest a wider variety of available foods [2, 4]. The gut microbiota also produces vitamins and other beneficial substances that the host cannot synthesize .
Some microbiota members can flourish within the mucus layer, a glycan-rich and anaerobic environment offered by its host , including mucus-degrading specialist Akkermansia muciniphila . A. muciniphila is the only representative of the Verrucomicrobia phylum in the human gut. Mucin utilization by A. muciniphila has been shown i) in vitro, where it grows on mucin as sole carbon and nitrogen source , ii) in vivo, where it scavenges mucin efficiently  and iii) in silico using a genome-scale model and omics analysis . Recent mouse and human studies have demonstrated that intake of A. muciniphila has a series of health benefits, including improved barrier function, increased insulin sensitivity and reduction of obesity [10, 11]. In addition, mouse experiments have shown that while A. muciniphila can degrade mucin, its presence increases mucus production, mucus layer thickness and tight junction protein production [10, 12, 13].
A. muciniphila is abundantly present in the human intestinal tract, varying from 1 to 4% of the bacterial population in the colon . Its abundance was found to be linked to a healthy status in humans . Interestingly, A. muciniphila has also been detected in other mammals, such as the brown bear . Furthermore, A. muciniphila has been detected in several small animals, such as the ground squirrel and Syrian hamster [17, 18]. In addition, Akkermansia spp. are widely spread in the GI tract throughout the animal kingdom, including mammals  and other vertebrates such as python [19, 20], zebra fish , chicken , and salmon . Next to bacteria belonging to the genus Akkermansia, other Verrucomicrobia are detected in the animal GI tract of invertebrates including termites , ants , earthworm  and nematodes . A few examples of Verrucomicrobia isolates from the gut environment are A. muciniphila , python isolate A. glycaniphila , and termite GI isolate Diplosphaera colotermitum .
A. muciniphila is able to colonize a broad range of hosts, despite differences in GI tract anatomy (simple, foregut, hindgut), diet (carnivore, omnivore, herbivore), host physiology and body temperature. This distribution might be an indication of co-evolution of this organism with its host. Therefore, we explored the presence and genomic divergence of Verrucomicrobia and Akkermansia spp. within different mammalian hosts.
A. muciniphila and other Verrucomicrobia within the GI tract of different mammals
Detailed analysis of Verrucomicrobia 16S rRNA gene sequences derived from SILVA database 138  (> 1100 bp, pintail > 75) revealed that A. muciniphila is not the sole representative species of the Verrucomicrobia phylum in the GI tract of mammals. The phylogenetic tree was constructed using both the neighbor joining method and RAxML. Both trees showed similar output and the Verrucomicrobia-derived sequences could be grouped into 12 clades of which 9 clades contained samples obtained from the mammalian intestine (Fig. 1 and Fig. S1).
The first clade contains 1352 A. muciniphila like 16S rRNA sequences, not solely mammalian. In this clade, eight mammalian orders were detected: Proboscidea (African elephant), primates (human, gorilla, lemur, chimpanzee, pygmy loris), Carnivora (cheetah), Sirenia (Dugong), Cingulata (armadillo), Rodentia (rat, mice and thirteen-lined ground squirrel), Artiodactyla (eland, pig, cow (rumen fluid)), Perissodactyla (horse). The non-mammalian sequence present in this clade was derived from a chicken. In addition, the 16S rRNA sequences of A. muciniphila strains isolated in this study were positioned in clade 1 (Fig. 1).
Five clades solely contained mammalian sequences, including clade 2, 3, 4, 5 and 6. Clade two contains 277 sequences, mainly derived from primates (human) and Rodentia (mice) and one sequence derived from Sirenia (dugong). Interestingly, sequences derived from the snub-nosed monkey formed a separate clade (3) within the phylogenetic tree, as well as 57 human-derived sequences in clade 4. Other clades solely consisting of mammalian derived sequences were clade 5 and 6. Clade 5 contained 201 sequences derived from 6 mammalian orders, including Proboscidea (African elephant), Artiodactyla (cow, okapi, buffalo, babirusa, warty pig, gazelle, takin, giraffe, przewalskii gazelle, springbok), Diprotodontia (kangaroo), Rodentia (capybara, Prevost’s squirrel), Perissodactyla (horse, wild ass, rhinoceros, zebra), Chiroptera (flying fox). Clade 6 contained 20 sequences from mammalian orders Lagomorpha (rabbit) and Proboscidea (elephant). Sequences derived from animals both living in captivity and in the wild were represented in all animal-containing clades excluding mice. To be able to compare the similarity (%) of type strain A. muciniphila MucT to the different clades, A. muciniphila MucT was compared to representative sequences of each clade. The similarities and amount of representatives per clade are shown in Table 1.
Other mammalian GI tract derived sequences that belonged to the Verrucomicrobia phylum were found within clade 8, 10 and 12. Clade 8 consisted of only five sequences in total, four derived from a python (A. glycaniphila) and one sequence from the mammalian order Sirenia (Dugong). Clade 10 contained two sequences that belong to the Prosthecobacter genus obtained from intestinal samples of UC patients (Fig. 1 and Table 1). Furthermore clade 12 contained sequences from mammalian orders Chiroptera (flying fox) and primates (hamadryas baboon and eastern black and white colobus and Sumatran orangutan), Diprotodontia (red kangaroo), Rodentia (capybara) and Lagomorpha (European rabbit) belonging to the class Opitutae. The remaining Verrucomicrobia clades (7, 9 and 11) solely contained 16S rRNA gene sequences derived from the non-mammal Animalia GI tract (Fig. 1 and Table 1).
Verrucomicrobia prevalence in fecal samples of different mammals
Fecal samples of 108 different animals belonging to 47 species of 16 mammalian orders were collected. The prevalence of A. muciniphila was determined by quantitative PCR (qPCR) on all samples. Amplicons were generated with A. muciniphila-specific primers in 50 out of 108 samples, with abundances up to 4% (Fig. 2). In addition, 16S rRNA gene sequencing of Verrucomicrobia resulted in Verrucomicrobia sequences derived from the Caribbean manatee, echidna, Western gorilla and otter. These sequences were added to the phylogenetic tree shown in Fig. 1, in which they were positioned in clade 1 (Caribbean manatee), clade 5 (Western gorilla and echidna) and clade 10 (otter).
New A. muciniphila isolates show low physiologic divergence
Ten new A. muciniphila isolates were obtained from fecal samples of the chimpanzee, siamang, mouse, pig, reindeer, horse, and elephant. The 16S rRNA gene sequence of these ten isolates was determined, and showed high similarity (> 99.9%) to the 16S rRNA gene sequence of A. muciniphila MucT (Fig. 3A). All new A. muciniphila strains had small and oval shaped cells of approximately 700 nm in length, as described previously for the type strain . All cells stained Gram-negative, grew in single cells, doublets, and aggregates in mucin medium. The cell growth (determined by OD600) and short chain fatty acid (SCFA) production (determined by high performance liquid chromatography (HPLC)) in a mucin-based medium was similar to the type strain A. muciniphila MucT (Fig. 4). In addition, the strains had similar growth rates, and produced similar amounts of acetate, propionate and 1,2-propanediol. Taking into consideration the morphologic, physiologic and 16S rRNA gene sequence similarity, all strains should belong to A. muciniphila species.
Low genomic divergence between 10 new A. muciniphila isolates
The genomic DNA of all newly obtained isolates were sequenced and assembled into draft genomes (Table 2) consisting of 25–215 contigs. Isolates had genome sizes in the range of A. muciniphila MucT (2.7 Mb), although the genomes of Chimpanzee 2_1, Chimpanzee 2_2, and Mouse 1 were slightly larger (2.9 Mb) (Table 2). This was also reflected in the total predicted gene count (Table 2). All isolates had comparable GC content ranging from 55.2 to 55.9 (Table 2). The average nucleotide identity (ANI) was > 99.7% for 7 isolates (from chimpanzee 1, siamang, mouse 1, pig, reindeer, horse, and elephant). The ANI was lower for the isolates from chimpanzee 2_1 (93.9%), chimpanzee 2_2 (97.4%), and mouse 2 (93.9%). The BLAST similarity and the number of SNPs of these three genomes were also in line with these results. This indicates that these three isolates are phylogenetically more distant. The ANI BLAST between all new isolates is shown in Table S1. The phylogeny based on the domain presence in the genomes was constructed and reflected the 16S rRNA gene phylogeny (Fig. 3B and Table S2). In addition the pan-genome has been determined as shown in Fig. S2.
Comparing all draft genomes of the new isolates to the complete genome of A. muciniphila MucT using BLAST Ring Image Generator (BRIG) showed the low genomic divergence (Fig. 5). Several potential phage remnants were identified that showed different GC content or GC skew and were not conserved among all isolates. Moreover, on 3 points in the draft genome there was a gap, likely because of the presence of one of the 3 rRNA operons that interfered with the draft genome sequence assembly.
Mucin degradation and utilization proteins of all isolates were individually analyzed for sequence divergence. In addition to the 61 mucin-degrading proteins predicted to be excreted , 18 other genes were predicted to be involved in mucin degradation and utilization. The seven isolates (from chimpanzee 1, siamang, mouse 1, pig, reindeer, horse, and elephant) that had the highest sequence similarity compared to strain MucT harbored the same 78 mucin degradation and utilization proteins (Additional Table S3), though the draft genome isolate obtained from a horse lacked one α-L-fucosidase (Amuc_1699). In comparison with the complete genome of the type strain MucT, we found that the draft genome of chimpanzee strain 2_2 lacked genes for two glycosyl hydrolases (Amuc_1637 and Amuc_1120) and the PfkB (Amuc_0075), which is annotated to be involved in fructose degradation. In comparison with the complete genome of the type strain MucT, we found that the draft genome of chimpanzee strain 2_1 lacked genes for two glycosyl hydrolases (Amuc_1120 and Amuc_0824) and one α-N-acetylglucosaminidase (Amuc_1220) as compared to the type strain. In comparison with the complete genome of the type strain MucT, we found that the draft genome of mouse strain 2 solely lacked the gene for one glycosyl hydrolase (Amuc_1637).
Akkermansia muciniphila is an abundant member of the healthy human intestine, and colonizes the mucus layer that lines the intestinal epithelial cells. Apart from human, Akkermansia 16S rRNA gene sequences can be detected in intestinal samples of many vertebrates. A. muciniphila is the only cultured representative of the Verrucomicrobia phylum obtained from the human GI tract and was isolated from a fecal sample of a healthy adult. We determined the prevalence of Verrucomicrobia using the SILVA database, the clone libraries of Verrucomicrobia specific amplicons, and performed qPCR on 108 fecal samples collected for this study. In addition, we obtained ten new A. muciniphila isolates from non-human mammals that were characterized by determining their draft genome.
Our results indicate that A. muciniphila is not the sole species belonging to the phylum Verrucomicrobia that colonizes the GI tract of mammals. Clade 5 and 6, that contained sequences derived from fecal samples of a wide variety of mammals excluding human and mice, were found to be closely related to A. muciniphila MucT. The 16S rRNA gene sequence similarities compared to MucT ranged between 85.15–90.69% and 86.25–89.22% for clade 5 and 6, respectively. In addition, more distantly related Verrucomicrobia sequences were detected in fecal samples. Verrucomicrobia sequences belonging to the Prosthecobacter family (85% similarity) were detected in GI samples of UC patients. Furthermore, Verrucomicrobia sequences belonging to the Opitutae class have been detected in the following mammals: flying fox, European rabbit, capybara, red kangaroo, Sumatran orangutan, hamadryas baboon and, eastern black and white colobus.
Clade 5 and 6 contain only non-human mammalian sequences. Our isolation effort did not result in any of these Akkermansia species. This could mean that mucin is not the main nutritional source for these organisms or that there are differences in host mucin and therefore in mucus degradation by Akkermansia. However, more research is needed to investigate the type of mucins produced and secreted by the gastrointestinal tract of different mammals to confirm this hypothesis. Recently, co-colonization with multiple A. muciniphila strains has been described to be possible inhealthy human individuals . It would be of interest to investigate whether A. muciniphila and a species from clade 5 or 6 could co-colonize in a mucin-dominated environment within the gut.
Our results confirm earlier reports that A. muciniphila is widely spread throughout mammals . Moreover, we provide a comparative genomic analysis of new A. muciniphila isolates. By combining 16S rRNA gene database sequences and quantitative PCR, we could detect Akkermansia in animals that belong to 15 out of 16 mammalian orders included in this study. The Akkermansia 16S rRNA gene sequences present within different animals were highly similar. This suggests that Akkermansia is highly conserved and only minor changes upon co-evolution with its different host species occurred.
The qPCR analysis did not confirm the presence of A. muciniphila in two of the samples from which new A. muciniphila isolates were obtained (chimpanzee 1 and reindeer). The Verrucomicrobia PCR did confirm the presence of Verrucomicrobia in these samples. Possible inhibition of qPCR amplification might be the reason that A. muciniphila was not detected in these samples.
Although for some orders only few animals were tested, there does not seem to be a correlation between the abundance of Akkermansia spp. and the host phylogeny. A possible explanation for this high degree of conservation could be horizontal spread of the microorganism among animals in captivity due to close contact with humans as most samples were taken in the Dutch zoos. However, we could confirm that Akkermansia spp. and other Verrucomicrobia can also be detected in fecal samples of animals that live in the wild, both from our and a previous study . The spread of Akkermansia spp. is not restricted to any geographical location, since highly similar 16S rRNA gene sequences were found in fecal samples taken from different parts of the world, more specific in rural Africa , rural Australia (this study), North-America [2, 32, 33] Asia  and Europe (this study) . Based on the observation that highly similar 16S rRNA gene sequences were found to be present within the different animals, we hypothesize that Akkermansia spp. are highly conserved within its different hosts. Potentially because only minor changes upon co-evolution with its different host species are needed to colonize the mucosal niche.
The presence of A. muciniphila in the mammalian GI-tract has also been assessed by reconstructing metagenome assembled genomes from datasets of human, mouse and pig gut microbiomes . In line with our findings, the presence of A. muciniphila in the mammalian gut was found to be globally distributed. In this study, we used culturing techniques to isolate new A. muciniphila strains from the mammalian GI-tract. However, metagenome assembled genomes of Akkermansia spp. that could not be obtained by culturing techniques in this study, may give more insight into the function of these strains and their ability to degrade mucus in the mammalian gut.
Mucin proteins are conserved among mammals, and even within the chordate phylum. Mucin glycoproteins are rich in proline, threonine and serine and are highly glycosylated . These properties enable A. muciniphila to use mucin as both nitrogen and carbon source . It is not known how the mucin proteins are glycosylated in the numerous vertebrates, but regardless of the potential glycosylation patterns the genomes of the A. muciniphila isolates encode many enzymes that can cleave a wide variety of glycan chains. Differences in the presence of mucus degradation genes in comparison to the type strain were only detected in three out of ten isolates. The genes lacking in these genomes may not directly have an effect on the mucus degrading capability, since other glycosyl hydrolases, α-L-fucosidases and N-acetylglucosaminidases are also present in the genome of A. muciniphila. Furthermore, it is important to note that the genomes of the new isolates were not closed. Therefore, it is a possibility that these genes are present in the genome but not detected in our analysis. Overall, this data suggests that A. muciniphila is a mucin-degrading specialist that has the potential to colonize different mammals regardless of their potential differences in mucin structure.
We did not observe indications for animal-species specific colonization when connecting the sequences of the Akkermansia clades with the hosts of origin. This contrasts what has been described for H. pylori, a well-studied mucosal pathogen that is mainly found in human and in very narrow range of other hosts . Testifying for the adaptation of H. pylori to the human host is the observation that its genome can be linked to human migration over our planet .
Our findings indicate that A. muciniphila is frequently colonizing the GI tract of mammals. In this study, we isolated 10 new A. muciniphila strains from feces of chimpanzee, siamang, mouse, pig, reindeer, horse and elephant. All new A. muciniphila isolates grew on mucin as sole carbon and nitrogen source suggesting that representatives of this species colonize the mucus layer of its host and therefore seem not be affected by the host diet or physiology. The low genomic divergence observed in the new strains may indicate that A. muciniphila favors mucosal colonization independent of the differences in hosts. In addition, the conserved mucus degradation capability points towards a similar beneficial role of the new strains in regulating host metabolic health.
Fecal samples from animals were obtained at three Dutch zoo’s: Burgers Zoo (Arnhem, The Netherlands), Dolfinarium (Hardewijk, The Netherlands), Natura Artis Magistra (Amsterdam, The Netherlands), Plankendael (Antwerp, Belgium), animal facilities of Wageningen UR (Wageningen, The Netherlands), animal facilities of Erasmus MC (Rotterdam, The Netherlands), animal facility of the Institute of Microbiology, (ETH Zurich, Switzerland), Animal facilities of Leiden University Medical Centre (Leiden, The Netherlands), from pets living at Dutch homes, but also from wild animals that live in either rural Africa or rural Australia. Mouse 2 has been deposited under the name A. muciniphila YL44, DSM26127 (https://www.dsmz.de/catalogues/dzif-sammlung-der-dsmz/maus-mikrobiomliste.html). Additional information about the samples can be found in Table S4.
Isolation and growth conditions
Approximately 0.2 g of fecal sample was taken and dissolved in anaerobic PBS (pH 7) containing 0.5 g/l of cysteine-HCL within 24 h of defecation. A fraction of all samples was used to prepare a glycerol (25%v/v) stock and stored at − 80 °C. The other fraction of all samples were 10-fold diluted in anaerobic mucin medium, composed of a bicarbonate-buffered basal medium  with a pH of 6.5–7.0, supplemented with 0.5% (vol/vol) purified and dialyzed hog gastric mucin (Type III, Sigma) as sole carbon and nitrogen source as described previously . All incubations were performed until growth was observed at 37 °C in 30 ml serum bottles, containing 10 ml mucin media, sealed with butyl rubber stoppers under anaerobic conditions provided by a gas phase of 1.5 atm N2/CO2 (80:20 vol/vol). Enrichment was achieved by repeated serial dilutions. After this primary enrichment, the strains were purified by repeated plating of single colonies onto anaerobic mucin medium agar (0.8% (w/v) agar (Bacto Agar, BD), only selecting Akkermansia-like colonies based on previously described morphology  and of which an Akkermansia-specific PCR  was found positive. The Short-chain fatty acids in cultures containing the purified strains were measured using a Thermo Electron spectrasystem HPLC equipped with an Agilent Metacarb 67H column. Purified strains were stored in mucin medium containing glycerol (25% v/v) at − 80 °C. DNA was extracted using the Masterpure™ Gram Positive DNA Purification Kit (Epicentre®).
A. muciniphila 16S rRNA gene abundance
Quantitative PCR amplification was performed as previously described  with minor modifications: samples were analyzed in a total volume of 10 μl consisting of 1 x iQ SYBR Green Supermix (BioRad), 200 nM forward primer AM1: CAGCACGTGAAGGTGGGGC  or 1369F: CGGTGAATACGTTCYCGG , 200 nM reverse primer AM2: CCTTGCGGTTGGCTTCAGAT  or 1492R: CGGCTACCTTGTTACGAC , 1 x VisiBlue Master Mix colorant (Tataabiocenter), and 0.2 ng/μl sample DNA, Nuclease-Free Water (Promega) was added to 10 μl. The primerset including AM1 and AM2 specifically amplifies A. muciniphila DNA and the primerset including 1369F and 1492R is a general 16S rRNA primerset to determine the total abundance. All reactions were performed in triplicates in a BioRad CFX-384 device (Veenendaal, The Netherlands). Standard curves of 16S rRNA from A. muciniphila cloned into pGMTeasy vector (Promega) were prepared, corresponding to a range from 108 to 100 cells. The quality of the standard curves were assessed using qPCR. The abundance of A. muciniphila 16S rRNA genes was determined by dividing the amount of Akkermansia 16S rRNA gene amplicon by that obtained from total 16S after correcting for the 16S rRNA gene copy of A. muciniphila (3 copies), and the average number of 16S rRNA genes (4.1 copies) in intestinal bacteria . The starting quantity (SQ) values used for the calculations are available in Table S5.
16S rRNA gene sequencing of Verrucomicrobia
DNA obtained from the fecal samples was amplified in a final volume of 25 μl consisting of 1 x Green GoTaq reaction buffer (Promega), 200 nM of each dNTPs (Promega), 200 nM forward primer VER_37: TGGCGGCGTGGWTAAGA , 200 nM of reverse primer VER_673: TGCTACACCGWGAATTC , 1 U GoTaq DNA polymerase (Promega), Nuclease-Free Water (Promega) was added to obtain a total volume of 25 μl. Samples were amplified with a Dinxperlo BV G Storm thermocycler (Somerton Biotechnology) with the following program: Denaturation at 95 °C for 5 min, followed by 35 cycles of denaturation at 95 °C for 30 s, annealing at 50 °C for 30 s, extension at 72 °C for 1 min, and a final extension step at 72 °C for 10 min. The amplicons were purified using a High pure PCR Cleanup micro kit following the manufacturer’s protocol (Roche, Woerden, the Netherlands). Ligation of these amplicons in pGEMTeasy vector system as described by the manufacturer (Promega) and subsequent transformation into E. coli XL1-blue competent cells (Agilent Technologies). Inserts were sequenced at GATC (Biotech, Konstanz, Germany) using the flanking binding sites for T7: TATTTAGGTGACACTATAG and SP6: TAATACGACTCACTATAGGG. Vector, primers and low quality ends of the sequences were trimmed using DNA-baser v.354. 16S rRNA gene sequences were aligned using the SINA online alignment services  and subsequently imported into ARB .
16S rRNA gene database mining and phylogenetic tree construction
All intestinal Verrucomicrobia sequences > 1100 bp, with pintails > 75 were downloaded from the SILVA database version 138. The isolation source and host organism, if lacking, were retrieved from the original publications if possible and added to the designated fields in the database. All analysis concerning 16S rRNA gene sequences used for data mining performed on this dataset. The selected outgroup for the phylogenetic analysis consisted of 13 sequences from three phyla: Lentisphaerae, Omnitrophica and Chlamydiae. All 16S rRNA based phylogenetic analysis were performed with a single trimmed alignment file. The phylogenetic tree was constructed in ARB (version 5.3-org-8209) using a randomized axelerated maximum likelihood (RAxML) method (version 7.0.3) and a 40% positional conservatory filter . Depending on the amount of sequences in each clade, up until 20 representatives of each clade were selected for sequence similarity comparisons of each clade to A. muciniphila MucT.
DNA isolation and genome sequencing
High molecular weight genomic DNA was extracted from overnight-grown cultures as previously described . DNA quality and concentrations were determined by spectrophotometric analysis using NanoDrop equipment (Thermo Scientific) and by electrophoresis on a 1% agarose gel. DNA was stored at − 20 °C until subsequent sequencing.
Genome sequencing was carried out at the Institute of Biotechnology, University of Helsinki (Finland). A MiSeq library was generated and sequenced on an Illumina MiSeq Personal Sequencer with 250 bp paired-end reads and an insert size of 500 bp. Reads were assembled using Ray (k-mer 101) .
Annotation was carried out with an in-house pipeline consisting of Prodigal v2.5 for prediction of protein coding DNA sequences , InterProScan 5RC7 for protein annotation , tRNAscan-SE v1.3.1 for prediction of tRNAs  and RNAmmer v1.2 for prediction of rRNAs . Additional protein function predictions were derived via BLAST identifications against the UniRef50  and Swissprot (UniProt, 2014) databases (download August 2013). Subsequently, the annotation was further enhanced by adding EC numbers via PRIAM version 2013-03-06 . Non-coding RNAs were identified using rfam_scan.pl v1.04, on release 11.0 of the RFAM database . CRISPRs were annotated using CRISPR Recognition Tool v1.1 . A further step of automatic curation was performed by weighing the annotation of the different associated domains, penalizing uninformative functions (e.g. “Domain of unknown function”), and prioritizing functions of interest (e.g. domains containing “virus”, “bacteriophage”, “integrase” for bacteriophage related elements; similar procedure for different other functions).
To determine the pan-genome, the genomes were annotated using Prokka . Subsequently, Roary was used to obtain the core-genome alignment . Based on this information, a maximum likelihood phylogenetic tree was constructed. The phylogenetic tree and the core-genome alignment were combined in Phandango to visualize the results .
16S rRNA gene sequence retrieval
For each organism the 16S rRNA reads were retrieved by filtering the FASTQ file through sortmeRNA using default settings while only using the 16S SILVA 118 database . The obtained reads where then assembled using IDBA_UD into a 16S rRNA gene contig and used for further analysis .
Availability of data and materials
The whole genome sequencing data of the A. muciniphila isolates obtained during the current study have been uploaded to the BioProject PRJEB21068 repository.
- GI tract:
Operational Taxonomic Units
Short Chain Fatty Acids
high performance liquid chromatography
Average Nucleotide Identity
BLAST Ring Image Generator
Randomized Axelerated maximum likelihood
Ley RE, Lozupone CA, Hamady M, Knight R, Gordon JI. Worlds within worlds: evolution of the vertebrate gut microbiota. Nat Rev Microbiol. 2008;6:776–88.
Ley RE, Hamady M, Lozupone C, Turnbaugh PJ, Ramey RR, Bircher JS, et al. Evolution of mammals and their gut microbes. Science (80- ). 2008;320:1647–51.
McFall-Ngai M, Hadfield MG, Bosch TCG, Carey HV, Domazet-Lošo T, Douglas AE, et al. Animals in a bacterial world, a new imperative for the life sciences. Proc Natl Acad Sci U S A. 2013;110:3229–36.
Muegge BD, Kuczynski J, Knights D, Clemente JC, González A, Fontana L, et al. Diet drives convergence in gut microbiome functions across mammalian phylogeny and within humans. Science (80- ). 2011;332:970–4.
Rajilić-Stojanović M, Shanahan F, Guarner F, de Vos WM. Phylogenetic analysis of dysbiosis in ulcerative colitis during remission. Inflamm Bowel Dis. 2013;19:481–8.
Van den Abbeele P, Van de Wiele T, Verstraete W, Possemiers S. The host selects mucosal and luminal associations of coevolved gut microorganisms: a novel concept. FEMS Microbiol Rev. 2011;35:681–704.
Derrien M, Vaughan EE, Plugge CM, de Vos WM. Akkermansia municiphila gen. nov., sp. nov., a human intestinal mucin-degrading bacterium. Int J Syst Evol Microbiol. 2004;54:1469–76.
Berry D, Stecher B, Schintlmeister A, Reichert J, Brugiroux S, Wild B, et al. Host-compound foraging by intestinal microbiota revealed by single-cell stable isotope probing. Proc Natl Acad Sci U S A. 2013;110:4720–5.
Ottman N, Davids M, Suarez-Diez M, Boeren S, Schaap PJ, VAP M, et al. Genome-scale model and omics analysis of metabolic capacities of Akkermansia muciniphila reveal a preferential mucin-degrading lifestyle. Appl Environ Microbiol. 2017;83:2020.
Plovier H, Everard A, Druart C, Depommier C, Van Hul M, Geurts L, et al. A purified membrane protein from Akkermansia muciniphila or the pasteurized bacterium improves metabolism in obese and diabetic mice. Nat Med. 2017;23:107–13.
Depommier C, Everard A, Druart C, Plovier H, Van Hul M, Vieira-Silva S, et al. Supplementation with Akkermansia muciniphila in overweight and obese human volunteers: a proof-of-concept exploratory study. Nat Med. 2019;25:1096–103.
Everard A, Belzer C, Geurts L, Ouwerkerk JP, Druart C, Bindels LB, et al. Cross-talk between Akkermansia muciniphila and intestinal epithelium controls diet-induced obesity. Proc Natl Acad Sci U S A. 2013;110:9066–71.
van der Lugt B, van Beek AA, Aalvink S, Meijer B, Sovran B, Vermeij WP, et al. Akkermansia muciniphila ameliorates the age-related decline in colonic mucus thickness and attenuates immune activation in accelerated aging Ercc1−/Δ7 mice. Immun Ageing. 2019;16:6.
Collado MC, Derrien M, Isolauri E, De Vos WM, Salminen S. Intestinal integrity and Akkermansia muciniphila, a mucin-degrading member of the intestinal microbiota present in infants, adults, and the elderly. Appl Environ Microbiol. 2007;73:7767–70.
Cani PD, de Vos WM. Next-generation beneficial microbes: the case of Akkermansia muciniphila. Front Microbiol. 2017;8:1765.
Sommer F, Ståhlman M, Ilkayeva O, Arnemo JM, Kindberg J, Josefsson J, et al. The gut microbiota modulates energy metabolism in the hibernating brown bear Ursus arctos. Cell Rep. 2016;14:1655–61.
Carey HV, Walters WA, Knight R. Seasonal restructuring of the ground squirrel gut microbiota over the annual hibernation cycle. Am J Physiol - Regul Integr Comp Physiol. 2013;304:R33.
Sonoyama K, Fujiwara R, Takemura N, Ogasawara T, Watanabe J, Ito H, et al. Response of gut microbiota to fasting and hibernation in Syrian hamsters. Appl Environ Microbiol. 2009;75:6451–6.
Costello EK, Gordon JI, Secor SM, Knight R. Postprandial remodeling of the gut microbiota in Burmese pythons. ISME J. 2010;4:1375–85.
Ouwerkerk JP, Aalvink S, Belzer C, de Vos WM. Akkermansia glycaniphila sp. nov., an anaerobic mucin-degrading bacterium isolated from reticulated python faeces. Int J Syst Evol Microbiol. 2016;66:4614–20.
Belzer C, De Vos WM. Microbes insidefrom diversity to function: the case of Akkermansia. ISME J. 2012;6:1449–58.
Green TJ, Smullen R, Barnes AC. Dietary soybean protein concentrate-induced intestinal disorder in marine farmed Atlantic salmon, Salmo salar is associated with alterations in gut microbiota. Vet Microbiol. 2013;166:286–92.
Isanapong J, Sealy Hambright W, Willis AG, Boonmee A, Callister SJ, Burnum KE, et al. Development of an ecophysiological model for Diplosphaera colotermitum TAV2, a termite hindgut Verrucomicrobium. ISME J. 2013;7:1803–13.
Sanders JG, Powell S, Kronauer DJC, Vasconcelos HL, Frederickson ME, Pierce NE. Stability and phylogenetic correlation in gut microbiota: lessons from ants and apes. Mol Ecol. 2014;23:1268–83.
Wüst PK, Horn MA, Drake HL. Clostridiaceae and Enterobacteriaceae as active fermenters in earthworm gut content. ISME J. 2011;5:92–106.
Vandekerckhove TTM, Coomans A, Cornelis K, Baert P, Gillis M. Use of the Verrucomicrobia-specific probe EUB338-III and fluorescent in situ hybridization for detection of “Candidatus xiphinematobacter” cells in nematode hosts. Appl Environ Microbiol. 2002;68:3121–5.
Wertz JT, Kim E, Breznak JA, Schmidt TM, Rodrigues JLM. Genomic and physiological characterization of the Verrucomicrobia isolate Diplosphaera colitermitum gen. nov., sp. nov., reveals microaerophily and nitrogen fixation genes. Appl Environ Microbiol. 2012;78:1544–55.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590.
Perelman P, Johnson WE, Roos C, Seuanez HN, Horvath JE, Moreira MA, Kessing B, Pontius J, Roelke M, Rumpler Y, Schneider MP, Silva A, O'Brien SJ, Pecon-Slattery J. A molecular phylogeny of living primates. PLoS Genet. 2011;7:e1001342.
van Passel MWJ, Kant R, Zoetendal EG, Plugge CM, Derrien M, Malfatti SA, et al. The genome of Akkermansia muciniphila, a dedicated intestinal mucin degrader, and its use in exploring intestinal metagenomes. PLoS One. 2011;6:e16876.
Guo X, Zhang J, Wu F, Zhang M, Yi M, Peng Y. Different subtype strains of Akkermansia muciniphila abundantly colonize in southern China. J Appl Microbiol. 2016;120:452–9.
Eckburg PB, Bik EM, Bernstein CN, Purdom E, Dethlefsen L, Sargent M, et al. Microbiology: diversity of the human intestinal microbial flora. Science (80- ). 2005;308:1635–8.
Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457:480–4.
Li M, Wang B, Zhang M, Rantalainen M, Wang S, Zhou H, et al. Symbiotic gut microbes modulate human metabolic phenotypes. Proc Natl Acad Sci U S A. 2008;105:2117–22.
Guo X, Li S, Zhang J, Wu F, Li X, Wu D, et al. Genome sequencing of 39 Akkermansia muciniphila isolates reveals its population structure, genomic and functional diverisity, and global distribution in mammalian gut microbiotas. BMC Genomics. 2017;18:800.
Moran AP, Gupta A, Joshi L. Sweet-talk: role of host glycosylation in bacterial pathogenesis of the gastrointestinal tract. Gut. 2011;60:1412–25.
Ottman NA. Host immunostimulation and substrate utilization of the gut symbiont Akkermansia muciniphila. PhD Thesis; 2015. p. 208.
Lee A, Fox J, Hazell S. Pathogenicity of helicobacter pylori: a perspective. Infect Immun. 1993;61:1601–10.
Falush D, Wirth T, Linz B, Pritchard JK, Stephens M, Kidd M, et al. Traces of human migrations in helicobacter pylori populations. Science (80- ). 2003;299:1582–5.
Miller RS, Hoskins LC. Mucin degradation in human colon ecosystems. Fecal population densities of mucin-degrading bacteria estimated by a “most probable number” method. Gastroenterology. 1981;81:759–65.
Suzuki MT, Taylor LT, DeLong EF. Quantitative analysis of small-subunit rRNA genes in mixed microbial populations via 5′-nuclease assays. Appl Environ Microbiol. 2000;66:4605–14.
Weisburg WG, Barns SM, Pelletier DA, Lane DJ. 16S ribosomal DNA amplification for phylogenetic study. J Bacteriol. 1991;173:697–703.
Case RJ, Boucher Y, Dahllöf I, Holmström C, Doolittle WF, Kjelleberg S. Use of 16S rRNA and rpoB genes as molecular markers for microbial ecology studies. Appl Environ Microbiol. 2007;73:278–88.
Ranjan K. Verrucomicrobia: a model phylum to study the effects of deforestation on microbial diversity in the Amazon forest. Arlington: The University of Texas; 2010.
Pruesse E, Peplies J, Glöckner FO. SINA: accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics. 2012;28:1823–9.
Ludwig W, Strunk O, Westram R, Richter L, Meier H, Yadhukumar A, et al. ARB: a software environment for sequence data. Nucleic Acids Res. 2004;32:1363–71.
Douillard FP, Ribbera A, Kant R, Pietilä TE, Järvinen HM, Messing M, et al. Comparative genomic and functional analysis of 100 lactobacillus rhamnosus strains and their comparison with strain GG. PLoS Genet. 2013;9:e1003683.
Boisvert S, Raymond F, Godzaridis É, Laviolette F, Corbeil J. Ray Meta: scalable de novo metagenome assembly and profiling. Genome Biol. 2012;13:1–13.
Hyatt D, Chen GL, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11:119.
Hunter S, Jones P, Mitchell A, Apweiler R, Attwood TK, Bateman A, et al. InterPro in 2011: new developments in the family and domain prediction database. Nucleic Acids Res. 2012;40:D306.
Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:955–64.
Lagesen K, Hallin P, Rødland EA, Stærfeldt HH, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35:3100–8.
Suzek BE, Huang H, McGarvey P, Mazumder R, Wu CH. UniRef: comprehensive and non-redundant UniProt reference clusters. Bioinformatics. 2007;23:1282–8.
Claudel-Renard C, Chevalet C, Faraut T, Kahn D. Enzyme-specific profiles for genome annotation: PRIAM. Nucleic Acids Res. 2003;31:6633–9.
Burge SW, Daub J, Eberhardt R, Tate J, Barquist L, Nawrocki EP, et al. Rfam 11.0: 10 years of RNA families. Nucleic Acids Res. 2013;41.D1:D226–D232.
Bland C, Ramsey TL, Sabree F, Lowe M, Brown K, Kyrpides NC, et al. CRISPR recognition tool (CRT): a tool for automatic detection of clustered regularly interspaced palindromic repeats. BMC Bioinformatics. 2007;8:209.
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30:2068–9.
Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MTG, et al. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics. 2015;31:3691–3.
Hadfield J, Croucher NJ, Goater RJ, Abudahab K, Aanensen DM, Harris SR. Phandango: an interactive viewer for bacterial population genomics. Bioinformatics. 2018;34:292–3.
Kopylova E, Noé L, Touzet H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics. 2012;28:3211–7.
Peng Y, Leung HCM, Yiu SM, Chin FYL. IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics. 2012;28:1420–8.
We thank Nicolette Snijders (NAM), Paulien Bunskoek (DF), Simone Kools (BZ), Plankendael for fecal sample collections, and Hans Heilig, Aylin Oymaci, Marloes Witte, Eline Stroobach for technical assistance. The authors gratefully appreciate Prof. Bert ‘t Hart (Rijswijk, the Netherlands) and Prof. Jon Lamans (Erasmus University) for availability of marmoset fecal samples, Prof. Els Urbanus for the gift of the initial fecal sample of which Mouse 1 was isolated, Tom van den Bogert for the gift of the cow rumen sample and Jing Zhang for the gift of the pig fecal samples. We are grateful to our colleagues at the laboratory of Microbiology for stimulating discussions.
J.P.O, C. B and J. R were funded by the ERC Advanced Grant 250172 - Microbes Inside from the European Research Council awarded to W.M.dV, whose work was further supported by an unrestricted Spinoza Award and the SIAM Gravity Grant 024.002.002 of the Netherlands Organization for Scientific Research (NWO), which also funded. SY.G.
Ethics approval and consent to participate
A subset of the human derived fecal samples were the baseline samples of a cohort study approved by the medical ethical board of Wageningen University (NL29070.081.09). In addition, consent was obtained from the zoos listed in Table S4 to collect dropped mammalian fecal samples.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Original detailed randomized axelerated maximum likelihood (RAxML) tree.
Table S1. ANI BLAST and aligned percentage between all new A. muciniphila isolates.
Table S2. Overview of the functional domains of the type strain and all new A. muciniphila isolates.
Figure S2. Pan-genome of the type strain and all new A. muciniphila isolates visualized using Phandango.
Table S3. Mucin degradation and utilization genes. Depicted are the sequence similarities of the genes from the newly obtained strains compared to the genes of strain MucT.
Table S4. Metadata of obtained fecal samples and isolates.
Table S5. Starting quantity (SQ) values of both Akkermansia and the general 16S rRNA qPCR runs.
About this article
Cite this article
Geerlings, S.Y., Ouwerkerk, J.P., Koehorst, J.J. et al. Genomic convergence between Akkermansia muciniphila in different mammalian hosts. BMC Microbiol 21, 298 (2021). https://doi.org/10.1186/s12866-021-02360-6
- Comparative genomics
- Phylogenetic analysis
- Gut bacteria
- Beneficial microbe