Exploration of Survival Traits, Probiotic Determinants, Host Interactions, and Functional Evolution of Bifidobacterial Genomes Using Comparative Genomics

Multi-omics of the gut microbial ecosystem in patients with microsatellite-instability-high gastrointestinal cancer resistant to immunotherapy

Despite the encouraging efficacy of anti-PD-1/PD-L1 immunotherapy in microsatellite-instability-high/deficient mismatch repair (MSI-H/dMMR) advanced gastrointestinal cancer, many patients exhibit primary or acquired resistance. Using multi-omics approaches, we interrogate gut microbiome, blood metabolome, and cytokines/chemokines of patients with MSI-H/dMMR gastrointestinal cancer (N = 77) at baseline and during the treatment. We identify a number of microbes (e.g., Porphyromonadaceae) and metabolites (e.g., arginine) highly associated with primary resistance to immunotherapy. An independent validation cohort (N = 39) and mouse model are used to further confirm our findings. A predictive machine learning model for primary resistance is also built and achieves an accuracy of 0.79 on the external validation set. Furthermore, several microbes are pinpointed that gradually changed during the process of acquired resistance. In summary, our study demonstrates the essential role of gut microbiome in drug resistance, and this can be utilized as a preventative diagnosis tool and therapeutic target in the future.

Meta-analysis reveals different functional characteristics of human gut Bifidobacteria associated with habitual diet

Dietary habits contribute to the composition and function of the gut microbiota. Different dietary structures, including vegan, vegetarian, and omnivorous diets, affect intestinal Bifidobacteria; however, the relationship between Bifidobacterial function and host metabolism in subjects with different dietary patterns is unclear. Here, we analyzed five metagenomics studies and six 16S sequencing studies, including 206 vegetarians (VG), 249 omnivores (O), and 270 vegans (V), through an unbiased theme-level meta-analysis framework and discovered that diet significantly affects the composition and functionality of intestinal Bifidobacteria. The relative abundance of Bifidobacterium pseudocatenulatum was significantly higher in V than in O and Bifidobacterium longum, Bifidobacterium adolescentis, and B. pseudocatenulatum differed significantly in carbohydrate transport and metabolism in subjects with different diet types. Diets high in fiber were associated with B. longum with increased capacity for carbohydrate catabolism and genes encoding GH29 and GH43_27 were significantly enriched in V. Bifidobacterium adolescentis and B. pseudocatenulatum, associated with O, had a higher prevalence of the genes related to carbohydrate transport and metabolism, which showed the enrichment of GH26 and GH27 families. The same Bifidobacterium species has different functions in subjects with different diet types, resulting in different physiological significance. The diversification and functionalities of Bifidobacterial species in the gut microbiome can be influenced by the host diet and this aspect should be considered when studying host-microbe associations.

Identification of conserved genomic signatures specific to Bifidobacterium species colonising the human gut

Bifidobacterium species are known for their ability to inhabit various habitats and are often regarded as the first colonisers of the human gut. In the present work, we have used comparative genomics to identify conserved genomic signatures specific to Bifidobacterium species associated with the human gut. Our approach discovered five genomic signatures with varying lengths and confidence. Among the predicted five signatures, a 1790 bp multi-drug resistance (MDR) signature was found to be remarkably specific to only those species that can colonise the human gut. The signature codes for a membrane transport protein belonging to the major facilitator superfamily (MFS) generally involved in MDR. Phylogenetic analyses of the MDR signature suggest a lineage-specific evolution of the MDR signature in bifidobacteria colonising the human gut. Functional annotation led to the discovery of two conserved domains in the protein; a catalytic MFS domain involved in the efflux of drugs and toxins, and a regulatory cystathionine-β-synthase (CBS) domain that can interact with adenosyl-carriers. Molecular docking simulation performed with the modelled tertiary structure of the MDR signature revealed the putative functional role of the covalently linked domains. The MFS domain displayed a high affinity towards various protein synthesis inhibitor antibiotics and human bile acids, whereas the C-terminally linked CBS domain exhibited favourable binding with molecular structures of ATP and AMP. Therefore, we believe that the predicted signature represents a niche-specific survival trait involved in bile and antibiotic resistance, imparting an adaptive advantage to the Bifidobacterium species colonising the human gut.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03492-4.

Comparative genomics reveals the evolution of antimicrobial resistance in Bacteroides nordii

