Genetic adaptation to amoxicillin in Escherichia coli: The limited role of dinB and katE

Bacteria can quickly adapt to sub-lethal concentrations of antibiotics. Several stress and DNA repair genes contribute to this adaptation process. However, the pathways leading to adaptation by acquisition of de novo mutations remain poorly understood. This study explored the roles of DNA polymerase IV (dinB) and catalase HP2 (katE) in E. coli's adaptation to amoxicillin. These genes are thought to play essential roles in beta-lactam resistance-dinB in increasing mutation rates and katE in managing oxidative stress. By comparing the adaptation rates, transcriptomic profiles, and genetic changes of wild-type and knockout strains, we aimed to clarify the contributions of these genes to beta-lactam resistance. While all strains exhibited similar adaptation rates and mutations in the frdD gene and ampC operon, several unique mutations were acquired in the ΔkatE and ΔdinB strains. Overall, this study distinguishes the contributions of general stress-related genes on the one hand, and dinB, and katE on the other hand, in development of beta-lactam resistance.

Holistic understanding of trimethoprim resistance in Streptococcus pneumoniae using an integrative approach of genome-wide association study, resistance reconstruction, and machine learning

Antimicrobial resistance (AMR) is a public health threat worldwide. Next-generation sequencing (NGS) has opened unprecedented opportunities to accelerate AMR mechanism discovery and diagnostics. Here, we present an integrative approach to investigate trimethoprim (TMP) resistance in the key pathogen Streptococcus pneumoniae. We explored a collection of 662 S. pneumoniae genomes by conducting a genome-wide association study (GWAS), followed by functional validation using resistance reconstruction experiments, combined with machine learning (ML) approaches to predict TMP minimum inhibitory concentration (MIC). Our study showed that multiple additive mutations in the folA and sulA loci are responsible for TMP non-susceptibility in S. pneumoniae and can be used as key features to build ML models for digital MIC prediction, reaching an average accuracy within ±1 twofold dilution factor of 86.3%. Our roadmap of in silico analysis-wet-lab validation-diagnostic tool building could be adapted to explore AMR in other combinations of bacteria-antibiotic.

IMPORTANCE

In the age of next-generation sequencing (NGS), while data-driven methods such as genome-wide association study (GWAS) and machine learning (ML) excel at finding patterns, functional validation can be challenging due to the high numbers of candidate variants. We designed an integrative approach combining a GWAS on S. pneumoniae clinical isolates, followed by whole-genome transformation coupled with NGS to functionally characterize a large set of GWAS candidates. Our study validated several phenotypic folA mutations beyond the standard Ile100Leu mutation, and showed that the overexpression of the sulA locus produces trimethoprim (TMP) resistance in Streptococcus pneumoniae. These validated loci, when used to build ML models, were found to be the best inputs for predicting TMP minimal inhibitory concentrations. Integrative approaches can bridge the genotype-phenotype gap by biological insights that can be incorporated in ML models for accurate prediction of drug susceptibility.

ProcaryaSV: structural variation detection pipeline for bacterial genomes using short-read sequencing

BACKGROUND

Structural variations play an important role in bacterial genomes. They can mediate genome adaptation quickly in response to the external environment and thus can also play a role in antibiotic resistance. The detection of structural variations in bacteria is challenging, and the recognition of even small rearrangements can be important. Even though most detection tools are aimed at and benchmarked on eukaryotic genomes, they can also be used on prokaryotic genomes. The key features of detection are the ability to detect small rearrangements and support haploid genomes. Because of the limiting performance of a single detection tool, combining the detection abilities of multiple tools can lead to more robust results. There are already available workflows for structural variation detection for long-reads technologies and for the detection of single-nucleotide variation and indels, both aimed at bacteria. Yet we are unaware of structural variations detection workflows for the short-reads sequencing platform. Motivated by this gap we created our workflow. Further, we were interested in increasing the detection performance and providing more robust results.

RESULTS

We developed an open-source bioinformatics pipeline, ProcaryaSV, for the detection of structural variations in bacterial isolates from paired-end short sequencing reads. Multiple tools, starting with quality control and trimming of sequencing data, alignment to the reference genome, and multiple structural variation detection tools, are integrated. All the partial results are then processed and merged with an in-house merging algorithm. Compared with a single detection approach, ProcaryaSV has improved detection performance and is a reproducible easy-to-use tool.

CONCLUSIONS

The ProcaryaSV pipeline provides an integrative approach to structural variation detection from paired-end next-generation sequencing of bacterial samples. It can be easily installed and used on Linux machines. It is publicly available on GitHub at https://github.com/robinjugas/ProcaryaSV .

Microevolution of Salmonella 4,[5],12:i:- derived from Salmonella enterica serovar Typhimurium through complicated transpositions

