Horizontally Acquired Genes Are Often Shared between Closely Related Bacterial Species

Genetic Transfer in Action: Uncovering DNA Flow in an Extremophilic Microbial Community

Horizontal genetic transfer (HGT) is a significant driver of genomic novelty in all domains of life. HGT has been investigated in many studies however, the focus has been on conspicuous protein-coding DNA transfers that often prove to be adaptive in recipient organisms and are therefore fixed longer-term in lineages. These results comprise a subclass of HGTs and do not represent exhaustive (coding and non-coding) DNA transfer and its impact on ecology. Uncovering exhaustive HGT can provide key insights into the connectivity of genomes in communities and how these transfers may occur. In this study, we use the term frequency-inverse document frequency (TF-IDF) technique, that has been used successfully to mine DNA transfers within real and simulated high-quality prokaryote genomes, to search for exhaustive HGTs within an extremophilic microbial community. We establish a pipeline for validating transfers identified using this approach. We find that most DNA transfers are within-domain and involve non-coding DNA. A relatively high proportion of the predicted protein-coding HGTs appear to encode transposase activity, restriction-modification system components, and biofilm formation functions. Our study demonstrates the utility of the TF-IDF approach for HGT detection and provides insights into the mechanisms of recent DNA transfer.

Horizontal transmission of a multidrug-resistant IncN plasmid isolated from urban wastewater

Wastewater treatment plants (WWTPs) are considered reservoirs of antibiotic resistance genes (ARGs). Given that plasmid-mediated horizontal gene transfer plays a critical role in disseminating ARGs in the environment, it is important to inspect the transfer potential of transmissible plasmids to have a better understanding of whether these mobile ARGs can be hosted by opportunistic pathogens and should be included in One Health's considerations. In this study, we used a fluorescent-reporter-gene based exogenous isolation approach to capture extended-spectrum beta-lactamases encoding mobile determinants from sewer microbiome samples that enter an urban water system (UWS) in Denmark. After screening and sequencing, we isolated a ∼73 Kbp IncN plasmid (pDK_DARWIN) that harboured and expressed multiple ARGs. Using a dual fluorescent reporter gene system, we showed that this plasmid can transfer into resident urban water communities. We demonstrated the transfer of pDK_DARWIN to microbiome members of both the sewer (in the upstream UWS compartment) and wastewater treatment (in the downstream UWS compartment) microbiomes. Sequence similarity search across curated plasmid repositories revealed that pDK_DARWIN derives from an IncN backbone harboured by environmental and nosocomial Enterobacterial isolates. Furthermore, we searched for pDK_DARWIN sequence matches in UWS metagenomes from three countries, revealing that this plasmid can be detected in all of them, with a higher relative abundance in hospital sewers compared to residential sewers. Overall, this study demonstrates that this IncN plasmid is prevalent across Europe and an efficient vector capable of disseminating multiple ARGs in the urban water systems.

Detecting horizontal gene transfer among microbiota: an innovative pipeline for identifying co-shared genes within the mobilome through advanced comparative analysis

Horizontal gene transfer (HGT) is a key driver in the evolution of bacterial genomes. The acquisition of genes mediated by HGT may enable bacteria to adapt to ever-changing environmental conditions. Long-term application of antibiotics in intensive agriculture is associated with the dissemination of antibiotic resistance genes among bacteria with the consequences causing public health concern. Commensal farm-animal-associated gut microbiota are considered the reservoir of the resistance genes. Therefore, in this study, we identified known and not-yet characterized mobilized genes originating from chicken and porcine fecal samples using our innovative pipeline followed by network analysis to provide appropriate visualization to support proper interpretation.

Molecular epidemiology and characterization of antimicrobial-resistant Staphylococcus haemolyticus strains isolated from dairy cattle milk in Northwest, China

Introduction

Non-aureus Staphylococcus (NAS) species are currently the most commonly identified microbial agents causing sub-clinical infections of the udder and are also deemed as opportunistic pathogens of clinical mastitis in dairy cattle. More than 10 NAS species have been identified and studied but little is known about S. haemolyticus in accordance with dairy mastitis. The present study focused on the molecular epidemiology and genotypic characterization of S. haemolyticus isolated from dairy cattle milk in Northwest, China.

