Next Article in Journal
The Bacterial Amyloid-Like Hfq Promotes In Vitro DNA Alignment
Next Article in Special Issue
Comparative Pangenomics of the Mammalian Gut Commensal Bifidobacterium longum
Previous Article in Journal
Paenibacillus lutrae sp. nov., A Chitinolytic Species Isolated from A River Otter in Castril Natural Park, Granada, Spain
Previous Article in Special Issue
Bifidobacterial Distribution Across Italian Cheeses Produced from Raw Milk
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 7
Issue 12
10.3390/microorganisms7120638
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mobilome and Resistome Reconstruction from Genomes Belonging to Members of the Bifidobacterium Genus

by
Walter Mancino
1,
Gabriele Andrea Lugli
1,
Douwe van Sinderen
2,
Marco Ventura
1,3 and
Francesca Turroni
1,3,*
1
Laboratory of Probiogenomics, Department of Chemistry, Life Sciences and Environmental Sustainability, University of Parma, 43124 Parma, Italy
2
School of Microbiology, APC Microbiome Institute, University College Cork, Cork T12 K8AF, Ireland
3
Microbiome Research Hub, University of Parma, 43124 Parma, Italy
*
Author to whom correspondence should be addressed.
Microorganisms 2019, 7(12), 638; https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms7120638
Submission received: 29 October 2019 / Revised: 27 November 2019 / Accepted: 29 November 2019 / Published: 2 December 2019
(This article belongs to the Special Issue Bifidobacteria: Insights from Ecology to Genomics of a Key Microbial Group of the Mammalian Gut Microbiota)
Browse Figures
Review Reports Versions Notes

Abstract

:
Specific members of the genus Bifidobacterium are among the first colonizers of the human/animal gut, where they act as important intestinal commensals associated with host health. As part of the gut microbiota, bifidobacteria may be exposed to antibiotics, used in particular for intrapartum prophylaxis, especially to prevent Streptococcus infections, or in the very early stages of life after the birth. In the current study, we reconstructed the in silico resistome of the Bifidobacterium genus, analyzing a database composed of 625 bifidobacterial genomes, including partial assembled strains with less than 100 genomic sequences. Furthermore, we screened bifidobacterial genomes for mobile genetic elements, such as transposases and prophage-like elements, in order to investigate the correlation between the bifido-mobilome and the bifido-resistome, also identifying genetic insertion hotspots that appear to be prone to horizontal gene transfer (HGT) events. These insertion hotspots were shown to be widely distributed among analyzed bifidobacterial genomes, and suggest the acquisition of antibiotic resistance genes through HGT events. These data were further corroborated by growth experiments directed to evaluate bacitracin A resistance in Bifidobacterium spp., a property that was predicted by in silico analyses to be part of the HGT-acquired resistome.

1. Introduction

Bifidobacteria are gram-positive, anaerobic, non-motile, and non-spore-forming bacteria with a high G + C genomic content [1]. They represent one of the dominant microbial groups inhabiting the gastrointestinal tract (GIT) of humans and animals, including mammals, birds, and social insects [2,3,4]. Members of the Bifidobacterium genus are believed to be crucial for the development of a healthy gut microbiota in early life, colonizing the GIT within the first days following birth [5]. Notably, exposure to antibiotic agents as a result of intrapartum prophylaxis, commonly applied to prevent group B Streptococcus infections [6,7] and in infants, presenting a highly dynamic microbiota [8], may disrupt the balance between microbial members of the gut microbiota [9,10,11]. It has previously been demonstrated that the relative abundance of bifidobacterial species in the gut microbiota of infants who had not been exposed to any antibiotic treatment is higher than that in children that had been subjected to antibiotic therapy [12,13]. The extensive use of antibiotics can promote the development of antibiotic resistance in members of the microbiota and consequently in the selection of antibiotic-resistant microorganisms [11]. In this context, the collective genetic arsenal responsible for conferring antibiotic resistance (AR) through inactivation and/or removal of antibiotics is commonly referred to as the resistome [14,15,16]. The occurrence of AR genes may increase the ecological fitness of a bacterium and thus its ability to colonize and persist in a specific environment [17]. Recent studies have indicated frequent occurrence of horizontal gene transfer (HGT) events among bacteria residing in the gut of humans and animals [18,19]. The presence of AR genes in mobile genetic elements (MGEs) or near transposable elements, in pathogenic and non-pathogenic microorganisms, may be the cause for the relatively frequent transfer of such elements to human/animal pathogens or to other non-resident microorganisms of the gastrointestinal tract [19,20].
The mobilome of a microorganism refers to the collection of all MGEs, including transposases, insertion elements (IS elements), and also plasmids and prophages [21,22,23]. A recent study investigating the presence of prophages in 48 members of the genus Bifidobacterium predicted the presence of 90 different prophages, called bifidoprophages [24].
In the current study, we reconstructed the entire mobilome of the genus Bifidobacterium, including prophages and transposases, based on 625 different bifidobacterial genomes belonging to 67 different (sub)species [25]. Moreover, we reconstructed the bifidobacterial resistome, i.e., the genes whose products are predicted to be responsible for resistance against antibiotic molecules. Combining the gathered bifido-mobilome with the bifido-resistome data, we identified genetic insertion signatures that may be involved in horizontal transfer of AR genes in bifidobacterial genomes.

2. Materials and Methods

2.1. Bacterial Strains

We retrieved the genome sequence of 625 public available Bifidobacterium genomes from the National Center for Biotechnology Information (NCBI) public database (Table S1). Collected genomes with more than 100 genomic sequences were discarded to analyze high quality genome sequences only. As reported in Table S1, strains that were not classified at the species level were validated using the average nucleotide identity (ANI) approach. Strains used for this analysis were compared with the 67 type strains of the Bifidobacterium genus. Notably, two bifidobacterial strains displaying an ANI value of <94% may be considered to belong to two different species [25,26,27,28].

2.2. IS Elements Identification

Predicted genes of 625 bifidobacterial strains used for this study were screened for the presence of IS elements. We used a custom database composed of 329,372 RefSeq sequences belonging to the Actinobacteria phylum, retrieved from the NCBI database. The alignment was performed through BLASTP analysis with an E-value cutoff of 1e−5 [29]. After the manual control of the sequences with an amino acid length less than 100 amino acids of the type strains of the species studied, we decided to discard these sequences because they were considered non-functional or truncated. Finally, the selected IS element sequences were validated and classified into IS families using the IS finder database [30].

2.3. Bifidoprophages Identification

The 625 bifidobacterial genomes were screened for prophage-like elements using a custom database based on already identified sequences through BLASTP analysis [29] (E-value cutoff of 1e−5). The custom database was constructed through previously bifidoprophage-validated sequences retrieved from 60 bifidoprophages identified by Lugli et al. [24], considering genetic islands presenting different genes encoding for phage functions. Following this, a manual examination of the DNA region surrounding a putative phage-encoding gene was performed. These manual screenings allowed us to identify complete prophage-like sequences while discarding incomplete or remnant phage sequences, as previously performed by Lugli et al. [24].

2.4. Prediction of the Antibiotic Resistance Genes

