Chromosomal Housekeeping Genes Research Articles

Three Distinct Annotation Platforms Differ in Detection of Antimicrobial Resistance Genes in Long-Read, Short-Read, and Hybrid Sequences Derived from Total Genomic DNA or from Purified Plasmid DNA.

Recent advances and lower costs in rapid high-throughput sequencing have engendered hope that whole genome sequencing (WGS) might afford complete resistome characterization in bacterial isolates. WGS is particularly useful for the clinical characterization of fastidious and slow-growing bacteria. Despite its potential, several challenges should be addressed before adopting WGS to detect antimicrobial resistance (AMR) genes in the clinical laboratory. Here, with three distinct ESKAPE bacteria (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter spp.), different approaches were compared to identify best practices for detecting AMR genes, including: total genomic DNA and plasmid DNA extractions, the solo assembly of Illumina short-reads and of Oxford Nanopore Technologies (ONT) long-reads, two hybrid assembly pipelines, and three in silico AMR databases. We also determined the susceptibility of each strain to 21 antimicrobials. We found that all AMR genes detected in pure plasmid DNA were also detectable in total genomic DNA, indicating that, at least in these three enterobacterial genera, the purification of plasmid DNA was not necessary to detect plasmid-borne AMR genes. Illumina short-reads used with ONT long-reads in either hybrid or polished assemblies of total genomic DNA enhanced the sensitivity and accuracy of AMR gene detection. Phenotypic susceptibility closely corresponded with genotypes identified by sequencing; however, the three AMR databases differed significantly in distinguishing mobile dedicated AMR genes from non-mobile chromosomal housekeeping genes in which rare spontaneous resistance mutations might occur. This study indicates that each method employed in a WGS workflow has an impact on the detection of AMR genes. A combination of short- and long-reads, followed by at least three different AMR databases, should be used for the consistent detection of such genes. Further, an additional step for plasmid DNA purification and sequencing may not be necessary. This study reveals the need for standardized biochemical and informatic procedures and database resources for consistent, reliable AMR genotyping to take full advantage of WGS in order to expedite patient treatment and track AMR genes within the hospital and community.

Genetic diversity of microsymbionts nodulating Trifolium pratense in subpolar and temperate climate regions

Rhizobia are soil-borne bacteria forming symbiotic associations with legumes and fixing atmospheric dinitrogen. The nitrogen-fixation potential depends on the type of host plants and microsymbionts as well as environmental factors that affect the distribution of rhizobia. In this study, we compared genetic diversity of bacteria isolated from root nodules of Trifolium pratense grown in two geographical regions (Tromsø, Norway and Lublin, Poland) located in distinct climatic (subpolar and temperate) zones. To characterize these isolates genetically, three PCR-based techniques (ERIC, BOX, and RFLP of the 16S-23S rRNA intergenic spacer), 16S rRNA sequencing, and multi-locus sequence analysis of chromosomal house-keeping genes (atpD, recA, rpoB, gyrB, and glnII) were done. Our results indicate that a great majority of the isolates are T. pratense microsymbionts belonging to Rhizobium leguminosarum sv. trifolii. A high diversity among these strains was detected. However, a lower diversity within the population derived from the subpolar region in comparison to that of the temperate region was found. Multi-locus sequence analysis showed that a majority of the strains formed distinct clusters characteristic for the individual climatic regions. The subpolar strains belonged to two (A and B) and the temperate strains to three R. leguminosarum genospecies (B, E, and K), respectively.

Specialization of small non-conjugative plasmids in Escherichia coli according to their family types.

We undertook a comprehensive comparative analysis of a collection of 30 small (<25 kb) non-conjugative Escherichia coli plasmids previously classified by the gene sharing approach into 10 families, as well as plasmids found in the National Center for Biotechnology Information (NCBI) nucleotide database sharing similar genomic sequences. In total, 302 mobilizable (belonging to 2 MOBrep and 5 MOBRNA families) and 106 non-transferable/relaxase-negative (belonging to three ReLRNA families) plasmids were explored. The most striking feature was the specialization of the plasmid family types that was not related to their transmission mode and replication system. We observed a range of host strain specificity, from narrow E. coli host specificity to broad host range specificity, including a wide spectrum of Enterobacteriaceae . We found a wide variety of toxin/antitoxin systems and colicin operons in the plasmids, whose numbers and types varied according to the plasmid family type. The plasmids carried genes conferring resistance spanning almost all of the antibiotic classes, from those to which resistance developed early, such as sulphonamides, to those for which resistance has only developed recently, such as colistin. However, the prevalence of the resistance genes varied greatly according to the family type, ranging from 0 to 100 %. The evolutionary history of the plasmids based on the family type core genes showed variability within family nucleotide divergences in the range of E. coli chromosomal housekeeping genes, indicating long-term co-evolution between plasmids and host strains. In rare cases, a low evolutionary divergence suggested the massive spread of an epidemic plasmid. Overall, the importance of these small non-conjugative plasmids in bacterial adaptation varied greatly according to the type of family they belonged to, with each plasmid family having specific hosts and genetic traits.

Rhizobium laguerreae is the main nitrogen-fixing symbiont of cultivated lentil (Lens culinaris) in Morocco

Genetic diversity and population structure of 268 Lens culinaris symbiotic rhizobia collected from 40 cultivated fields in the main lentil production regions in Morocco were estimated. Three chromosomal housekeeping genes (recA, glnII and atpD) and one common symbiotic gene (nodC) were sequenced and analyzed in order to identify the local symbionts of lentil. The molecular phylogeny of the concatenated housekeeping genes clustered more than 95% of the isolates in one main clade together with Rhizobium laguerreae species. R. laguerreae represents the main symbiont of cultivated lentil in Morocco and, for the first time, a large sample of individuals is obtained for this species. There is a significant and high genetic differentiation of bacterial populations among the four regions for their symbiotic gene, and much lower for their housekeeping genes. The reasons why R. laguerreae is so frequently recovered in our study is discussed.

Phylogenomic analysis of the genus Ralstonia based on 686 single-copy genes.

The genus Ralstonia contains species that are devastating plant pathogens, opportunistic human pathogens, and/or important degraders of xenobiotic and recalcitrant compounds. However, significant nomenclature problems exist, especially for the Ralstonia solanacearum species complex which consists of four phylotypes. Phylogenomics of the Ralstonia genus was investigated via a comprehensive analysis of 39 Ralstonia genomes as well as four genomes of Cupriavidus necator (more commonly known by its previous name Ralstonia eutropha). These data revealed 686 single-copy orthologs that could be extracted from the Ralstonia core-genome and used to reconstruct the phylogeny of the genus Ralstonia. The generated tree has strong bootstrap support for almost all branches. We also estimated the in silico DNA-DNA hybridization (isDDH) and the average nucleotide identity (ANI) values between each genome. Our data confirmed that whole genome sequence data provides a powerful tool to resolve the complex taxonomic questions of the genus Ralstonia, e.g. strains of Ralstonia solanacearum phylotype IIA and IIB may represent two subspecies of R. solanacearum, and strains of R. solanacearum phylotype I and III may be classified into two subspecies of Ralstonia pseudosolanacearum. Recently, strains of R. solanacearum phylotype IV were proposed to be reclassified into different subspecies of Ralstonia syzygii; our study, however, showed that phylotype IV strains had high isDDH values (83.8-96.1%), indicating it may be not appropriate to classify these closely related strains into different subspecies. We also evaluated the performance of six chromosomal housekeeping genes (gdhA, mutS, adk, leuS, rplB and gyrB) used in Ralstonia phylogenetic inference. The multilocus sequence analysis of these six marker genes was able to reliably infer the phylogenetic relationships of the genus Ralstonia.

Assessment of the genetic and phenotypic diversity among rhizogenic Agrobacterium biovar 1 strains infecting solanaceous and cucurbit crops.

Rhizogenic Agrobacterium biovar 1 strains have been found to cause extensive root proliferation on hydroponically grown Cucurbitaceae and Solanaceae crops, resulting in substantial economic losses. As these agrobacteria live under similar ecological conditions, infecting a limited number of crops, it may be hypothesized that genetic and phenotypic variation among such strains is relatively low. In this study we assessed the phenotypic diversity as well as the phylogenetic and evolutionary relationships of several rhizogenic Agrobacterium biovar 1 strains from cucurbit and solanaceous crops. A collection of 41 isolates was subjected to a number of phenotypic assays and characterized by MLSA targeting four housekeeping genes (16S rRNA gene, recA, rpoB and trpE) and two loci from the root-inducing Ri-plasmid (part of rolB and virD2). Besides phenotypic variation, remarkable genotypic diversity was observed, especially for some chromosomal loci such as trpE. In contrast, genetic diversity was lower for the plasmid-borne loci, indicating that the studied chromosomal housekeeping genes and Ri-plasmid-borne loci might not exhibit the same evolutionary history. Furthermore, phylogenetic and network analyses and several recombination tests suggested that recombination could be contributing in some extent to the evolutionary dynamics of rhizogenic Agrobacterium populations. Finally, a genomospecies-level identification analysis revealed that at least four genomospecies may occur on cucurbit and tomato crops (G1, G3, G8 and G9). Together, this study gives a first glimpse at the genetic and phenotypic diversity within this economically important plant pathogenic bacterium.

Population structure of root nodulating Rhizobium leguminosarum in Vicia cracca populations at local to regional geographic scales

The genetic diversity and population structure of about 350 Rhizobium leguminosarum biovar viciae isolates from Vicia cracca were analysed. A hierarchical sampling design was used covering three regions, one region in Belgium and two in France, in which multiple local V. cracca populations were sampled. Rhizobium isolates were genotyped using RAPD and by sequencing two chromosomal housekeeping genes (glnII and recA) and one plasmid-borne gene (nodC). Twenty-six nodC types and sixty-seven chromosomal types were identified, many of which appeared to be regional or local endemics. We found strong genetic differentiation both among V. cracca populations that are separated by only a few kilometres, and among regions that are 50 to 350km apart. Despite significant plasmid exchange, chromosomal and nod types were similarly structured among host populations and regions. We found two lineages of which one prevailed in the Belgian region while the other dominated the French regions. Although a significant correlation between genetic differentiation and geographic distance was found, it is deemed more likely that the observed biogeographic patterns are rather due to coevolutionary interactions and environmental pressures. Furthermore, the impact of recombination on the chromosomal differentiation was found to be considerable.

Review and International Recommendation of Methods for Typing Neisseria gonorrhoeae Isolates and Their Implications for Improved Knowledge of Gonococcal Epidemiology, Treatment, and Biology

Gonorrhea, which may become untreatable due to multiple resistance to available antibiotics, remains a public health problem worldwide. Precise methods for typing Neisseria gonorrhoeae, together with epidemiological information, are crucial for an enhanced understanding regarding issues involving epidemiology, test of cure and contact tracing, identifying core groups and risk behaviors, and recommending effective antimicrobial treatment, control, and preventive measures. This review evaluates methods for typing N. gonorrhoeae isolates and recommends various methods for different situations. Phenotypic typing methods, as well as some now-outdated DNA-based methods, have limited usefulness in differentiating between strains of N. gonorrhoeae. Genotypic methods based on DNA sequencing are preferred, and the selection of the appropriate genotypic method should be guided by its performance characteristics and whether short-term epidemiology (microepidemiology) or long-term and/or global epidemiology (macroepidemiology) matters are being investigated. Currently, for microepidemiological questions, the best methods for fast, objective, portable, highly discriminatory, reproducible, typeable, and high-throughput characterization are N. gonorrhoeae multiantigen sequence typing (NG-MAST) or full- or extended-length porB gene sequencing. However, pulsed-field gel electrophoresis (PFGE) and Opa typing can be valuable in specific situations, i.e., extreme microepidemiology, despite their limitations. For macroepidemiological studies and phylogenetic studies, DNA sequencing of chromosomal housekeeping genes, such as multilocus sequence typing (MLST), provides a more nuanced understanding.

Evolutionary Dynamics of Insertion Sequences in Relation to the Evolutionary Histories of the Chromosome and Symbiotic Plasmid Genes of Rhizobium etli Populations

Insertion sequences (IS) are mobile genetic elements that are distributed in many prokaryotes. In particular, in the genomes of the symbiotic nitrogen-fixing bacteria collectively known as rhizobia, IS are fairly abundant in plasmids or chromosomal islands that carry the genes needed for symbiosis. Here, we report an analysis of the distribution and genetic conservation of the IS found in the genome of Rhizobium etli CFN42 in a collection of 87 Rhizobium strains belonging to populations with different geographical origins. We used PCR to generate presence/absence profiles of the 39 IS found in R. etli CFN42 and evaluated whether the IS were located in consistent genomic contexts. We found that the IS from the symbiotic plasmid were frequently present in the analyzed strains, whereas the chromosomal IS were observed less frequently. We then examined the evolutionary dynamics of these strains based on a population genetic analysis of two chromosomal housekeeping genes (glyA and dnaB) and three symbiotic sequences (nodC and the two IS elements). Our results indicate that the IS contained within the symbiotic plasmid have a higher degree of genomic context conservation, lower nucleotide diversity and genetic differentiation, and fewer recombination events than the chromosomal housekeeping genes. These results suggest that the R. etli populations diverged recently in Mexico, that the symbiotic plasmid also had a recent origin, and that the IS elements have undergone a process of cyclic infection and expansion.

Fine-scale phylogeographic structure of Borrelia lusitaniae revealed by multilocus sequence typing.