Methods

In this study, a total of 356 milk samples were collected from large dairy farms in three provinces in Northwest, China. The bacterial isolation and presumptive identification were done by microbiological and biochemical methods following the molecular confirmation by 16S rRNA gene sequencing. The antimicrobial susceptibility testing (AST) was done by Kirby-Bauer disk diffusion assay and antibiotic-resistance genes (ARGs) were identified by PCR. The phylogenetic grouping and sequence typing was done by Pulsed Field Gel Electrophoresis (PFGE) and Multi-Locus Sequence Typing (MLST) respectively.

Results

In total, 39/356 (11.0%) were identified as positive for S. haemolyticus. The overall prevalence of other Staphylococcus species was noted to be 39.6% (141/356), while the species distribution was as follows: S. aureus 14.9%, S. sciuri 10.4%, S. saprophyticus 7.6%, S. chromogenes 4.2%, S. simulans 1.4%, and S. epidermidis 1.1%. The antimicrobial susceptibility of 39 S. haemolyticus strains exhibited higher resistance to erythromycin (92.3%) followed by trimethoprim-sulfamethoxazole (51.3%), ciprofloxacin (43.6%), florfenicol (30.8%), cefoxitin (28.2%), and gentamicin (23.1%). All of the S. haemolyticus strains were susceptible to tetracycline, vancomycin, and linezolid. The overall percentage of multi-drug resistant (MDR) S. haemolyticus strains was noted to be 46.15% (18/39). Among ARGs, mphC was identified as predominant (82.05%), followed by ermB (33.33%), floR (30.77%), gyrA (30.77%), sul1 (28.21%), ermA (23.08%), aadD (12.82%), grlA (12.82%), aacA-aphD (10.26%), sul2 (10.26%), dfrA (7.69%), and dfrG (5.13%). The PFGE categorized 39 S. haemolyticus strains into A-H phylogenetic groups while the MLST categorized strains into eight STs with ST8 being the most predominant while other STs identified were ST3, ST11, ST22, ST32, ST19, ST16, and ST7.

Conclusion

These findings provided new insights into our understanding of the epidemiology and genetic characteristics of S. haemolyticus in dairy farms to inform interventions limiting the spread of AMR in dairy production.

Molecular characterization of antimicrobial resistance related genes in E. coli, Salmonella and Klebsiella isolates from broilers in the West Region of Cameroon

BACKGROUND

Antibiotic resistance has become an enduring threat to human health. This has prompted extensive research to identify the determinants responsible in a bid to fight the spread of resistance and also develop new antibiotics. However, routine procedures focus on identifying genetic determinants of resistance only on phenotypically resistant isolates. We aimed to characterise plasmid mediated resistance determinants in key Enterobacteriaceae isolates with differential phenotypic susceptibility profiles and evaluated the contribution of resistance genes on phenotypic expression of susceptibility.

METHODS

The study was carried out on 200 Enterobacteriaceae isolates belonging to the genera E. coli, Salmonella, and Klebsiella; 100 resistant and 100 susceptible to quinolones, aminoglycosides, and ESBL-producing as determined by disk diffusion. Reduced susceptibility in susceptible isolates was determined as an increased MIC by broth microdilution. Plasmid-borne resistance genes were sought in all isolates by endpoint PCR. We performed correlations tests to determine the relationship between the occurrence of resistance genes and increased MIC in susceptible isolates. We then used the notion of penetrance to show adequacy between resistance gene carriage and phenotypic resistance as well as diagnostic odds ratio to evaluate how predictable phenotypic susceptibility profile could determine the presence of resistant genes in the isolates.

RESULTS

Reduced susceptibility was detected in 30% (9/30) ESBL negative, 50% (20/40) quinolone-susceptible and 53.33% (16/30) aminoglycoside-susceptible isolates. Plasmid-borne resistance genes were detected in 50% (15/30) of ESBL negative, 65% (26/40) quinolone susceptible and 66.67% (20/30) aminoglycoside susceptible isolates. Reduced susceptibility increased the risk of susceptible isolates carrying resistance genes (ORs 4.125, 8.36, and 8.89 respectively for ESBL, quinolone, and aminoglycoside resistance genes). Resistance gene carriage correlated significantly to reduced susceptibility for quinolone and aminoglycoside resistance genes (0.002 and 0.015 at CI95). Gene carriage correlated with phenotypic resistance at an estimated 64.28% for ESBL, 56.90% for quinolone, and 58.33% for aminoglycoside resistance genes.

CONCLUSIONS

A high carriage of plasmid-mediated genes for ESBL, quinolone, and aminoglycoside resistance was found among the Enterobacteriaceae tested. However, gene carriage was not always correlated with phenotypic expression. This allows us to suggest that assessing genetic determinants of resistance should not be based on AST profile only. Further studies, including assessing the role of chromosomal determinants will shed light on other factors that undermine antimicrobial susceptibility locally.

Selection of Multidrug-Resistant Enterobacteria in Weaned Pigs and Its Association With In-feed Subtherapeutic Combination of Colistin and Tylosin

In-feed antibiotics are administered to piglets to improve performance and production efficiency. However, the use of growth promoters in the swine industry can select for multidrug-resistant (MDR) bacteria. Here, we evaluate the resistance profile of enterobacteria isolated from fecal samples of weaned pigs (21-35 days) fed or not with antibiotics (colistin and tylosin) and investigated the piglets gut microbiota in both groups. Six hundred and eighteen bacterial cultures were isolated from the control group (CON; n = 384) and antibiotic-fed pigs (ATB; n = 234). All isolates were tested for resistance to 12 antibiotics belonging to six distinct antibiotic classes. Isolates were highly resistant to ampicillin (90%; n = 553), amoxicillin (85%; n = 525), and tetracycline (81%; n = 498). A significant increase (P < 0.05) in resistance to cephalexin, kanamycin, doxycycline, and colistin was observed for bacteria from the ATB group. Piglets allocated in the ATB and CON groups shared similar intestinal microbiota, as revealed by alpha- and beta-diversity analyses. Our findings demonstrate that colistin and tylosin contribute to select MDR enterobacteria in weaned piglets. The high frequency of antibiotic resistance among isolates from the CON group suggests that environmental sources (e.g., fecal contents, aerosols, soil, water, food) also represent a potential reservoir of multidrug-resistant enterobacteria in pig production systems.

Phylogenomic disentangling of the Bifidobacterium longum subsp. infantis taxon

Members of the Bifidobacterium longum species have been shown to possess adaptive abilities to allow colonization of different mammalian hosts, including humans, primates and domesticated mammalian species, such as dogs, horses, cattle and pigs. To date, three subspecies have formally been recognized to belong to this bifidobacterial taxon, i.e. B. longum subsp. longum, B. longum subsp. infantis and B. longum subsp. suis. Although B. longum subsp. longum is widely distributed in the human gut irrespective of host age, B. longum subsp. infantis appears to play a significant role as a prominent member of the gut microbiota of breast-fed infants. Nevertheless, despite the considerable scientific relevance of these taxa and the vast body of genomic data now available, an accurate dissection of the genetic features that comprehensively characterize the B. longum species and its subspecies is still missing. In the current study, we employed 261 publicly available B. longum genome sequences, combined with those of 11 new isolates, to investigate genomic diversity of this taxon through comparative genomic and phylogenomic approaches. These analyses allowed us to highlight a remarkable intra-species genetic and physiological diversity. Notably, characterization of the genome content of members of B. longum subsp. infantis subspecies suggested that this taxon may have acquired genetic features for increased competitiveness in the gut environment of suckling hosts. Furthermore, specific B. longum subsp. infantis genomic features appear to be responsible for enhanced horizontal gene transfer (HGT) occurrences, underpinning an intriguing dedication toward acquisition of foreign DNA by HGT events.

Evolutionary responses to codon usage of horizontally transferred genes in Pseudomonas aeruginosa: gene retention, amelioration and compensatory evolution

Prokaryote genome evolution is characterized by the frequent gain of genes through horizontal gene transfer (HGT). For a gene, being horizontally transferred can represent a strong change in its genomic and physiological context. If the codon usage of a transferred gene deviates from that of the receiving organism, the fitness benefits it provides can be reduced due to a mismatch with the expression machinery. Consequently, transferred genes with a deviating codon usage can be selected against or elicit evolutionary responses that enhance their integration, such as gene amelioration and compensatory evolution. Within bacterial species, the extent and relative importance of these different mechanisms has never been considered altogether. In this study, a phylogeny-based method was used to investigate the occurrence of these different evolutionary responses in Pseudomonas aeruginosa. Selection on codon usage of genes acquired through HGT was observed over evolutionary time, with the overall codon usage converging towards that of the core genome. Gene amelioration, through the accumulation of synonymous mutations after HGT, did not seem to systematically affect transferred genes. This pattern therefore seemed to be mainly driven by selective retention of transferred genes with an initial codon usage similar to that of the core genes. Additionally, variation in the copy number of tRNA genes was often associated with the acquisition of genes for which the observed variation could enhance their expression. This provides evidence that compensatory evolution might be an important mechanism for the integration of horizontally transferred genes.

Genomic GC content drifts downward in most bacterial genomes

In every kingdom of life, GC->AT transitions occur more frequently than any other type of mutation due to the spontaneous deamination of cytidine. In eukaryotic genomes, this slow loss of GC base pairs is counteracted by biased gene conversion which increases genomic GC content as part of the recombination process. However, this type of biased gene conversion has not been observed in bacterial genomes, so we hypothesized that GC->AT transitions cause a reduction of genomic GC content in prokaryotic genomes on an evolutionary time scale. To test this hypothesis, we used a phylogenetic approach to analyze triplets of closely related genomes representing a wide range of the bacterial kingdom. The resulting data indicate that genomic GC content is drifting downward in bacterial genomes where GC base pairs comprise 40% or more of the total genome. In contrast, genomes containing less than 40% GC base pairs have fewer opportunities for GC->AT transitions to occur so genomic GC content is relatively stable or actually increasing. It should be noted that this observed change in genomic GC content is the net change in shared parts of the genome and does not apply to parts of the genome that have been lost or acquired since the genomes being compared shared common ancestor. However, a more detailed analysis of two Caulobacter genomes revealed that the acquisition of mobile elements by the two genomes actually reduced the total genomic GC content as well.

Insights on the Evolutionary Genomics of the Blautia Genus: Potential New Species and Genetic Content Among Lineages

Blautia, a genus established in 2008, is a relevantly abundant taxonomic group present in the microbiome of human and other mammalian gastrointestinal (GI) tracts. Several described (or proposed) Blautia species are available at this date. However, despite the increasing level of knowledge about Blautia, its diversity is still poorly understood. The increasing availability of Blautia genomic sequences in the public databases opens the possibility to study this genus from a genomic perspective. Here we report the pangenome analysis and the phylogenomic study of 225 Blautia genomes available in RefSeq. We found 33 different potential species at the genomic level, 17 of them previously undescribed; we also confirmed by genomic standards the status of 4 previously proposed new Blautia species. Comparative genomic analyses suggest that the Blautia pangenome is open, with a relatively small core genome (∼ 700-800 gene families). Utilizing a set of representative genomes, we performed a gene family gain/loss model for the genus, showing that despite terminal nodes suffered more massive gene gain events than internal nodes (i.e., predicted ancestors), some ancestors were predicted to have gained an important number of gene families, some of them associated with the possible acquisition of metabolic abilities. Gene loss events remained lower than gain events in most cases. General aspects regarding pangenome composition and gene gain/loss events are discussed, as well as the proposition of changes in the taxonomic assignment of B. coccoides TY and the proposition of a new species, "B. pseudococcoides.".

Time-calibrated genomic evolution of a monomorphic bacterium during its establishment as an endemic crop pathogen

