Pangenomic analysis identifies correlations between Akkermansia species and subspecies and human health outcomes

Abstract

Aim:Akkermansia are common members of the human gastrointestinal microbiota. The prevalence of these mucophilic bacteria, especially Akkermansia muciniphila (A. muciniphila), correlates with immunological and metabolic health. The genus Akkermansia in humans includes species with significantly larger genomes than A. muciniphila, leading us to postulate that this added genetic content may influence how they impact human metabolic and immunological health.

Methods: We conducted a pangenomic analysis of 234 Akkermansia complete or near-complete genomes. We also used high-resolution species and subspecies assignments to reanalyze publicly available metagenomic datasets to determine if there are relationships between Akkermansia species and A. muciniphila clades with various disease outcomes.

Results: Analysis of genome-wide average nucleotide identity, 16S rRNA gene identity, conservation of core Akkermansia genes, and analysis of the fatty acid composition of representative isolates support the partitioning of the genus Akkermansia into several species. In addition, A. muciniphila sensu stricto, the most prevalent Akkermansia species in humans, should be subdivided into two subspecies. For a pediatric cohort, we observed species-specific correlations between Akkermansia abundance with baseline obesity or after various interventions. For inflammatory bowel disease cohorts, we identified a decreased abundance of Akkermansia in patients with ulcerative colitis or Crohn’s disease, which was species and subspecies-dependent. In patients undergoing immune checkpoint inhibitor therapies for non-small cell lung carcinoma, we observed a significant association between one A. muciniphila subspecies and survival outcomes.

Conclusion: Our findings suggest that the prevalence of specific Akkermansia species and/or subspecies can be crucial in evaluating their association with human health, particularly in different disease contexts, and is an important consideration for their use as probiotics.

Keywords

Akkermansia / phylogeny / microbiome / IBD / obesity / immunotherapy

References

[1] Derrien M,Plugge CM.Akkermansia muciniphila gen. nov., sp. nov., a human intestinal mucin-degrading bacterium.Int J Syst Evol Microbiol2004;54:1469-76

[2] Derrien M,Ben-Amor K,de Vos WM.The Mucin degrader Akkermansia muciniphila is an abundant resident of the human intestinal tract.Appl Environ Microbiol2008;74:1646-8 PMCID:PMC2258631

[3] Everard A,Geurts L.Cross-talk between Akkermansia muciniphila and intestinal epithelium controls diet-induced obesity.Proc Natl Acad Sci U S A2013;110:9066-71 PMCID:PMC3670398

[4] Lukovac S,Pellis L.Differential modulation by Akkermansia muciniphila and Faecalibacterium prausnitzii of host peripheral lipid metabolism and histone acetylation in mouse gut organoids.mBio2014;5:e01438-14 PMCID:PMC4145684

[5] Plovier H,Druart C.A purified membrane protein from Akkermansia muciniphila or the pasteurized bacterium improves metabolism in obese and diabetic mice.Nat Med2017;23:107-13

[6] Routy B,Derosa L.Gut microbiome influences efficacy of PD-1-based immunotherapy against epithelial tumors.Science2018;359:91-7

[7] Schneeberger M,Gómez-Valadés AG.Akkermansia muciniphila inversely correlates with the onset of inflammation, altered adipose tissue metabolism and metabolic disorders during obesity in mice.Sci Rep2015;5:16643 PMCID:PMC4643218

[8] Earley H,Balfe Á,Winter DC.The abundance of Akkermansia muciniphila and its relationship with sulphated colonic mucins in health and ulcerative colitis.Sci Rep2019;9:15683 PMCID:PMC6821857

[9] Olson CA,Yano JM,Nusbaum DJ.The gut microbiota mediates the anti-seizure effects of the ketogenic diet.Cell2018;174:497 PMCID:PMC6062008

[10] Blacher E,Shapiro H.Potential roles of gut microbiome and metabolites in modulating ALS in mice.Nature2019;572:474-80

[11] Cox LM,Liu S.Gut microbiome in progressive multiple sclerosis.Ann Neurol2021;89:1195-211 PMCID:PMC8132291

[12] Salazar N,Fernández-Navarro T,Gonzalez S.Age-associated changes in gut microbiota and dietary components related with the immune system in adulthood and old age: a cross-sectional study.Nutrients2019;11:1765 PMCID:PMC6722604

[13] Perraudeau F,Bullard J.Improvements to postprandial glucose control in subjects with type 2 diabetes: a multicenter, double blind, randomized placebo-controlled trial of a novel probiotic formulation.BMJ Open Diabetes Res Care2020;8:e001319 PMCID:PMC7368581

[14] Cani PD,Derrien M,de Vos WM.Akkermansia muciniphila: paradigm for next-generation beneficial microorganisms.Nat Rev Gastroenterol Hepatol2022;19:625-37

[15] Depommier C,Druart C.Supplementation with Akkermansia muciniphila in overweight and obese human volunteers: a proof-of-concept exploratory study.Nat Med2019;25:1096-103 PMCID:PMC6699990

[16] Spreafico A,Araujo DV.First-in-class Microbial Ecosystem Therapeutic 4 (MET4) in combination with immune checkpoint inhibitors in patients with advanced solid tumors (MET4-IO trial).Ann Oncol2023;34:520-30

[17] Derrien M,Hooiveld G,Müller M.Modulation of mucosal immune response, tolerance, and proliferation in mice colonized by the mucin-degrader Akkermansia muciniphila.Front Microbiol2011;2:166 PMCID:PMC3153965

[18] Collado MC,Isolauri E,Salminen S.Intestinal integrity and Akkermansia muciniphila, a mucin-degrading member of the intestinal microbiota present in infants, adults, and the elderly.Appl Environ Microbiol2007;73:7767-70 PMCID:PMC2168041

[19] Guo X,Zhang J.Genome sequencing of 39 Akkermansia muciniphila isolates reveals its population structure, genomic and functional diverisity, and global distribution in mammalian gut microbiotas.BMC Genomics2017;18:800 PMCID:PMC5648452

[20] Becken B,Middleton DR.Genotypic and phenotypic diversity among human isolates of Akkermansia muciniphila.mBio2021;12:e00478-21 PMCID:PMC8262928

[21] Kelly C,Davey L.Spontaneous episodic inflammation in the intestines of mice lacking HNF4A is driven by microbiota and associated with early life microbiota alterations.mBio2023;14:e0150423 PMCID:PMC10470520

[22] Kirmiz N,Cross KL.Comparative genomics guides elucidation of vitamin B₁₂ biosynthesis in novel human-associated akkermansia strains.Appl Environ Microbiol2020;86:e02117-19 PMCID:PMC6974653

[23] Padilla L,Luna E.Mechanism of 2’-fucosyllactose degradation by human-associated Akkermansia.J Bacteriol2024;206:e0033423 PMCID:PMC10886448

[24] Xing J,Sun Y.Comparative genomic and functional analysis of Akkermansia muciniphila and closely related species.Genes Genomics2019;41:1253-64 PMCID:PMC6828834

[25] Karcher N,Punčochář M.Genomic diversity and ecology of human-associated Akkermansia species in the gut microbiome revealed by extensive metagenomic assembly.Genome Biol2021;22:209 PMCID:PMC8278651

[26] Kim JS,Lee JH,Lee JS.The evolution and competitive strategies of Akkermansia muciniphila in gut.Gut Microbes2022;14:2025017 PMCID:PMC8920140

[27] Li W,Jing Y.Comparative genomics revealed wide intra-species genetic heterogeneity and lineage-specific genes of Akkermansia muciniphila.Microbiol Spectr2022;10:e0243921 PMCID:PMC9241678

[28] Bukhari SAR,Ahmad I.Comparative genomics and pan-genome driven prediction of a reduced genome of Akkermansia muciniphila.Microorganisms2022;10:1350 PMCID:PMC9315967

[29] Kobayashi Y,Inoue S.Akkermansia biwaensis sp. nov., an anaerobic mucin-degrading bacterium isolated from human faeces.Int J Syst Evol Microbiol2023;73:005697

[30] Ndongo S,Raoult D.Reclassification of eight Akkermansia muciniphila strains and description of Akkermansia massiliensis sp. nov. and Candidatus Akkermansia timonensis, isolated from human feces.Sci Rep2022;12:21747 PMCID:PMC9758162

[31] Zhai R,Zhang L,Zhao L.Strain-specific anti-inflammatory properties of two Akkermansia muciniphila strains on chronic colitis in mice.Front Cell Infect Microbiol2019;9:239 PMCID:PMC6624636

