- Volume 8, Issue 4, 2022

The large-scale and high-intensity application of   Streptococcus thermophilus   species in milk fermentation processes is associated with a persistent threat of (bacterio)phage infection. Phage infection of starter cultures may cause inconsistent, slow or even failed fermentations with consequent diminished product quality and/or output. The phage life cycle commences with the recognition of, and binding to, a specific host-encoded and surface-exposed receptor, which in the case of   S. thermophilus   can be the rhamnose-glucose polysaccharide (RGP; specified by the rgp gene cluster) or exopolysaccharide (EPS; specified by the eps gene cluster). The genomic diversity of 23   S  .   thermophilus   strains isolated from unpasteurized dairy products was evaluated, including a detailed analysis of the rgp and eps loci. In the present study, five novel eps genotypes were identified while variations of currently recognized rgp gene cluster types were also observed. Furthermore, the diversity of rgp genotypes amongst retrieved isolates positively correlated with phage diversity based on phageome analysis of eight representative dairy products. Our findings therefore substantially expand our knowledge on S. thermophilus’ strain and phage diversity in (artisanal) dairy products and highlight the merit of phageome analysis of artisanal and traditional fermented foods as a sensitive marker of dominant microbiota involved in the fermentation.

  Wolbachia   is a genus of obligate bacterial endosymbionts that infect a diverse range of arthropod species as well as filarial nematodes, with its single described species,   Wolbachia pipientis  , divided into several ‘supergroups’ based on multilocus sequence typing.   Wolbachia   strains in mosquitoes have been shown to inhibit the transmission of human pathogens, including Plasmodium malaria parasites and arboviruses. Despite their large host range,   Wolbachia   strains within the major malaria vectors of the Anopheles gambiae and Anopheles funestus complexes appear at low density, established solely on PCR-based methods. Questions have been raised as to whether this represents a true endosymbiotic relationship. However, recent definitive evidence for two distinct, high-density strains of supergroup B   Wolbachia   within Anopheles demeilloni and Anopheles moucheti has opened exciting possibilities to explore naturally occurring   Wolbachia   endosymbionts in Anopheles for biocontrol strategies to block Plasmodium transmission. Here, we utilize genomic analyses to demonstrate that both   Wolbachia   strains have retained all key metabolic and transport pathways despite their smaller genome size, with this reduction potentially attributable to degenerated prophage regions. Even with this reduction, we confirmed the presence of cytoplasmic incompatibility (CI) factor genes within both strains, with wAnD maintaining intact copies of these genes while the cifB gene was interrupted in wAnM, so functional analysis is required to determine whether wAnM can induce CI. Additionally, phylogenetic analysis indicates that these   Wolbachia   strains may have been introduced into these two Anopheles species via horizontal transmission events, rather than by ancestral acquisition and subsequent loss events in the Anopheles gambiae species complex. These are the first   Wolbachia   genomes, to our knowledge, that enable us to study the relationship between natural strain Plasmodium malaria parasites and their anopheline hosts.

A hospital outbreak of carbapenem-resistant Enterobacterales was detected by routine surveillance. Whole genome sequencing and subsequent analysis revealed a conserved promiscuous bla OXA-48 carrying plasmid as the defining factor within this outbreak. Four different species of Enterobacterales were involved in the outbreak. Escherichia coli ST399 accounted for 35 of all the 55 isolates. Comparative genomics analysis using publicly available E. coli ST399 genomes showed that the outbreak E. coli ST399 isolates formed a unique clade. We developed a mathematical model of pOXA-48-like plasmid transmission between host lineages and used it to estimate its conjugation rate, giving a lower bound of 0.23 conjugation events per lineage per year. Our analysis suggests that co-evolution between the pOXA-48-like plasmid and E. coli ST399 could have played a role in the outbreak. This is the first study to report carbapenem-resistant E. coli ST399 carrying blaOXA-48 as the main cause of a plasmid-borne outbreak within a hospital setting. Our findings suggest complementary roles for both plasmid conjugation and clonal expansion in the emergence of this outbreak.

Symbioses between bacteria and their insect hosts can range from loose associations through to obligate interdependence. While fundamental evolutionary insights have been gained from the in-depth study of obligate mutualisms, there is increasing interest in the evolutionary potential of flexible symbiotic associations between hosts and their gut microbiomes. Understanding relationships between microbes and hosts also offers the potential for exploitation for insect control. Here, we investigate the gut microbiome of a global agricultural pest, the Mediterranean fruit fly (Ceratitis capitata). We used 16S rRNA profiling to compare the gut microbiomes of laboratory and wild strains raised on different diets and from flies collected from various natural plant hosts. The results showed that medfly guts harbour a simple microbiome that is primarily determined by the larval diet. However, regardless of the laboratory diet or natural plant host on which flies were raised, Klebsiella spp. dominated medfly microbiomes and were resistant to removal by antibiotic treatment. We sequenced the genome of the dominant putative Klebsiella spp. (‘Medkleb’) isolated from the gut of the Toliman wild-type strain. Genome-wide ANI analysis placed Medkleb within the K. oxytoca / michiganensis group. Species level taxonomy for Medkleb was resolved using a mutli-locus phylogenetic approach - and molecular, sequence and phenotypic analyses all supported its identity as K. michiganensis . Medkleb has a genome size (5825435 bp) which is 1.6 standard deviations smaller than the mean genome size of free-living Klebsiella spp. Medkleb also lacks some genes involved in environmental sensing. Moreover, the Medkleb genome contains at least two recently acquired unique genomic islands as well as genes that encode pectinolytic enzymes capable of degrading plant cell walls. This may be advantageous given that the medfly diet includes unripe fruits containing high proportions of pectin. The results suggest that the medfly harbours a commensal gut bacterium that may have developed a mutualistic association with its host and provide nutritional benefits.

The transient upsurge of G2P[4] group A rotavirus (RVA) after Rotarix vaccine introduction in several countries has been a matter of concern. To gain insight into the diversity and evolution of G2P[4] strains in South Africa pre- and post-RVA vaccination introduction, whole-genome sequencing was performed for RVA positive faecal specimens collected between 2003 and 2017 and samples previously sequenced were obtained from GenBank (n=103; 56 pre- and 47 post-vaccine). Pre-vaccine G2 sequences predominantly clustered within sub-lineage IVa-1. In contrast, post-vaccine G2 sequences clustered mainly within sub-lineage IVa-3, whereby a radical amino acid (AA) substitution, S15F, was observed between the two sub-lineages. Pre-vaccine P[4] sequences predominantly segregated within sub-lineage IVa while post-vaccine sequences clustered mostly within sub-lineage IVb, with a radical AA substitution R162G. Both S15F and R162G occurred outside recognised antigenic sites. The AA residue at position 15 is found within the signal sequence domain of Viral Protein 7 (VP7) involved in translocation of VP7 into endoplasmic reticulum during infection process. The 162 AA residue lies within the hemagglutination domain of Viral Protein 4 (VP4) engaged in interaction with sialic acid-containing structure during attachment to the target cell. Free energy change analysis on VP7 indicated accumulation of stable point mutations in both antigenic and non-antigenic regions. The segregation of South African G2P[4] strains into pre- and post-vaccination sub-lineages is likely due to erstwhile hypothesized stepwise lineage/sub-lineage evolution of G2P[4] strains rather than RVA vaccine introduction. Our findings reinforce the need for continuous whole-genome RVA surveillance and investigation of contribution of AA substitutions in understanding the dynamic G2P[4] epidemiology.

Comparative analysis of Clostridioides difficile whole-genome sequencing (WGS) data enables fine scaled investigation of transmission and is increasingly becoming part of routine surveillance. However, these analyses are constrained by the computational requirements of the large volumes of data involved. By decomposing WGS reads or assemblies into k-mers and using the dimensionality reduction technique MinHash, it is possible to rapidly approximate genomic distances without alignment. Here we assessed the performance of MinHash, as implemented by sourmash, in predicting single nucleotide differences between genomes (SNPs) and C. difficile ribotypes (RTs). For a set of 1905 diverse C. difficile genomes (differing by 0–168 519 SNPs), using sourmash to screen for closely related genomes, at a sensitivity of 100 % for pairs ≤10 SNPs, sourmash reduced the number of pairs from 1 813 560 overall to 161 934, i.e. by 91 %, with a positive predictive value of 32 % to correctly identify pairs ≤10 SNPs (maximum SNP distance 4144). At a sensitivity of 95 %, pairs were reduced by 94 % to 108 266 and PPV increased to 45 % (maximum SNP distance 1009). Increasing the MinHash sketch size above 2000 produced minimal performance improvement. We also explored a MinHash similarity-based ribotype prediction method. Genomes with known ribotypes (n=3937) were split into a training set (2937) and test set (1000) randomly. The training set was used to construct a sourmash index against which genomes from the test set were compared. If the closest five genomes in the index had the same ribotype this was taken to predict the searched genome’s ribotype. Using our MinHash ribotype index, predicted ribotypes were correct in 780/1000 (78 %) genomes, incorrect in 20 (2 %), and indeterminant in 200 (20 %). Relaxing the classifier to 4/5 closest matches with the same RT improved the correct predictions to 87 %. Using MinHash it is possible to subsample C. difficile genome k-mer hashes and use them to approximate small genomic differences within minutes, significantly reducing the search space for further analysis.

Streptococcus pneumoniae is a major human pathogen that can cause severe invasive diseases such as pneumonia, septicaemia and meningitis. Young children are at a particularly high risk, with an estimated 3–4 million cases of severe disease and between 300 000 and 500 000 deaths attributable to pneumococcal disease each year. The haemolytic toxin pneumolysin (Ply) is a primary virulence factor for this bacterium, yet despite its key role in pathogenesis, immune evasion and transmission, the regulation of Ply production is not well defined. Using a genome-wide association approach, we identified a large number of potential affectors of Ply activity, including a gene acquired horizontally on the antibiotic resistance-conferring Integrative and Conjugative Element (ICE) ICESp23FST81. This gene encodes a novel modular protein, ZomB, which has an N-terminal UvrD-like helicase domain followed by two Cas4-like domains with potent ATP-dependent nuclease activity. We found the regulatory effect of ZomB to be specific for the ply operon, potentially mediated by its high affinity for the BOX repeats encoded therein. Using a murine model of pneumococcal colonization, we further demonstrate that a ZomB mutant strain colonizes both the upper respiratory tract and lungs at higher levels when compared to the wild-type strain. While the antibiotic resistance-conferring aspects of ICESp23FST81 are often credited with contributing to the success of the S. pneumoniae lineages that acquire it, its ability to control the expression of a major virulence factor implicated in bacterial transmission is also likely to have played an important role.

The clustered regularly interspaced short palindromic repeat (CRISPR)-CRISPR-associated protein (Cas) system is an important adaptive immune system for bacteria to resist foreign DNA infection, which has been widely used in genotyping and gene editing. To provide a theoretical basis for the application of the CRISPR-Cas system in Bifidobacterium breve , the occurrence and diversity of CRISPR-Cas systems were analysed in 150 B. breve strains. Specifically, 47 % (71/150) of B. breve genomes possessed the CRISPR-Cas system, and type I-C CRISPR-Cas system was the most widely distributed among those strains. The spacer sequences present in B. breve can be used as a genotyping marker. Additionally, the phage assembly-related proteins were important targets of the type I-C CRISPR-Cas system in B. breve , and the protospacer adjacent motif sequences were further characterized in B. breve type I-C system as 5′-TTC-3′. All these results might provide a molecular basis for the development of endogenous genome editing tools in B. breve .

The arbuscular mycorrhizal fungi (AMFs) are obligate root symbionts in the subphylum Glomeromycotina that can benefit land plants by increasing their soil nutrient uptake in exchange for photosynthetically fixed carbon sources. To date, annotated genome data from representatives of the AMF orders Glomerales, Diversisporales and Archaeosporales have shown that these organisms have large and highly repeated genomes, and no genes to produce sugars and fatty acids. This led to the hypothesis that the most recent common ancestor (MRCA) of Glomeromycotina was fully dependent on plants for nutrition. Here, we aimed to further test this hypothesis by obtaining annotated genome data from a member of the early diverging order Paraglomerales (Paraglomus occultum). Genome analyses showed this species carries a 39.6 Mb genome and considerably fewer genes and repeats compared to most AMF relatives with annotated genomes. Consistent with phylogenies based on ribosomal genes, our phylogenetic analyses suggest P. occultum as the earliest diverged branch within Glomeromycotina. Overall, our analyses support the view that the MRCA of Glomeromycotina carried hallmarks of obligate plant biotrophy. The small genome size and content of P. occultum could either reflect adaptive reductive processes affecting some early AMF lineages, or indicate that the high gene and repeat family diversity thought to drive AMF adaptability to host and environmental change was not an ancestral feature of these prominent plant symbionts.

Although Salmonella Derby ST71 strains have been recognized as poultry-specific by previous studies, multiple swine-associated S. Derby ST71 strains were identified in this long-term, multi-site epidemic study. Here, 15 representative swine-associated S. Derby ST71 strains were sequenced and compared with 65 (one swine-associated and 64 poultry-associated) S. Derby ST71 strains available in the NCBI database at a pangenomic level through comparative genomics analysis to identify genomic features related to the differentiation of swine-associated strains and previously reported poultry-associated strains. The distribution patterns of known Salmonella pathogenicity islands (SPIs) and virulence factor (VF) encoding genes were not capable of differentiating between the two strain groups. The results demonstrated that the S. Derby ST71 population harbours an open pan-genome, and swine-associated ST71 strains contain many more genes than the poultry-associated strains, mainly attributed to the prophage sequence contents in the genomes. The numbers of prophage sequences identified in the swine-associated strains were higher than those in the poultry-associated strains. Prophages specifically harboured by the swine-associated strains were found to contain genes that facilitate niche adaptation for the bacterial hosts. Gene deletion experiments revealed that the dam gene specifically present in the prophage of the swine-associated strains is important for S. Derby to adhere onto the host cells. This study provides novel insights into the roles of prophages during the genome differentiation of Salmonella .

Pneumococcal serotype 35B is an important non-conjugate vaccine (non-PCV) serotype. Its continued emergence, post-PCV7 in the USA, was associated with expansion of a pre-existing 35B clone (clonal complex [CC] 558) along with post-PCV13 emergence of a non-35B clone previously associated with PCV serotypes (CC156). This study describes lineages circulating among 35B isolates in South Africa before and after PCV introduction. We also compared 35B isolates belonging to a predominant 35B lineage in South Africa (GPSC5), with isolates belonging to the same lineage in other parts of the world. Serotype 35B isolates that caused invasive pneumococcal disease in South Africa in 2005–2014 were characterized by whole-genome sequencing (WGS). Multi-locus sequence types and global pneumococcal sequence clusters (GPSCs) were derived from WGS data of 63 35B isolates obtained in 2005–2014. A total of 262 isolates that belong to GPSC5 (115 isolates from South Africa and 147 from other countries) that were sequenced as part of the global pneumococcal sequencing (GPS) project were included for comparison. Serotype 35B isolates from South Africa were differentiated into seven GPSCs and GPSC5 was most common (49 %, 31/63). While 35B was the most common serotype among GPSC5/CC172 isolates in South Africa during the PCV13 period (66 %, 29/44), 23F was the most common serotype during both the pre-PCV (80 %, 37/46) and PCV7 period (32 %, 8/25). Serotype 35B represented 15 % (40/262) of GPSC5 isolates within the global GPS database and 75 % (31/40) were from South Africa. The predominance of the GPSC5 lineage within non-vaccine serotype 35B, is possibly unique to South Africa and warrants further molecular surveillance of pneumococci.