Bacteroides nordii, is an understudied member of the pathogenic B. fragilis group which comprises several multidrug-resistant (MDR) strains. Thus, it is of great interest to study the genome biology of Bacteroides nordii. However, no detailed study is available that characterized B. nordii at the genetic level and explored its role as a potential pathogen. We isolated an MDR strain viz., B. nordii PGMM4098 from the pus sample and subjected it to whole genome sequencing using Illumina technology. The draft genome was de-novo assembled and annotated, followed by comprehensive comparative genomics analyses using the publicly available genome dataset of B. nordii. The pan-genome analysis revealed the open nature of B. nordii, indicating the continuous accumulation of novel genes in non-core components leading to the emergence of new strains of this species. The thirteen antimicrobial resistance (AMR) genes identified in the genomes of all B. nordii strains were part of the non-core component of the pan-genome. Of these, four AMR genes, nimE, aadS, mef(En2), and ermB/F/G were found to be acquired via the process of horizontal gene transfer (HGT) from anaerobic Bacteroidetes. Importantly, the nimE gene conferring metronidazole resistance was found to be present only in B. nordii PGMM4098, which harbors five other AMR genes encoded in its genome. Of these, nimE (metronidazole resistance), ermB/F/G (macrolide-lincosamide-streptogramin B resistance), and cfxA2/A3 (class A β-lactam resistance) genes were further validated using targeted polymerase chain reaction assay. Notably, these three genes were also found to be under the operation of positive selective pressure suggesting the diversification of these genes, which might lead to the emergence of new MDR strains of B. nordii in the near future. Our study reported and characterized the genome of the first MDR strain of B. nordii and revealed the AMR evolution in this species using a comprehensive comparative genomics approach.

Genome-Wide Assessment of Stress-Associated Genes in Bifidobacteria

Over the last decade, the genomes of several Bifidobacterium strains have been sequenced, delivering valuable insights into their genetic makeup. However, bifidobacterial genomes have not yet been systematically mined for genes associated with stress response functions and their regulation. In this work, a list of 76 genes related to stress response in bifidobacteria was compiled from previous studies. The prevalence of the genes was evaluated among the genome sequences of 171 Bifidobacterium strains. Although genes of the protein quality control and DNA repair systems appeared to be highly conserved, genome-wide in silico screening for consensus sequences of putative regulators suggested that the regulation of these systems differs among phylogenetic groups. Homologs of multiple oxidative stress-associated genes are shared across species, albeit at low sequence similarity. Bee isolates were confirmed to harbor unique genetic features linked to oxygen tolerance. Moreover, most studied Bifidobacterium adolescentis and all Bifidobacterium angulatum strains lacked a set of reactive oxygen species-detoxifying enzymes, which might explain their high sensitivity to oxygen. Furthermore, the presence of some putative transcriptional regulators of stress responses was found to vary across species and strains, indicating that different regulation strategies of stress-associated gene transcription contribute to the diverse stress tolerance. The presented stress response gene profiles of Bifidobacterium strains provide a valuable knowledge base for guiding future studies by enabling hypothesis generation and the identification of key genes for further analyses.

IMPORTANCE Bifidobacteria are Gram-positive bacteria that naturally inhabit diverse ecological niches, including the gastrointestinal tract of humans and animals. Strains of the genus Bifidobacterium are widely used as probiotics, since they have been associated with health benefits. In the course of their production and administration, probiotic bifidobacteria are exposed to several stressors that can challenge their survival. The stress tolerance of probiotic bifidobacteria is, therefore, an important selection criterion for their commercial application, since strains must maintain their viability to exert their beneficial health effects. As the ability to cope with stressors varies among Bifidobacterium strains, comprehensive understanding of the underlying stress physiology is required for enabling knowledge-driven strain selection and optimization of industrial-scale production processes.

Complete Genome Sequencing and Comparative Genomics of Three Potential Probiotic Strains, Lacticaseibacillus casei FBL6, Lacticaseibacillus chiayiensis FBL7, and Lacticaseibacillus zeae FBL8

Lacticaseibacillus casei, Lacticaseibacillus chiayiensis, and Lacticaseibacillus zeae are very closely related Lacticaseibacillus species. L. casei has long been proposed as a probiotic, whereas studies on functional characterization for L. chiayiensis and L. zeae are some compared to L. casei. In this study, L. casei FBL6, L. chiayiensis FBL7, and L. zeae FBL8 were isolated from raw milk, and their probiotic properties were investigated. Genomic analysis demonstrated the role of L. chiayiensis and L. zeae as probiotic candidates. The three strains were tolerant to acid and bile salt, with inhibitory action against pathogenic bacterial strains and capacity of antioxidants. Complete genome sequences of the three strains were analyzed to highlight the probiotic properties at the genetic level, which results in the discovery of genes corresponding to phenotypic characterization. Moreover, genes known to confer probiotic characteristics were identified, including genes related to biosynthesis, defense machinery, adhesion, and stress adaptation. The comparative genomic analysis with other available genomes revealed 256, 214, and 32 unique genes for FBL6, FBL7, and FBL8, respectively. These genomes contained individual genes encoding proteins that are putatively involved in carbohydrate transport and metabolism, prokaryotic immune system for antiviral defense, and physiological control processes. In particular, L. casei FBL6 had a bacteriocin gene cluster that was not present in other genomes of L. casei, resulting in this strain may exhibit a wide range of antimicrobial activity compared to other L. casei strains. Our data can help us understand the probiotic functionalities of the three strains and suggest that L. chiayiensis and L. zeae species, which are closely related to L. casei, can also be considered as novel potential probiotic candidate strains.

The bifidobacterial distribution in the microbiome of captive primates reflects parvorder and feed specialization of the host

Bifidobacteria, which commonly inhabit the primate gut, are beneficial contributors to host wellbeing. Anatomical differences and natural habitat allow an arrangement of primates into two main parvorders; New World monkeys (NWM) and Old World monkeys (OWM). The number of newly described bifidobacterial species is clearly elevated in NWM. This corresponds to our finding that bifidobacteria were the dominant group of cultivated gut anaerobes in NWM, while their numbers halved in OWM and were often replaced by Clostridiaceae with sarcina morphology. We examined an extended MALDI-TOF MS database as a potential identification tool for rapid screening of bifidobacterial distribution in captive primates. Bifidobacterial isolates of NWM were assigned mainly to species of primate origin, while OWM possessed typically multi-host bifidobacteria. Moreover, bifidobacterial counts reflected the feed specialization of captive primates decreasing from frugivore-insectivores, gummivore-insectivores, frugivore-folivores to frugivore-omnivores. Amplicon sequencing analysis supported this trend with regards to the inverse ratio of Actinobacteria and Firmicutes. In addition, a significantly higher diversity of the bacterial population in OWM was found. The evolution specialization of primates seems to be responsible for Bifidobacterium abundance and species occurrence. Balanced microbiota of captive primates could be supported by optimized prebiotic and probiotic stimulation based on the primate host.

Invited review: Characterization of new probiotics from dairy and nondairy products-Insights into acid tolerance, bile metabolism and tolerance, and adhesion capability

The selection of potential probiotic strains that possess the physiological capacity of performing successfully in the gastrointestinal tract (GIT) is a critical challenge. Probiotic microorganisms must tolerate the deleterious effects of various stresses to survive passage and function in the human GIT. Adhesion to the intestinal mucosa is also an important aspect. Recently, numerous studies have been performed concerning the selection and evaluation of novel probiotic microorganisms, mainly probiotic bacteria isolated from dairy and nondairy products. Therefore, it would be crucial to critically review the assessment methods employed to select the potential probiotics. This article aims to review and discuss the recent approaches, methods used for the selection, and outcomes of the evaluation of novel probiotic strains with the main purpose of supporting future probiotic microbial assessment studies. The findings and approaches used for assessing acid tolerance, bile metabolism and tolerance, and adhesion capability are the focus of this review. In addition, probiotic bile deconjugation and bile salt hydrolysis are explored. The selection of a new probiotic strain has mainly been based on the in vitro tolerance of physiologically related stresses including low pH and bile, to ensure that the potential probiotic microorganism can survive the harsh conditions of the GIT. However, the varied experimental conditions used in these studies (different types of media, bile, pH, and incubation time) hamper the comparison of the results of these investigations. Therefore, standardization of experimental conditions for characterizing and selecting probiotics is warranted.

Comparative gut microbiome analysis of the Prakriti and Sasang systems reveals functional level similarities in constitutionally similar classes

The traditional medicinal systems (TMS) of India (Prakriti) and Korea (Sasang) classify human individuals based on their constitution determined by the physiological and psychological traits of individuals. Similarities in the constitutions are already found between the classes of Prakriti (Vata, Pitta, and Kapha) and Sasang (TE: Taeeumin, SE: Soeumin, and SY: Soyangin) systems. Gut health is an important aspect of this constitution based classification in TMS. To determine the role of gut microbes in such classifications, we have analyzed the gut microbiome (taxa and imputed functions) in the constitutionally similar Prakriti and Sasang classes. An enrichment of Bacteroides and Prevotella enterotypes is observed in the Sasang and Prakriti samples, respectively. The impact of the constitution is found to be more prominent with respect to the taxa and predicted-functions within the Prakriti classes. Gut microbiome functional-level similarities are found to correlate well with the host phenotypes of the constitutionally similar Prakriti and Sasang classes. An enrichment of carbohydrate and amino-acid metabolism is observed in the Vata and SE classes which may be responsible for meeting with their high energy demands and lean phenotype. The Pitta and SY classes exhibit the high capacity to metabolize toxins. An enrichment of functions responsible for predisposition to obesity and high drug metabolism is observed in the Kapha and TE classes. The contribution of gut adaptive functions is found to correlate with the constitution-based classification in both Prakriti and Sasang systems. The TE class harboured the highest number of biofilm-forming and stress-tolerant microbes thus exhibiting the maximum tolerance of environmental stress. Similarities in the gut microbiota and the resulting disease predisposition patterns are found to exist between the constitutionally matching Prakriti and Sasang classes.

Bifidobacterium β-Glucosidase Activity and Fermentation of Dietary Plant Glucosides Is Species and Strain Specific

Dietary plant glucosides are phytochemicals whose bioactivity and bioavailability can be modified by glucoside hydrolase activity of intestinal microbiota through the release of acylglycones. Bifidobacteria are gut commensals whose genomic potential indicates host-adaption as they possess a diverse set of glycosyl hydrolases giving access to a variety of dietary glycans. We hypothesized bifidobacteria with β-glucosidase activity could use plant glucosides as fermentation substrate and tested 115 strains assigned to eight different species and from different hosts for their potential to express β-glucosidases and ability to grow in the presence of esculin, amygdalin, and arbutin. Concurrently, the antibacterial activity of arbutin and its acylglycone hydroquinone was investigated. Beta-glucosidase activity of bifidobacteria was species specific and most prevalent in species occurring in human adults and animal hosts. Utilization and fermentation profiles of plant glucosides differed between strains and might provide a competitive benefit enabling the intestinal use of dietary plant glucosides as energy sources. Bifidobacterial β-glucosidase activity can increase the bioactivity of plant glucosides through the release of acylglycone.

Evolutionary relationships among bifidobacteria and their hosts and environments

BACKGROUND

The assembly of animal microbiomes is influenced by multiple environmental factors and host genetics, although the relative importance of these factors remains unclear. Bifidobacteria (genus Bifidobacterium, phylum Actinobacteria) are common first colonizers of gut microbiomes in humans and inhabit other mammals, social insects, food, and sewages. In humans, the presence of bifidobacteria in the gut has been correlated with health-promoting benefits. Here, we compared the genome sequences of a subset of the over 400 Bifidobacterium strains publicly available to investigate the adaptation of bifidobacteria diversity. We tested 1) whether bifidobacteria show a phylogenetic signal with their isolation sources (hosts and environments) and 2) whether key traits encoded by the bifidobacteria genomes depend on the host or environment from which they were isolated. We analyzed Bifidobacterium genomes available in the PATRIC and NCBI repositories and identified the hosts and/or environment from which they were isolated. A multilocus phylogenetic analysis was conducted to compare the genetic relatedness the strains harbored by different hosts and environments. Furthermore, we examined differences in genomic traits and genes related to amino acid biosynthesis and degradation of carbohydrates.

RESULTS

We found that bifidobacteria diversity appears to have evolved with their hosts as strains isolated from the same host were non-randomly associated with their phylogenetic relatedness. Moreover, bifidobacteria isolated from different sources displayed differences in genomic traits such as genome size and accessory gene composition and on particular traits related to amino acid production and degradation of carbohydrates. In contrast, when analyzing diversity within human-derived bifidobacteria, we observed no phylogenetic signal or differences on specific traits (amino acid biosynthesis genes and CAZymes).

CONCLUSIONS

Overall, our study shows that bifidobacteria diversity is strongly adapted to specific hosts and environments and that several genomic traits were associated with their isolation sources. However, this signal is not observed in human-derived strains alone. Looking into the genomic signatures of bifidobacteria strains in different environments can give insights into how this bacterial group adapts to their environment and what types of traits are important for these adaptations.

Enterotype Variations of the Healthy Human Gut Microbiome in Different Geographical Regions

Enterotypes are used for classifying individuals based on the gut microbiome. A number of studies are available to find the Enterotypes in healthy individuals; however, most of them lack comparisons at the world level. We analyzed the healthy human gut microbiomes of 495 datasets available in the European Nucleotide Archive (ENA) database derived from fifteen countries from four continents. Firmicutes and Bacteroidetes were the two most abundant phyla in the healthy human gut, worldwide. A high ratio of Proteobacteriato Actinobacteria and a low abundance of Prevotella were identified as the indicators of IBD. Prevotella, Bacteroides, and Bifidobacterium were identified as the Enterotypes in the inter-continental comparisons. At the intra-continental level, two (Bacteroides and Ruminococcaceae), four (Faecalibacterium, Bacteroides, Prevotella, and Clostridiales), and two (Prevotella, Bacteroides/Bifidobacterium) Enterotypes were identified in the American, European, and Asian continents, respectively. In addition, a high abundance of the unknown genus of Ruminococcaeae was observed in the Colombian human gut microbiome. A substantial impact of the geographical distance was observed on human gut microbiome variations, demonstrating a cumulative effect of factors, including dietary habits, genetics, lifestyle, environment, and climate, etc.