[32] Quast C,Yilmaz P.The SILVA ribosomal RNA gene database project: improved data processing and web-based tools.Nucleic Acids Res2013;41:D590-6 PMCID:PMC3531112

[33] McIver LJ,Franzosa EA.bioBakery: a meta’omic analysis environment.Bioinformatics2018;34:1235-7 PMCID:PMC6030947

[34] Oxford Nanopore Technologies. FAQs. Available from: https://nanoporetech.com/support. [Last accessed on 6 Jun 2024]

[35] Kolmogorov M,Lin Y.Assembly of long, error-prone reads using repeat graphs.Nat Biotechnol2019;37:540-6

[36] BCL convert. Illumina, Inc. Available from: https://support.illumina.com/sequencing/sequencing_software/bcl-convert.html. [Last accessed on 6 Jun 2024]

[37] Walker BJ,Shea T.Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement.PLoS One2014;9:e112963 PMCID:PMC4237348

[38] Gillespie JJ,Cammer SA.PATRIC: the comprehensive bacterial bioinformatics resource with a focus on human pathogenic species.Infect Immun2011;79:4286-98 PMCID:PMC3257917

[39] McMurdie PJ.phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data.PLoS One2013;8:e61217 PMCID:PMC3632530

[40] McCann JR,Roche K.The pediatric obesity microbiome and metabolism study (POMMS): methods, baseline data, and early insights.Obesity2021;29:569-78 PMCID:PMC7927749

[41] Derosa L,Thomas AM.Intestinal Akkermansia muciniphila predicts clinical response to PD-1 blockade in patients with advanced non-small-cell lung cancer.Nat Med2022;28:315-24 PMCID:PMC9330544

[42] Mills RH,Vázquez-Baeza Y.Multi-omics analyses of the ulcerative colitis gut microbiome link Bacteroides vulgatus proteases with disease severity.Nat Microbiol2022;7:262-76 PMCID:PMC8852248

[43] Lloyd-Price J, Arze C, Ananthakrishnan AN, et al; IBDMDB Investigators. Multi-omics of the gut microbial ecosystem in inflammatory bowel diseases.Nature2019;569:655-62 PMCID:PMC6650278

[44] Franzosa EA,Avila-Pacheco J.Gut microbiome structure and metabolic activity in inflammatory bowel disease.Nat Microbiol2019;4:293-305 PMCID:PMC6342642

[45] Eren AM,Shaiber A.Community-led, integrated, reproducible multi-omics with anvi’o.Nat Microbiol2021;6:3-6 PMCID:PMC8116326

[46] Eddy SR.Accelerated profile HMM searches.PLoS Comput Biol2011;7:e1002195 PMCID:PMC3197634

[47] Delmont TO.Linking pangenomes and metagenomes: the Prochlorococcus metapangenome.PeerJ2018;6:e4320 PMCID:PMC5804319

[48] Pritchard L,Humphris S,Toth IK.Genomics and taxonomy in diagnostics for food security: soft-rotting enterobacterial plant pathogens.Anal Methods2016;8:12-24

[49] Price MN,Arkin AP.FastTree: computing large minimum evolution trees with profiles instead of a distance matrix.Mol Biol Evol2009;26:1641-50 PMCID:PMC2693737

[50] Veseli I,Schechter MS.Microbes with higher metabolic independence are enriched in human gut microbiomes under stress.eLife2023;12:RP89862

[51] Sievers F,Dineen D.Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega.Mol Syst Biol2011;7:539 PMCID:PMC3261699

[52] RCoreTeam. The R project for statistical computing. Available from: https://www.r-project.org/. [Last accessed on 6 Jun 2024]

[53] Beghini F,Blanco-Míguez A.Integrating taxonomic, functional, and strain-level profiling of diverse microbial communities with bioBakery 3.Elife2021;10:e65088 PMCID:PMC8096432

[54] Truong DT,Pasolli E,Segata N.Microbial strain-level population structure and genetic diversity from metagenomes.Genome Res2017;27:626-38 PMCID:PMC5378180

[55] Bisanz JE,Noecker C.A genomic toolkit for the mechanistic dissection of intractable human gut bacteria.Cell Host Microbe2020;27:1001-13.e9 PMCID:PMC7292766