Resistance to piperacillin/tazobactam (TZP) in Escherichia coli has predominantly been associated with mechanisms that confer resistance to third-generation cephalosporins. Recent reports have identified E. coli strains with phenotypic resistance to piperacillin/tazobactam but susceptibility to third-generation cephalosporins (TZP-R/3GC-S). In this study we sought to determine the genetic diversity of this phenotype in E. coli (n=58) isolated between 2014–2017 at a single tertiary hospital in Liverpool, UK, as well as the associated resistance mechanisms. We compare our findings to a UK-wide collection of invasive E. coli isolates (n=1509) with publicly available phenotypic and genotypic data. These data sets included the TZP-R/3GC-S phenotype (n=68), and piperacillin/tazobactam and third-generation cephalosporin-susceptible (TZP-S/3GC-S, n=1271) phenotypes. The TZP-R/3GC-S phenotype was displayed in a broad range of sequence types, which was mirrored in the same phenotype from the UK-wide collection, and the overall diversity of invasive E. coli isolates. The TZP-R/3GC-S isolates contained a diverse range of plasmids, indicating multiple acquisition events of TZP resistance mechanisms rather than clonal expansion of a particular plasmid or sequence type. The putative resistance mechanisms were equally diverse, including hyperproduction of TEM-1, either via strong promoters or gene amplification, carriage of inhibitor-resistant β-lactamases, and an S133G bla CTX-M-15 mutation detected for the first time in clinical isolates. Several of these mechanisms were present at a lower abundance in the TZP-S/3GC-S isolates from the UK-wide collection, but without the associated phenotypic resistance to TZP. Eleven (19%) of the isolates had no putative mechanism identified from the genomic data. Our findings highlight the complexity of this cryptic phenotype and the need for continued phenotypic monitoring, as well as further investigation to improve detection and prediction of the TZP-R/3GC-S phenotype from genomic data.

Infections with globally disseminated Shiga toxin-producing Escherichia coli (STEC) of the O113:H21 serotype can progress to severe clinical complications, such as hemolytic uremic syndrome (HUS). Two phylogeographically distinct clonal complexes have been established by multi locus sequence typing (MLST). Infections with ST-820 isolates circulating exclusively in Australia have caused severe human disease, such as HUS. Conversely, ST-223 isolates prevalent in the US and outside Australia seem to rarely cause severe human disease but are frequent contaminants. Following a genomic epidemiology approach, we wanted to gain insights into the underlying cause for this disparity. We examined the plasticity in the genome make-up and Shiga toxin production in a collection of 20 ST-820 and ST-223 strains isolated from produce, the bovine reservoir, and clinical cases. STEC are notorious for assembly into fragmented draft sequences when using short-read sequencing technologies due to the extensive and partly homologous phage complement. The application of long-read technology (LRT) sequencing yielded closed reference chromosomes and plasmids for two representative ST-820 and ST-223 strains. The established high-resolution framework, based on whole genome alignments, single nucleotide polymorphism (SNP)-typing and MLST, includes the chromosomes and plasmids of other publicly available O113:H21 sequences and allowed us to refine the phylogeographical boundaries of ST-820 and ST-223 complex isolates and to further identify a historic non-shigatoxigenic strain from Mexico as a quasi-intermediate. Plasmid comparison revealed strong correlations between the strains’ featured pO113 plasmid genotypes and chromosomally inferred ST, which suggests coevolution of the chromosome and virulence plasmids. Our pathogenicity assessment revealed statistically significant differences in the Stx2a-production capabilities of ST-820 as compared to ST-223 strains under RecA-induced Stx phage mobilization, a condition that mimics Stx-phage induction. These observations suggest that ST-820 strains may confer an increased pathogenic potential in line with the strain-associated epidemiological metadata. Still, some of the tested ST-223 cultures sourced from contaminated produce or the bovine reservoir also produced Stx at levels comparable to those of ST-820 isolates, which calls for awareness and for continued surveillance of this lineage.