The in silico proteome of the 625 Bifidobacterium genomes used in this study was screened for proteins that can act as antibiotic resistance proteins through inactivation and/or removal of antibiotic molecules. The screening was carried out using the MEGAREs database through BLASTP analysis (E-value cutoff of 1e−18) [29,31]. The E-value cutoff was chosen based on a manual editing performed to identify false positive sequences. The core database was obtained by non-redundant compilation of sequences contained in Resfinder, ARG-ANNOT, the Comprehensive Antibiotic Resistance Database (CARD), and the NCBI Lahey Clinic beta-lactamase archive [32,33,34,35]. Following this, a manual examination of the sequence with an E-value less than 1e−18 was performed in order to explore all the biodiversity of the AR genes of the Bifidobacterium species. We excluded the putative AR genes encoding for transporters for low accuracy in their prediction [16]. The predicted AR genes were classified according to the presumed mechanism of action and the antibiotic molecules they counteract.
Moreover, for the 625 bifidobacterial genomes analyzed, we manually evaluated the genes flanked by the predicted AR genes, forming the Bifidobacterium resistome, in order to identify mobile genetic hotspots that may promote HGT events.

2.5. Phylogenomic Analyses

The nucleotide similarity of each obtained bifidoprophage sequence was calculated using the software package LAST [36]. Results were employed to build a matrix representing the genome similarity among different prophage and to generate a clustering tree. The bifidophage sequences were aligned using Mafft software [37] and the clustering tree was constructed using ClustalW [38]. The constructed clustering tree was visualized using the FIGTREE software (http://tree.bio.ed.ac.uk/software/figtree/).

2.6. Bacitracin A Antibiotic Susceptibility Tests

The minimal inhibitory concentration (MIC) breakpoints (micrograms per milliliter) of bacitracin A were determined using the broth microdilution method (MDIL) according to the ISO standard guidelines [39]. Bacitracin A antibiotic was purchased from Merck (Germany). Microplates were incubated under anaerobic conditions for 48 h at 37 °C. Cell density was monitored by optical density measurements at 600 nm (OD600) using a plate reader (BioTek, Vermont, USA). The MIC breakpoint represents the highest concentration of a given antibiotic to which a particular bacterial strain is resistant.

2.7. Statistical Analyses

SPSS software (IBM, Italy) was used to perform statistical analysis between the BacA strains group and control group by T-student test.

3. Results and Discussion

3.1. The Putative Resistome of the Genus Bifidobacterium

In order to investigate the genetic AR arsenal carried by members of the Bifidobacterium genus, we investigated the resistome of 625 bifidobacterial genomes. We enlarged the previously published database on the resistome of the Bifidobacterium genus, which were based on 91 different genomes [16]. Putative AR genes encoding transporters were excluded from this analysis due to the inaccuracy of their bioinformatic prediction [16]. The overall number of putative antibiotic-resistance genes identified among these 625 genomes was 13,870, representing less than 1% of the total Bifidobacterium genes analyzed (Table S1). According to the predicted mechanism of action and the antibiotic molecules that could be counteracted, seven different AR gene classes were identified (Figure 1). The AR class with the highest number of representatives was the one conferring glycopeptide resistance, which corresponds to 5999 putative enzymes acting against glycopeptide antibiotics, such as vancomycin, teicoplanin, and telavancin (Figure 1) [40,41,42]. Notably, Bifidobacterium bifidum 791, Bifidobacterium longum subsp. infantis 1888B and B. bifidum AM42-15AC were strains containing the highest number of genes predicted to belong to this glycopeptide-resistance class, each encoding 29 distinct enzymes predicted to confer such resistance. Moreover, we identified 2178 genes putatively belonging to the tetracycline-resistance class (Figure 1) [43,44,45,46]. Members of the B. bifidum species, isolated from fecal samples of healthy Chinese individuals [47], i.e., strains TM05-15, TF07-22, TM02-15, TM02-17, TM06-10, and TM07-4AC, were shown to contain the highest number of genes encoding proteins predicted to counteract tetracycline antibiotics, ranging from 29 genes of B. bifidum TM05-15 to 28 genes for the other B. bifidum strains. Notably, 484 analyzed strains did not appear to contain genes encoding tetracycline-resistance proteins, representing 77.5% of the total Bifidobacterium strains analyzed. Furthermore, 2437 genes were found to belong to the beta-lactamase class and Bifidobacterium animalis subsp. animalis ATCC 25527 was shown to be the strain with the highest number (i.e., 32) of predicted beta-lactamase-encoding genes, while 469 of the 625 analyzed genomes did not appear to encompass genes belonging to this AR class (Figure 1).
Moreover, 2618 genes were predicted to belong to the methyltransferase AR class, including 23S rRNA methyltransferase, which may confer resistance toward erythromycin and clindamycin, as demonstrated in a previous study [48] (Figure 1). In addition, we identified 500 genes predicted to belong to the sulfonamide-resistance class, which includes genes encoding enzymes counteracting sulfonamide antimicrobial agents, also known as sul genes [49] (Figure 1). The sul gene appears to be present as three variants in the investigated genomes, i.e., sul1, sul2, and sul3, all encoding a dihydropteroate synthase [49,50]. Interestingly, in the assessed genomes of the Bifidobacterium genus, the most prevalent gene variant was sul3, found in 90.4% of all identified sulfonamide-resistance genes.
Finally, the aminoglycoside class and the metronidazole class were the two least represented classes of AR genes in bifidobacteria, with just 73 and 64 identified genes predicted to be members of these two respective classes (Figure 1). Notably, B. bifidum AF11-25B was predicted to contain the highest number of genes encoding enzymes that counteract aminoglycoside antibiotics, such as streptomycin, kanamycin, and gentamicin.
Moreover, B. bifidum TF05-1 was the only strain whose chromosome contains a gene encoding a putative quinolone-resistance protein. This gene encodes a pentapeptide repeat protein, which is predicted to be involved in fluoroquinolone resistance [51,52,53].
Interestingly, comparing the identified Bifidobacterium resistome with AR determinants of other gut commensal, such as members of the Lactobacillus genus, we observed a lesser complexity of the resistome [54]. In fact, the Lactobacillus genus included different genes that could confer resistance toward a wide range of antibiotic molecules, such as vancomycin, erythromycin, and penicillin, but also tetracycline, chloramphenicol, and aminoglycoside antibiotics [55,56,57,58,59]. Furthermore, different Escherichia coli strains presented in their genomes AR genes that counteracted carbapenem antibiotics [60,61], whereas bifidobacteria seemed to be very sensitive to this antibiotic class, and their genomes do not encompass any genes that could confer resistance toward this antibiotic. Notably, a recent study based on metagenomics analyses of the human gut microbiota revealed that Entrococcus and Enterobacter genera possessed a very high antibiotic resistance load [62]. These genera presented AR genes that could counteract different antibiotics, such as trimethoprim/sulfamethoxazole, metronizadole, cycloserine, and cefixime [62]. Moreover, different studies have demonstrated the presence of AR genes in the genomes of the members of Bacillus genus, used as probiotic bacteria in functional food and for animal feed [54,63]. In the latter genus, macrolide-resistance genes have been identified on extra-chromosomal elements, tetracycline resistance genes, but also cfr-like genes (i.e., conferring resistance toward several classes of antibiotics, including phenicols, oxalozidinone, lincosamides, pleuromutilinis, and streptogramin) that have not been identified in the genomes of the members of the Bifidobacterium genus [64,65,66]. Our resistome analyses revealed a lack of specific Bifidobacterium AR genes, corroborating the safer behavior of the Bifidobacterium genus compared to other human gut commensals.
Although we do acknowledge the limitations of the in silico analysis in assigning antibiotic resistance functions to these identified genes, they are nonetheless considered to represent a potential arsenal to counter antimicrobial molecules.

3.2. The Predicted Mobilome of the Bifidobacterium Genus

The mobilome is defined as genetic elements that can move within a genome and between different genomes, including transposable elements, bacteriophages, and plasmids [21,22,23]. Similar to other members of the gut microbiota, it has been demonstrated that bifidobacteria possess genetic elements whose action is responsible for shaping their genomes [23,24,67,68]. In order to explore the mobile element repertoire of the Bifidobacterium genus, we analyzed the same 625 Bifidobacterium genomes as indicated above (Table S1). Our analyses updated previously published data based on the reconstruction of 60 different bifidoprophage-like elements of 48 species of the Bifidobacterium genus.
A screening among the analyzed bifidobacterial genomes revealed 16,065 different genes encoding transposases, excluding genes that were truncated at the start or end codon (Table S1). The sequence of each IS element was classified according to the ISFinder database [69], showing that members of the IS3 family are the most widespread among the Bifidobacterium genus (Table S2). Notably, members of the Bifidobacterium breve species showed the highest number of IS elements, i.e., strains BR-06, BR-H29, BR-21, BR-L29, and BR-C29, ranging from 174 to 102 (Table S1). Moreover, 16.5% of the analyzed genomes were predicted to contain less than 10 genes encoding transposases in their chromosomes while Bifidobacterium commune LMG 28292 does not appear to encompass any IS element at all (Table S1).
Recently, Lugli et. al. recognized and classified all prophage-like elements (referred to as bifidoprophages) present in 48 genomes of type strains belonging to different bifidobacterial species [24] and Mavrich et al. characterized three of these identified groups of prophages integrated in members of B. breve and B. longum species by means of induction experiments [70]. In the current study, the screening for bifidoprophages was further extended to 625 different bifidobacterial genomes, resulting in the identification of 598 putative and apparently complete prophage sequences (Table S1). Notably, the genomes of Bifidobacterium biavatii DSM 23969, Bifidobacterium imperatoris LMG 30297, and Bifidobacterium cuniculi LMG 10738 were predicted to contain the highest number of prophages in their genomes, i.e., seven, six, and five prophage-like elements, respectively (Table S1). In order to evaluate the homology among the identified prophage-like elements, a genomic-based alignment clustering was performed. We observed the presence of four main homology clusters, in which the taxonomic origin of the corresponding Bifidobacterium hosts was highly heterogeneous. Each identified cluster showed several sub-clusters consisting of different prophage-like elements belonging to bifidobacterial strains of the same species, highlighting a sub-cluster phage specificity that appears to be host related (Figure 2). As reported in previous studies, prophages contribute to the genetic individuality of bacterial strains, containing many unique genes that in some instances may confer a fitness advantage to the host, such as a gene related to antibiotic resistance [71,72,73]. The mobile nature of phages may then allow the transfer of such advantageous genes to human/animal pathogens or to other non-resident microorganisms of the gut microbiota.

3.3. Identification of the Putative Mobile Resistome of the Bifidobacterium Genus

In order to evaluate the insurgence of AR genes located on or close to mobile elements, such as transposases and bifidoprophages, we investigated the flanking genes of the predicted resistome of the 625 bifidobacterial genomes. These regions may represent mobile genetic hotspots (MGHs) that promote HGT events, thereby transferring antibiotic resistance to other bacteria. We identified 201 putative MGHs distributed in 120 of the 625 Bifidobacterium strains studied. The number of AR genes involved in MGHs were very small compared to the total number of resistome genes (i.e., 13,870), representing less than 1.5% of the total Bifidobacterium resistome. Interestingly, we could not observe a correlation between a specific type of IS element and a class of AR genes.
As already noted in previous studies, 37 of the 41 strains of Bifidobacterium animalis subsp. lactis species contain a tetW gene flanked by a putative conjugative transposon (Figure 3) [74,75,76,77]. The tetW gene encodes a protein belonging to the Guanosine-5’-triphosphate (GTP)-binding elongation factor family that protects ribosomes from the translation inhibition activity of tetracycline [78]. Notably, this MGH is also present in 15 other genomes belonging to members of the Bifidobacterium adolescentis, B. animalis subsp. animalis, B. breve, Bifidobacterium longum spp., Bifidobacterium pseudolongum subsp. pseudolongum, and Bifidobacterium pullorum species. Remarkably, tetW appears to be well conserved among different species (Figure 3), suggesting the involvement of HGT events that could have transferred this tetracycline resistance gene to different bifidobacterial strains.
Interestingly, 67 MGHs involved prophage-like elements, which harbor a gene encoding for an UDP pyrophosphate phosphatase within their sequence (Figure 3), revealing a domain in the amino acid sequence that resembles a bacitracin-resistance protein (BacA) [79,80]. These 67 MGHs were present in members of three different Bifidobacterium species, i.e., B. breve, B. longum spp., and Bifidobacterium pseudocatenulatum, putatively conferring resistance to bacitracin through the phosphorylation of undecaprenol [79,80]. Prophages influence the biodiversity and abundance of bacteria in the human/animal intestinal tract, conferring new capabilities to their host [72]. The acquisition of a prophage-like element may thus confer a fitness advantage [72,73], in this particular case by conferring bacitracin resistance to these Bifidobacterium strains.
Moreover, a gene encoding a 23S rRNA methyltransferase flanked by a transposase was identified in 53 putative MGHs (Figure 3). In a recent study, Martinez et al. demonstrated the existence of a 23S rRNA methylase that confers erythromycin and clindamycin resistance to B. breve CECT7263 [48]. We found these MGHs in 10 different Bifidobacterium species, including B. adolescentis, B. bifidum, B. breve, Bifidobacterium choerinum, Bifidobacterium kashiwanohense, B. longum spp., B. pseudocatenulatum, B. pseudolongum subsp. pseudolongum, and B. pullorum. The highest occurrence of this genetic hotspot was in B. breve strains, where this hotspot was present in 15 out of 88 B. breve genomes analyzed. This methyltransferase is responsible for the enzymatic modification of the nucleotide sequence of the 23S rRNA gene, adding a methyl group, and preventing the linking of macrolide molecules [48]. Notably, the transposases that encompass these MGHs are predicted to be replicative transposons that may cause a rearrangement within bifidobacterial genomes, indicating that these MGHs rarely transfer to other genomes.
Remarkably, B. longum subsp. longum E18, isolated from healthy adult feces samples [81], is the only strain whose chromosome contains a prophage-like element, including a gene predicted to encode a protein with a complete beta-lactamase domain (Figure 3). Furthermore, the genome of strain Bifidobacterium parmae LMG 30295 contains a vanZ homolog flanked by a predicted transposase-encoding gene, belonging to the transposon family IS256. The vanZ gene is predicted to confer low-level resistance to the glycopeptide antibiotic, teicoplanin (Te), which prevents incorporation of D-alanine into peptidoglycan precursors [40]. This hotspot did not include a conjugative transposon, decreasing possible transfer events and bringing possible genomic rearrangements [82,83]. Therefore, more than 50% of putative MGHs identified encompassed transposons that cannot be classified as conjugative transposons, reducing possible HGT events involving AR genes, and corroborating previously published data [16].
The distribution of putative AR genes among analyzed bifidobacteria could be due to selective pressure imposed by intensive antibiotic use in their animal/human hosts, similar to what has been observed for lactic acid bacteria (LAB) [16,54]. These findings underline the safety of this genus and the very low frequency by which these AR genes may transfer to other members of the gut microbiota.

3.4. Assessment of Bacitracin A Resistance of Bifidobacterium spp

In order to validate our in silico predictions, we further investigated the antibiotic resistance of bifidobacterial strains whose genomes were shown to contain a bacA gene located in the sequence of a prophage-like element. Thus, in vitro measurements of MIC breakpoints for the bacitracin A antibiotic were monitored, including three Bifidobacterium strains, i.e., B. breve 1891B, B. longum subsp. longum 35B, and B. longum subsp. infantis ATCC 15697, whose genomes encompass a predicted bacA gene and three additional strains as a control, i.e., B. breve LMG 13208, B. longum subsp. longum LMG 13197, and B. longum subsp. infantis 1888B, whose chromosomes do not include a predicted bacA gene.
As indicated by this in silico analysis, those strains containing the bacA gene in their genomes exhibit a higher resistance level to bacitracin A (ranging from 16-fold to 32-fold) when compared to control strains (Figure S1). In this context, the bacitracin A breakpoints MIC values of B. breve 1891B, B. longum subsp. longum 35B, and B. longum subsp. infantis ATCC 15697 were, respectively, 16, 32, and 64 µg/mL, whereas the MIC values of the members of the control group were 2 µg/mL for B. breve LMG 13208 and 1 µg/mL for B. longum subsp. longum LMG 13197 and B. longum subsp. infantis 1888B (Figure S1). Statistical analyses were performed to corroborate the observed MIC differences, resulting in a significant growth difference between the two groups analyzed (p-value < 0.001) (Figure S1). These results confirmed the in silico-predicted resistance to bacitracin of those bifidobacterial strains possessing MGHs related to the bacA gene. The fact that the in silico analyses matched with the in vitro data highlighted the validity of an in silico resistome prediction [84].

4. Conclusions

Bifidobacteria are dominant members of the human/animal GIT, especially during the early stage of life. It has previously been demonstrated that the presence of this genus in the microbiota is associated with health-promoting effects [4]. In the current study, we reconstructed the resistome and the mobilome of members of the Bifidobacterium genus, evaluating genetic hotspots that could be involved in HGT events. The reconstructed putative resistome revealed that only a limited number of bifidobacterial genes are likely to be involved in putative AR spread. Moreover, the AR genetic arsenal of the Bifidobacterium genus seems to be less complex compared to the resistome of other gram-positive bacteria, such as members of the Lactobacillus genus, or other species included in food supplements and used as probiotics, such as members of the Bacillus genus [54,56,65,66,85,86]. Identified MGHs were restricted to less than 20% of the analyzed strains, of which most were isolated from the human GIT, suggesting the occurrence of AR in members of the human microbiota because of intense antibiotic therapies. Remarkably, the acquisition of phages encompassing AR genes in their sequence could confer ecological advantages, increasing the biological fitness of their host [72,73]. Nevertheless, the vast majority of identified MGHs in the Bifidobacterium genus are unlikely to be transferred to other microorganisms, due to the transposition mechanisms of the identified IS elements flanking putative AR genes. Moreover, in vitro bacitracin A antibiotic resistance tests based on bifidobacterial strains containing bacA located in an MGH confirmed our in silico prediction. Finally, these findings underpin the safety of the Bifidobacterium genus compared to other taxa, such as Escherichia coli and members of the Gammaproteobacteria class, which were shown to contribute to a high antibiotic resistance load in the human microbiota [87].

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2076-2607/7/12/638/s1, Figure S1: Bacitracin A antibiotic breakpoint values, Table S1: The resistome and the mobilome of the Bifidobacterium genus, Table S2: Distribution of IS family in Bifidobacterium genus.

Author Contributions

For research articles with several authors, a short paragraph specifying their individual contributions must be provided. The following statements should be used “conceptualization, F.T. and M.V.; methodology, W.M. and G.A.L.; software, W.M. and G.A.L.; validation, F.T., M.V. and D.v.S.; formal analysis, W.M. and G.A.L.; investigation, W.M. and G.A.L.; resources, F.T., M.V. and D.v.S.; data curation, W.M. and G.A.L.; writing—original draft preparation, W.M. and G.A.L.; writing—review and editing, F.T., M.V. and D.v.S.; visualization, F.T. and M.V.; supervision, F.T. and M.V.; project administration, F.T. and M.V.; funding acquisition, F.T. and M.V., please turn to the CRediT taxonomy for the term explanation. Authorship must be limited to those who have contributed substantially to the work reported.

Funding

This work was primarily funded by the EU Joint Programming Initiative – A Healthy Diet for a Healthy Life (JPI HDHL, http://www.healthydietforhealthylife.eu/) to DVS (in conjunction with Science Foundation Ireland [SFI], Grant number 15/JP-HDHL/3280) and to MV (in conjunction with MIUR, Italy). DVS is member of APC microbiome Ireland which is funded by SFI through the Irish Government’s National Development Plan (Grant Numbers SFI/12/RC/2273-P1 and SFI/12/RC/2273-P2). The study is supported by Fondazione Cariparma, under TeachInParma Project (DVS).

Acknowledgments

This work was primarily funded by the EU Joint Programming Initiative—A Healthy Diet for a Healthy Life (JPI HDHL, http://www.healthydietforhealthylife.eu/) to DvS (in conjunction with Science Foundation Ireland [SFI], Grant number 15/JP-HDHL/3280) and to MV (in conjunction with MIUR, Italy). D.v.S. is member of APC microbiome Ireland which is funded by SFI through the Irish Government’s National Development Plan (Grant Numbers SFI/12/RC/2273-P1 and SFI/12/RC/2273-P2). The study is supported by Fondazione Cariparma, under TeachInParma Project (DV). GA is supported by Fondazione Cariparma, Parma, Italy. We furthermore thank GenProbio srl for financial support of the Laboratory of Probiogenomics. Part of this research is conducted using the High Performance Computing (HPC) facility of the University of Parma.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Milani, C.; Lugli, G.A.; Duranti, S.; Turroni, F.; Bottacini, F.; Mangifesta, M.; Sanchez, B.; Viappiani, A.; Mancabelli, L.; Taminiau, B.; et al. Genomic encyclopedia of type strains of the genus Bifidobacterium. Appl. Environ. Microbiol. 2014, 80, 6290–6302. [Google Scholar] [CrossRef] [PubMed]
  2. Ventura, M.; van Sinderen, D.; Fitzgerald, G.F.; Zink, R. Insights into the taxonomy, genetics and physiology of bifidobacteria. Antonie Van Leeuwenhoek 2004, 86, 205–223. [Google Scholar] [CrossRef] [PubMed]
  3. Ventura, M.; Canchaya, C.; Tauch, A.; Chandra, G.; Fitzgerald, G.F.; Chater, K.F.; van Sinderen, D. Genomics of Actinobacteria: Tracing the evolutionary history of an ancient phylum. Microbiol. Mol. Biol. Rev. 2007, 71, 495–548. [Google Scholar] [CrossRef] [PubMed]
  4. Ventura, M.; Turroni, F.; Lugli, G.A.; van Sinderen, D. Bifidobacteria and humans: Our special friends, from ecological to genomics perspectives. J. Sci. Food Agric. 2014, 94, 163–168. [Google Scholar] [CrossRef]
  5. Turroni, F.; Peano, C.; Pass, D.A.; Foroni, E.; Severgnini, M.; Claesson, M.J.; Kerr, C.; Hourihane, J.; Murray, D.; Fuligni, F.; et al. Diversity of bifidobacteria within the infant gut microbiota. PLoS ONE 2012, 7, e36957. [Google Scholar] [CrossRef]
  6. Verani, J.R.; McGee, L.; Schrag, S.J. Division of Bacterial Diseases, National Center for Immunization and Respiratory Diseases Prevention of perinatal group B streptococcal disease—revised guidelines from CDC, 2010. MMWR Recomm. Rep. 2010, 59, 1–36. [Google Scholar]
  7. Cox, L.M.; Blaser, M.J. Antibiotics in early life and obesity. Nat. Rev. Endocrinol. 2015, 11, 182–190. [Google Scholar] [CrossRef]
  8. Gibson, M.K.; Crofts, T.S.; Dantas, G. Antibiotics and the developing infant gut microbiota and resistome. Curr. Opin. Microbiol. 2015, 27, 51–56. [Google Scholar] [CrossRef]
  9. Zou, Z.H.; Liu, D.; Li, H.D.; Zhu, D.P.; He, Y.; Hou, T.; Yu, J.L. Prenatal and postnatal antibiotic exposure influences the gut microbiota of preterm infants in neonatal intensive care units. Ann. Clin. Microbiol. Antimicrob. 2018, 17, 9. [Google Scholar] [CrossRef]
  10. Marshall, B.M.; Levy, S.B. Food animals and antimicrobials: Impacts on human health. Clin. Microbiol. Rev. 2011, 24, 718–733. [Google Scholar] [CrossRef]
  11. Ouwehand, A.C.; Forssten, S.; Hibberd, A.A.; Lyra, A.; Stahl, B. Probiotic approach to prevent antibiotic resistance. Ann. Med. 2016, 48, 246–255. [Google Scholar] [CrossRef]
  12. Arboleya, S.; Sanchez, B.; Milani, C.; Duranti, S.; Solis, G.; Fernandez, N.; de los Reyes-Gavilan, C.G.; Ventura, M.; Margolles, A.; Gueimonde, M. Intestinal microbiota development in preterm neonates and effect of perinatal antibiotics. J. Pediatr. 2015, 166, 538–544. [Google Scholar] [CrossRef] [PubMed]
  13. Fouhy, F.; Guinane, C.M.; Hussey, S.; Wall, R.; Ryan, C.A.; Dempsey, E.M.; Murphy, B.; Ross, R.P.; Fitzgerald, G.F.; Stanton, C.; et al. High-throughput sequencing reveals the incomplete, short-term recovery of infant gut microbiota following parenteral antibiotic treatment with ampicillin and gentamicin. Antimicrob. Agents Chemother. 2012, 56, 5811–5820. [Google Scholar] [CrossRef]
  14. D’Costa, V.M.; McGrann, K.M.; Hughes, D.W.; Wright, G.D. Sampling the antibiotic resistome. Science 2006, 311, 374–377. [Google Scholar] [CrossRef] [PubMed]
  15. Wright, G.D. The antibiotic resistome: The nexus of chemical and genetic diversity. Nat. Rev. Microbiol. 2007, 5, 175–186. [Google Scholar] [CrossRef]
  16. Duranti, S.; Lugli, G.A.; Mancabelli, L.; Turroni, F.; Milani, C.; Mangifesta, M.; Ferrario, C.; Anzalone, R.; Viappiani, A.; van Sinderen, D.; et al. Prevalence of Antibiotic Resistance Genes among Human Gut-Derived Bifidobacteria. Appl. Environ. Microbiol. 2017, 83. [Google Scholar] [CrossRef] [PubMed]
  17. Wiedenbeck, J.; Cohan, F.M. Origins of bacterial diversity through horizontal genetic transfer and adaptation to new ecological niches. FEMS Microbiol. Rev. 2011, 35, 957–976. [Google Scholar] [CrossRef]
  18. Brito, I.L.; Yilmaz, S.; Huang, K.; Xu, L.; Jupiter, S.D.; Jenkins, A.P.; Naisilisili, W.; Tamminen, M.; Smillie, C.S.; Wortman, J.R.; et al. Mobile genes in the human microbiome are structured from global to individual scales. Nature 2016, 535, 435–439. [Google Scholar] [CrossRef]
  19. Hagbo, M.; Ravi, A.; Angell, I.L.; Sunde, M.; Ludvigsen, J.; Diep, D.B.; Foley, S.L.; Vento, M.; Collado, M.C.; Perez-Martinez, G.; et al. Experimental support for multidrug resistance transfer potential in the preterm infant gut microbiota. Pediatr. Res. 2019, 231. [Google Scholar] [CrossRef]
  20. Martinez, J.L.; Coque, T.M.; Baquero, F. What is a resistance gene? Ranking risk in resistomes. Nat. Rev. Microbiol. 2015, 13, 116–123. [Google Scholar] [CrossRef]
  21. Siefert, J.L. Defining the mobilome. Methods Mol. Biol. 2009, 532, 13–27. [Google Scholar] [CrossRef] [PubMed]
  22. Smets, B.F.; Barkay, T. Horizontal gene transfer: Perspectives at a crossroads of scientific disciplines. Nat. Rev. Microbiol. 2005, 3, 675–678. [Google Scholar] [CrossRef] [PubMed]
  23. Guglielmetti, S.; Mayo, B.; Alvarez-Martin, P. Mobilome and genetic modification of bifidobacteria. Benef. Microbes. 2013, 4, 143–166. [Google Scholar] [CrossRef] [PubMed]
  24. Lugli, G.A.; Milani, C.; Turroni, F.; Tremblay, D.; Ferrario, C.; Mancabelli, L.; Duranti, S.; Ward, D.V.; Ossiprandi, M.C.; Moineau, S.; et al. Prophages of the genus Bifidobacterium as modulating agents of the infant gut microbiota. Environ. Microbiol. 2016, 18, 2196–2213. [Google Scholar] [CrossRef]
  25. Lugli, G.A.; Milani, C.; Duranti, S.; Mancabelli, L.; Mangifesta, M.; Turroni, F.; Viappiani, A.; van Sinderen, D.; Ventura, M. Tracking the Taxonomy of the Genus Bifidobacterium Based on a Phylogenomic Approach. Appl. Environ. Microbiol. 2018, 84. [Google Scholar] [CrossRef]
  26. Lugli, G.A.; Milani, C.; Turroni, F.; Duranti, S.; Ferrario, C.; Viappiani, A.; Mancabelli, L.; Mangifesta, M.; Taminiau, B.; Delcenserie, V.; et al. Investigation of the evolutionary development of the genus Bifidobacterium by comparative genomics. Appl. Environ. Microbiol. 2014, 80, 6383–6394. [Google Scholar] [CrossRef]
  27. Lugli, G.A.; Milani, C.; Turroni, F.; Duranti, S.; Mancabelli, L.; Mangifesta, M.; Ferrario, C.; Modesto, M.; Mattarelli, P.; Jiri, K.; et al. Comparative genomic and phylogenomic analyses of the Bifidobacteriaceae family. BMC Genom. 2017, 18, 568. [Google Scholar] [CrossRef]
  28. Richter, M.; Rossello-Mora, R. Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. USA 2009, 106, 19126–19131. [Google Scholar] [CrossRef]
  29. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  30. Siguier, P.; Perochon, J.; Lestrade, L.; Mahillon, J.; Chandler, M. ISfinder: The reference centre for bacterial insertion sequences. Nucleic Acids Res. 2006, 34, D32–D36. [Google Scholar] [CrossRef]
  31. Lakin, S.M.; Dean, C.; Noyes, N.R.; Dettenwanger, A.; Ross, A.S.; Doster, E.; Rovira, P.; Abdo, Z.; Jones, K.L.; Ruiz, J.; et al. MEGARes: An antimicrobial resistance database for high throughput sequencing. Nucleic Acids Res. 2017, 45, D574–D580. [Google Scholar] [CrossRef] [PubMed]
  32. Gupta, S.K.; Padmanabhan, B.R.; Diene, S.M.; Lopez-Rojas, R.; Kempf, M.; Landraud, L.; Rolain, J.M. ARG-ANNOT, a new bioinformatic tool to discover antibiotic resistance genes in bacterial genomes. Antimicrob. Agents Chemother. 2014, 58, 212–220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Zankari, E.; Hasman, H.; Cosentino, S.; Vestergaard, M.; Rasmussen, S.; Lund, O.; Aarestrup, F.M.; Larsen, M.V. Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother. 2012, 67, 2640–2644. [Google Scholar] [CrossRef] [PubMed]
  34. McArthur, A.G.; Waglechner, N.; Nizam, F.; Yan, A.; Azad, M.A.; Baylay, A.J.; Bhullar, K.; Canova, M.J.; De Pascale, G.; Ejim, L.; et al. The comprehensive antibiotic resistance database. Antimicrob. Agents Chemother. 2013, 57, 3348–3357. [Google Scholar] [CrossRef] [Green Version]
  35. Bush, K.; Jacoby, G.A. Updated functional classification of beta-lactamases. Antimicrob. Agents Chemother. 2010, 54, 969–976. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Kielbasa, S.M.; Wan, R.; Sato, K.; Horton, P.; Frith, M.C. Adaptive seeds tame genomic sequence comparison. Genome Res. 2011, 21, 487–493. [Google Scholar] [CrossRef] [Green Version]
  37. Katoh, K.; Misawa, K.; Kuma, K.; Miyata, T. MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30, 3059–3066. [Google Scholar] [CrossRef] [Green Version]
  38. Chenna, R.; Sugawara, H.; Koike, T.; Lopez, R.; Gibson, T.J.; Higgins, D.G.; Thompson, J.D. Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003, 31, 3497–3500. [Google Scholar] [CrossRef] [Green Version]
  39. Authority, E.F.S. Guidance on the assessment of bacterial susceptibility to antimicrobials of human and veterinary importance. EFSA J. 2012, 10, 2740. [Google Scholar]
  40. Arthur, M.; Depardieu, F.; Molinas, C.; Reynolds, P.; Courvalin, P. The vanZ gene of Tn1546 from Enterococcus faecium BM4147 confers resistance to teicoplanin. Gene 1995, 154, 87–92. [Google Scholar] [CrossRef]
  41. Bugg, T.D.; Wright, G.D.; Dutka-Malen, S.; Arthur, M.; Courvalin, P.; Walsh, C.T. Molecular basis for vancomycin resistance in Enterococcus faecium BM4147: Biosynthesis of a depsipeptide peptidoglycan precursor by vancomycin resistance proteins VanH and VanA. Biochemistry 1991, 30, 10408–10415. [Google Scholar] [CrossRef] [PubMed]
  42. Evers, S.; Courvalin, P. Regulation of VanB-type vancomycin resistance gene expression by the VanS(B)-VanR (B) two-component regulatory system in Enterococcus faecalis V583. J. Bacteriol. 1996, 178, 1302–1309. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Giovanetti, E.; Brenciani, A.; Lupidi, R.; Roberts, M.C.; Varaldo, P.E. Presence of the tet(O) gene in erythromycin- and tetracycline-resistant strains of Streptococcus pyogenes and linkage with either the mef(A) or the erm(A) gene. Antimicrob. Agents Chemother. 2003, 47, 2844–2849. [Google Scholar] [CrossRef] [Green Version]
  44. Roberts, M.C. Tetracycline resistance determinants: Mechanisms of action, regulation of expression, genetic mobility, and distribution. FEMS Microbiol. Rev. 1996, 19, 1–24. [Google Scholar] [CrossRef] [PubMed]
  45. Hedayatianfard, K.; Akhlaghi, M.; Sharifiyazdi, H. Detection of tetracycline resistance genes in bacteria isolated from fish farms using polymerase chain reaction. Vet. Res. Forum Int. Q. J. 2014, 5, 269–275. [Google Scholar]
  46. Wang, N.; Hang, X.; Zhang, M.; Liu, X.; Yang, H. Analysis of newly detected tetracycline resistance genes and their flanking sequences in human intestinal bifidobacteria. Sci. Rep. 2017, 7, 6267. [Google Scholar] [CrossRef]
  47. Zou, Y.; Xue, W.; Luo, G.; Deng, Z.; Qin, P.; Guo, R.; Sun, H.; Xia, Y.; Liang, S.; Dai, Y.; et al. 1,520 reference genomes from cultivated human gut bacteria enable functional microbiome analyses. Nat. Biotechnol. 2019, 37, 179–185. [Google Scholar] [CrossRef]
  48. Martinez, N.; Luque, R.; Milani, C.; Ventura, M.; Banuelos, O.; Margolles, A. A Gene Homologous to rRNA Methylase Genes Confers Erythromycin and Clindamycin Resistance in Bifidobacterium breve. Appl. Environ. Microbiol. 2018, 84. [Google Scholar] [CrossRef] [Green Version]
  49. Skold, O. Resistance to trimethoprim and sulfonamides. Vet. Res. 2001, 32, 261–273. [Google Scholar] [CrossRef] [Green Version]
  50. Phuong Hoa, P.T.; Nonaka, L.; Hung Viet, P.; Suzuki, S. Detection of the sul1, sul2, and sul3 genes in sulfonamide-resistant bacteria from wastewater and shrimp ponds of north Vietnam. Sci. Total Environ. 2008, 405, 377–384. [Google Scholar] [CrossRef]
  51. Vetting, M.W.; Hegde, S.S.; Fajardo, J.E.; Fiser, A.; Roderick, S.L.; Takiff, H.E.; Blanchard, J.S. Pentapeptide repeat proteins. Biochemistry 2006, 45, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Merens, A.; Matrat, S.; Aubry, A.; Lascols, C.; Jarlier, V.; Soussy, C.J.; Cavallo, J.D.; Cambau, E. The pentapeptide repeat proteins MfpAMt and QnrB4 exhibit opposite effects on DNA gyrase catalytic reactions and on the ternary gyrase-DNA-quinolone complex. J. Bacteriol. 2009, 191, 1587–1594. [Google Scholar] [CrossRef] [Green Version]
  53. Park, K.S.; Lee, J.H.; Jeong, D.U.; Lee, J.J.; Wu, X.; Jeong, B.C.; Kang, C.M.; Lee, S.H. Determination of pentapeptide repeat units in Qnr proteins by the structure-based alignment approach. Antimicrob. Agents Chemother. 2011, 55, 4475–4478. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Gueimonde, M.; Sanchez, B.; de Los Reyes-Gavilán, C.G.; Margolles, A. Antibiotic resistance in probiotic bacteria. Front. Microbiol. 2013, 4, 202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Delcour, J.; Ferain, T.; Deghorain, M.; Palumbo, E.; Hols, P. The biosynthesis and functionality of the cell-wall of lactic acid bacteria. Antonie Van Leeuwenhoek 1999, 76, 159–184. [Google Scholar] [CrossRef]
  56. Florez, A.B.; Ladero, V.; Alvarez-Martin, P.; Ammor, M.S.; Alvarez, M.A.; Mayo, B. Acquired macrolide resistance in the human intestinal strain Lactobacillus rhamnosus E41 associated with a transition mutation in 23S rRNA genes. Int. J. Antimicrob. Agents 2007, 30, 341–344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Hummel, A.S.; Hertel, C.; Holzapfel, W.H.; Franz, C.M. Antibiotic resistances of starter and probiotic strains of lactic acid bacteria. Appl. Environ. Microbiol. 2007, 73, 730–739. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Ammor, M.S.; Gueimonde, M.; Danielsen, M.; Zagorec, M.; van Hoek, A.H.; de Los Reyes-Gavilan, C.G.; Mayo, B.; Margolles, A. Two different tetracycline resistance mechanisms, plasmid-carried tet(L) and chromosomally located transposon-associated tet(M), coexist in Lactobacillus sakei Rits 9. Appl. Environ. Microbiol. 2008, 74, 1394–1401. [Google Scholar] [CrossRef] [Green Version]
  59. Rojo-Bezares, B.; Saenz, Y.; Poeta, P.; Zarazaga, M.; Ruiz-Larrea, F.; Torres, C. Assessment of antibiotic susceptibility within lactic acid bacteria strains isolated from wine. Int. J. Food Microbiol. 2006, 111, 234–240. [Google Scholar] [CrossRef]
  60. Johnning, A.; Karami, N.; Tang Hallback, E.; Muller, V.; Nyberg, L.; Buongermino Pereira, M.; Stewart, C.; Ambjornsson, T.; Westerlund, F.; Adlerberth, I.; et al. The resistomes of six carbapenem-resistant pathogens—a critical genotype-phenotype analysis. Microb. Genom. 2018, 4. [Google Scholar] [CrossRef]
  61. Dagher, C.; Salloum, T.; Alousi, S.; Arabaghian, H.; Araj, G.F.; Tokajian, S. Molecular characterization of Carbapenem resistant Escherichia coli recovered from a tertiary hospital in Lebanon. PLoS ONE 2018, 13, e0203323. [Google Scholar] [CrossRef] [PubMed]
  62. Khan, I.; Yasir, M.; Farman, M.; Kumosani, T.; AlBasri, S.F.; Bajouh, O.S.; Azhar, E.I. Evaluation of gut bacterial community composition and antimicrobial resistome in pregnant and non-pregnant women from Saudi population. Infect. Drug Resist. 2019, 12, 1749–1761. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Hong, H.A.; Duc le, H.; Cutting, S.M. The use of bacterial spore formers as probiotics. Fems Microbiol. Rev. 2005, 29, 813–835. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Monod, M.; Denoya, C.; Dubnau, D. Sequence and properties of pIM13, a macrolide-lincosamide-streptogramin B resistance plasmid from Bacillus subtilis. J. Bacteriol. 1986, 167, 138–147. [Google Scholar] [CrossRef] [Green Version]
  65. Phelan, R.W.; Clarke, C.; Morrissey, J.P.; Dobson, A.D.; O’Gara, F.; Barbosa, T.M. Tetracycline resistance-encoding plasmid from Bacillus sp. strain #24, isolated from the marine sponge Haliclona simulans. Appl. Environ. Microbiol. 2011, 77, 327–329. [Google Scholar] [CrossRef] [Green Version]
  66. Dai, L.; Wu, C.M.; Wang, M.G.; Wang, Y.; Wang, Y.; Huang, S.Y.; Xia, L.N.; Li, B.B.; Shen, J.Z. First report of the multidrug resistance gene cfr and the phenicol resistance gene fexA in a Bacillus strain from swine feces. Antimicrob. Agents Chemother. 2010, 54, 3953–3955. [Google Scholar] [CrossRef] [Green Version]
  67. Mahony, J.; Lugli, G.A.; van Sinderen, D.; Ventura, M. Impact of gut-associated bifidobacteria and their phages on health: Two sides of the same coin? Appl. Microbiol. Biotechnol. 2018, 102, 2091–2099. [Google Scholar] [CrossRef]
  68. Ventura, M.; Canchaya, C.; Fitzgerald, G.F.; Gupta, R.S.; van Sinderen, D. Genomics as a means to understand bacterial phylogeny and ecological adaptation: The case of bifidobacteria. Antonie Van Leeuwenhoek 2007, 91, 351–372. [Google Scholar] [CrossRef]
  69. Siguier, P.; Varani, A.; Perochon, J.; Chandler, M. Exploring bacterial insertion sequences with ISfinder: Objectives, uses, and future developments. Methods Mol. Biol. 2012, 859, 91–103. [Google Scholar] [CrossRef] [PubMed]
  70. Mavrich, T.N.; Casey, E.; Oliveira, J.; Bottacini, F.; James, K.; Franz, C.; Lugli, G.A.; Neve, H.; Ventura, M.; Hatfull, G.F.; et al. Characterization and induction of prophages in human gut-associated Bifidobacterium hosts. Sci. Rep. 2018, 8, 12772. [Google Scholar] [CrossRef]
  71. Bottacini, F.; Medini, D.; Pavesi, A.; Turroni, F.; Foroni, E.; Riley, D.; Giubellini, V.; Tettelin, H.; van Sinderen, D.; Ventura, M. Comparative genomics of the genus Bifidobacterium. Microbiology 2010, 156, 3243–3254. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Menouni, R.; Hutinet, G.; Petit, M.A.; Ansaldi, M. Bacterial genome remodeling through bacteriophage recombination. Fems Microbiol. Lett. 2015, 362, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Bondy-Denomy, J.; Davidson, A.R. When a virus is not a parasite: The beneficial effects of prophages on bacterial fitness. J. Microbiol. 2014, 52, 235–242. [Google Scholar] [CrossRef] [PubMed]
  74. Milani, C.; Duranti, S.; Lugli, G.A.; Bottacini, F.; Strati, F.; Arioli, S.; Foroni, E.; Turroni, F.; van Sinderen, D.; Ventura, M. Comparative genomics of Bifidobacterium animalis subsp. lactis reveals a strict monophyletic bifidobacterial taxon. Appl. Environ. Microbiol. 2013, 79, 4304–4315. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Ammor, M.S.; Florez, A.B.; Alvarez-Martin, P.; Margolles, A.; Mayo, B. Analysis of tetracycline resistance tet(W) genes and their flanking sequences in intestinal Bifidobacterium species. J. Antimicrob. Chemother. 2008, 62, 688–693. [Google Scholar] [CrossRef] [Green Version]
  76. Florez, A.B.; Ammor, M.S.; Alvarez-Martin, P.; Margolles, A.; Mayo, B. Molecular analysis of tet(W) gene-mediated tetracycline resistance in dominant intestinal Bifidobacterium species from healthy humans. Appl. Environ. Microbiol. 2006, 72, 7377–7379. [Google Scholar] [CrossRef] [Green Version]
  77. Gueimonde, M.; Florez, A.B.; van Hoek, A.H.; Stuer-Lauridsen, B.; Stroman, P.; de los Reyes-Gavilan, C.G.; Margolles, A. Genetic basis of tetracycline resistance in Bifidobacterium animalis subsp. lactis. Appl. Environ. Microbiol. 2010, 76, 3364–3369. [Google Scholar] [CrossRef] [Green Version]
  78. Scott, K.P.; Barbosa, T.M.; Forbes, K.J.; Flint, H.J. High-frequency transfer of a naturally occurring chromosomal tetracycline resistance element in the ruminal anaerobe Butyrivibrio fibrisolvens. Appl. Environ. Microbiol. 1997, 63, 3405–3411. [Google Scholar]
  79. Cain, B.D.; Norton, P.J.; Eubanks, W.; Nick, H.S.; Allen, C.M. Amplification of the bacA gene confers bacitracin resistance to Escherichia coli. J. Bacteriol. 1993, 175, 3784–3789. [Google Scholar] [CrossRef] [Green Version]
  80. El Ghachi, M.; Bouhss, A.; Blanot, D.; Mengin-Lecreulx, D. The bacA gene of Escherichia coli encodes an undecaprenyl pyrophosphate phosphatase activity. J. Biol. Chem. 2004, 279, 30106–30113. [Google Scholar] [CrossRef] [Green Version]
  81. Zhurina, D.; Dudnik, A.; Waidmann, M.S.; Grimm, V.; Westermann, C.; Breitinger, K.J.; Yuan, J.; van Sinderen, D.; Riedel, C.U. High-Quality Draft Genome Sequence of Bifidobacterium longum E18, Isolated from a Healthy Adult. Genome Announc. 2013, 1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Wright, L.D.; Grossman, A.D. Autonomous Replication of the Conjugative Transposon Tn916. J. Bacteriol. 2016, 198, 3355–3366. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Cury, J.; Touchon, M.; Rocha, E.P.C. Integrative and conjugative elements and their hosts: Composition, distribution and organization. Nucleic Acids Res. 2017, 45, 8943–8956. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Moradigaravand, D.; Palm, M.; Farewell, A.; Mustonen, V.; Warringer, J.; Parts, L. Prediction of antibiotic resistance in Escherichia coli from large-scale pan-genome data. PLoS Comput. Biol. 2018, 14, e1006258. [Google Scholar] [CrossRef] [Green Version]
  85. Mayrhofer, S.; van Hoek, A.H.; Mair, C.; Huys, G.; Aarts, H.J.; Kneifel, W.; Domig, K.J. Antibiotic susceptibility of members of the Lactobacillus acidophilus group using broth microdilution and molecular identification of their resistance determinants. Int. J. Food Microbiol. 2010, 144, 81–87. [Google Scholar] [CrossRef]
  86. Lin, C.F.; Fung, Z.F.; Wu, C.L.; Chung, T.C. Molecular characterization of a plasmid-borne (pTC82) chloramphenicol resistance determinant (cat-TC) from Lactobacillus reuteri G4. Plasmid 1996, 36, 116–124. [Google Scholar] [CrossRef]
  87. Parnanen, K.; Karkman, A.; Hultman, J.; Lyra, C.; Bengtsson-Palme, J.; Larsson, D.G.J.; Rautava, S.; Isolauri, E.; Salminen, S.; Kumar, H.; et al. Maternal gut and breast milk microbiota affect infant gut antibiotic resistome and mobile genetic elements. Nat. Commun. 2018, 9, 3891. [Google Scholar] [CrossRef] [Green Version]
Microorganisms 07 00638 g001 550
Figure 1. Predicted resistome of the Bifidobacterium genus. Abundance of different predicted antibiotic resistance gene classes identified among the 625 analyzed Bifidobacterium genomes.
Figure 1. Predicted resistome of the Bifidobacterium genus. Abundance of different predicted antibiotic resistance gene classes identified among the 625 analyzed Bifidobacterium genomes.
Microorganisms 07 00638 g001
Microorganisms 07 00638 g002 550
Figure 2. Phylogenetic tree of identified bifidoprophages. Genomic alignment-based clustering of 598 prophages identified within bifidobacterial strain genomes. Each colored dot represents the Bifidobacterium host species origin of a given bifidoprophage. The dot size refers to the number of prophage-like sequences identified within the same branch tree. The four different clusters are highlighted with different colors.
Figure 2. Phylogenetic tree of identified bifidoprophages. Genomic alignment-based clustering of 598 prophages identified within bifidobacterial strain genomes. Each colored dot represents the Bifidobacterium host species origin of a given bifidoprophage. The dot size refers to the number of prophage-like sequences identified within the same branch tree. The four different clusters are highlighted with different colors.
Microorganisms 07 00638 g002
Microorganisms 07 00638 g003 550
Figure 3. Mobile genetic hotspots identified in the Bifidobacterium genus. Bifidobacterial genomic regions containing putative Mobile Genetic Hotspots (MGHs). Different species and gene names are reported next to each genomic region. Panels (a–c) show the genomic regions conserved among different Bifidobacterium species. Panels (d,e) display unique mobile genetic hotspots identified in B. longum subsp. longum E18 and B. parmae LMG 30,295 strains. Each arrow indicates a gene and the different colors indicate the function of the gene product.
Figure 3. Mobile genetic hotspots identified in the Bifidobacterium genus. Bifidobacterial genomic regions containing putative Mobile Genetic Hotspots (MGHs). Different species and gene names are reported next to each genomic region. Panels (a–c) show the genomic regions conserved among different Bifidobacterium species. Panels (d,e) display unique mobile genetic hotspots identified in B. longum subsp. longum E18 and B. parmae LMG 30,295 strains. Each arrow indicates a gene and the different colors indicate the function of the gene product.
Microorganisms 07 00638 g003

Share and Cite

MDPI and ACS Style

Mancino, W.; Lugli, G.A.; van Sinderen, D.; Ventura, M.; Turroni, F. Mobilome and Resistome Reconstruction from Genomes Belonging to Members of the Bifidobacterium Genus. Microorganisms 2019, 7, 638. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms7120638

AMA Style

Mancino W, Lugli GA, van Sinderen D, Ventura M, Turroni F. Mobilome and Resistome Reconstruction from Genomes Belonging to Members of the Bifidobacterium Genus. Microorganisms. 2019; 7(12):638. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms7120638

Chicago/Turabian Style

Mancino, Walter, Gabriele Andrea Lugli, Douwe van Sinderen, Marco Ventura, and Francesca Turroni. 2019. "Mobilome and Resistome Reconstruction from Genomes Belonging to Members of the Bifidobacterium Genus" Microorganisms 7, no. 12: 638. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms7120638

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop