Genome Analysis for AR Genes
Based on the phenotype and genotype results (;
Supplementary Figure S1), six strains were subjected to genome sequencing:
L. crispatus VA50-4AN (resistant to kanamycin, ampicillin, and trimethoprim),
L. jensenii VA04-2AN (resistant to trimethoprim),
L. salivarius VA40-10 (resistant to gentamicin, kanamycin, streptomycin, neomycin, and vancomycin),
L. paracasei VA02-1AN (resistant to chloramphenicol and vancomycin),
L. reuteri VA24-5 (resistant to vancomycin and trimethoprim), and
B. bifidum VA07-1AN (resistant to streptomycin).
Supplementary Table S2 shows the general features of their genomes. Their size was, in all cases, around 2.2 Mbp but the number of contigs obtained after assembly ranged from 17 to 300.
Supplementary Table S3 summarizes some of the key genetic features of the genomes of the sequenced strains. Genes coding for penicillin binding proteins (PBP) were found in all the genomes, although with different numbers and types for the distinct species. Mutations in PBPs encoding-genes known to confer AR were not identified. One gene coding for a D-alanine-D-alanine ligase (Ddl) was detected in each of the strains. In several LAB species, the presence of phenylalanine at the enzyme active site in Ddl is correlated with intrinsic resistance to vancomycin [
18]. In addition, in each of the strains, 9-32 genes were classified by the RAST server as belonging to the category “Virulence, Disease, and Defence”, subcategory “Resistance to Antibiotic and Toxic Compounds”. The majority of these genes encoded components dedicated to homeostasis or resistance to heavy metals, such as copper, mercury, and the cobalt-zinc-cadmium triad. Genes encoding elongation factors, efflux pumps, DNA gyrases, and topoisomerases were also included by RAST in this subcategory.
By comparing the genome sequences against the databases CARD, ResFinder, and ARG-ANNOT, no genes known to be involved in AR in
L. jensenii VA04-2AN (resistant to trimethoprim),
L. paracasei VA02-1AN (resistant to chloramphenicol and vancomycin), and
L. reuteri VA24-5 (resistant to trimethoprim and vancomycin) were detected. The only positive correlation between phenotype and genotype was the presence of a conserved phenylalanine (F) residue in the active site of the Ddl ligase, corresponding to amino acid 261 of the
Leuconostoc mesenteroides enzyme [
18], in the deduced sequence of all vancomycin-resistant (Vm
r) strains, while the susceptible (Vm
s) strains were characterized by the presence of a tyrosine (Y) residue at this position ().
Figure 2. Alignment of amino acid sequences around the active site of D-Ala-D-Ala ligases of the five Lactobacillus spp. strains sequenced. Strains with phenylalanine (F) at the enzyme active site (green) show a vancomycin-resistant phenotype, while those having a tyrosine (Y) (pale blue) display a vancomycin-susceptible phenotype.
Genome analysis of L. crispatus VA50-4AN, L. salivarius VA40-10, and B. bifidum VA07-1AN identified no genes known to be involved in aminoglycoside resistance. Therefore, mutations in key genes, such as those coding for the ribosomal S12 protein and others acting on the 16S rRNA molecule, were therefore sought by comparing the DNA and deduced protein sequences from our strains with those in databases. No amino acid differences were observed in the sequences of the ribosomal protein S12 for L. crispatus VA50-4AN and L. salivarius VA40-10 from those belonging to susceptible strains of the same species. Further, alignment of the deduced amino acid sequences for the 16S rRNA guanine(527)-N(7)-methyltranferase (RsmG) proteins of the sequenced strains, showed heterogeneity at several positions between themselves and with respect to sequences in databases. In particular, the RsmG sequence of L. crispatus VA50-4AN showed one amino acid change at position 38 (N→H), while that of B. bifidum VA07-1AN showed three amino acid changes at positions 105 (E→A), 150 (G→D), and 206 (R→G), and that of L. salivarius VA40-10 showed six exclusive amino acid changes at positions 12 (G→E), 67 (D→N), 186 (N→D), 199 (I→V), 208 (Q→K), and 209 (V→I). However, by comparing RsmG sequences from resistant and susceptible strains, none of the changes considered could be associated with streptomycin resistance.
As expected, the genome analysis confirmed the presence of
tet(W) in
B. bifidum VA07-1AN; this gene was also unequivocally identified by searches in the three AR databases used. The
tet(W) gene in
B. bifidum VA07-1AN was located in a contig of 76,748 bp. shows the genetic organization of the 40-kbp left extreme of the contig that included the
tet(W) gene. The
tet(W) sequence of VA07-1AN (1922 bp) was almost identical to that described for
Bifidobacterium longum LTBL16 (CP034089.1). Similar
tet(W) sequences have also been found in the chromosome of strains belonging to other species such as
B. bifidum L22 (NG_048301.1),
Lachnospiraceae bacterium KGMB03038 (CP041667.1), and
Ruminococcus sp. JE7A12 (CP039381.1); and in plasmids, such as pTZC1 from
Cutibacterium acnes TP-CU389 (LC473083.1). Compared to
tet(W) in
B. longum LTBL16, the
tet(W) in
B. bifidum VA07-1AN contained an insertion of two extra guanine residues (GG) after nucleotide 731 in the ORF resulting in a frameshift, which produced only a short peptide—289 amino acids long compared to 639 residues for the functional Tet(W). This likely explains the susceptibility of VA07-1AN to tetracycline. The
tet(W) gene was flanked by ORFs coding for proteins showing the greatest homology to others from
B. longum in the upstream region, and proteins typical of
B. bifidum in the downstream region (
Supplementary Table S4).
Figure 3. Diagram showing the genetic organization of ORFs in the contig harboring the tet(W) gene of Bifidobacterium bifidum VA07-1AN. Color key: purple, tet(W) gene (the position of the GG insertion disrupting the ORF is indicated); yellow, conjugation-associated gene; pale blue, gene encoding a transcription regulator; white, genes involved in other processes. The broken line symbol indicates the contig extends beyond this point.
The CARD database further identified in the genome of
B. bifidum VA07-1AN a single nucleotide polymorphism (SNP) point mutation in the
rpsL gene (encoding the ribosomal S12 protein), a variation causing an amino acid substitution (K→R) at position 43 of the protein (). This amino acid change has been associated with streptomycin resistance in many species [
19].
Figure 4. Alignment of the deduced amino acid sequence of S12 ribosomal proteins encoded by the rpsL gene from streptomycin-resistant (Smr) and -susceptible (Sms) Bifidobacterium bifidum strains. The amino acid replacement K→R at position 43 in the resistant strains is highlighted in pale blue. In bold, the strain of this study (VA07-1AN).