Carbapenem-resistant   Enterobacterales   (CRE) are a serious public health threat because of their rapid dissemination. To determine the epidemiological and genetic characteristics of CRE infections in Thailand, we performed whole-genome sequencing of 577 carbapenem-resistant   Klebsiella pneumoniae   isolates and 170 carbapenem-resistant   Escherichia coli   isolates from hospitals across the nation. The four most prevalent carbapenemase genes harboured by these bacteria were bla NDM-1, bla NDM-5, bla OXA-181 and bla OXA-232. The gene bla NDM-1 was identified in diverse sequence types. The gene bla NDM-5 was identified almost exclusively in   E. coli  . The genes bla OXA-181, bla OXA-232, and co-carriage of bla NDM-1 and bla OXA-232 were found in specific sequence types from certain provinces. Replicon typing revealed the diverse backbones of bla NDM-1- and bla NDM-5-harbouring plasmids and successful expansion of bla NDM-1-harbouring IncN2-type plasmids. Core-genome single-nucleotide polymorphism analysis suggested that bla OXA-181-, bla OXA-232-, bla NDM-5-, and co-carriage of bla NDM-1 and bla OXA-232-associated sub-clonal lineages have recently predominated in the provinces from where these isolates were isolated. Thus, we demonstrate genotype-dependent dissemination of CRE in Thailand, which is helpful for establishing infection-control strategies in CRE-endemic areas.

Extensively drug-resistant tuberculosis (XDR-TB), defined as resistance to at least isoniazid (INH), rifampicin (RIF), a fluoroquinolone (FQ) and a second-line injectable drug (SLID), is difficult to treat and poses a major threat to TB control. The transmission dynamics and distribution of XDR Mycobacterium tuberculosis (Mtb) strains have not been thoroughly investigated. Using whole genome sequencing data on 461 XDR-Mtb strains, we aimed to investigate the geographical distribution of XDR-Mtb strains in the Western Cape Province of South Africa over a 10 year period (2006–2017) and assess the association between Mtb sub-lineage, age, gender, geographical patient location and membership or size of XDR-TB clusters. First, we identified transmission clusters by excluding drug resistance-conferring mutations and using the 5 SNP cutoff, followed by merging clusters based on their most recent common ancestor. We then consecutively included variants conferring resistance to INH, RIF, ethambutol (EMB), pyrazinamide (PZA), SLIDs and FQs in the cluster definition. Cluster sizes were classified as small (2–4 isolates), medium (5–20 isolates), large (21–100 isolates) or very large (>100 isolates) to reflect the success of individual strains. We found that most XDR-TB strains were clustered and that including variants conferring resistance to INH, RIF, EMB, PZA and SLIDs in the cluster definition did not significantly reduce the proportion of clustered isolates (85.5–82.2 %) but increased the number of patients belonging to small clusters (4.3–12.4 %, P=0.56). Inclusion of FQ resistance-conferring variants had the greatest effect, with 11 clustered isolates reclassified as unique while the number of clusters increased from 17 to 37. Lineage 2 strains (lineage 2.2.1 typical Beijing or lineage 2.2.2 atypical Beijing) showed the large clusters which were spread across all health districts of the Western Cape Province. We identified a significant association between residence in the Cape Town metropole and cluster membership (P=0.016) but no association between gender, age and cluster membership or cluster size (P=0.39). Our data suggest that the XDR-TB epidemic in South Africa probably has its origin in the endemic spread of MDR Mtb and pre-XDR Mtb strains followed by acquisition of FQ resistance, with more limited transmission of XDR Mtb strains. This only became apparent with the inclusion of drug resistance-conferring variants in the definition of a cluster. In addition to the prevention of amplification of resistance, rapid diagnosis of MDR, pre-XDR and XDR-TB and timely initiation of appropriate treatment is needed to reduce transmission of difficult-to-treat TB.