Salmonella enterica subsp. enterica serovar 4,[5],12:i:- (Salmonella 4,[5],12:i:-), derived from S. Typhimurium, has become the dominant serotype causing human salmonellosis. In this study, we define the genetic mechanism of the generation of Salmonella 4,[5],12:i:- from S. Typhimurium through complicated transpositions and demonstrate that Salmonella 4,[5],12:i:- displays more efficient colonization and survival abilities in mice than its parent S. Typhimurium strain. We identified intermediate strains carrying both resistance regions (RRs) and the fljAB operon for the generation of Salmonella 4,[5],12:i:-. The insertion of RR3 into the chromosomal hin-iroB site of S. Typhimurium produced RR3-S. Typhimurium as a primary intermediate. Salmonella 4,[5],12:i:- was then produced by replacing the fljAB operon and/or its flanking sequences through intramolecular transpositions mediated by IS26 and/or IS1R elements in RR3-S. Typhimurium, which was further confirmed both in vitro and in vivo. Overall, we demonstrate the molecular mechanism underlying the origin, generation, and advantage of RRs-Salmonella 4,[5],12:i:- from S. Typhimurium.

In Vitro Selection Identifies Staphylococcus aureus Genes Influencing Biofilm Formation

Staphylococcus aureus generates biofilms during many chronic human infections, which contributes to its growth and persistence in the host. Multiple genes and pathways necessary for S. aureus biofilm production have been identified, but knowledge is incomplete, and little is known about spontaneous mutations that increase biofilm formation as infection progresses. Here, we performed in vitro selection of four S. aureus laboratory strains (ATCC 29213, JE2, N315, and Newman) to identify mutations associated with enhanced biofilm production. Biofilm formation increased in passaged isolates from all strains, exhibiting from 1.2- to 5-fold the capacity of parental lines. Whole-genome sequencing identified nonsynonymous mutations affecting 23 candidate genes and a genomic duplication encompassing sigB. Six candidate genes significantly impacted biofilm formation as isogenic transposon knockouts: three were previously reported to impact S. aureus biofilm formation (icaR, spdC, and codY), while the remaining three (manA, narH, and fruB) were newly implicated by this study. Plasmid-mediated genetic complementation of manA, narH, and fruB transposon mutants corrected biofilm deficiencies, with high-level expression of manA and fruB further enhancing biofilm formation over basal levels. This work recognizes genes not previously identified as contributing to biofilm formation in S. aureus and reveals genetic changes able to augment biofilm production by that organism.

Carbapenemase-Encoding Gene Copy Number Estimator (CCNE): a Tool for Carbapenemase Gene Copy Number Estimation

Carbapenemase production is one of the leading mechanisms of carbapenem resistance in Gram-negative bacteria. An increase in carbapenemase gene (blaCarb) copies is an important mechanism of carbapenem resistance. No currently available bioinformatics tools allow for reliable detection and reporting of carbapenemase gene copy numbers. Here, we describe the carbapenemase-encoding gene copy number estimator (CCNE), a ready-to-use bioinformatics tool that was developed to estimate blaCarb copy numbers from whole-genome sequencing data. Its performance on Klebsiella pneumoniae carbapenemase gene (blaKPC) copy number estimation was evaluated by simulation and quantitative PCR (qPCR), and the results were compared with available algorithms. CCNE has two components, CCNE-acc and CCNE-fast. CCNE-acc detects blaCarb copy number in a comprehensive and high-accuracy way, while CCNE-fast rapidly screens blaCarb copy numbers. CCNE-acc achieved the best accuracy (100%) and the lowest root mean squared error (RMSE; 0.07) in simulated noise data sets, compared to the assembly-based method (23.4% accuracy, 1.697 RMSE) and the OrthologsBased method (78.9% accuracy, 0.395 RMSE). In the qPCR validation, a high consistency was observed between the blaKPC copy number determined by qPCR and that determined with CCNE. Reverse transcription-qPCR transcriptional analysis of 40 isolates showed that blaKPC expression was positively correlated with the blaKPC copy numbers detected by CCNE (P < 0.001). An association study of 357 KPC-producing K. pneumoniae isolates and their antimicrobial susceptibility identified a significant association between the estimated blaKPC copy number and MICs of imipenem (P < 0.001) and ceftazidime-avibactam (P < 0.001). Overall, CCNE is a useful genomic tool for the analysis of antimicrobial resistance genes copy number; it is available at https://github.com/biojiang/ccne. IMPORTANCE Globally disseminated carbapenem-resistant Enterobacterales is an urgent threat to public health. The most common carbapenem resistance mechanism is the production of carbapenemases. Carbapenemase-producing isolates often exhibit a wide range of carbapenem MICs. Higher carbapenem MICs have been associated with treatment failure. The increase of carbapenemase gene (blaCarb) copy numbers contributes to increased carbapenem MICs. However, blaCarb gene copy number detection is not routinely conducted during a genomic analysis, in part due to the lack of optimal bioinformatics tools. In this study, we describe a ready-to-use tool we developed and designated the carbapenemase-encoding gene copy number estimator (CCNE) that can be used to estimate the blaCarb copy number directly from whole-genome sequencing data, and we extended the data to support the analysis of all known blaCarb genes and some other antimicrobial resistance genes. Furthermore, CCNE can be used to interrogate the correlations between genotypes and susceptibility phenotypes and to improve our understanding of antimicrobial resistance mechanisms.

Pre-existing chromosomal polymorphisms in pathogenic E. coli potentiate the evolution of resistance to a last-resort antibiotic

Bacterial pathogens show high levels of chromosomal genetic diversity, but the influence of this diversity on the evolution of antibiotic resistance by plasmid acquisition remains unclear. Here, we address this problem in the context of colistin, a 'last line of defence' antibiotic. Using experimental evolution, we show that a plasmid carrying the MCR-1 colistin resistance gene dramatically increases the ability of Escherichia coli to evolve high-level colistin resistance by acquiring mutations in lpxC, an essential chromosomal gene involved in lipopolysaccharide biosynthesis. Crucially, lpxC mutations increase colistin resistance in the presence of the MCR-1 gene, but decrease the resistance of wild-type cells, revealing positive sign epistasis for antibiotic resistance between the chromosomal mutations and a mobile resistance gene. Analysis of public genomic datasets shows that lpxC polymorphisms are common in pathogenic E. coli, including those carrying MCR-1, highlighting the clinical relevance of this interaction. Importantly, lpxC diversity is high in pathogenic E. coli from regions with no history of MCR-1 acquisition, suggesting that pre-existing lpxC polymorphisms potentiated the evolution of high-level colistin resistance by MCR-1 acquisition. More broadly, these findings highlight the importance of standing genetic variation and plasmid/chromosomal interactions in the evolutionary dynamics of antibiotic resistance.

Niche-specific genome degradation and convergent evolution shaping Staphylococcus aureus adaptation during severe infections

During severe infections, Staphylococcus aureus moves from its colonising sites to blood and tissues and is exposed to new selective pressures, thus, potentially driving adaptive evolution. Previous studies have shown the key role of the agr locus in S. aureus pathoadaptation; however, a more comprehensive characterisation of genetic signatures of bacterial adaptation may enable prediction of clinical outcomes and reveal new targets for treatment and prevention of these infections. Here, we measured adaptation using within-host evolution analysis of 2590 S. aureus genomes from 396 independent episodes of infection. By capturing a comprehensive repertoire of single nucleotide and structural genome variations, we found evidence of a distinctive evolutionary pattern within the infecting populations compared to colonising bacteria. These invasive strains had up to 20-fold enrichments for genome degradation signatures and displayed significantly convergent mutations in a distinctive set of genes, linked to antibiotic response and pathogenesis. In addition to agr-mediated adaptation, we identified non-canonical, genome-wide significant loci including sucA-sucB and stp1. The prevalence of adaptive changes increased with infection extent, emphasising the clinical significance of these signatures. These findings provide a high-resolution picture of the molecular changes when S. aureus transitions from colonisation to severe infection and may inform correlation of infection outcomes with adaptation signatures.

Localized pmrB hypermutation drives the evolution of colistin heteroresistance

Colistin has emerged as an important last line of defense for the treatment of infections caused by antibiotic-resistant gram-negative pathogens, but colistin resistance remains poorly understood. Here, we investigate the responses of ≈1,000 populations of a multi-drug-resistant (MDR) strain of P. aeruginosa to a high dose of colistin. Colistin exposure causes rapid cell death, but some populations eventually recover due to the growth of sub-populations of heteroresistant cells. Heteroresistance is unstable, and resistance is rapidly lost under culture in colistin-free medium. The evolution of heteroresistance is primarily driven by selection for heteroresistance at two hotspot sites in the PmrAB regulatory system. Localized hypermutation of pmrB generates colistin resistance at 10³–10⁴ times the background resistance mutation rate (≈2 × 10^-5 per cell division). PmrAB provides resistance to antimicrobial peptides that are involved in host immunity, suggesting that this pathogen may have evolved a highly mutable pmrB as an adaptation to host immunity.

•

Pseudomonas populations recover from colistin due to the growth of heteroresistant cells

•

Heteroresistance is driven by pre-existing mutations in the PmrAB regulatory system

•

pmrB mutations arise at 10³–10⁴ times the background mutation rate

•

Heteroresistance is unstable and is rapidly lost in the absence of colistin

OUP accepted manuscript

Objectives

Cefiderocol is a siderophore cephalosporin active against MDR Gram-negatives including Stenotrophomonas maltophilia. Cefiderocol resistance remains uncommon and incompletely understood. We selected for cefiderocol-resistant S. maltophilia in vitro and characterized the genetic mechanisms and potential for cross-resistance to other antimicrobials.

Methods

We selected cefiderocol resistance in three clinical strains of S. maltophilia by serial passage in escalating concentrations of cefiderocol. Emergent cefiderocol-resistant isolates were subjected to repeat susceptibility testing against a panel of relevant antimicrobials. Isolates with confirmed MIC changes were whole genome sequenced.

Results