Horizontal gene transfer is of major evolutionary importance as it allows for the redistribution of phenotypically important genes among lineages. Such genes with essential functions include those involved in resistance to antimicrobial compounds and virulence factors in pathogenic bacteria. Understanding gene turnover at microevolutionary scales is critical to assess the pace of this evolutionary process. Here, we characterized and quantified gene turnover for the epidemic lineage of a bacterial plant pathogen of major agricultural importance worldwide. Relying on a dense geographic sampling spanning 39 years of evolution, we estimated both the dynamics of single nucleotide polymorphism accumulation and gene content turnover. We identified extensive gene content variation among lineages even at the smallest phylogenetic and geographic scales. Gene turnover rate exceeded nucleotide substitution rate by three orders of magnitude. Accessory genes were found preferentially located on plasmids, but we identified a highly plastic chromosomal region hosting ecologically important genes such as transcription activator-like effectors. Whereas most changes in the gene content are probably transient, the rapid spread of a mobile element conferring resistance to copper compounds widely used for the management of plant bacterial pathogens illustrates how some accessory genes can become ubiquitous within a population over short timeframes.

Bacterial Interactions as Detected by Pooled Antibiotic Susceptibility Testing (P-AST) in Polymicrobial Urine Specimens

INTRODUCTION

Antimicrobial susceptibility is well characterized in monomicrobial infections, but bacterial species often coexist with other bacterial species. Antimicrobial susceptibility is often tested against single bacterial isolates; this approach ignores interactions between cohabiting bacteria that could impact susceptibility. Here, we use Pooled Antibiotic Susceptibility Testing to compare antimicrobial susceptibility patterns exhibited by polymicrobial and monomicrobial urine specimens obtained from patients with urinary tract infection symptoms.

METHODS

Urine samples were collected from patients who had symptoms consistent with a urinary tract infection. Multiplex polymerase chain reaction testing was performed to identify and quantify 31 bacterial species. Antibiotic susceptibility was determined using a novel Pooled Antibiotic Susceptibility Testing method. Antibiotic resistance rates in polymicrobial specimens were compared with those in monomicrobial infections. Using a logistic model, resistance rates were estimated when specific bacterial species were present. To assess interactions between pairs of bacteria, the predicted resistance rates were compared when a pair of bacterial species were present versus when just one bacterial species was present.

RESULTS

Urine specimens were collected from 3,124 patients with symptoms of urinary tract infection. Of these, multiplex polymerase chain reaction testing detected bacteria in 61.1% (1910) of specimens. Pooled Antibiotic Susceptibility Testing results were available for 70.8% (1352) of these positive specimens. Of these positive specimens, 43.9% (594) were monomicrobial, while 56.1% (758) were polymicrobial. The odds of resistance to ampicillin (p = 0.005), amoxicillin/clavulanate (p = 0.008), five different cephalosporins, vancomycin (p = <0.0001), and tetracycline (p = 0.010) increased with each additional species present in a polymicrobial specimen. In contrast, the odds of resistance to piperacillin/tazobactam decreased by 75% for each additional species present (95% CI 0.61, 0.94, p = 0.010). For one or more antibiotics tested, thirteen pairs of bacterial species exhibited statistically significant interactions compared with the expected resistance rate obtained with the Highest Single Agent Principle and Union Principle.

CONCLUSION

Bacterial interactions in polymicrobial specimens can result in antimicrobial susceptibility patterns that are not detected when bacterial isolates are tested by themselves. Optimizing an effective treatment regimen for patients with polymicrobial infections may depend on accurate identification of the constituent species, as well as results obtained by Pooled Antibiotic Susceptibility Testing.

PPanGGOLiN: Depicting microbial diversity via a partitioned pangenome graph

The use of comparative genomics for functional, evolutionary, and epidemiological studies requires methods to classify gene families in terms of occurrence in a given species. These methods usually lack multivariate statistical models to infer the partitions and the optimal number of classes and don't account for genome organization. We introduce a graph structure to model pangenomes in which nodes represent gene families and edges represent genomic neighborhood. Our method, named PPanGGOLiN, partitions nodes using an Expectation-Maximization algorithm based on multivariate Bernoulli Mixture Model coupled with a Markov Random Field. This approach takes into account the topology of the graph and the presence/absence of genes in pangenomes to classify gene families into persistent, cloud, and one or several shell partitions. By analyzing the partitioned pangenome graphs of isolate genomes from 439 species and metagenome-assembled genomes from 78 species, we demonstrate that our method is effective in estimating the persistent genome. Interestingly, it shows that the shell genome is a key element to understand genome dynamics, presumably because it reflects how genes present at intermediate frequencies drive adaptation of species, and its proportion in genomes is independent of genome size. The graph-based approach proposed by PPanGGOLiN is useful to depict the overall genomic diversity of thousands of strains in a compact structure and provides an effective basis for very large scale comparative genomics. The software is freely available at https://github.com/labgem/PPanGGOLiN.

