Abstract
Background:
Numerous studies have demonstrated that Akkermansia is closely associated with human health. These bacteria colonize the mucus layer of the gastrointestinal tract and utilize mucin as their sole source of carbon and nitrogen. Akkermansia spp. exhibit potential as probiotics under specific conditions. However, the gene accumulation curve derived from pan-genome analysis suggests that the genome of Akkermansia strains remains open. Consequently, current genome mining efforts are insufficient to fully capture the intraspecific and interspecific characteristics of Akkermansia, necessitating continuous exploration of the genomic and phenotypic diversity of new isolates.
Methods:
Based on this finding, we sequenced, assembled, and functionally annotated the whole genomes of four new human isolates from our laboratory: AKK-HX001, AKK-HX002, AKK-HX003, and AKK-HX004.
Results:
Phylogenetic analysis revealed that all four isolates belonged to the AmII phylogroup, whereas the type strain DSM 22959 is classified within the AmI phylogroup. Moreover, 2,184 shared homologous genes were identified among the four isolates. Functional annotation using the COG, KEGG, and CAZy databases indicated that the functional genes of the four isolates were primarily associated with metabolism. Two antibiotic resistance genes were identified in AKK-HX001 and AKK-HX002, while three resistance genes were detected in AKK-HX003 and AKK-HX004. Additionally, each of the four isolates possessed two virulence genes and three pathogenicity genes, none of which were associated with pathogenicity. The prediction of mobile genetic elements indicated unequal distributions of GIs among the isolates, and a complete CRISPR system was identified in all isolates except AKK-HX003. Two annotated regions of secondary metabolite biosynthesis genes, both belonging to Terpene, were detected using the antiSMASH online tool.
Conclusion:
These findings indicate that the four Akkermansia isolates, which belong to a phylogroup distinct from the model strain DSM 22959, exhibit lower genetic risk and may serve as potential probiotic resources for future research.
1 Introduction
The intestinal tract possesses the richest microbiota in the body: the intestinal microbiota, which is closely related to the development of human diseases (; ). concluded that the composition of the gut microbiota is unique to each individual and varies at all stages of life. With the rise of large-scale sequencing technologies, the composition and function of the gut microbiota are gradually being decoded, and more and more links between diseases and the gut microbiota are being characterized. Hence, it is possible to design personalized drugs for specific diseases that are targeted to an individual’s specific gut microbiota (). A subset of microorganisms derived from human commensals have potential to be the next-generation probiotics (NGPs), which attract increasing attention. With the deepening of multi-faceted research, it has been found that the adaptability of these strains to the intestinal environment allows them to produce bioactive compounds and have a positive impact on gut health, immune function, and metabolism, which may not be possible with traditional probiotics (). Their ability to target specific health problems may include the production of bioactive compounds, the modulation of immune responses, and competitive interactions with pathogens, all of which may give NGPs an advantage in targeted therapies for chronic diseases (; ).
Akkermansia muciniphila (A. muciniphila), a representative of the Verrucomicrobia phylum, was first isolated from human feces in 2004 by , is a focal point of such studies due to its ability to colonize the intestinal mucus layer and influence host physiology (; ). Numerous studies have linked A. muciniphila to beneficial health outcomes, including its role in reducing metabolic disorders, regulating the immune response, and maintaining gut barrier integrity (; ). For example, its abundance correlates positively with lower obesity risk, reduced inflammation in inflammatory bowel disease (), and improved glucose metabolism in diabetic patients (). Additionally, novel associations with familial Mediterranean fever and tumor modulation have been reported (; ; ). However, recent studies present a contrasting perspective. Research by and suggests that A. muciniphila may exacerbate immune nephropathy and induce inflammation in the absence of dietary fiber. These findings emphasize the importance of context and strain specificity in evaluating its probiotic potential. Such conflicting results underscore the urgent need for strain-level functional studies to clarify safety profiles and therapeutic applications.
The ability to isolate and characterize Akkermansia strains has significantly expanded in recent years (; ). Improved culturing methods have revealed substantial inter- and intra-species genomic variability, which influences host adaptation and functional traits (; ; ; ). Pan-genome analyses indicate that Akkermansia has an open genome structure (), with new genetic features discovered in every additional strain isolated. This highlights the need to continue isolating strains and analyzing their phenotypic diversity to fully understand their roles in the microbiota. Developing an Akkermansia gene library could aid in identifying functional traits, enabling targeted applications in precision medicine. Furthermore, strain-specific adaptations to hosts and geographic regions suggest the potential for designing probiotics tailored to specific populations.
In this study, four Akkermansia isolates from the intestinal tract of Chinese infants and young children, previously isolated in the laboratory, were subjected to whole genome sequencing to analyze and characterize their potential functional features, safety, and horizontal gene transfer ability at the gene level. This study aims to provide a theoretical basis for subsequent research on the functional and beneficial effects of Akkermansia, contribute new insights into the regional and population diversity of Akkermansia strains, and use genomic data to guide the selection of practical strains.
2 Materials and methods
2.1 Microorganisms and culture conditions
Akkermansia muciniphila (A. muciniphila, DSM22959) was generously provided by Prof. Liu Li, Nanjing Agricultural University. The 4 Akkermansia strains were isolated from the feces of infants provided by Sichuan Provincial Maternity and Child Health Care Hospital using a conventional isolation process. Briefly, the fecal glycerol tube was thawed in a 37°C water bath, diluted tenfold with saline, and cultured under anaerobic conditions for 5 days at 37°C in mucin liquid medium supplemented with vancomycin (6 μg/mL) and kanamycin (12 μg/mL) for bacterial enrichment. During the growth process, 100 μL of a uniformly turbid bacterial solution was taken from the anaerobic tube, diluted tenfold to 10−5, and 100 μL of the diluted solution was spread onto BHI agar plates containing vancomycin and kanamycin for isolation and culture. The plates were incubated anaerobically at 37°C for 1–2 weeks, with PCR detection conducted in parallel. Single colonies suspected of containing the target bacteria were isolated and purified through plate streaking and BHI liquid medium enrichment. The DNA of the isolates was then extracted, and 16S rRNA sequencing was performed after PCR amplification using universal primers 1492R and 27F. The sequencing results were uploaded to GenBank for comparison, confirming that the 16S rRNA sequences of the 4 isolates showed 100% similarity to those of Akkermansia, verifying their classification within the genus. The strains were designated as AKK-HX001, AKK-HX002, AKK-HX003, and AKK-HX004. The bacterial solutions were preserved at −80°C with the addition of 30% glycerol in a 1:1 ratio.
The preserved strains were recovered by streaking on BHI agar plates (1 week), activated, and multiplied in BHI liquid medium (48 h), with all processes carried out anaerobically at 37°C.
2.2 Extraction of DNA, library construction, and whole-genome sequencing
Bacterial cells were harvested by centrifugation of the logarithmic growth phase of the bacterial solution at 12,000 r/min for 5 min. DNA extraction was performed using the TIANamp Bacteria DNA Kit, according to the manufacturer’s instructions. DNA quantity and purity were measured by the NanoDrop 2000 (Thermo Fisher Scientific).
Library preparation was performed using the NEBNext®Ultra™ II DNA Library Prep Kit for Illumina (NEB, United States, Catalog #: E7370L), with 0.2 μg of total DNA per sample. Genomic DNA samples were fragmented by sonication, and 350 bp fragments were selected for end-polishing, A-tailing, and ligation to full-length adapters, followed by PCR amplification. Agencourt AMPure XP (Beverly, United States) purified the PCR products. Sequencing was performed on the Illumina Miseq PE300 (Novogene, China, Peking) platform after the evaluation of library construction.
2.3 Genome assembly and prediction
Quality control and filtering of whole genome sequencing results were performed using Fastp (). De novo assembly was performed using SPAdes (), followed by genome annotation using Prokka (), both of which are available sources. Predictions and mapped genomic circular maps were created using the CGView online database (). Homology analysis was performed using OrthoFinder () to obtain the gene family classification of the 4 Akkermansia isolates and their respective numbers.
2.4 Phylogenomic analysis and comparison of average nucleotide identity
Fifteen human Akkermansia isolates of known phylogroups as well as the model strain DSM22959 were selected and their genome sequences were downloaded from NCBI. The 16S rRNA sequences were extracted from the 20 whole-genome sequences using Barrnap, imported into MEGA11 () for comparison and clipping, and selected the “Maximum Likelihood Tree” approach, with Bootstrap set to 1,000 for phylogenetic tree construction. Phylogenomic analysis based on single-copy core genes was performed using GToTree (), and the evolutionary tree was landscaped using the ITOL online website (). FastANI () was used to calculate the average nucleotide identity (ANI) between the 16 sequences and the 4 isolate genomes at two intervals.
2.5 Functional annotation
In order to recognize the function of the predicted genes, the web annotation tool eggNOG-mapper () was used to blast with the Clusters of Orthologous Groups of proteins (COG) database. In parallel, an analysis was conducted using KofamKOALA () to compare the function and composition of the 4 isolates’ biological systems by Kyoto Encyclopedia of Genes and Genomes (KEGG) database annotation. Carbohydrate-Active enZYmes Database (CAZy) database () annotation was performed by importing the amino acid sequences of predicted genes into the DbCAN3 website (). Secondary metabolite gene clusters were identified by antiSMASH ().
2.6 Genomic safety assessment based on CARD, VFDB, and PHI databases
To assess the safety of the 4 Akkermansia isolates at the genetic level, their antibiotic resistance, virulence factors, and pathogenicity genes were analyzed separately by blasting against The Comprehensive Antibiotic Resistance Database (CARD) (), Virulence Factor Database (VFDB) (), and Pathogen-Host Interactions database (PHI) (), where CARD is based on “the perfect and strict hits” select criteria, VFDB and PHI employed the results with the cut-off value at >70% identity and optimal comparison results for further analysis.
2.7 Analysis of mobile genetic elements
The 4 isolates were predicted to mobile genetic elements (MGEs) after genome rearrangement using MAVUE software (Darling et al., 2004). In this study, the MGEs we analyzed included genomic islands (GI) and the Clustered Regularly Interspaced Short Palindromic Repeats-Cas (CRISPR-Cas) system, which were predicted by IslandView 4 () and CRISPRCas Finder (), respectively.
2.8 Statistical analysis
If not specifically indicated, the visualization of the data was achieved through Origin (Version 2021, OriginLab Corporation, Northampton, MA, United States).
3 Results
3.1 General genomic characteristics of the isolates
The genome size, GC content, protein-coding sequences (CDS), and repeat_regions of the 4 Akkermansia isolates generally exceeded those of the model strain, and rRNA counts were lower than the type strain. In contrast, the amounts of tRNA were 52, 54, 51, and 53, respectively, which were close to the DSM22959 (54), and all of them encompassed 1 tmRNA, distinctively. The details were displayed in Table 1 below and the genome circular diagram was shown in Figure 1A. Homology analysis revealed 2,184 homologous genes in the genomes of the 4 Akkermansia isolates, with 75, 22, 105, and 26 unique genes, respectively (Figure 1B). COG functional annotation of genes specific to the 4 isolates revealed that only 10.53% (24/228) of the genes were assigned to the corresponding functional categories and were mainly related to information storage and processing (8/24), with the majority of the genes encoding hypothetical proteins for which the specific functions they exercise are not known (11/24).
Table 1
| Class definition | AKK-HX001 | AKK-HX002 | AKK-HX003 | AKK-HX004 | DSM22959 |
|---|---|---|---|---|---|
| Genome size (Mb) | 3.21 | 3.28 | 3.28 | 3.51 | 2.67 |
| GC percent (%) | 58.18 | 57.82 | 57.78 | 57.72 | 55.80 |
| CDS | 2,474 | 2,653 | 2,658 | 2,587 | 2,184 |
| Number of repeat_region | 4 | 5 | 4 | 4 | 1 |
| Number of rRNA | 3 | 3 | 3 | 3 | 9 |
| Number of tRNA | 52 | 54 | 51 | 53 | 54 |
| Number of tmRNA | 1 | 1 | 1 | 1 | 1 |
General characterization of the genomes of 4 Akkermansia isolates.
Figure 1
3.2 Phylogenetic analysis
3.2.1 16S rRNA phylogenetic tree analysis and phylogenomic analysis
Previously, classified the genus A. muciniphila into four phylogenetic groups based on the research of . To clarify the phylogroup classification of the 4 Akkermansia isolates from our laboratory, 15 human Akkermansia isolates with known phylogroup classifications and DSM22959 were selected. Their whole-genome sequences were downloaded from the NCBI, and their 16S rRNAs were extracted to construct the phylogenetic tree (Figure 2A). The 4 Akkermansia isolates were subsumed into the AmII phylogroup, while DSM22959 belonged to the AmI phylogroup by contrast. The results showed that AKK-HX001 belonged to a different evolutionary branch from AKK-HX002, AKK-HX003, and AKK-HX004, suggesting that AKK-HX001 may be more distantly related to the other three strains. Phylogenomic analysis based on single-copy core genes also yielded consistent results (Supplementary Figure S1). This was confirmed by a phylogenetic tree built based on ANI (Figure 2B). The ANI values between the 4 isolates and DSM22959 were 87.49, 87.48, 87.32, and 87.41%, respectively, while the ANI values among the 4 isolates were as high as 98.53–99.93%.
Figure 2
3.3 Genome functional classification
3.3.1 COG annotation
COG includes 4 major classifications, information storage and processing, cellular processes and signaling, metabolism, and category unknown, which are also subdivided into 26 functional categories. Four Akkermansia isolates were predicted to have 1,772, 1,788, 1,782, and 1,788 functional protein genes annotated into 20 functional categories by COG analysis, respectively. Among them, the higher number of functional genes was category M (Cell wall/membrane/envelope biogenesis), category E (Amino acid transport and metabolism), category J (Translation, ribosomal structure and biogenesis), category G (Carbohydrate transport and metabolism), category L (Replication, recombination and repair), and category C (Energy production and conversion) (Figure 3). Notably, AKK-HX001 lacked the category Z (Cytoskeleton) functional gene, while the other three strains annotated two genes of category Z.
Figure 3
3.3.2 KEGG annotation
KEGG pathway annotation of the predicted protein-coding genes of the 4 Akkermansia isolates resulted in 1,129, 1,132, 1,135, and 1,132 genes being included in the annotation of 203, 204, 204, and 204 metabolic pathways, respectively. These genes belonged to six major categories in the KEGG pathway first-level classification: metabolism, genetic information processing, environmental information processing, cellular processes, organismal systems, and human diseases, and were mainly annotated to the first two major categories. Each of the 37, 37, 38, and 37 metabolic pathways was annotated in the second-level classification (Figure 4). AKK-HX003 had more genes annotated to “metabolism,” especially “carbohydrate metabolism,” which suggested that AKK-HX003 may have a stronger ability to utilize and degrade carbohydrates. AKK-HX003 had one gene annotated to the Chromosome, while the other three strains did not have. The coding genes that were annotated to metabolism were the most numerous and were mainly related to carbohydrate metabolism, amino acid metabolism, metabolism of cofactors and vitamins. Next up was genetic information processing, and the coding genes annotated to translation, replication and repair processes were mainly involved. In terms of environmental information processing, the coding genes were mostly related to membrane transport and signal transduction. As for human diseases, the coding genes were significantly associated with antibiotic resistance.
Figure 4
3.3.3 CAZy annotation
The 4 Akkermansia isolates, with 146, 149, 147, and 149 coding genes, were annotated to six carbohydrate-active enzyme families, respectively. These included 2 AAs (Auxiliary Activities), 4 CBMs (Carbohydrate-Binding Modules) (AKK-HX003 contained only three and lacked CBM9), which has only been found in xylanases with the ability to bind cellulose so far (), 6 CEs (Carbohydrate Esterases), 29 GHs (Glycoside Hydrolases), 14 GTs (GlycosylTransferases), and 1 PLs (Polysaccharide Lyases) (Figure 5). The largest number of genes was assigned to the GHs and GTs families, and the results of the annotations for the AAs and PLs were consistent, as cataloged in Table 2.
Figure 5
Table 2
| CAZy-categories | AKK-HX001 | AKK-HX002 | AKK-HX003 | AKK-HX004 |
|---|---|---|---|---|
| Number of AAs | 2 | 2 | 2 | 2 |
| Number of CBMs | 10 | 10 | 9 | 10 |
| Number of CEs | 9 | 9 | 9 | 9 |
| Number of GHs | 66 | 67 | 65 | 67 |
| Number of GTs | 67 | 68 | 69 | 68 |
| Number of PLs | 1 | 1 | 1 | 1 |
CAZy classification statistics of the genomes of 4 Akkermansia isolates.
3.4 Genomic safety assessment
Safety is a basic, essential characteristic of probiotics. The safety of A. muciniphila, a potential next-generation probiotic, is currently controversial. To explore the safety of the 4 Akkermansia isolates at the genetic level, they were analyzed for antibiotic resistance and mechanisms, virulence genes, and pathogenicity genes for prediction, respectively.
3.4.1 Antibiotic resistance and mechanism
By comparing the predicted amino acid sequences with the CARD database, the results of AKK-HX001 and AKK-HX002 were consistent, and both of them were annotated with two antibiotic resistance genes (qacG, adeF), and the antibiotic resistance mechanism was antibiotic efflux pump (resistance-nodulation-cell division, RND; and small multidrug resistance, SMR). On the other hand, AKK-HX003 and AKK-HX004 annotated one more antibiotic resistance gene: ErmB, encoding 23S rRNA methyltransferase, whose antibiotic resistance mechanism is related to the change of antibiotic targets () (Table 3). qacG is an SMP subclass protein of the SMR family encoded by a plasmid or class I integrator that confers resistance to benzalkonium chloride in bacteria (). adeF is involved in the membrane fusion protein encoding the RND efflux pump system, AdeFGH, which equips the bacteria with resistance to fluoroquinolone and tetracycline antibiotics (). Expression of ErmB is inducible by erythromycin. The leader peptide causes attenuation of the mRNA and stabilizes the structure, preventing further translation. When erythromycin is present, it binds the leader peptide causing a change in conformation allowing for the expression of ErmB (). It renders bacteria resistant to erythromycin, roxithromycin, and lincomycin antibiotics ().
Table 3
| Class definition | Resistance genes | Resistance mechanism | AMR gene family | Antibiotic |
|---|---|---|---|---|
| AKK-HX001 | qacG | Antibiotic efflux | SMR | Benzalkonium chloride |
| adeF | Antibiotic efflux | RND | Fluoroquinolone; tetracycline | |
| AKK-HX002 | qacG | Antibiotic efflux | SMR | Benzalkonium chloride |
| adeF | Antibiotic efflux | RND | Fluoroquinolone; tetracycline | |
| AKK-HX003 | ErmB | Antibiotic target alteration | Erm 23S ribosomal RNA methyltransferase | Streptogramin; streptogramin B; streptogramin A; lincosamide; macrolide |
| qacG | Antibiotic efflux | SMR | Benzalkonium chloride | |
| adeF | Antibiotic efflux | RND | Fluoroquinolone; tetracycline | |
| AKK-HX004 | ErmB | Antibiotic target alteration | Erm 23S ribosomal RNA methyltransferase | Streptogramin; streptogramin B; streptogramin A; lincosamide; macrolide |
| qacG | Antibiotic efflux | SMR | Benzalkonium chloride | |
| adeF | Antibiotic efflux | RND | Fluoroquinolone; tetracycline |
Prediction results of antibiotic resistance genes in 4 Akkermansia isolates.
3.4.2 Virulent factors
The coding genes of the 4 isolates were annotated in the VFDB database, and the results were consistent, with two virulence factor-coding genes, katA and tufA, annotated in each strain, accounting for only 0.075–0.080% of the CDS.
3.4.3 Pathogenicity
The PHI annotation results basically corresponded to the virulence gene annotation results, and all 4 isolates were annotated with three pathogenic genes, katA, tufA, and fumB, which accounted for only 0.113–0.121% of the CDS.
3.5 Mobile genetic elements analysis
Prokaryotes, such as bacteria, evolve mainly through horizontal gene transfer (HGT) to obtain new genes from other individuals and often produce abundant MGEs in the process (). MGEs can reflect the stability of the strain and its ability to adapt to the environment (), and the 4 isolates were analyzed and found to have MGEs, including genomic islands and the CRISPR-Cas system.
3.5.1 Genomic islands
The 4 isolates had 39, 41, 129, and 91 coding genes involved in constituting 2, 3, 7, and 5 GIs, respectively. COG functional annotation of the above coding genes showed that the coding genes in the GI of the 4 isolates were mainly involved in transcription, replication, recombination and repair processes (Figure 6). To understand the security of the above coding genes at the gene level, they were annotated with CARD, VFDB, and PHI databases, and no other antibiotic resistance genes, virulence factor-coding genes, or pathogenic genes were found, except for one antibiotic resistance gene (ErmB) predicted on one GI of AKK-HX004.
Figure 6
3.5.2 CRISPR-Cas system
The results of CRISPR-Cas system prediction showed that AKK-HX001, AKK-HX002, and AKK-HX004 each had a complete IC-based CRISPR system, including cas2, cas1, cas4, cas3, cas5c, cas8c, and cas7c proteins. AKK-HX003 was predicted to have only 4 CRISPRs; no cas proteins were found, and no complete CRISPR-Cas system (Table 4).
Table 4
| Class definition | AKK-HX001 | AKK-HX002 | AKK-HX003 | AKK-HX004 |
|---|---|---|---|---|
| Number of CRISPR | 4 | 5 | 4 | 4 |
| Number of Cas cluster | 1 | 1 | – | 1 |
| Cas type | CAS-TypeIC | CAS-TypeIC | – | CAS-TypeIC |
| Number of cas gene | 7 | 7 | – | 7 |
| Cas gene | cas2_TypeI-II-III cas1_TypeIC cas4_TypeI-II cas3_TypeI cas5c_TypeIC cas8c_TypeIC cas7c_TypeIC | cas2_TypeI-II-III cas1_TypeIC cas4_TypeI-II cas3_Type cas5c_TypeIC cas8c_TypeIC cas7c_TypeIC | – | cas2_TypeI-II-III cas1_TypeIC cas4_TypeI-II cas3_Type cas5c_TypeIC cas8c_TypeIC cas7c_TypeIC |
Prediction results of the CRISPR-Cas system of 4 Akkermansia isolates.
3.6 Secondary metabolite gene cluster prediction
Probiotic strains have to be competitively colonized and survive in the gut, and bacteria with the ability to produce bacteriocins are more advantageous ().
By comparison with the antiSMASH database, two clusters of secondary metabolite genes were predicted for all 4 isolates, both representing terpenoids (Table 5).
Table 5
| Class definition | AKK-HX001 | AKK-HX002 | AKK-HX003 | AKK-HX004 |
|---|---|---|---|---|
| Region | Region 1.1 | Region 3.1 | Region 7.1 | Region 4.1 |
| Region 7.1 | Region 6.1 | Region 13.1 | Region 6.1 | |
| Region position | 569,022–589,942 | 309,741–330,634 | 229,105–250,025 | 401,871–422,764 |
| 1–20,664 | 9,058–29,978 | 90,731–111,624 | 31,644–52,564 | |
| Type | Terpene | Terpene | Terpene | Terpene |
Characteristics of gene clusters in 4 Akkermansia isolates.
4 Discussion
NGPs are suitable for conventional food or dietary supplements and pharmaceutical applications. It offers a variety of advantages over traditional probiotics, including opportunities for personalized probiotic therapies, involvement in synthetic biology and gene editing, participation in combination therapies, targeted delivery methods, and application in therapeutic settings (; ). A. muciniphila, a member of the next generation of probiotics, has received much attention in vivo. However, due to its harsh requirements for the growth environment, the acquisition of its isolates has become a formidable challenge, which has hindered the study of the genomes of Akkermansia isolates. Currently, studies on the genome of Akkermansia have mainly come from the analysis of gastrointestinal macro-genome assembly data. Therefore, in order to fill part of the gap of the genomic data of Akkermansia isolates as much as possible, the present study was carried out to analyze the genomes of 4 Akkermansia isolates obtained from the preliminary isolation in our laboratory, and revealed their genomic features, potential genomic functions, and also evaluated their safety and security at the genetic level.
Whole genome sequencing results showed that the genome size and CDS counts of the 4 Akkermansia isolates were much higher than the type strain DSM22959, suggesting that there might be a substantial difference between the genomes of our 4 isolates and the type strain. The GC content usually reflects the stability of the genome, and a high GC content may lead to the organisms consuming more energy during the replication process (). Whereas the results indicated that the GC content of these 4 isolates was higher than that of the model strain, suggesting that there may be some differences between them and the model strain in terms of metabolism. Homology analysis of the 4 isolates displayed that each possessed a different number of specific genes, given that most of these specific genes were hypothetical proteins that were not annotated to their corresponding functional classifications, we could not analyze them further. However, the presence of specific genes corroborates the diversity of the Akkermansia genome to a certain extent.
Phylogenetic analysis based on 16S rRNA and phylogenomic analysis based on single-copy core genes disclosed that these 4 isolates did not belong to the same phylogroup as the type strain. Phylogenetic analysis based on 16S rRNA revealed that AKK-HX001 belonged to a different branch from the other three strains. In contrast, the phylogenomic analysis based on single-copy core genes more finely classified AKK-HX002 and AKK-HX004 into the same branch, while AKK-HX003 was classified in the same branch as AKK-HX001. As for the reason why AKK-HX001 and AKK-HX003 were in a different evolutionary branch from the rest of the strains, we hypothesized that it is host-associated. The 4 isolates came from different human hosts, whose habits, such as living and eating, may have influenced their genotypes and phenotypes, which in turn induced their evolution. Functional annotation of the genomes of the four isolates also revealed that AKK-HX001 and AKK-HX003 were different from the other two strains to a different extent, and the presence of a higher number of unique genes in AKK-HX001 and AKK-HX003 may also contribute to their evolutionary differences. This further confirmed that phylogenomic analysis may be more rigorous than based on 16S rRNA phylogenetic analysis when full genomes are available. This classification result was also confirmed in the genome average nucleotide concordance. The ANI values between the 4 isolates and DSM22959 were 87.49, 87.48, 87.32, and 87.41%, respectively. Based on the threshold ANI > 96% between prokaryotic species (), it can be assumed that the 4 isolates are not of the same species of Akkermansia as DSM22959. Rather, the ANI values between the 4 isolates were as high as 98.53–99.93%, and it can be assumed that they belong to the same species within the genus Akkermansia. This further confirmed that differences exist with respect to the genomes of the 4 isolates and the model strain.
COG functional annotation is a relatively simple and reliable way to identify potential orthologs and paralogs based on complete genomic protein sequences (). The KEGG database is a large-scale integrated repository of systematic, genomic, chemical, and health information (). KEGG is based on the annotation of the functional ortholog list, which connects genomic information with higher-order functional information, stored in the GENES and KEGG pathway databases, providing rich pathway information and helping us understand the biological functions of genes systematically (). COG analysis of the predicted coding genes of the 4 isolates indicated that most of the genes were assigned to the functional categories of information storage and processing, cellular signaling, and metabolism, which are frequently associated with bacterial life activities and energy metabolism, suggesting that they possessed a strong ability to metabolize amino acids and carbohydrates. Immediately following the KEGG database annotation of the coding genes, similar results were obtained, with the vast majority of genes assigned to metabolism (especially carbohydrate metabolism) and genetic information processing pathways, further corroborating the importance of carbohydrates as an energy source in the biological functions of Akkermansia strains. The carbohydrate-active enzyme studies of the 4 isolates will hopefully deepen our insights into their growth and metabolism mechanisms. CAZy is a database resource on enzymes capable of synthesizing or breaking down complex carbohydrates and sugar complexes, categorizing carbohydrate-active enzymes into different protein families based on amino acid sequence similarities in the structural domains of the proteins (). The database supplies classification and related information on enzymes for synthesizing, metabolizing, and transporting carbon compounds (). Annotation of the CAZy database identified the 4 isolates as being rich in genes coding for carbohydrate-active enzymes (CAZymes), mainly GHs and GTs. GHs hydrolyze glycosidic bonds between two or more carbohydrates, or between carbohydrates and non-carbohydrate fractions, and the genes of the GHs annotated this time belonged to the GH2, GH20, and GH29 families, which are the enzymes with β-galactosidase activity, β-N-acetylamino glucosidase and β-N-acetyl galactosidase activity, and Α-L-fucosidase activity, respectively. These enzymes are sufficiently required for A. muciniphila to fully utilize 2′-fucosyl lactose (). GTs are enzymes that catalyze the formation of glycosides from glycosidic bonds, and the 4 isolates annotated with GTs were mainly derived from the GT2 and GT4 families, all of which possessed multiple glycosyltransferase activities (). Carbohydrate-active enzymes digest most of our complex pool of dietary polysaccharides, ultimately producing short-chain fatty acids, which are involved in numerous physiological and biochemical processes in our body (). Several clinical studies have suggested that short-chain fatty acids may bridge the involvement of A. muciniphila in glucose metabolism (), a role that our CAZy annotation results seem to echo, providing a gene-level explanation for.
A review described the role of A. muciniphila in reducing intestinal inflammation, modulating immune response, and enhancing intestinal barrier function, supporting its potential as a probiotic (). A recent study also validated the safety of A. muciniphila for use as a probiotic (). Since Akkermansia is a little-studied genus and its mechanism of action is not well understood, it is difficult to analyze it directly at the species level, which should be studied at the strain level (). A significant need exists to explore the antibiotic resistance genes, virulence genes, and pathogenicity genes of these 4 isolates of Akkermansia at the genetic level to assess their safety. Antibiotic resistance mechanisms include intrinsic resistance, acquired resistance, and adaptive resistance, among which intrinsic and acquired resistance are closely related to the strain’s antibiotic resistance genes. Therefore, the comparison between intrinsic antibiotic resistance genes and adaptive resistance genes should be emphasized (). Intrinsic antibiotic resistance genes may mediate nonspecific resistance mechanisms, including nonspecific efflux pumps and inactivating enzymes, whereas acquired antibiotic resistance genes can arise from mutations and/or horizontal transfer of genes that mediate specific resistance mechanisms, including enzymes modifying the antibiotic or the antibiotic target (). In the present study, a total of three antibiotic resistance genes were annotated in the 4 isolates, among which the adeF and ermB genes conferred multiple antibiotic resistances to the bacteria, and were hypothesized to mediate nonspecific resistance mechanisms, potentially belonging to intrinsic resistance genes. This is consistent with the report of Filardi et al. that the adeF gene is widely present in the genome of Akkermansia strains (). The qacG gene confers bacterial benzalkonium chloride (quaternary ammonium disinfectant) resistance and may mediate specific resistance mechanisms; however, available reports on Akkermansia have not identified such a gene for this resistance. It is worth mentioning that in our unpublished study on antibiotic susceptibility, the 4 isolates showed resistance to glycopeptides, aminoglycosides, and other antibiotics, but the corresponding resistance genes were not found in this genomic analysis. The reasons for this may involve several factors: firstly, A. muciniphila is a Gram-negative bacterium, which makes it intrinsically resistant to glycopeptides due to its outer membrane (). In addition, aminoglycosides are known to be less active against Gram-negative anaerobes because energy is required for this process, which is mostly unavailable due to their metabolism (). Therefore, it is not surprising that there is genotypic and phenotypic inconsistency in antibiotic resistance in Akkermansia strains. The predictions of virulence and pathogenicity genes were more aligned, with all 4 isolates predicting the katA gene, which is associated with stress survival, and the tufA gene, which is associated with adhesion ability. Among these genes, the katA gene encodes catalase, which helps bacteria resist reactive oxygen species (), and belongs to the category of unaffected pathogenicity, which is associated with host opportunistic pathogenic infections (). Su et al. showed that Pseudomonas aeruginosa (PA) exhibited higher katA expression under anaerobic conditions compared to aerobic conditions and acted as a protector during anaerobic respiration in PA (), suggesting that the ability of A. muciniphila to tolerate trace amounts of oxygen may be attributed to katA. tufA encodes the elongation factor Tu (EF-Tu), which is involved in the binding of aminoacyl-tRNA to the ribosomal A site (). Annotation results suggest that its mutant phenotype belongs to unaffected pathogenicity and/or increased virulence, associated with skin infection, food poisoning, and respiratory disease (). A fumB pathogenic gene, encoding fumarate hydratase class I, the absence of which causes Salmonella typhimurium to exhibit lower virulence, was also annotated (). The presence of several of these virulence factors may result in a pathogenic risk for the strain, but several previous animal experiments (; ) and clinical trials (; ) have demonstrated that A. muciniphila is safe. showed that it is already present within the first year of life, accounting for 3–5% of the gut microbiota, and gradually increasing with age until adulthood (). Since its safety has been controversial, we favor the presence of these virulence factors to aid in its adaptation to the environment and facilitate colonization.
In addition, to assess the intergenomic transfer of antibiotic resistance genes, virulence genes, and pathogenic genes, as well as resistance to invasion by exogenous genetic material, 4 isolates were predicted to have MGEs (GIs, CRISPR-Cas system). GIs are the most significant class of HGT, carrying genes affecting the host pathobiology and providing a selective advantage for host adaptation (). Four isolates were each predicted to have variable numbers of GIs, and no virulence factors were identified in GI-encoding genes except for one EmrB resistance gene. They were mostly associated with genetic information storage and processing. All isolates except AKK-HX003 were matched with the detection of an IC-type CRISPR-Cas system, which is most typical of Akkermansia strains (). The absence of CRISPR sites may enhance the stability of the strain’s genome and improve its adaptability to the environment (). Therefore, it can be assumed that the 4 isolates have a low risk of potential transfer and some resistance to invasion by exogenous genetic material.
Terpenes, a general term for a range of terpenoids, many of which have important biological activities, were identified in a previous study as terpene gene clusters in the genus Pseudovibrio with the potential to produce novel active compounds (). Two terpene gene clusters were identified in each of the 4 isolates in this study, but were not compared to a particular class of terpene compounds known in databases, probably because little research has been done on Akkermansia secondary metabolite gene clusters and due to the high percentage of genes encoding hypothetical proteins in the genome of this genus, suggesting a high potential for mining the functions of its genome.
5 Conclusion
This study, involving whole genome sequencing and genome analysis of 4 isolates of Akkermansia, resulted in a striking similarity between their genomes, with only about 87% ANI similarity to the model strain DSM22959. Additionally, the 4 isolates did not belong to the same phylogroup as DSM22959. The 4 isolates also showed homogeneity in terms of genome function, being mainly involved in the processing of genetic information, metabolic processes, and carrying abundant carbohydrate-active enzymes. Although a small number of antibiotic resistance genes, virulence genes, and pathogenicity genes were annotated, it was almost proven that these genes were not associated with pathogenicity. No significant mobile genetic elements were detected, indicating a low potential risk of gene transfer and a certain level of safety at the genetic level. In contrast to the relatively well-studied Lactobacillus and Bifidobacterium species, Akkermansia exhibited a certain number of annotated metabolic genes, albeit slightly fewer than those of the former two, suggesting that its probiotic effects may not be entirely metabolism-dependent. Additionally, a significant proportion of the genes in the four Akkermansia isolates were annotated as “Function Unknown,” highlighting substantial gaps in our understanding of the Akkermansia genome. This observation indirectly suggests that Akkermansia possesses considerable potential for further exploration, with these uncharacterized genes potentially holding the key to elucidating its probiotic role. Therefore, it is imperative to isolate and analyze additional strains to enhance genomic data and facilitate the functional characterization of these “Function Unknown” genes.
Statements
Data availability statement
The whole genome sequence datasets analyzed for this study can be found in the NCBI GenBank with the accession numbers JBGKAZ000000000, JBGKBA000000000, JBGKBB000000000, and JBGKBC000000000.
Author contributions
WL: Conceptualization, Data curation, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. BZ: Conceptualization, Data curation, Software, Visualization, Writing – review & editing. JL: Project administration, Supervision, Writing – review & editing. CL: Investigation, Resources, Writing – review & editing. JC: Investigation, Resources, Writing – review & editing. SD: Investigation, Resources, Writing – review & editing. XZ: Investigation, Resources, Writing – review & editing. LL: Project administration, Resources, Supervision, Writing – review & editing. GW: Conceptualization, Project administration, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
Appreciate Sichuan Provincial Maternal and Child Health Care Hospital for providing the source of strain isolation samples.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1500886/full#supplementary-material
References
1
AbouelelaM. E.HelmyY. A. (2024). Next-generation probiotics as novel therapeutics for improving human health: current trends and future perspectives. Microorganisms12:430. doi: 10.3390/microorganisms12030430
2
AlcockB. P.HuynhW.ChalilR.SmithK. W.RaphenyaA. R.WlodarskiM. A.et al. (2023). CARD 2023: expanded curation, support for machine learning, and resistome prediction at the comprehensive antibiotic resistance database. Nucleic Acids Res.51, D690–D699. doi: 10.1093/nar/gkac920
3
Al-EmranH. M.MoonJ. F.MiahM. L.MeghlaN. S.ReubenR. C.UddinM. J.et al. (2022). Genomic analysis and in vivo efficacy of Pediococcus acidilactici as a potential probiotic to prevent hyperglycemia, hypercholesterolemia and gastrointestinal infections. Sci. Rep.12:20429. doi: 10.1038/s41598-022-24791-5
4
AramakiT.Blanc-MathieuR.EndoH.OhkuboK.KanehisaM.GotoS.et al. (2020). KofamKOALA: KEGG Ortholog assignment based on profile HMM and adaptive score threshold. Bioinformatics36, 2251–2252. doi: 10.1093/bioinformatics/btz859
5
BeckenB.DaveyL.MiddletonD. R.MuellerK. D.SharmaA.HolmesZ. C.et al. (2021). Genotypic and phenotypic diversity among human isolates of Akkermansia muciniphila. MBio12, e00478–e00421. doi: 10.1128/mBio.00478-21
6
BertelliC.LairdM. R.WilliamsK. P.LauB. Y.HoadG.WinsorG. L.et al. (2017). IslandViewer 4: expanded prediction of genomic islands for larger-scale datasets. Nucleic Acids Res.45, W30–W35. doi: 10.1093/nar/gkx343
7
BlinK.ShawS.AugustijnH. E.ReitzZ. L.BiermannF.AlanjaryM.et al. (2023). antiSMASH 7.0: new and improved predictions for detection, regulation, chemical structures and visualisation. Nucleic Acids Res.51, W46–W50. doi: 10.1093/nar/gkad344
8
CantalapiedraC. P.Hernández-PlazaA.LetunicI.BorkP.Huerta-CepasJ. (2021). eggNOG-mapper v2: functional annotation, Orthology assignments, and domain prediction at the metagenomic scale. Mol. Biol. Evol.38, 5825–5829. doi: 10.1093/molbev/msab293
9
CantarelB. L.CoutinhoP. M.RancurelC.BernardT.LombardV.HenrissatB. (2009). The carbohydrate-active EnZymes database (CAZy): an expert resource for Glycogenomics. Nucleic Acids Res.37, D233–D238. doi: 10.1093/nar/gkn663
10
ChenS.ZhouY.ChenY.GuJ. (2018). fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics34, i884–i890. doi: 10.1093/bioinformatics/bty560
11
ColladoM. C.DerrienM.IsolauriE.De VosW. M.SalminenS. (2007). Intestinal integrity and Akkermansia muciniphila, a mucin-degrading member of the intestinal microbiota present in infants, adults, and the elderly. Appl. Environ. Microbiol.73, 7767–7770. doi: 10.1128/aem.01477-07
12
CouvinD.BernheimA.Toffano-NiocheC.TouchonM.MichalikJ.NéronB.et al. (2018). CRISPRCasFinder, an update of CRISRFinder, includes a portable version, enhanced performance and integrates search for Cas proteins. Nucleic Acids Res.46, W246–W251. doi: 10.1093/nar/gky425
13
CoyneS.RosenfeldN.LambertT.CourvalinP.PérichonB. (2010). Overexpression of resistance-nodulation-cell division pump AdeFGH confers multidrug resistance in Acinetobacter baumannii. Antimicrob. Agents Chemother.54, 4389–4393. doi: 10.1128/aac.00155-10
14
DepommierC.EverardA.DruartC.PlovierH.Van HulM.Vieira-SilvaS.et al. (2019). Supplementation with Akkermansia muciniphila in overweight and obese human volunteers: a proof-of-concept exploratory study. Nat. Med.25, 1096–1103. doi: 10.1038/s41591-019-0495-2
15
DepommierC.Van HulM.EverardA.DelzenneN. M.De VosW. M.CaniP. D. (2020). Pasteurized Akkermansia muciniphila increases whole-body energy expenditure and fecal energy excretion in diet-induced obese mice. Gut Microbes11, 1231–1245. doi: 10.1080/19490976.2020.1737307
16
DerrienM.BelzerC.De VosW. M. (2017). Akkermansia muciniphila and its role in regulating host functions. Microb. Pathog.106, 171–181. doi: 10.1016/j.micpath.2016.02.005
17
DerrienM.VaughanE. E.PluggeC. M.De VosW. M. (2004). Akkermansia muciniphila gen. Nov., sp. nov., a human intestinal mucin-degrading bacterium. Int. J. Syst. Evol. Microbiol.54, 1469–1476. doi: 10.1099/ijs.0.02873-0
18
DrulaE.GarronM. L.DoganS.LombardV.HenrissatB.TerraponN. (2022). The carbohydrate-active enzyme database: functions and literature. Nucleic Acids Res.50, D571–D577. doi: 10.1093/nar/gkab1045
19
DarlingA. C.MauB.BlattnerF. R.PernaN. T. (2004). Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res, 14:1394–1403. doi: 10.1101/gr.2289704
20
EmmsD. M.KellyS. (2019). OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol.20:238. doi: 10.1186/s13059-019-1832-y
21
EverardA.BelzerC.GeurtsL.OuwerkerkJ. P.DruartC.BindelsL. B.et al. (2013). Cross-talk between Akkermansia muciniphila and intestinal epithelium controls diet-induced obesity. Proc. Natl. Acad. Sci.110, 9066–9071. doi: 10.1073/pnas.1219451110
22
FilardiR.GargariG.MoraD.ArioliS. (2022). Characterization of antibiotic-resistance traits in Akkermansia muciniphila strains of human origin. Sci. Rep.12:19426. doi: 10.1038/s41598-022-23980-6
23
GalperinM. Y.KristensenD. M.MakarovaK. S.WolfY. I.KooninE. V. (2019). Microbial genome analysis: the COG approach. Brief. Bioinform.20, 1063–1070. doi: 10.1093/bib/bbx117
24
GaubaA.RahmanK. M. (2023). Evaluation of antibiotic resistance mechanisms in gram-negative Bacteria. Antibiotics (Basel)12:1590. doi: 10.3390/antibiotics12111590
25
GeerlingsS. Y.KostopoulosI.De VosW. M.BelzerC. (2018). Akkermansia muciniphila in the human gastrointestinal tract: when, where, and how?Microorganisms6:75. doi: 10.3390/microorganisms6030075
26
GeerlingsS. Y.OuwerkerkJ. P.KoehorstJ. J.RitariJ.AalvinkS.StecherB.et al. (2021). Genomic convergence between Akkermansia muciniphila in different mammalian hosts. BMC Microbiol.21:298. doi: 10.1186/s12866-021-02360-6
27
GleesonP. J.BenechN.ChemounyJ.MetallinouE.BerthelotL.DaS. J.et al. (2024). The gut microbiota posttranslationally modifies IgA1 in autoimmune glomerulonephritis. Sci. Transl. Med.16:eadl6149. doi: 10.1126/scitranslmed.adl6149
28
GrantJ. R.EnnsE.MarinierE.MandalA.HermanE. K.ChenC.-Y.et al. (2023). Proksee: in-depth characterization and visualization of bacterial genomes. Nucleic Acids Res.51, W484–W492. doi: 10.1093/nar/gkad326
29
GuoX.LiS.ZhangJ.WuF.LiX.WuD.et al. (2017). Genome sequencing of 39 Akkermansia muciniphila isolates reveals its population structure, genomic and functional diverisity, and global distribution in mammalian gut microbiotas. BMC Genomics18:800. doi: 10.1186/s12864-017-4195-3
30
HajdukP. J.DingesJ.SchkeryantzJ. M.JanowickD.KaminskiM.TufanoM.et al. (1999). Novel inhibitors of Erm methyltransferases from NMR and parallel synthesis. J. Med. Chem.42, 3852–3859. doi: 10.1021/jm990293a
31
HeirE.SundheimG.HolckA. L. (1999). The qacG gene on plasmid pST94 confers resistance to quaternary ammonium compounds in staphylococci isolated from the food industry. J. Appl. Microbiol.86, 378–388. doi: 10.1046/j.1365-2672.1999.00672.x
32
HeoS.KimJ. H.KwakM. S.SungM. H.JeongD. W. (2021). Functional annotation genome unravels potential probiotic Bacillus velezensis strain KMU01 from traditional Korean fermented kimchi. Food Secur.10:563. doi: 10.3390/foods10030563
33
HillC.GuarnerF.ReidG.GibsonG. R.MerensteinD. J.PotB.et al. (2014). The international scientific Association for Probiotics and Prebiotics consensus statement on the scope and appropriate use of the term probiotic. Nat. Rev. Gastroenterol. Hepatol.11, 506–514. doi: 10.1038/nrgastro.2014.66
34
HouK.WuZ.-X.ChenX.-Y.WangJ.-Q.ZhangD.XiaoC.et al. (2022). Microbiota in health and diseases. Signal Transduct. Target. Ther.7:135. doi: 10.1038/s41392-022-00974-4
35
IoannouA.BerkhoutM. D.GeerlingsS. Y.BelzerC. (2024). Akkermansia muciniphila: biology, microbial ecology, host interactions and therapeutic potential. Nat. Rev. Microbiol. doi: 10.1038/s41579-024-01106-1
36
JainC.RodriguezR. L.PhillippyA. M.KonstantinidisK. T.AluruS. (2018). High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun.9:5114. doi: 10.1038/s41467-018-07641-9
37
KandasamyS.YooJ.YunJ.LeeK. H.KangH. B.KimJ. E.et al. (2022). Probiogenomic in-silico analysis and safety assessment of Lactiplantibacillus plantarum DJF10 strain isolated from Korean raw Milk. Int. J. Mol. Sci.23:14494. doi: 10.3390/ijms232214494
38
KanehisaM.FurumichiM.SatoY.KawashimaM.Ishiguro-WatanabeM. (2023). KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res.51, D587–D592. doi: 10.1093/nar/gkac963
39
KanehisaM.GotoS. (2000). KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res.28, 27–30. doi: 10.1093/nar/28.1.27
40
KarcherN.NigroE.PunčochářM.Blanco-MíguezA.CicianiM.ManghiP.et al. (2021). Genomic diversity and ecology of human-associated Akkermansia species in the gut microbiome revealed by extensive metagenomic assembly. Genome Biol.22:209. doi: 10.1186/s13059-021-02427-7
41
KirmizN.GalindoK.CrossK. L.LunaE.RhoadesN.PodarM.et al. (2020). Comparative genomics guides elucidation of vitamin B12 biosynthesis in novel human-associated Akkermansia strains. Appl. Environ. Microbiol.86, e02117–e02119. doi: 10.1128/AEM.02117-19
42
LalowskiP.ZielińskaD. (2024). The Most promising next-generation probiotic candidates—impact on human health and potential application in food technology. Fermentation10:444. doi: 10.3390/fermentation10090444
43
LawsonM.O’NeillI. J.KujawskaM.Gowrinadh JavvadiS.WijeyesekeraA.FleggZ.et al. (2020). Breast milk-derived human milk oligosaccharides promote Bifidobacterium interactions within a single ecosystem. ISME J.14, 635–648. doi: 10.1038/s41396-019-0553-2
44
LeeM. D. (2019). GToTree: a user-friendly workflow for phylogenomics. Bioinformatics35, 4162–4164. doi: 10.1093/bioinformatics/btz188
45
LetunicI.BorkP. (2024). Interactive tree of life (iTOL) v6: recent updates to the phylogenetic tree display and annotation tool. Nucleic Acids Res.52, W78–W82. doi: 10.1093/nar/gkae268
46
LiuB.ZhengD.ZhouS.ChenL.YangJ. (2022). VFDB 2022: a general classification scheme for bacterial virulence factors. Nucleic Acids Res.50, D912–D917. doi: 10.1093/nar/gkab1107
47
LombardV.Golaconda RamuluH.DrulaE.CoutinhoP. M.HenrissatB. (2014). The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res.42, D490–D495. doi: 10.1093/nar/gkt1178
48
LvN.WangC.ZhouH.MaX.YuX.RenD. (2024). Genome- and toxicology-based safety assessment of probiotic Akkermansia muciniphila ONE isolated from humans. Food Secur.13:1979. doi: 10.3390/foods13131979
49
McleanK.HolmesE. A.PenewitK.LeeD. K.HardyS. R.RenM.et al. (2019). Artificial selection for pathogenicity mutations in Staphylococcus aureus identifies novel factors relevant to chronic infection. Infect. Immun.87, e00884–e00818. doi: 10.1128/iai.00884-18
50
McmurdieP. J.StoevaM. K.JusticeN.NemchekM.SieberC. M. K.TyagiS.et al. (2022). Increased circulating butyrate and ursodeoxycholate during probiotic intervention in humans with type 2 diabetes. BMC Microbiol.22:19. doi: 10.1186/s12866-021-02415-8
51
MinY. H.KwonA. R.YoonE. J.ShimM. J.ChoiE. C. (2008). Translational attenuation and mRNA stabilization as mechanisms of erm(B) induction by erythromycin. Antimicrob. Agents Chemother.52, 1782–1789. doi: 10.1128/aac.01376-07
52
NaughtonL. M.RomanoS.O’GaraF.DobsonA. D. W. (2017). Identification of secondary metabolite gene clusters in the Pseudovibrio genus reveals encouraging biosynthetic potential toward the production of novel bioactive compounds. Front. Microbiol.8:1494. doi: 10.3389/fmicb.2017.01494
53
NosterJ.HansmeierN.PersickeM.ChaoT. C.KurreR.PoppJ.et al. (2019). Blocks in tricarboxylic acid cycle of Salmonella enterica cause global perturbation of carbon storage, motility, and host-pathogen interaction. mSphere4:e00796-19. doi: 10.1128/mSphere.00796-19
54
NovickR. P.RamG. (2016). The floating (pathogenicity) island: a genomic dessert. Trends Genet.32, 114–126. doi: 10.1016/j.tig.2015.11.005
55
OgawaraH. (2019). Comparison of antibiotic resistance mechanisms in antibiotic-producing and pathogenic Bacteria. Molecules24:3430. doi: 10.3390/molecules24193430
56
Palmnäs-BédardM. S. A.CostabileG.VetraniC.ÅbergS.HjalmarssonY.DicksvedJ.et al. (2022). The human gut microbiota and glucose metabolism: a scoping review of key bacteria and the potential role of SCFAs. Am. J. Clin. Nutr.116, 862–874. doi: 10.1093/ajcn/nqac217
57
PepoyanA. (2024). Gut Akkermansia muciniphila, Prevotellaceae, and Enterobacteriaceae spp. as possible markers in women-related nutritional and clinical trials: familial Mediterranean fever disease. Womens Health Rep. (New Rochelle)5, 785–793. doi: 10.1089/whr.2024.0076
58
PepoyanE.MarottaF.ManvelyanA.GalstyanA.StepanyanL.GrigoryanH.et al. (2024). Placebo-resistant gut bacteria: Akkermansia muciniphila spp. and familial Mediterranean fever disease. Front. Cell. Infect. Microbiol.14:1336752. doi: 10.3389/fcimb.2024.1336752
59
PereiraG. V.BoudaudM.WolterM.AlexanderC.De SciscioA.GrantE. T.et al. (2024). Opposing diet, microbiome, and metabolite mechanisms regulate inflammatory bowel disease in a genetically susceptible host. Cell Host Microbe32, 527–542.e9. doi: 10.1016/j.chom.2024.03.001
60
PetersonE.KaurP. (2018). Antibiotic resistance mechanisms in Bacteria: relationships between resistance determinants of antibiotic producers, environmental Bacteria, and clinical pathogens. Front. Microbiol.9:2928. doi: 10.3389/fmicb.2018.02928
61
PittayanonR.LauJ. T.LeontiadisG. I.TseF.YuanY.SuretteM.et al. (2020). Differences in gut microbiota in patients with vs without inflammatory bowel diseases: a systematic review. Gastroenterology158, 930–946.e1. doi: 10.1053/j.gastro.2019.11.294
62
PrjibelskiA.AntipovD.MeleshkoD.LapidusA.KorobeynikovA. (2020). Using SPAdes De novo assembler. Curr. Protoc. Bioinformatics70:e102. doi: 10.1002/cpbi.102
63
RichterM.Rossello-MoraR. (2009). Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. USA106, 19126–19131. doi: 10.1073/pnas.0906412106
64
RinninellaE.RaoulP.CintoniM.FranceschiF.MiggianoG.GasbarriniA.et al. (2019). What is the healthy gut microbiota composition? A changing ecosystem across age, environment, diet, and diseases. Microorganisms7:14. doi: 10.3390/microorganisms7010014
65
RodninaM. V.FischerN.MaracciC.StarkH. (2017). Ribosome dynamics during decoding. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci.372:20160182. doi: 10.1098/rstb.2016.0182
66
RuigrokR.WeersmaR. K.Vich VilaA. (2023). The emerging role of the small intestinal microbiota in human health and disease. Gut Microbes15:2201155. doi: 10.1080/19490976.2023.2201155
67
SelvarajT.KimS. K.KimY. H.JeongY. S.KimY. J.PhuongN. D.et al. (2010). The role of carbohydrate-binding module (CBM) repeat of a multimodular xylanase (XynX) from Clostridium thermocellum in cellulose and xylan binding. J. Microbiol.48, 856–861. doi: 10.1007/s12275-010-0285-5
68
SinghT. P.NatrajB. H. (2021). Next-generation probiotics: a promising approach towards designing personalized medicine. Crit. Rev. Microbiol.47, 479–498. doi: 10.1080/1040841x.2021.1902940
69
SoodA.RayP.AngrupA. (2023). Anaerobic gram-negative Bacteria: role as a reservoir of antibiotic resistance. Antibiotics (Basel)12:942. doi: 10.3390/antibiotics12050942
70
StothardP.GrantJ. R.Van DomselaarG. (2019). Visualizing and comparing circular genomes using the CGView family of tools. Brief. Bioinform.20, 1576–1582. doi: 10.1093/bib/bbx081
71
SuS.PanmaneeW.WilsonJ. J.MahtaniH. K.LiQ.VanderwielenB. D.et al. (2014). Catalase (KatA) plays a role in protection against anaerobic nitric oxide in Pseudomonas aeruginosa. PLoS One9:e91813. doi: 10.1371/journal.pone.0091813
72
TamuraK.StecherG.KumarS.BattistuzziF. U. (2021). MEGA11: molecular evolutionary genetics analysis version 11. Mol. Biol. Evol.38, 3022–3027. doi: 10.1093/molbev/msab120
73
UrbanM.CuzickA.SeagerJ.WoodV.RutherfordK.VenkateshS. Y.et al. (2022). PHI-base in 2022: a multi-species phenotype database for pathogen-host interactions. Nucleic Acids Res.50, D837–D847. doi: 10.1093/nar/gkab1037
74
van DijkB. (2020). Can mobile genetic elements rescue genes from extinction?Curr. Genet.66, 1069–1071. doi: 10.1007/s00294-020-01104-9
75
WangL.TangL.FengY.ZhaoS.HanM.ZhangC.et al. (2020). A purified membrane protein from Akkermansia muciniphila or the pasteurised bacterium blunts colitis associated tumourigenesis by modulation of CD8+ T cells in mice. Gut69, 1988–1997. doi: 10.1136/gutjnl-2019-320105
76
WeberH.EngelmannS.BecherD.HeckerM. (2004). Oxidative stress triggers thiol oxidation in the glyceraldehyde-3-phosphate dehydrogenase of Staphylococcus aureus. Mol. Microbiol.52, 133–140. doi: 10.1111/j.1365-2958.2004.03971.x
77
WongJ. M.De SouzaR.KendallC. W.EmamA.JenkinsD. J. (2006). Colonic health: fermentation and short chain fatty acids. J. Clin. Gastroenterol.40, 235–243. doi: 10.1097/00004836-200603000-00015
78
XiaY.WengY.XuC.WangD.PanX.TianZ.et al. (2020). Endoribonuclease YbeY is essential for RNA processing and virulence in Pseudomonas aeruginosa. MBio11, e00659–e00620. doi: 10.1128/mBio.00659-20
79
YangM.BoseS.LimS.SeoJ.ShinJ.LeeD.et al. (2020). Beneficial effects of newly isolated Akkermansia muciniphila strains from the human gut on obesity and metabolic dysregulation. Microorganisms8:1413. doi: 10.3390/microorganisms8091413
80
ZhaoX.ZhaoJ.LiD.YangH.ChenC.QinM.et al. (2023). Akkermansia muciniphila: a potential target and pending issues for oncotherapy. Pharmacol. Res.196:106916. doi: 10.1016/j.phrs.2023.106916
81
ZhengJ.GeQ.YanY.ZhangX.HuangL.YinY. (2023). dbCAN3: automated carbohydrate-active enzyme and substrate annotation. Nucleic Acids Res.51, W115–W121. doi: 10.1093/nar/gkad328
Summary
Keywords
Akkermansia, whole-genome sequencing, genomic analysis, genomic diversity, probiotics
Citation
Lu W, Zha B, Lyu J, LingHu C, Chen J, Deng S, Zhang X, Li L and Wang G (2024) Whole-genome sequencing and genomic analysis of four Akkermansia strains newly isolated from human feces. Front. Microbiol. 15:1500886. doi: 10.3389/fmicb.2024.1500886
Received
24 September 2024
Accepted
25 November 2024
Published
16 December 2024
Volume
15 - 2024
Edited by
Astghik Zaveni Pepoyan, Armenian National Agrarian University, Armenia
Reviewed by
Joana Cristina Barbosa, Centre for Biotechnology and Fine Chemistry - Universidade Católica Portuguesa, Portugal
Jéssica Assis Pereira, Federal University of Alfenas, Brazil
Updates
Copyright
© 2024 Lu, Zha, Lyu, LingHu, Chen, Deng, Zhang, Li and Wang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Liang Li, lil16@by-health.comGuoqing Wang, huaxiwgq@163.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.