Each parent strain was initially susceptible to cefiderocol (MICs of 0.03125, 0.03125 and 0.125 mg/L), and one initially tested susceptible to ceftazidime/avibactam (MIC 4 mg/L). We recovered evolved isolates achieving cefiderocol resistance at MICs of 8–32 mg/L from each parental strain. Some cefiderocol resistant isolates reverted following one to four drug-free passages. Ceftazidime/avibactam MICs of passaged isolates repeatedly increased to ≥256 mg/L, and while other MICs were largely unchanged, trimethoprim/sulfamethoxazole MICs declined 4-fold in two strains. WGS revealed one evolved isolate carrying six coding mutations, while four were isogenic mutants of tonB, tolQ, smf-1 and the smeT promoter. Mutation of the smeT promoter downregulated the smeDEF efflux pump and reduced susceptibility to penicillins but increased susceptibility to several other classes including sulphonamides. Other mutations occurred in genes putatively involved in iron metabolism including smlt1148 and cirA.

Conclusions

S. maltophilia strains evolved cefiderocol resistance through different genetic pathways, but often involved iron transport. Future work is required to fully understand the role(s) of other genes in cefiderocol resistance.

Introduction and adaptation of an emerging pathogen to olive trees in Italy

The invasive plant pathogen Xylella fastidiosa currently threatens European flora through the loss of economically and culturally important host plants. This emerging vector-borne bacterium, native to the Americas, causes several important diseases in a wide range of plants including crops, ornamentals, and trees. Previously absent from Europe, and considered a quarantine pathogen, X. fastidiosa was first detected in Apulia, Italy in 2013 associated with a devastating disease of olive trees (Olive Quick Decline Syndrome, OQDS). OQDS has led to significant economic, environmental, cultural, as well as political crises. Although the biology of X. fastidiosa diseases have been studied for over a century, there is still no information on the determinants of specificity between bacterial genotypes and host plant species, which is particularly relevant today as X. fastidiosa is expanding in the naive European landscape. We analysed the genomes of 79 X. fastidiosa samples from diseased olive trees across the affected area in Italy as well as genomes of the most genetically closely related strains from Central America. We provided insights into the ecological and evolutionary emergence of this pathogen in Italy. We first showed that the outbreak in Apulia is due to a single introduction from Central America that we estimated to have occurred in 2008 [95 % HPD: 1930–2016]. By using a combination of population genomic approaches and evolutionary genomics methods, we further identified a short list of genes that could play a major role in the adaptation of X. fastidiosa to this new environment. We finally provided experimental evidence for the adaptation of the strain to this new environment.

Alterations in chromosomal genes nfsA, nfsB, and ribE are associated with nitrofurantoin resistance in Escherichia coli from the United Kingdom

Antimicrobial resistance in enteric or urinary Escherichia coli is a risk factor for invasive E. coli infections. Due to widespread trimethoprim resistance amongst urinary E. coli and increased bacteraemia incidence, a national recommendation to prescribe nitrofurantoin for uncomplicated urinary tract infection was made in 2014. Nitrofurantoin resistance is reported in <6% urinary E. coli isolates in the UK, however, mechanisms underpinning nitrofurantoin resistance in these isolates remain unknown. This study aimed to identify the genetic basis of nitrofurantoin resistance in urinary E. coli isolates collected from north west London and then elucidate resistance-associated genetic alterations in available UK E. coli genomes. As a result, an algorithm was developed to predict nitrofurantoin susceptibility. Deleterious mutations and gene-inactivating insertion sequences in chromosomal nitroreductase genes nfsA and/or nfsB were identified in genomes of nine confirmed nitrofurantoin-resistant urinary E. coli isolates and additional 11 E. coli isolates that were highlighted by the prediction algorithm and subsequently validated to be nitrofurantoin-resistant. Eight categories of allelic changes in nfsA , nfsB , and the associated gene ribE were detected in 12412 E. coli genomes from the UK. Evolutionary analysis of these three genes revealed homoplasic mutations and explained the previously reported order of stepwise mutations. The mobile gene complex oqxAB , which is associated with reduced nitrofurantoin susceptibility, was identified in only one of the 12412 genomes. In conclusion, mutations and insertion sequences in nfsA and nfsB were leading causes of nitrofurantoin resistance in UK E. coli . As nitrofurantoin exposure increases in human populations, the prevalence of nitrofurantoin resistance in carriage E. coli isolates and those from urinary and bloodstream infections should be monitored.

CNproScan: Hybrid CNV detection for bacterial genomes

Discovering copy number variation (CNV) in bacteria is not in the spotlight compared to the attention focused on CNV detection in eukaryotes. However, challenges arising from bacterial drug resistance bring further interest to the topic of CNV and its role in drug resistance. General CNV detection methods do not consider bacteria's features and there is space to improve detection accuracy. Here, we present a CNV detection method called CNproScan focused on bacterial genomes. CNproScan implements a hybrid approach and other bacteria-focused features and depends only on NGS data. We benchmarked our method and compared it to the previously published methods and we can resolve to achieve a higher detection rate together with providing other beneficial features, such as CNV classification. Compared with other methods, CNproScan can detect much shorter CNV events.

Shared and Unique Evolutionary Trajectories to Ciprofloxacin Resistance in Gram-Negative Bacterial Pathogens