Recombination and gene loss occur simultaneously during bacterial horizontal gene transfer

Bacteria can acquire new genes by incorporating environmental DNA into their genomes, yet genome sizes stay relatively constant. In nature, gene acquisition is a rare event so it is difficult to observe. However, the Caulobacter crescentus CB2A genome contains 114 insertions of genetic material from the closely-related NA1000 strain, providing a unique opportunity to analyze the horizontal transfer of genetic material. Analyses of these insertions led to a new model that involves preferential recombination at non-homologous regions that are flanked by regions of homology and does not involve any mutational processes. The net result is the replacement of segments of the recipient genome instead of the simple addition of genetic material during horizontal gene transfer. Analyses of the genomes of closely related strains of other bacterial and archaea genera, suggested that horizontal gene transfer occurs preferentially in non-homologous regions in these organisms as well. Thus, it appears to be a general phenomenon that prokaryotic horizontal gene transfer occurs preferentially at sites where the incoming DNA contains a non-homologous region that is flanked by regions of homology. Therefore, gene replacement is a common phenomenon during horizontal gene transfer.

MetaCHIP: community-level horizontal gene transfer identification through the combination of best-match and phylogenetic approaches

BACKGROUND

Metagenomic datasets provide an opportunity to study horizontal gene transfer (HGT) on the level of a microbial community. However, current HGT detection methods cannot be applied to community-level datasets or require reference genomes. Here, we present MetaCHIP, a pipeline for reference-independent HGT identification at the community level.

RESULTS

Assessment of MetaCHIP's performance on simulated datasets revealed that it can predict HGTs with various degrees of genetic divergence from metagenomic datasets. The results also indicated that the detection of very recent gene transfers (i.e. those with low levels of genetic divergence) from metagenomics datasets is largely affected by the read assembly step. Comparison of MetaCHIP with a previous analysis on soil bacteria showed a high level of consistency for the prediction of recent HGTs and revealed a large number of additional non-recent gene transfers, which can provide new biological and ecological insight. Assessment of MetaCHIP's performance on real metagenomic datasets confirmed the role of HGT in the spread of genes related to antibiotic resistance in the human gut microbiome. Further testing also showed that functions related to energy production and conversion as well as carbohydrate transport and metabolism are frequently transferred among free-living microorganisms.

CONCLUSION

MetaCHIP provides an opportunity to study HGTs among members of a microbial community and therefore has several applications in the field of microbial ecology and evolution. MetaCHIP is implemented in Python and freely available at https://github.com/songweizhi/MetaCHIP .

Genome Comparisons of Wild Isolates of Caulobacter crescentus Reveal Rates of Inversion and Horizontal Gene Transfer

Since previous interspecies comparisons of Caulobacter genomes have revealed extensive genome rearrangements, we decided to compare the nucleotide sequences of four C. crescentus genomes, NA1000, CB1, CB2, and CB13. To accomplish this goal, we used PacBio sequencing technology to determine the nucleotide sequence of the CB1, CB2, and CB13 genomes, and obtained each genome sequence as a single contig. To correct for possible sequencing errors, each genome was sequenced twice. The only differences we observed between the two sets of independently determined sequences were random omissions of a single base in a small percentage of the homopolymer regions where a single base is repeated multiple times. Comparisons of these four genomes indicated that horizontal gene transfer events that included small numbers of genes occurred at frequencies in the range of 10^-3 to 10^-4 insertions per generation. Large insertions were about 100 times less frequent. Also, in contrast to previous interspecies comparisons, we found no genome rearrangements when the closely related NA1000, CB1, and CB2 genomes were compared, and only eight inversions and one translocation when the more distantly related CB13 genome was compared to the other genomes. Thus, we estimate that inversions occur at a rate of one per 10 to 12 million generations in Caulobacter genomes. The inversions seem to be complex events that include the simultaneous creation of indels.