[56] Emiola A,Oh J.Metagenomic growth rate inferences of strains in situ.Sci Adv2020;6:eaaz2299 PMCID:PMC7176420

[57] Inkscape. Available from: https://inkscape.org/. [Last accessed on 6 Jun 2024]

[58] Letunic I.Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation.Nucleic Acids Res2021;49:W293-6 PMCID:PMC8265157

[59] GraphPad Prism. Available from: https://www.graphpad.com/. [Last accessed on 6 Jun 2024]

[60] Luna E,Kirmiz N.Utilization efficiency of human milk oligosaccharides by human-associated Akkermansia is strain dependent.Appl Environ Microbiol2022;88:e0148721 PMCID:PMC8752153

[61] Goris J,Klappenbach JA,Vandamme P.DNA-DNA hybridization values and their relationship to whole-genome sequence similarities.Int J Syst Evol Microbiol2007;57:81-91

[62] Richter M.Shifting the genomic gold standard for the prokaryotic species definition.Proc Natl Acad Sci U S A2009;106:19126-31 PMCID:PMC2776425

[63] Pearce ME,Lauer AC,Maiden MCJ.An evaluation of the species and subspecies of the genus Salmonella with whole genome sequence data: Proposal of type strains and epithets for novel S. enterica subspecies VII, VIII, IX, X and XI.Genomics2021;113:3152-62 PMCID:PMC8426187

[64] Kim M,Park SC.Towards a taxonomic coherence between average nucleotide identity and 16S rRNA gene sequence similarity for species demarcation of prokaryotes.Int J Syst Evol Microbiol2014;64:346-51

[65] Stackebrandt E.Taxonomic note: a place for DNA-DNA reassociation and 16S rRNA sequence analysis in the present species definition in bacteriology.Int J Syst Evol Microbiol1994;44:846-9

[66] Lv QB,Zhang Y.A thousand metagenome-assembled genomes of Akkermansia reveal phylogroups and geographical and functional variations in the human gut.Front Cell Infect Microbiol2022;12:957439 PMCID:PMC9378777

[67] Tindall BJ,Busse HJ,Kämpfer P.Notes on the characterization of prokaryote strains for taxonomic purposes.Int J Syst Evol Microbiol2010;60:249-66

[68] Diogo A,Nobre MF.Usefulness of fatty acid composition for differentiation of Legionella species.J Clin Microbiol1999;37:2248-54 PMCID:PMC85129

[69] Hall AB,Sauk J.A novel Ruminococcus gnavus clade enriched in inflammatory bowel disease patients.Genome Med2017;9:103 PMCID:PMC5704459

[70] Rice P,Bleasby A.EMBOSS: the european molecular biology open software suite.Trends Genet2000;16:276-7

[71] Choi Y,Seo J.Effects of live and pasteurized forms of akkermansia from the human gut on obesity and metabolic dysregulation.Microorganisms2021;9:2039 PMCID:PMC8538271

[72] Zhang T,Wu X.Alterations of Akkermansia muciniphila in the inflammatory bowel disease patients with washed microbiota transplantation.Appl Microbiol Biotechnol2020;104:10203-15

[73] Lopez-Siles M,Aldeguer X.Alterations in the abundance and co-occurrence of Akkermansia muciniphila and Faecalibacterium prausnitzii in the colonic mucosa of inflammatory bowel disease subjects.Front Cell Infect Microbiol2018;8:281 PMCID:PMC6137959

[74] Shaw KA,Hofmekler T.Dysbiosis, inflammation, and response to treatment: a longitudinal study of pediatric subjects with newly diagnosed inflammatory bowel disease.Genome Med2016;8:75 PMCID:PMC4944441

[75] Jin Y,Zhou J.Genome-based classification of Burkholderia cepacia complex provides new insight into its taxonomic status.Biol Direct2020;15:6 PMCID:PMC7057466

[76] Caudill MT.The use and limitations of the 16S rRNA sequence for species classification of anaplasma samples.Microorganisms2022;10:605 PMCID:PMC8949108

[77] González D,Vidal-Veuthey B.Insights into early evolutionary adaptations of the Akkermansia genus to the vertebrate gut.Front Microbiol2023;14:1238580 PMCID:PMC10540074