Resistance to the broad-spectrum antibiotic ciprofloxacin is detected at high rates for a wide range of bacterial pathogens. To investigate the dynamics of ciprofloxacin resistance development, we applied a comparative resistomics workflow for three clinically relevant species of Gram-negative bacteria: Escherichia coli, Acinetobacter baumannii, and Pseudomonas aeruginosa. We combined experimental evolution in a morbidostat with deep sequencing of evolving bacterial populations in time series to reveal both shared and unique aspects of evolutionary trajectories. Representative clone characterization by sequencing and MIC measurements enabled direct assessment of the impact of mutations on the extent of acquired drug resistance. In all three species, we observed a two-stage evolution: (i) early ciprofloxacin resistance reaching 4- to 16-fold the MIC for the wild type, commonly as a result of single mutations in DNA gyrase target genes (gyrA or gyrB), and (ii) additional genetic alterations affecting the transcriptional control of the drug efflux machinery or secondary target genes (DNA topoisomerase parC or parE).

New Genotype of Yersinia pestis Found in Live Rodents in Yunnan Province, China

Yunnan Province, China is thought to be the original source of biovar Orientalis of Yersinia pestis, the causative agent of the third plague pandemic that has spread globally since the end of the 19th century. Although encompassing a large area of natural plague foci, Y. pestis strains have rarely been found in live rodents during surveillance in Yunnan, and most isolates are from rodent corpses and their fleas. In 2017, 10 Y. pestis strains were isolated from seven live rodents and three fleas in Heqing County of Yunnan. These strains were supposed to have low virulence to local rodents Eothenomys miletus and Apodemus chevrieri because the rodents were healthy and no dead animals were found in surrounding areas, as had occurred in previous epizootic disease. We performed microscopic and biochemical examinations of the isolates, and compared their whole-genome sequences and transcriptome with those of 10 high virulence Y. pestis strains that were isolated from nine rodents and one parasitic flea in adjacent city (Lijiang). We analyzed the phenotypic, genomic, and transcriptomic characteristics of live rodent isolates. The isolates formed a previously undefined monophyletic branch of Y. pestis that was named 1.IN5. Six SNPs, two indels, and one copy number variation were detected between live rodent isolates and the high virulence neighbors. No obvious functional consequence of these variations was found according to the known annotation information. Among genes which expression differential in the live rodent isolates compared to their high virulent neighbors, we found five iron transfer related ones that were significant up-regulated (| log₂ (FC) | > 1, p.adjust < 0.05), indicating these genes may be related to the low-virulence phenotype. The novel genotype of Y. pestis reported here provides further insights into the evolution and spread of plague as well as clues that may help to decipher the virulence mechanism of this notorious pathogen.

Genomic analyses of Staphylococcus aureus clonal complex 45 isolates does not distinguish nasal carriage from bacteraemia

Staphylococcus aureus is a colonizing opportunistic pathogen and a leading cause of bloodstream infection with high morbidity and mortality. S. aureus carriage frequency is reportedly between 20 and 40 % among healthy adults, with S. aureus colonization considered to be a risk factor for S. aureus bacteraemia. It is unknown whether a genetic component of the bacterium is associated with S. aureus bacteraemia in comparison to nasal carriage strains. Previous association studies primarily focusing on the clinical outcome of an S. aureus infection have produced conflicting results, often limited by study design challenged by sample collections and the clonal diversity of S. aureus. To date, no study has investigated whether genomic features separate nasal carriage isolates from S. aureus bacteraemia isolates within a single clonal lineage. Here we have investigated whether genomic features, including single-nucleotide polymorphisms (SNPs), genes, or kmers, distinguish S. aureus nasal carriage isolates from bacteraemia isolates that all belong to the same clonal lineage [clonal complex 45 (CC45)] using whole-genome sequencing (WGS) and a genome-wide association (GWA) approach. From CC45, 100 isolates (50 bacteraemia and 50 nasal carriage, geographically and temporally matched) from Denmark were whole-genome sequenced and subjected to GWA analyses involving gene copy number variation, SNPs, gene content, kmers and gene combinations, while correcting for lineage effects. No statistically significant association involving SNPs, specific genes, gene variants, gene copy number variation, or a combination of genes was identified that could distinguish bacteraemia isolates from nasal carriage isolates. The presented results suggest that all S. aureus nasal CC45 isolates carry the potential to cause invasive disease, as no core or accessory genome content or variations were statistically associated with invasiveness.

CNV-BAC: Copy number Variation Detection in Bacterial Circular Genome

Motivation

Whole-genome sequencing (WGS) is widely used for copy number variation (CNV) detection. However, for most bacteria, their circular genome structure and high replication rate make reads more enriched near the replication origin. CNV detection based on read depth could be seriously influenced by such replication bias.

Results

We show that the replication bias is widespread using ∼200 bacterial WGS data. We develop CNV-BAC (CNV-Bacteria) that can properly normalize the replication bias and other known biases in bacterial WGS data and can accurately detect CNVs. Simulation and real data analysis show that CNV-BAC achieves the best performance in CNV detection compared with available algorithms.

Availability and implementation

CNV-BAC is available at https://github.com/XiDsLab/CNV-BAC.

Supplementary information