Cladosporium fulvum is a fungal pathogen that causes leaf mould of tomato. The reference genome of this pathogen was released in 2012 but its high repetitive DNA content prevented a contiguous assembly and further prohibited the analysis of its genome architecture. In this study, we combined third generation sequencing technology with the Hi-C chromatin conformation capture technique, to produce a high-quality and near complete genome assembly and gene annotation of a Race 5 isolate of C. fulvum. The resulting genome assembly contained 67.17 Mb organized into 14 chromosomes (Chr1-to-Chr14), all of which were assembled telomere-to-telomere. The smallest of the chromosomes, Chr14, is only 460 kb in size and contains 25 genes that all encode hypothetical proteins. Notably, PCR assays revealed that Chr14 was absent in 19 out of 24 isolates of a world-wide collection of C. fulvum, indicating that Chr14 is dispensable. Thus, C. fulvum is currently the second species of Capnodiales shown to harbour dispensable chromosomes. The genome of C. fulvum Race 5 is 49.7 % repetitive and contains 14 690 predicted genes with an estimated completeness of 98.9%, currently one of the highest among the Capnodiales. Genome structure analysis revealed a compartmentalized architecture composed of gene-dense and repeat-poor regions interspersed with gene-sparse and repeat-rich regions. Nearly 39.2 % of the C. fulvum Race 5 genome is affected by Repeat-Induced Point (RIP) mutations and evidence of RIP leakage toward non-repetitive regions was observed in all chromosomes, indicating the RIP plays an important role in the evolution of this pathogen. Finally, 345 genes encoding candidate effectors were identified in C. fulvum Race 5, with a significant enrichment of their location in gene-sparse regions, in accordance with the ‘two-speed genome’ model of evolution. Overall, the new reference genome of C. fulvum presents several notable features and is a valuable resource for studies in plant pathogens.

Soybean is the most important legume cropped worldwide and can highly benefit from the biological nitrogen fixation (BNF) process. Brazil is recognized for its leadership in the use of inoculants and two strains, Bradyrhizobium japonicum CPAC 15 (=SEMIA 5079) and Bradyrhizobium diazoefficiens CPAC 7 (=SEMIA 5080) compose the majority of the 70 million doses of soybean inoculants commercialized yearly in the country. We studied a collection of natural variants of these two strains, differing in properties of competitiveness and efficiency of BNF. We sequenced the genomes of the parental strain SEMIA 566 of B. japonicum , of three natural variants of this strain (S 204, S 340 and S 370), and compared with another variant of this group, strain CPAC 15. We also sequenced the genome of the parental strain SEMIA 586 of B. diazoefficiens , of three natural variants of this strain (CPAC 390, CPAC 392 and CPAC 394) and compared with the genome of another natural variant, strain CPAC 7. As the main genes responsible for nodulation (nod, noe, nol) and BNF (nif, fix) in soybean Bradyrhizobium are located in symbiotic islands, our objective was to identify genetic variations located in this region, including single nucleotide polymorphisms (SNPs) and insertions and deletions (indels), that could be potentially related to their different symbiotic phenotypes. We detected 44 genetic variations in the B. japonicum strains and three in B. diazoefficiens . As the B. japonicum strains have gone through a longer period of adaptation to the soil, the higher number of genetic variations could be explained by survival strategies under the harsh environmental conditions of the Brazilian Cerrado biome. Genetic variations were detected in genes enconding proteins such as a dephospho-CoA kinase, related to the CoA biosynthesis; a glucosamine-fructose-6-phosphate aminotransferase, key regulator of the hexosamine biosynthetic pathway; a LysR family transcriptional regulator related to nodulation genes; and NifE and NifS proteins, directly related to the BNF process. We suggest potential genetic variations related to differences in the symbiotic phenotypes.