Borrelia lusitaniae is an Old World species of the Lyme borreliosis (LB) group of tick-borne spirochetes and prevails mainly in countries around the Mediterranean Basin. Lizards of the family Lacertidae have been identified as reservoir hosts of B. lusitaniae. These reptiles are highly structured geographically, indicating limited migration. In order to examine whether host geographic structure shapes the evolution and epidemiology of B. lusitaniae, we analyzed the phylogeographic population structure of this tick-borne bacterium using a recently developed multilocus sequence typing (MLST) scheme based on chromosomal housekeeping genes. A total of 2,099 questing nymphal and adult Ixodes ricinus ticks were collected in two climatically different regions of Portugal, being ∼130 km apart. All ticks were screened for spirochetes by direct PCR. Attempts to isolate strains yielded 16 cultures of B. lusitaniae in total. Uncontaminated cultures as well as infected ticks were included in this study. The results using MLST show that the regional B. lusitaniae populations constitute genetically distinct populations. In contrast, no clear phylogeographic signals were detected in sequences of the commonly used molecular markers ospA and ospC. The pronounced population structure of B. lusitaniae over a short geographic distance as captured by MLST of the housekeeping genes suggests that the migration rates of B. lusitaniae are rather low, most likely because the distribution of mediterranean lizard populations is highly parapatric. The study underlines the importance of vertebrate hosts in the geographic spread of tick-borne microparasites.

MLST of housekeeping genes captures geographic population structure and suggests a European origin ofBorrelia burgdorferi

Lyme borreliosis, caused by the tick-borne bacterium Borrelia burgdorferi, has become the most common vector-borne disease in North America over the last three decades. To understand the dynamics of the epizootic spread and to predict the evolutionary trajectories of B. burgdorferi, accurate information on the population structure and the evolutionary relationships of the pathogen is crucial. We, therefore, developed a multilocus sequence typing (MLST) scheme for B. burgdorferi based on eight chromosomal housekeeping genes. We validated the MLST scheme on B. burgdorferi specimens from North America and Europe, comprising both cultured isolates and infected ticks. These data were compared with sequences for the commonly used genetic markers rrs-rrlA intergenic spacer (IGS) and the gene encoding the outer surface protein C (ospC). The study demonstrates that the concatenated sequences of the housekeeping genes of B. burgdorferi provide highly resolved phylogenetic signals and that the housekeeping genes evolve differently compared with the IGS locus and ospC. Using sequence data, the study reveals that North American and European populations of B. burgdorferi correspond to genetically distinct populations. Importantly, the MLST data suggest that B. burgdorferi originated in Europe rather than in North America as proposed previously.

The Bacillus cereus group: novel aspects of population structure and genome dynamics

To provide new insights into the population and genomic structure of the Bacillus cereus group of bacteria. The genetic relatedness among B. cereus group strains was assessed by multilocus sequence typing (MLST) using an optimized scheme based on seven chromosomal housekeeping genes. A set of 48 strains from different clinical sources was included, and six clonal complexes containing several genetically similar isolates from unrelated patients were identified. Interestingly, several clonal groups contained strains that were isolated from similar human sources. Furthermore, comparative whole genome sequence analysis of 16 strains led to the discovery of novel ubiquitous genome features of the B. cereus group, such as atypical group II introns, IStrons, and hitherto uncharacterized repeated elements. The B. cereus group constitutes a coherent population unified by the presence of ubiquitous and specific genetic elements which do not show any pattern, either in their sequences or genomic locations, which allows to differentiate between the member species of the group. Nevertheless, the population is very dynamic, as particular lineages of clinical origin can evolve to form clonal complexes. At the genome level, the dynamic behaviour is indicated by the presence of numerous mobile and repeated elements. The B. cereus group of bacteria comprises species that are of medical and economic importance. The MLST data, along with the primers and protocols used, will be available in a public, web-accessible database (http://mlstoslo.uio.no).

Intergeneric transfer and recombination of the 6-phosphogluconate dehydrogenase gene (gnd) in enteric bacteria.

The gnd gene, encoding 6-phosphogluconate dehydrogenase (EC 1.1.1.44), was sequenced in 87 strains of 15 species assigned to five nominal genera of the Enterobacteriaceae, including 36 isolates of Salmonella enterica and 32 strains of Escherichia coli. In S. enterica, the effective (realized) rate of recombination of horizontally transferred gnd sequences is only moderately higher than the rates for other chromosomal housekeeping genes. In contrast, recombination at gnd has occurred with such high frequency in Escherichia coli that the indicated evolutionary relationships among strains are not congruent with those estimated by sequence analysis of other genes and by multilocus enzyme electrophoresis. E. coli and S. enterica apparently have not exchanged gnd sequences, but those of several strains of E. coli have been imported from species of Citrobacter and Klebsiella. The relatively frequent exchange of gnd within and among taxonomic groups of the Enterobacteriaceae, compared with other housekeeping genes, apparently results from its close linkage with genes that are subject to diversifying selection, including those of the rfb region determining the structure of the O antigen polysaccharide.