Supplementary data are available at Bioinformatics online.

Current Affairs of Microbial Genome-Wide Association Studies: Approaches, Bottlenecks and Analytical Pitfalls

Microbial genome-wide association studies (mGWAS) are a new and exciting research field that is adapting human GWAS methods to understand how variations in microbial genomes affect host or pathogen phenotypes, such as drug resistance, virulence, host specificity and prognosis. Several computational tools and methods have been developed or adapted from human GWAS to facilitate the discovery of novel mutations and structural variations that are associated with the phenotypes of interest. However, no comprehensive, end-to-end, user-friendly tool is currently available. The development of a broadly applicable pipeline presents a real opportunity among computational biologists. Here, (i) we review the prominent and promising tools, (ii) discuss analytical pitfalls and bottlenecks in mGWAS, (iii) provide insights into the selection of appropriate tools, (iv) highlight the gaps that still need to be filled and how users and developers can work together to overcome these bottlenecks. Use of mGWAS research can inform drug repositioning decisions as well as accelerate the discovery and development of more effective vaccines and antimicrobials for pressing infectious diseases of global health significance, such as HIV, TB, influenza, and malaria.

The Copy Number of the spoVA 2mob Operon Determines Pressure Resistance of Bacillus Endospores

The spoVA ^2mob operon confers heat resistance to Bacillus spp., and the resistance correlates to the copy number of the operon. Bacillus endospores also exhibit a strong variation in resistance to pressure, but the underlying mechanisms of endospore resistance to pressure are not fully understood. We determined the effects of multiple spoVA ^2mob operons on high-pressure resistance in Bacillus endospores. The copy numbers of the spoVA ^2mob operon in 17 strains of Bacillus subtilis, Bacillus amyloliquefaciens, Bacillus cereus, Bacillus velezensis, and Bacillus pumilus were determined via droplet digital PCR (ddPCR) and genome sequencing. These strains contained between 0 and 3 copies of the spoVA ^2mob operon; the quantification of the gene copy number by ddPCR was as accurate as whole-genome sequencing. We further tested the pressure resistance of 17 Bacillus endospores at 600 MPa and 80°C. Strains with one or no spoVA ^2mob operon had significantly lower pressure resistance than strains with two or three copies of the operons (P < 0.001), indicating that redundant spoVA ^2mob operons in Bacillus contributed to higher pressure resistance of endospores. The copy number of the spoVA ^2mob operon was not related to the dipicolinic acid (DPA) content of endospores. Overall, the copy number of the spoVA ^2mob operon contributes to pressure resistance of Bacillus endospores. This improves our understanding of the pressure resistance mechanisms in Bacillus spp. and may inform the development of high-pressure sterilization in food processing.IMPORTANCE Bacillus spp. are considered pressure-resistant microorganisms, but the resistance mechanisms remain unknown. The spoVA ^2mob operon is a mobile genetic element, and it can transfer to pathogenic or spoilage organisms by horizontal gene transfer. Results in this study indicate that multiple copies of the spoVA ^2mob operon mediate high-pressure resistance of Bacillus endospores, and it might contribute to the identification of the source of pressure-resistant pathogens and spoilage organisms that may contaminate the food supply. The droplet digital PCR (ddPCR) system is well suited for analysis in some human diseases due to its high efficiency and capability to provide high precision; however, no relevant studies in food microbiology have been reported so far. This study demonstrates a novel application of ddPCR in food microbiology.

A large-scale whole-genome comparison shows that experimental evolution in response to antibiotics predicts changes in naturally evolved clinical Pseudomonas aeruginosa

Pseudomonas aeruginosa is an opportunistic pathogen that causes a wide range of acute and chronic infections. An increasing number of isolates have mutations that make them antibiotic resistant, making treatment difficult. To identify resistance-associated mutations we experimentally evolved the antibiotic sensitive strain P. aeruginosa PAO1 to become resistant to three widely used anti-pseudomonal antibiotics, ciprofloxacin, meropenem and tobramycin. Mutants could tolerate up to 2048-fold higher concentrations of antibiotic than strain PAO1. Genome sequences were determined for thirteen mutants for each antibiotic. Each mutant had between 2 and 8 mutations. For each antibiotic at least 8 genes were mutated in multiple mutants, demonstrating the genetic complexity of resistance. For all three antibiotics mutations arose in genes known to be associated with resistance, but also in genes not previously associated with resistance. To determine the clinical relevance of mutations uncovered in this study we analysed the corresponding genes in 558 isolates of P. aeruginosa from patients with chronic lung disease and in 172 isolates from the general environment. Many genes identified through experimental evolution had predicted function-altering changes in clinical isolates but not in environmental isolates, showing that mutated genes in experimentally evolved bacteria can predict those that undergo mutation during infection. Additionally, large deletions of up to 479kb arose in experimentally evolved meropenem resistant mutants and large deletions were present in 87 of the clinical isolates. These findings significantly advance understanding of antibiotic resistance in P. aeruginosa and demonstrate the validity of experimental evolution in identifying clinically-relevant resistance-associated mutations.

Loss and Gain in the Evolution of the Salmonella enterica Serovar Gallinarum Biovar Pullorum Genome

Pullorum disease, an acute poultry septicemia caused by Salmonella Gallinarum biovar Pullorum, is fatal for young chickens and is a heavy burden on poultry industry. The pathogen is rare in most developed countries but still extremely difficult to eliminate in China. Efficient epidemiological surveillance necessitates clarifying the origin of the isolates from different regions and their phylogenic relationships. Genomic epidemiological analysis of 97 S. Pullorum strains was carried out to reconstruct the phylogeny and transmission history of S. Pullorum. Further analysis demonstrated that functional gene loss and acquisition occurred simultaneously throughout the evolution of S. Pullorum, both of which reflected adaptation to the changing environment. The result of our study will be helpful in surveillance and prevention of pullorum disease.

Salmonella enterica subspecies enterica serovar Gallinarum biovar Pullorum (S. Pullorum) is the etiological agent of pullorum disease, causing white diarrhea with high mortality in chickens. There are many unsolved issues surrounding the epidemiology of S. Pullorum, including its origin and transmission history as well as the discordance between its phenotypic heterogeneity and genetic monomorphism. In this paper, we report the results of whole-genome sequencing of a panel of 97 S. Pullorum strains isolated between 1962 and 2014 from four countries across three continents. We utilized 6,795 core genome single nucleotide polymorphisms (SNPs) to reconstruct a phylogenetic tree within a spatiotemporal Bayesian framework, estimating that the most recent common ancestor of S. Pullorum emerged in ∼914 CE (95% confidence interval [95%CI], 565 to 1273 CE). The extant S. Pullorum strains can be divided into four distinct lineages, each of which is significantly associated with geographical distribution. The intercontinental transmissions of lineages III and IV can be traced to the mid-19th century and are probably related to the “Hen Fever” prevalent at that time. Further genomic analysis indicated that the loss or pseudogenization of functional genes involved in metabolism and virulence in S. Pullorum has been ongoing since before and after divergence from the ancestor. In contrast, multiple prophages and plasmids have been acquired by S. Pullorum, and these have endowed it with new characteristics, especially the multidrug resistance conferred by two large plasmids in lineage I. The results of this study provide insight into the evolution of S. Pullorum and prove the efficiency of whole-genome sequencing in epidemiological surveillance of pullorum disease.

IMPORTANCE Pullorum disease, an acute poultry septicemia caused by Salmonella Gallinarum biovar Pullorum, is fatal for young chickens and is a heavy burden on poultry industry. The pathogen is rare in most developed countries but still extremely difficult to eliminate in China. Efficient epidemiological surveillance necessitates clarifying the origin of the isolates from different regions and their phylogenic relationships. Genomic epidemiological analysis of 97 S. Pullorum strains was carried out to reconstruct the phylogeny and transmission history of S. Pullorum. Further analysis demonstrated that functional gene loss and acquisition occurred simultaneously throughout the evolution of S. Pullorum, both of which reflected adaptation to the changing environment. The result of our study will be helpful in surveillance and prevention of pullorum disease.

Sequana coverage: detection and characterization of genomic variations using running median and mixture models

Background

In addition to mapping quality information, the Genome coverage contains valuable biological information such as the presence of repetitive regions, deleted genes, or copy number variations (CNVs). It is essential to take into consideration atypical regions, trends (e.g., origin of replication), or known and unknown biases that influence coverage. It is also important that reported events have robust statistics (e.g. z-score) associated with their detections as well as precise location.

Results

We provide a stand-alone application, sequana_coverage, that reports genomic regions of interest (ROIs) that are significantly over- or underrepresented in high-throughput sequencing data. Significance is associated with the events as well as characteristics such as length of the regions. The algorithm first detrends the data using an efficient running median algorithm. It then estimates the distribution of the normalized genome coverage with a Gaussian mixture model. Finally, a z-score statistic is assigned to each base position and used to separate the central distribution from the ROIs (i.e., under- and overcovered regions). A double thresholds mechanism is used to cluster the genomic ROIs. HTML reports provide a summary with interactive visual representations of the genomic ROIs with standard plots and metrics. Genomic variations such as single-nucleotide variants or CNVs can be effectively identified at the same time.

Read Depth Analysis to Identify CNV in Bacteria Using CNOGpro

Whole-genome sequencing with short-read technologies is well suited for calling single nucleotide polymorphisms, but has major problems with the detection of structural variants larger than the read length. One such type of variation is copy number variation (CNV), which entails deletion or duplication of genomic regions, and the expansion or contraction of repeated elements. Duplicated and deleted regions will typically be collapsed during de novo assembly of sequence data, or ignored when mapping reads toward a reference. However, signatures of the copy number variation can be detected in the resultant read depth at each position in the genome. We here provide instructions on how to analyze this read depth signal with the R package CNOGpro, allowing for estimation of copy numbers with uncertainty for each feature in a genome.

Parallel Evolution of Group B Streptococcus Hypervirulent Clonal Complex 17 Unveils New Pathoadaptive Mutations

The incidence of group B Streptococcus (GBS) neonatal disease continues to be a significant cause of concern worldwide. Strains belonging to clonal complex 17 (CC17) are the most frequently responsible for GBS infections in neonates, especially among late-onset disease cases. Therefore, we undertook the largest genomic study of GBS CC17 strains to date to decipher the genetic bases of their remarkable colonization and infection ability. We show that crucial functions involved in different steps of the colonization or infection process of GBS are distinctly mutated during the adaptation of CC17 to the human host. In particular, our results implicate the CovRS two-component regulator of virulence in the differentiation between carriage- and disease-associated isolates. Not only does this work raise important implications for the ongoing development of a vaccine against GBS but might also drive the discovery of key functions for GBS adaptation and pathogenesis that have been overlooked until now.

Group B Streptococcus (GBS) is a commensal of the gastrointestinal and genitourinary tracts, while a prevailing cause of neonatal disease worldwide. Of the various clonal complexes (CCs), CC17 is overrepresented in GBS-infected newborns for reasons that are still largely unknown. Here, we report a comprehensive genomic analysis of 626 CC17 isolates collected worldwide, identifying the genetic traits behind their successful adaptation to humans and the underlying differences between carriage and clinical strains. Comparative analysis with 923 GBS genomes belonging to CC1, CC19, and CC23 revealed that the evolution of CC17 is distinct from that of other human-adapted lineages and recurrently targets functions related to nucleotide and amino acid metabolism, cell adhesion, regulation, and immune evasion. We show that the most distinctive features of disease-specific CC17 isolates were frequent mutations in the virulence-associated CovS and Stk1 kinases, underscoring the crucial role of the entire CovRS regulatory pathway in modulating the pathogenicity of GBS. Importantly, parallel and convergent evolution of major components of the bacterial cell envelope, such as the capsule biosynthesis operon, the pilus, and Rib, reflects adaptation to host immune pressures and should be taken into account in the ongoing development of a GBS vaccine. The presence of recurrent targets of evolution not previously implicated in virulence also opens the way for uncovering new functions involved in host colonization and GBS pathogenesis.

IMPORTANCE The incidence of group B Streptococcus (GBS) neonatal disease continues to be a significant cause of concern worldwide. Strains belonging to clonal complex 17 (CC17) are the most frequently responsible for GBS infections in neonates, especially among late-onset disease cases. Therefore, we undertook the largest genomic study of GBS CC17 strains to date to decipher the genetic bases of their remarkable colonization and infection ability. We show that crucial functions involved in different steps of the colonization or infection process of GBS are distinctly mutated during the adaptation of CC17 to the human host. In particular, our results implicate the CovRS two-component regulator of virulence in the differentiation between carriage- and disease-associated isolates. Not only does this work raise important implications for the ongoing development of a vaccine against GBS but might also drive the discovery of key functions for GBS adaptation and pathogenesis that have been overlooked until now.

Author Video: An author video summary of this article is available.

Identifying copy number variation of the dominant virulence factors msa and p22 within genomes of the fish pathogen Renibacterium salmoninarum

Renibacterium salmoninarum is the causative agent of bacterial kidney disease, an important disease of farmed and wild salmonid fish worldwide. Despite the wide spatiotemporal distribution of this disease and habitat pressures ranging from the natural environment to aquaculture and rivers to marine environments, little variation has been observed in the R. salmoninarum genome. Here we use the coverage depth from genomic sequencing corroborated by real-time quantitative PCR to detect copy number variation (CNV) among the genes of R. salmoninarum. CNV was primarily limited to the known dominant virulence factors msa and p22. Among 68 isolates representing the UK, Norway and North America, the msa gene ranged from two to five identical copies and the p22 gene ranged from one to five copies. CNV for these two genes co-occurred, suggesting they may be functionally linked. Isolates carrying CNV were phylogenetically restricted and originated predominantly from sites in North America, rather than the UK or Norway. Although both phylogenetic relationship and geographical origin were found to correlate with CNV status, geographical origin was a much stronger predictor than phylogeny, suggesting a role for local selection pressures in the repeated emergence and maintenance of this trait.

16Stimator: statistical estimation of ribosomal gene copy numbers from draft genome assemblies

The 16S rRNA gene (16S) is an accepted marker of bacterial taxonomic diversity, even though differences in copy number obscure the relationship between amplicon and organismal abundances. Ancestral state reconstruction methods can predict 16S copy numbers through comparisons with closely related reference genomes; however, the database of closed genomes is limited. Here, we extend the reference database of 16S copy numbers to de novo assembled draft genomes by developing 16Stimator, a method to estimate 16S copy numbers when these repetitive regions collapse during assembly. Using a read depth approach, we estimate 16S copy numbers for 12 endophytic isolates from Arabidopsis thaliana and confirm estimates by qPCR. We further apply this approach to draft genomes deposited in NCBI and demonstrate accurate copy number estimation regardless of sequencing platform, with an overall median deviation of 14%. The expanded database of isolates with 16S copy number estimates increases the power of phylogenetic correction methods for determining organismal abundances from 16S amplicon surveys.