Search

The adaptive immune system protects against infection via selection of specific antigen receptors on B-cells and T-cells. We studied the prevalent CD8 ‘killer’ T-cell response mounted against SARS-CoV-2 Spike269-277 epitope YLQPRTFLL via the most frequent Human Leukocyte Antigen (HLA) class I worldwide, HLA A*02. The widespread Spike P272L mutation has arisen in at least 14 different SARS-CoV-2 lineages to date, including in lineages identified as variants of concern. P272L was common in the B.1.177 lineage associated with establishing the ‘second wave’ in Europe. The large CD8 T-cell response seen across a cohort of HLA A*02+ convalescent patients, comprising of over 120 different TCRs, failed to respond to the P272L. Sizable populations (0.01%-0.2%) of total CD8 T-cells from individuals vaccinated against SARS-CoV-2 stained with HLA A*02-YLQPRTFLL multimers but failed to bind to the P272L reagent. Viral escape at prevalent T-cell epitopes restricted by high frequency HLA may be particularly problematic when vaccine immunity is focussed on a single protein such as SARS-CoV-2 Spike and provides a strong argument for inclusion of multiple viral proteins in next generation vaccines and highlights the urgent need for monitoring T-cell escape in new SARS-CoV-2 variants.New in Version 2Updated as P272L now seen in 14 different SARS-CoV-2 lineages including the B.1.1.7/Alpha (UK variant) lineageAtomic structures of HLA A*02-YLQPRTFLL and HLA A*02-YLQLRTFLL addedPierre Rizkallah added as an author and author list changed to reflect the contributions of Aaron Wall and Anna Fuller to the newly added datasets

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections were rising during early summer 2021 in many countries as a result of the Delta variant. We assessed reverse transcription polymerase chain reaction swab positivity in the Real-time Assessment of Community Transmission–1 (REACT-1) study in England. During June and July 2021, we observed sustained exponential growth with an average doubling time of 25 days, driven by complete replacement of the Alpha variant by Delta and by high prevalence at younger, less-vaccinated ages. Prevalence among unvaccinated people [1.21% (95% credible interval 1.03%, 1.41%)] was three times that among double-vaccinated people [0.40% (95% credible interval 0.34%, 0.48%)]. However, after adjusting for age and other variables, vaccine effectiveness for double-vaccinated people was estimated at between ~50% and ~60% during this period in England. Increased social mixing in the presence of Delta had the potential to generate sustained growth in infections, even at high levels of vaccination.

In 1998 a nomenclature for the growing list of pesticidal proteins from Bacillus thuringiensis (Bt) was derived based solely on protein sequence comparisons. This nomenclature was widely adopted and provided a robust framework for the naming and classification of the proteins. The success of these proteins in integrated pest management schemes prompted an increased effort to find others with improved or more diverse activities. These discovery activities led to the characterization of proteins from a wider range of bacteria and with a variety of different protein folds. Since most of these new proteins were grouped together as Cry proteins it became apparent that the existing nomenclature had limitations in representing the diverse range of proteins that had been identified. This revised nomenclature retains the basic principles of the 1998 version but provides specific mnemonics to represent different structural groups. For the purposes of consistency, the vast majority of the proteins have either retained their name or have a new name that clearly references the previous one. Other pesticidal proteins not previously included in the nomenclature have been incorporated into this version.

We identify amino acid variants within dominant SARS-CoV-2 T-cell epitopes by interrogating global sequence data. Several variants within nucleocapsid and ORF3a epitopes have arisen independently in multiple lineages and result in loss of recognition by epitope-specific T-cells assessed by IFN-γ and cytotoxic killing assays. Complete loss of T-cell responsiveness was seen due to Q213K in the A*01:01-restricted CD8+ ORF3a epitope FTSDYYQLY207-215, due to P13L, P13S and P13T in the B*27:05-restricted CD8+ nucleocapsid epitope QRNAPRITF9-17, and due to T362I and P365S in the A*03:01/A*11:01-restricted CD8+ nucleocapsid epitope KTFPPTEPK361-369. CD8+ T-cell lines unable to recognise variant epitopes have diverse T-cell receptor repertoires. These data demonstrate the potential for T-cell evasion and highlight the need for ongoing surveillance for variants capable of escaping T-cell as well as humoral immunity., Graphical Abstract

We reviewed all genomic epidemiology studies on COVID-19 in long term care facilities (LTCF) that had been published to date. We found that staff and residents were usually infected with identical, or near identical, SARS-CoV-2 genomes. Outbreaks usually involved one predominant cluster, and the same lineages persisted in LTCFs despite infection control measures. Outbreaks were most commonly due to single or few introductions followed by spread rather than a series of seeding events from the community into LTCFs. Sequencing of samples taken consecutively from the same cases showed persistence of the same genome sequence indicating that the sequencing technique was robust over time. When combined with local epidemiology, genomics facilitated likely transmission sources to be better characterised. Transmission between LTCFs was detected in multiple studies. The mortality rate amongst residents was high in all cases, regardless of the lineage. Bioinformatics methods were inadequate in one third of the studies reviewed, and reproducing the analyses was difficult as sequencing data were not available in many cases., BBSRC, Research England E3 Fund, Wellcome Trust, Academy of Medical Sciences and the Health Foundation, MRC, NIHR

Natural products that possess alkyne or polyyne moieties have been\ud isolated from a variety of biological sources and possess a broad a range of bioactivities. In bacteria, the basic biosynthesis of polyynes is known, but their biosynthetic\ud gene cluster (BGC) distribution and evolutionary relationship to alkyne biosynthesis\ud have not been addressed. Through comprehensive genomic and phylogenetic analyses, the distribution of alkyne biosynthesis gene cassettes throughout bacteria was\ud explored, revealing evidence of multiple horizontal gene transfer events. After investigation of the evolutionary connection between alkyne and polyyne biosynthesis, a\ud monophyletic clade was identified that possessed a conserved seven-gene cassette\ud for polyyne biosynthesis that built upon the conserved three-gene cassette for\ud alkyne biosynthesis. Further diversity mapping of the conserved polyyne gene cassette revealed a phylogenetic subclade for an uncharacterized polyyne BGC present\ud in several Pseudomonas species, designated pgn. Pathway mutagenesis and high-resolution analytical chemistry showed the Pseudomonas protegens pgn BGC directed\ud the biosynthesis of a novel polyyne, protegencin. Exploration of the biosynthetic\ud logic behind polyyne production, through BGC mutagenesis and analytical chemistry,\ud highlighted the essentiality of a triad of desaturase proteins and a thioesterase in\ud both the P. protegens pgn and Trinickia caryophylli (formerly Burkholderia caryophylli)\ud caryoynencin pathways. We have unified and expanded knowledge of polyyne diversity and uniquely demonstrated that alkyne and polyyne biosynthetic gene clusters\ud are evolutionarily related and widely distributed within bacteria. The systematic mapping of conserved biosynthetic genes across the available bacterial genomic diversity\ud proved to be a fruitful method for discovering new natural products and better\ud understanding polyyne biosynthesis.\ud IMPORTANCE Natural products bearing alkyne (triple carbon bond) or polyyne (multiple alternating single and triple carbon bonds) moieties exhibit a broad range of important biological activities. Polyyne metabolites have been implicated in important\ud ecological roles such as cepacin mediating biological control of plant pathogens and caryoynencin protecting Lagriinae beetle eggs against pathogenic fungi. After further\ud phylogenetic exploration of polyyne diversity, we identified a novel gene cluster in\ud Pseudomonas bacteria with known biological control abilities and proved it was responsible for synthesizing a new polyyne metabolite, protegencin. The evolutionary\ud analysis of polyyne pathways showed that multiple biosynthetic genes were conserved, and using mutagenesis, their essentiality was demonstrated. Our research\ud provides a foundation for the future modification of polyyne metabolites and has\ud identified a novel polyyne, protegencin, with potential bioactive roles of ecological\ud and agricultural importance.

We report severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike ΔH69/V70 in multiple independent lineages, often occurring after acquisition of receptor binding motif replacements such as N439K and Y453F, known to increase binding affinity to the ACE2 receptor and confer antibody escape. In vitro, we show that, although ΔH69/V70 itself is not an antibody evasion mechanism, it increases infectivity associated with enhanced incorporation of cleaved spike into virions. ΔH69/V70 is able to partially rescue infectivity of spike proteins that have acquired N439K and Y453F escape mutations by increased spike incorporation. In addition, replacement of the H69 and V70 residues in the Alpha variant B.1.1.7 spike (where ΔH69/V70 occurs naturally) impairs spike incorporation and entry efficiency of the B.1.1.7 spike pseudotyped virus. Alpha variant B.1.1.7 spike mediates faster kinetics of cell-cell fusion than wild-type Wuhan-1 D614G, dependent on ΔH69/V70. Therefore, as ΔH69/V70 compensates for immune escape mutations that impair infectivity, continued surveillance for deletions with functional effects is warranted., Graphical abstract, Meng et al. report that the SARS-CoV-2 spike ΔH69/V70 has arisen multiple times. The deletion increases entry efficiency associated with increased cleaved spike in virions and can compensate for loss of infectivity. The B.1.1.7 spike requires ΔH69/V70 for efficient cell entry and cell-cell fusion activity.

SummaryWe present evidence for multiple independent origins of recombinant SARS-CoV-2 viruses sampled from late 2020 and early 2021 in the United Kingdom. Their genomes carry single nucleotide polymorphisms and deletions that are characteristic of the B.1.1.7 variant of concern, but lack the full complement of lineage-defining mutations. Instead, the remainder of their genomes share contiguous genetic variation with non-B.1.1.7 viruses circulating in the same geographic area at the same time as the recombinants. In four instances there was evidence for onward transmission of a recombinant-origin virus, including one transmission cluster of 45 sequenced cases over the course of two months. The inferred genomic locations of recombination breakpoints suggest that every community-transmitted recombinant virus inherited its spike region from a B.1.1.7 parental virus, consistent with a transmission advantage for B.1.1.7’s set of mutations.

The Burkholderia cepacia complex (Bcc) is a closely related group of bacteria, composed of at least 20 different species, the accurate identification of which is essential in the context of infectious diseases. In industry, they can contaminate non-food products, including home and personal care products and cosmetics. The Bcc are problematic contaminants due to their ubiquitous presence and intrinsic antimicrobial resistance, which enables them to occasionally overcome preservation systems in non-sterile products. Burkholderia lata and Burkholderia contaminans are amongst the Bcc bacteria encountered most frequently as industrial contaminants, but their identification is not straightforward. Both species were historically established as a part of a group known collectively as taxon K, based upon analysis of the recA gene and multilocus sequence typing (MLST). Here, we deploy a straightforward genomics-based workflow for accurate Bcc classification using average nucleotide identity (ANI) and core-gene analysis. The workflow was used to examine a panel of 23 Burkholderia taxon K industrial strains, which, based on MLST, comprised 13 B. lata, 4 B. contaminans and 6 unclassified Bcc strains. Our genomic identification showed that the B. contaminans strains retained their classification, whilst the remaining strains were reclassified as Burkholderia aenigmatica sp. nov. Incorrect taxonomic identification of industrial contaminants is a problematic issue. Application and testing of our genomic workflow allowed the correct classification of 23 Bcc industrial strains, and also indicated that B. aenigmatica sp. nov. may have greater importance than B. lata as a contaminant species. Our study illustrates how the non-food manufacturing industry can harness whole-genome sequencing to better understand antimicrobial-resistant bacteria affecting their products.

Beneficial microorganisms are widely used in agriculture for control of plant pathogens, but a lack of efficacy and safety information has limited the exploitation of multiple promising biopesticides. We applied phylogeny-led genome mining, metabolite analyses and biological control assays to define the efficacy of Burkholderia ambifaria, a naturally beneficial bacterium with proven biocontrol properties but potential pathogenic risk. A panel of 64 B. ambifaria strains demonstrated significant antimicrobial activity against priority plant pathogens. Genome sequencing, specialized metabolite biosynthetic gene cluster mining and metabolite analysis revealed an armoury of known and unknown pathways within B. ambifaria. The biosynthetic gene cluster responsible for the production of the metabolite cepacin was identified and directly shown to mediate protection of germinating crops against Pythium damping-off disease. B. ambifaria maintained biopesticidal protection and overall fitness in the soil after deletion of its third replicon, a non-essential plasmid associated with virulence in Burkholderia cepacia complex bacteria. Removal of the third replicon reduced B. ambifaria persistence in a murine respiratory infection model. Here, we show that by using interdisciplinary phylogenomic, metabolomic and functional approaches, the mode of action of natural biological control agents related to pathogens can be systematically established to facilitate their future exploitation.

SummaryCurrently the primary method for confirming acute SARS-CoV-2 infection is through the use of molecular assays that target highly conserved regions within the viral genome. Many, if not most of the diagnostic targets currently in use were produced early in the pandemic, using genomes sequenced and shared in early 2020. As viral diversity increases, mutations may arise in diagnostic target sites that have an impact on the performance of diagnostic tests. Here, we report on a local outbreak of SARS-CoV-2 which had gained an additional mutation at position 28890 of the nucleocapsid protein, on a background of pre-existing mutations at positions 28881, 28882, 28883 in one of the main circulating viral lineages in Wales at that time. The impact of this additional mutation had a statistically significant impact on the Ct value reported for the N gene target designed by the Chinese CDC and used in a number of commercial diagnostic products. Further investigation identified that, in viral genomes sequenced from Wales over the summer of 2020, the N gene had a higher rate of mutations in diagnostic target sites than other targets, with 115 issues identified affecting over 10% of all cases sequenced between February and the end of August 2020. In comparison an issue was identified for ORFab, the next most affected target, in less than 1.4% of cases over the same time period. This work emphasises the potential impact that mutations in diagnostic target sites can have on tracking local outbreaks, as well as demonstrating the value of genomics as a routine tool for identifying and explaining potential diagnostic primer issues as part of a laboratory quality management system. This work also indicates that with increasing genomic sequencing data availability, there is a need to re-evaluate the diagnostic targets that are in use for SARS-CoV-2 testing, to better target regions that are now demonstrated to be of lower variability.

The SARS-CoV-2 virus, the causative agent of COVID-19, is undergoing constant mutation. Here, we utilized an integrative approach combining epidemiology, virus genome sequencing, clinical phenotyping, and experimental validation to locate mutations of clinical importance. We identified 35 recurrent variants, some of which are associated with clinical phenotypes related to severity. One variant, containing a deletion in the Nsp1-coding region (Δ500-532), was found in more than 20% of our sequenced samples and associates with higher RT-PCR cycle thresholds and lower serum IFN-β levels of infected patients. Deletion variants in this locus were found in 37 countries worldwide, and viruses isolated from clinical samples or engineered by reverse genetics with related deletions in Nsp1 also induce lower IFN-β responses in infected Calu-3 cells. Taken together, our virologic surveillance characterizes recurrent genetic diversity and identified mutations in Nsp1 of biological and clinical importance, which collectively may aid molecular diagnostics and drug design., Graphical Abstract, Bringing together epidemiological, genetic, and clinical data, Lin et al. identify viral mutations that associate with clinical phenotypes. The 500-532 locus in the Nsp1 coding region is a deletion hotspot in the SARS-CoV-2 genome, and deletion variants lead to lower IFN-I response. Prevalence of 500-532 deletion variants is accumulating worldwide.

Natural products that possess alkyne or polyyne moieties have been isolated from a variety of biological sources. In bacteria their biosynthesis has been defined, however, the distribution of polyyne biosynthetic gene clusters (BGCs), and their evolutionary relationship to alkyne biosynthesis, have not been addressed. We explored the distribution of alkyne biosynthesis gene cassettes throughout bacteria, revealing evidence of multiple horizontal gene transfer events. Investigating the evolutionary connection between alkyne and polyyne biosynthesis identified a monophyletic clade possessing a conserved seven-gene cassette for polyyne biosynthesis. Mapping the diversity of these conserved genes revealed a phylogenetic clade representing a polyyne BGC in Pseudomonas, pgn, and subsequent pathway mutagenesis and analytical chemistry characterised the associated metabolite, protegencin. In addition to unifying and expanding our knowledge of polyyne diversity, our results show that alkyne and polyyne biosynthetic gene clusters are promiscuous within bacteria. Systematic mapping of conserved biosynthetic genes across bacterial genomic diversity has proven to be a successful method for discovering natural products.

Burkholderia gladioli is a bacterium with a broad ecology spanning disease in humans, animals and plants, but also encompassing multiple beneficial interactions. It is a plant pathogen, a toxin-producing food-poisoning agent, and causes lung infections in people with cystic fibrosis (CF). Contrasting beneficial traits include antifungal production exploited by insects to protect their eggs, plant protective abilities and antibiotic biosynthesis. We explored the genomic diversity and specialized metabolic potential of 206 B. gladioli strains, phylogenomically defining 5 clades. Historical disease pathovars (pv.) B. gladioli pv. allicola and B. gladioli pv. cocovenenans were distinct, while B. gladioli pv. gladioli and B. gladioli pv. agaricicola were indistinguishable; soft-rot disease and CF infection were conserved across all pathovars. Biosynthetic gene clusters (BGCs) for toxoflavin, caryoynencin and enacyloxin were dispersed across B. gladioli , but bongkrekic acid and gladiolin production were clade-specific. Strikingly, 13 % of CF infection strains characterized were bongkrekic acid-positive, uniquely linking this food-poisoning toxin to this aspect of B. gladioli disease. Mapping the population biology and metabolite production of B. gladioli has shed light on its diverse ecology, and by demonstrating that the antibiotic trimethoprim suppresses bongkrekic acid production, a potential therapeutic strategy to minimize poisoning risk in CF has been identified.

Global dispersal and increasing frequency of the SARS-CoV-2 spike protein variant D614G are suggestive of a selective advantage but may also be due to a random founder effect. We investigate the hypothesis for positive selection of spike D614G in the United Kingdom using more than 25,000 whole genome SARS-CoV-2 sequences. Despite the availability of a large dataset, well represented by both spike 614 variants, not all approaches showed a conclusive signal of positive selection. Population genetic analysis indicates that 614G increases in frequency relative to 614D in a manner consistent with a selective advantage. We do not find any indication that patients infected with the spike 614G variant have higher COVID-19 mortality or clinical severity, but 614G is associated with higher viral load and younger age of patients. Significant differences in growth and size of 614G phylogenetic clusters indicate a need for continued study of this variant., Graphical Abstract, Highlights • Increasing frequency of SARS-CoV-2 D614G is consistent with a selective advantage • Phylodynamic analyses do not show significantly different growth of D614G clusters • There is no association of D614G replacement with greater severity of infection • The D614G replacement is associated with higher viral loads and younger patient age, Analysis of the spread and frequency of SARS-CoV-2 D614G in the United Kingdom suggests a selective advantage for this strain that is associated with higher viral loads in younger patients but not higher COVID-19 clinical severity or mortality.

The genomes of 450 members of Burkholderiaceae , isolated from clinical and environmental sources, were sequenced and assembled as a resource for genome mining. Genomic analysis of the collection has enabled the identification of multiple metabolites and their biosynthetic gene clusters, including the antibiotics gladiolin, icosalide A, enacyloxin, and cepacin A.

The UK’s COVID-19 epidemic during early 2020 was one of world’s largest and unusually well represented by virus genomic sampling. Here we reveal the fine-scale genetic lineage structure of this epidemic through analysis of 50,887 SARS-CoV-2 genomes, including 26,181 from the UK sampled throughout the country’s first wave of infection. Using large-scale phylogenetic analyses, combined with epidemiological and travel data, we quantify the size, spatio-temporal origins and persistence of genetically-distinct UK transmission lineages. Rapid fluctuations in virus importation rates resulted in >1000 lineages; those introduced prior to national lockdown were larger and more dispersed. Lineage importation and regional lineage diversity declined after lockdown, whilst lineage elimination was size-dependent. We discuss the implications of our genetic perspective on transmission dynamics for COVID-19 epidemiology and control.

Burkholderia gladioli is one of few bacteria with a broad ecology spanning disease in humans, animals, and plants, and encompassing beneficial interactions with multiple eukaryotic hosts. It is a plant pathogen, a bongkrekic acid toxin producing food-poisoning agent, and a lung pathogen in people with cystic fibrosis (CF). Contrasting beneficial traits include antifungal production exploited by insects to protect their eggs, plant protective abilities and antibiotic biosynthesis. We explored the ecological diversity and specialized metabolite biosynthesis of 206 B. gladioli strains, phylogenomically defining 5 evolutionary clades. Historical disease pathovars (pv) B. gladioli pv. allicola and B. gladioli pv. cocovenenans were phylogenetically distinct, while B. gladioli pv. gladioli and B. gladioli pv. agaricicola were indistinguishable. Soft-rot disease and CF infection pathogenicity traits were conserved across all pathovars. Biosynthetic gene clusters for toxoflavin, caryoynencin and enacyloxin were dispersed across B. gladioli, but bongkrekic acid and gladiolin production were clade specific. Strikingly, 13% of CF-infection strains characterised (n=194) were bongkrekic acid toxin positive, uniquely linking this food-poisoning risk factor to chronic lung disease. Toxin production was suppressed by exposing strains to the antibiotic trimethoprim, providing a potential therapeutic strategy to minimise poisoning risk in CF.

In order to expand the limited understanding of the genomics of antimicrobial-resistant industrial bacteria, we report the genome sequence of Pluralibacter gergoviae ECO77, a historical contaminant strain of industrial origin. The multireplicon 6.16-Mbp genome of ECO77 consists of a 5.37-Mbp main chromosome, a megaplasmid (605,666 bp), and a large plasmid (182,007 bp).

Sequencing data from host-associated microbes can often be contaminated by the body of the investigator or research subject. Human DNA is typically removed from microbial reads either by subtractive alignment (dropping all reads that map to the human genome) or using a read classification tool to predict those of human origin, and then discarding them. To inform best practice guidelines, we benchmarked 8 alignment-based and 2 classification-based methods of human read detection using simulated data from 10 clinically prevalent bacteria and 3 viruses, into which contaminating human reads had been added.While the majority of methods successfully detected > 99% of the human reads, they were distinguishable by variance. The most precise methods, with negligible variance, were Bowtie2 and SNAP, both of which misclassified few, if any, bacterial reads (and no viral reads) as human. While correctly detecting a similar number of human reads, methods based on taxonomic classification, such as Kraken2 and Centrifuge, often misclassified bacterial reads as human, the extent of which was species-specific. Among the most sensitive methods of human read detection was BWA, although this also made the greatest number of false positive classifications. Across all methods, the set of human reads not identified as such, although often representing < 0.1% of the total reads, were non-randomly distributed along the human genome with many originating from the repeat-rich sex chromosomes.For viral reads and longer (> 300bp) bacterial reads, the highest performing approaches were classification-based, using Kraken2 or Centrifuge. For shorter (150-300bp) bacterial reads, combining multiple methods of human read detection maximised the recovery of human reads from contaminated short read datasets without being compromised by false positives. The highest-performing approach with shorter bacterial reads was a two-stage classification using Bowtie2 followed by SNAP. Using this approach, we re-examined 11,577 publicly archived bacterial readsets for hitherto undetected human contamination. We were able to extract a sufficient number of reads to call known human SNPs, including those with clinical significance, in 6% of the samples. These results show that phenotypically-distinct human sequence is widespread in publicly-archived (and nominally pure) bacterial datasets.

Objectives Shigella sonnei is a globally important diarrhoeal pathogen tracked through the surveillance network PulseNet Latin America and Caribbean (PNLA&C), which participates in PulseNet International. PNLA&C laboratories use common molecular techniques to track pathogens causing foodborne illness. We aimed to demonstrate the possibility and advantages of transitioning to whole genome sequencing (WGS) for surveillance within existing networks across a continent where S. sonnei is endemic. Methods We applied WGS to representative archive isolates of S. sonnei (n = 323) from laboratories in nine PNLA&C countries to generate a regional phylogenomic reference for S. sonnei and put this in the global context. We used this reference to contextualise 16 S. sonnei from three Argentinian outbreaks, using locally generated sequence data. Assembled genome sequences were used to predict antimicrobial resistance (AMR) phenotypes and identify AMR determinants. Results S. sonnei isolates clustered in five Latin American sublineages in the global phylogeny, with many (46%, 149 of 323) belonging to previously undescribed sublineages. Predicted multidrug resistance was common (77%, 249 of 323), and clinically relevant differences in AMR were found among sublineages. The regional overview showed that Argentinian outbreak isolates belonged to distinct sublineages and had different epidemiologic origins. Conclusions Latin America contains novel genetic diversity of S. sonnei that is relevant on a global scale and commonly exhibits multidrug resistance. Retrospective passive surveillance with WGS has utility for informing treatment, identifying regionally epidemic sublineages and providing a framework for interpretation of prospective, locally sequenced outbreaks., Graphical abstract Image 1

SARS-CoV-2 and the resulting COVID-19 pandemic represents one of the greatest recent threats to human health, wellbeing and economic growth. Wastewater-based epidemiology (WBE) of human viruses can be a useful tool for population-scale monitoring of SARS-CoV-2 prevalence and epidemiology to help prevent further spread of the disease, particularly within urban centres. Here we present a longitudinal analysis (March-July, 2020) of SARS-CoV-2 RNA prevalence in sewage across six major urban centres in the UK (total population equivalent 3 million) by q(RT-)PCR and viral genome sequencing. Our results demonstrate that levels of SARS-CoV-2 RNA generally correlated with the abundance of clinical cases recorded within the community in large urban centres, with a marked decline in SARS-CoV-2 RNA abundance following the implementation of lockdown measures. The strength of this association was weaker in areas with lower confirmed COVID-19 case numbers. Further sequencing analysis of SARS-CoV-2 from wastewater suggested that multiple genetically distinct clusters were co-circulating in the local populations covered by our sample sites, and that the genetic variants observed in wastewater reflected similar SNPs observed in contemporaneous samples from cases tested in clinical diagnostic laboratories. We demonstrate how WBE can be used for both community-level detection and tracking of SARS-CoV-2 and other virus’ prevalence, and can inform public health policy decisions. Although, greater understanding of the factors that affect SARS-CoV-2 RNA concentration in wastewater are needed for the full integration of WBE data into outbreak surveillance. In conclusion, our results lend support to the use of routine WBE for monitoring of SARS-CoV-2 and other human pathogenic viruses circulating in the population and assessment of the effectiveness of disease control measures., Graphical abstract Image, graphical abstract

Pseudomonas aeruginosa is a highly versatile, antibiotic-resistant Gram-negative bacterium known for causing opportunistic infections and contamination of industrial products. Despite extensive genomic analysis of clinical P. aeruginosa strains, no genomes exist for preservative-tolerant industrial strains. A unique collection of 69 industrial isolates was assembled and compared to clinical and environmental strains; 16 genetically distinct industrial strains were subjected to array tube genotyping, multilocus sequence typing and whole-genome sequencing. The industrial strains possessed high preservative tolerance and were dispersed widely across P. aeruginosa as a species, but recurrence of strains from the same lineage within specific industrial products and locations was identified. The industrial P. aeruginosa genomes (mean=7.0 Mb) were significantly larger than those of previously sequenced environmental (mean=6.5 Mb; n=19) and clinical (mean=6.6 Mb; n=66) strains. Complete sequencing of the P. aeruginosa industrial strain RW109, which encoded the largest genome (7.75 Mb), revealed a multireplicon structure including a megaplasmid (555 265 bp) and large plasmid (151 612 bp). The RW109 megaplasmid represented an emerging plasmid family conserved in seven industrial and two clinical P. aeruginosa strains, and associated with extremely stress-resilient phenotypes, including antimicrobial resistance and solvent tolerance. Here, by defining the detailed phylogenomics of P. aeruginosa industrial strains, we show that they uniquely possess multireplicon, megaplasmid-bearing genomes, and significantly greater genomic content worthy of further study.

BackgroundInfluenza viruses are associated with a significant global public health burden. The segmented RNA genome of influenza changes continually due to mutation, and the accumulation of these changes within the antigenic recognition sites of haemagglutinin (HA) and neuraminidase (NA) in turn leads to annual epidemics. Influenza A is also zoonotic, allowing for exchange of segments between human and non-human viruses, resulting in new strains with pandemic potential. These processes necessitate a global surveillance system for influenza monitoring. To this end, whole-genome sequencing (WGS) has begun to emerge as a useful tool. However, due to the diversity and mutability of the influenza genome, and noise in short-read data, bioinformatics processing can present challenges.ResultsConventional mapping approaches can be insufficient when a sub-optimal reference strain is chosen. For short-read datasets simulated from influenza H1N1 HA sequences, read recovery after single-reference mapping was routinely as low as 90% for human-origin influenza sequences, and often lower than 10% for those from avian hosts. To this end, we developed adeBruijn Graph (DBG)-based classifier of influenza WGS datasets: VAPOR. In real data benchmarking using 257 WGS read sets with correspondingde novoassemblies, VAPOR provided classifications for all samples with a mean of >99.8% identity to assembled contigs. This resulted in an increase in the number of mapped reads by 6.8% on average, up to a maximum of 13.3%. Additionally, using simulations, we demonstrate that classification from reads may be applied to detection of reassorted strains.ConclusionsVAPOR has potential to simplify bioinformatics pipelines for surveillance, providing a novel method for detection of influenza strains of human and non-human origin directly from reads, minimization of potential data loss and bias associated with conventional mapping, and allowing visualization of alignments that would otherwise require slowde novoassembly. Whilst with expertise and time these pitfalls can largely be avoided, with pre-classification they are remedied in a single step. Furthermore, our algorithm could be adapted in future to surveillance of other RNA viruses. VAPOR is available athttps://github.com/connor-lab/vapor.

Background Burkholderia bacteria are renowned for the biosynthesis of an extensive repertoire of specialised metabolites encoded by biosynthetic gene clusters (BGCs). Burkholderia ambifaria is a historical biopesticide and produces multiple antimicrobials with activity against nematodes, fungi and bacteria. We evaluated the distribution of antimicrobial activity and BGCs across the species which, despite their biopesticidal potency, is poorly understood. Results Genomic analysis of 64-B. ambifaria whole-genomes revealed considerable diversity. Pathway prediction identified 1272 BGCs across 64 strains that were de-replicated to 38 distinct BGCs. Known BGCs accounted for 13 of the BGCs; 7 of which corresponded to known antimicrobials, and included enacyloxins and bactobolins. The antimicrobial activity of the strains was assessed in a classical overlay assay against animal and plant pathogens. The combined genome mining and in vitro activity screens identified a potent anti-Gram-positive and anti-Oomycete compound, cepacin. A damping-off disease model with the pathogen Pythium ultimum and legume Pisum sativum (garden peas) highlighted cepacin as a key specialised metabolite in crop protection. B. ambifaria strain persistence in a mouse inhalation infection model was reduced by removal of the third replicon, whilst retaining the biopesticidal properties of the B. ambifaria strain. Conclusion This work highlights the diversity and prevalence of specialised metabolites in B. ambifaria, and the importance of the potent antimicrobial cepacin in protection against damping-off disease in agriculture. Construction of a B. ambifaria strain that retained biopesticidal properties yet was attenuated for persistence within a murine infection model provides a potential route towards developing safe biopesticidal agent.

The recent detection of SARS-CoV-2 RNA in feces has led to speculation that it can be transmitted via the fecal-oral/ocular route. This review aims to critically evaluate the incidence of gastrointestinal (GI) symptoms, the quantity and infectivity of SARS-CoV-2 in feces and urine, and whether these pose an infection risk in sanitary settings, sewage networks, wastewater treatment plants, and the wider environment (e.g. rivers, lakes and marine waters). A review of 48 independent studies revealed that severe GI dysfunction is only evident in a small number of COVID-19 cases, with 11 ± 2% exhibiting diarrhea and 12 ± 3% exhibiting vomiting and nausea. In addition to these cases, SARS-CoV-2 RNA can be detected in feces from some asymptomatic, mildly- and pre-symptomatic individuals. Fecal shedding of the virus peaks in the symptomatic period and can persist for several weeks, but with declining abundances in the post-symptomatic phase. SARS-CoV-2 RNA is occasionally detected in urine, but reports in fecal samples are more frequent. The abundance of the virus genetic material in both urine (ca. 102–105 gc/ml) and feces (ca. 102–107 gc/ml) is much lower than in nasopharyngeal fluids (ca. 105–1011 gc/ml). There is strong evidence of multiplication of SARS-CoV-2 in the gut and infectious virus has occasionally been recovered from both urine and stool samples. The level and infectious capability of SARS-CoV-2 in vomit remain unknown. In comparison to enteric viruses transmitted via the fecal-oral route (e.g. norovirus, adenovirus), the likelihood of SARS-CoV-2 being transmitted via feces or urine appears much lower due to the lower relative amounts of virus present in feces/urine. The biggest risk of transmission will occur in clinical and care home settings where secondary handling of people and urine/fecal matter occurs. In addition, while SARS-CoV-2 RNA genetic material can be detected by in wastewater, this signal is greatly reduced by conventional treatment. Our analysis also suggests the likelihood of infection due to contact with sewage-contaminated water (e.g. swimming, surfing, angling) or food (e.g. salads, shellfish) is extremely low or negligible based on very low predicted abundances and limited environmental survival of SARS-CoV-2. These conclusions are corroborated by the fact that tens of million cases of COVID-19 have occurred globally, but exposure to feces or wastewater has never been implicated as a transmission vector., Graphical abstract Unlabelled Image, Highlights • SARS-CoV-2 RNA can be readily detected in feces and occasionally urine. • Severe GI dysfunction only occurs in a small number of cases (11 ± 2%). • Likelihood of SARS-CoV-2 being transmitted via feces appears very low. • Likelihood of infection from sewage-contaminated water or food is extremely low.

Microevolution resulted in considerable genotypic variation., Microevolution associated with emergence and expansion of new epidemic clones of bacterial pathogens holds the key to epidemiologic success. To determine microevolution associated with monophasic Salmonella Typhimurium during an epidemic, we performed comparative whole-genome sequencing and phylogenomic analysis of isolates from the United Kingdom and Italy during 2005–2012. These isolates formed a single clade distinct from recent monophasic epidemic clones previously described from North America and Spain. The UK monophasic epidemic clones showed a novel genomic island encoding resistance to heavy metals and a composite transposon encoding antimicrobial drug resistance genes not present in other Salmonella Typhimurium isolates, which may have contributed to epidemiologic success. A remarkable amount of genotypic variation accumulated during clonal expansion that occurred during the epidemic, including multiple independent acquisitions of a novel prophage carrying the sopE gene and multiple deletion events affecting the phase II flagellin locus. This high level of microevolution may affect antigenicity, pathogenicity, and transmission.

Pseudomonas aeruginosa is a common contaminant associated with product recalls in the home and personal care industry. Preservation systems are used to prevent spoilage and protect consumers, but greater knowledge is needed of preservative resistance mechanisms used by P. aeruginosa contaminants. We aimed to identify genetic pathways associated with preservative exposure by using an industrial P. aeruginosa strain and implementing RNA-Seq to understand gene expression changes in response to industry relevant conditions. The consistent differential expression of five genetic pathways during exposure to multiple industrial growth conditions associated with benzisothiazolone (BIT) and phenoxyethanol (POE) preservatives, and a laundry detergent (LD) formulation, was observed. A MexPQ-OpmE Resistance Nodulation Division efflux pump system was commonly upregulated in response to POE, a combination of BIT and POE, and LD together with BIT. In response to all industry conditions, a putative sialic acid transporter and isoprenoid biosynthesis gnyRDBHAL operon demonstrated consistent upregulation. Two operons phnBA and pqsEDCBA involved in Pseudomonas quinolone signaling production and quorum-sensing were also consistently downregulated during exposure to all the industry conditions. The ability to identify consistently differentially expressed genetic pathways in P. aeruginosa can inform the development of future targeted preservation systems that maintain product safety and minimise resistance development., Exposing an industrial Pseudomonas aeruginosa strain to multiple industry relevant conditions associated with the preservatives benzisothiazolone and phenoxyethanol, and a laundry detergent, results in the consistent differential expression of five genetic pathways.

The emergence of Clostridium difficile as a significant human diarrheal pathogen is associated with the production of highly transmissible spores and the acquisition of antimicrobial resistance genes (ARGs) and virulence factors. Unlike the hospital-associated C. difficile RT027 lineage, the community-associated C. difficile RT078 lineage is isolated from both humans and farm animals; however, the geographical population structure and transmission networks remain unknown. Here, we applied whole-genome phylogenetic analysis of 248 C. difficile RT078 strains from 22 countries. Our results demonstrate limited geographical clustering for C. difficile RT078 and extensive coclustering of human and animal strains, thereby revealing a highly linked intercontinental transmission network between humans and animals. Comparative whole-genome analysis reveals indistinguishable accessory genomes between human and animal strains and a variety of antimicrobial resistance genes in the pangenome of C. difficile RT078. Thus, bidirectional spread of C. difficile RT078 between farm animals and humans may represent an unappreciated route disseminating antimicrobial resistance genes between humans and animals. These results highlight the importance of the “One Health” concept to monitor infectious disease emergence and the dissemination of antimicrobial resistance genes.

The genus Yersinia has been used as a model system to study pathogen evolution. Using whole-genome sequencing of all Yersinia species, we delineate the gene complement of the whole genus and define patterns of virulence evolution. Multiple distinct ecological specializations appear to have split pathogenic strains from environmental, nonpathogenic lineages. This split demonstrates that contrary to hypotheses that all pathogenic Yersinia species share a recent common pathogenic ancestor, they have evolved independently but followed parallel evolutionary paths in acquiring the same virulence determinants as well as becoming progressively more limited metabolically. Shared virulence determinants are limited to the virulence plasmid pYV and the attachment invasion locus ail. These acquisitions, together with genomic variations in metabolic pathways, have resulted in the parallel emergence of related pathogens displaying an increasingly specialized lifestyle with a spectrum of virulence potential, an emerging theme in the evolution of other important human pathogens.

Sourcing Antibiotic Resistance It is widely assumed that antibiotic resistance in farm animals contributes in a major way to antibiotic resistance in humans. Mather et al. (p. 1514 , published online 12 September; see the Perspective by Woolhouse and Ward ) analyzed hundreds of genome sequences from Salmonella isolates collected from both livestock and patients in Scotland between 1990 and 2004. The relative contributions of animal-derived and human-derived sources of infection were quantified and the phylogenetic diversity of resistance profiles was matched with bacterial phylogenies. The results suggest that most human infections are caught from other humans rather than from livestock and that humans harbor a greater diversity of antibiotic resistance.

A highly invasive form of non-typhoidal Salmonella (iNTS) disease has recently been documented in many countries in sub-Saharan Africa. The most common Salmonella enterica serovar causing this disease is Typhimurium (Salmonella Typhimurium). We applied whole-genome sequence-based phylogenetic methods to define the population structure of sub-Saharan African invasive Salmonella Typhimurium isolates and compared these to global Salmonella Typhimurium populations. Notably, the vast majority of sub-Saharan invasive Salmonella Typhimurium isolates fell within two closely related, highly clustered phylogenetic lineages that we estimate emerged independently ∼52 and ∼35 years ago in close temporal association with the current HIV pandemic. Clonal replacement of isolates from lineage I by those from lineage II was potentially influenced by the use of chloramphenicol for the treatment of iNTS disease. Our analysis suggests that iNTS disease is in part an epidemic in sub-Saharan Africa caused by highly related Salmonella Typhimurium lineages that may have occupied new niches associated with a compromised human population and antibiotic treatment.

Summary The Tasmanian devil (Sarcophilus harrisii), the largest marsupial carnivore, is endangered due to a transmissible facial cancer spread by direct transfer of living cancer cells through biting. Here we describe the sequencing, assembly, and annotation of the Tasmanian devil genome and whole-genome sequences for two geographically distant subclones of the cancer. Genomic analysis suggests that the cancer first arose from a female Tasmanian devil and that the clone has subsequently genetically diverged during its spread across Tasmania. The devil cancer genome contains more than 17,000 somatic base substitution mutations and bears the imprint of a distinct mutational process. Genotyping of somatic mutations in 104 geographically and temporally distributed Tasmanian devil tumors reveals the pattern of evolution and spread of this parasitic clonal lineage, with evidence of a selective sweep in one geographical area and persistence of parallel lineages in other populations. PaperClip, Graphical Abstract Highlights ► Whole-genome sequences of the Tasmanian devil and two distant cancer subclones ► The Tasmanian devil cancer lineage originated recently in a female devil ► The devil cancer genome is relatively stable despite ongoing evolution ► Clonal divergence and geographic spread elucidated through patterns of mutation, Whole-genome sequences of the Tasmanian devil and two devil cancer subclones suggest that the cancer first arose from a female devil and that the clone has subsequently genetically diverged during its spread across Tasmania.

Recombining Resistance Homologous recombination is frequent in many bacteria, but few studies have addressed whether subpopulations within a species are more or less likely to undergo this process and whether it has consequences for their evolution. Taking a large data set from the pathogen Streptococcus pneumoniae , Hanage et al. (p. 1454 ) discovered a group of strains characterized by an anomalous sequence of housekeeping genes. This sequence appeared to have been horizontally acquired from other pneumococci and related species and was associated with resistance to all classes of antibiotics for which data are available. Thus, hyper-recombination (in contrast to hypermutation) is important in the evolution and spread of antibiotic resistance and may play a role in determining the emergence of species clusters and the phenotypes associated with them.

Shigella flexneri is the most common cause of bacterial dysentery in low-income countries. Despite this, S. flexneri remains largely unexplored from a genomic standpoint and is still described using a vocabulary based on serotyping reactions developed over half-a-century ago. Here we combine whole genome sequencing with geographical and temporal data to examine the natural history of the species. Our analysis subdivides S. flexneri into seven phylogenetic groups (PGs); each containing two-or-more serotypes and characterised by distinct virulence gene complement and geographic range. Within the S. flexneri PGs we identify geographically restricted sub-lineages that appear to have persistently colonised regions for many decades to over 100 years. Although we found abundant evidence of antimicrobial resistance (AMR) determinant acquisition, our dataset shows no evidence of subsequent intercontinental spread of antimicrobial resistant strains. The pattern of colonisation and AMR gene acquisition suggest that S. flexneri has a distinct life-cycle involving local persistence. DOI: http://dx.doi.org/10.7554/eLife.07335.001, eLife digest Dysentery is a disease in which the intestine becomes inflamed due to infection by bacteria, viruses or other microbes. Of the bacteria that can cause dysentery, bacteria called Shigella are most often responsible. Humans can acquire Shigella through contaminated food or water. Over the last century, improvements to sanitation combined with access to clean drinking water and better food hygiene have decreased the number of cases of dysentery in many countries. However, the disease continues to be common in low-income countries, especially in young children. One species of Shigella bacteria, called S. flexneri, causes far more cases of dysentry than other species of Shigella. Across the world, there are many different strains of S. flexneri, but it is not clear how these strains are related to each other, or how variable the genes that they carry are—known as genetic diversity. Here, Connor, Barker, Baker et al. used a technique called whole genome sequencing to map the evolutionary relationships of over 300 S. flexneri samples collected from around the globe over the past 100 years. This revealed that the bacterial strains can be split into seven groups that each have distinct geographic ranges and combinations of genes that enable the bacteria to infect humans. Many of the strains of bacteria within these groups seem to have colonized, and remained in, quite small geographic areas over long periods of time. This is different to other Shigella species, which appear to have spread between continents far more frequently over much shorter timescales. Connor, Barker, Baker et al.'s findings reveal that S. flexneri is more diverse than other Shigella bacteria, and suggest that the ability of strains to persist in local areas may have contributed to the species' long-term success. These results point towards the importance of the provision of clean water in the fight against S. flexneri, and underline the need for a greater understanding of how disease-causing bacteria colonize and interact with the local environment. DOI: http://dx.doi.org/10.7554/eLife.07335.002

© 2015 Elsevier Ltd. Background: Shigellosis is an acute, severe bacterial colitis that, in high-income countries, is typically associated with travel to high-risk regions (Africa, Asia, and Latin America). Since the 1970s, shigellosis has also been reported as a sexually transmitted infection in men who have sex with men (MSM), in whom transmission is an important component of shigellosis epidemiology in high-income nations. We aimed to use sophisticated subtyping and international sampling to determine factors driving shigellosis emergence in MSM linked to an outbreak in the UK. Methods: We did a large-scale, cross-sectional genomic epidemiological study of shigellosis cases collected from 29 countries between December, 1995, and June 8, 2014. Focusing on an ongoing epidemic in the UK, we collected and whole-genome sequenced clinical isolates of Shigella flexneri serotype 3a from high-risk and low-risk regions, including cases associated with travel and sex between men. We examined relationships between geographical, demographic, and clinical patient data with the isolate antimicrobial susceptibility, genetic data, and inferred evolutionary relationships. Findings: We obtained 331 clinical isolates of S flexneri serotype 3a, including 275 from low-risk regions (44 from individuals who travelled to high-risk regions), 52 from high-risk regions, and four outgroup samples (ie, closely related, but genetically distinct isolates used to determine the root of the phylogenetic tree). We identified a recently emerged lineage of S flexneri 3a that has spread intercontinentally in less than 20 years throughout regions traditionally at low risk for shigellosis via sexual transmission in MSM. The lineage had acquired multiple antimicrobial resistance determinants, and prevailing sublineages were strongly associated with resistance to the macrolide azithromycin. Eight (4%) of 206 isolates from the MSM-associated lineage were obtained from patients who had previously provided an isolate; these serial isolations indicated atypical infection patterns (eg, reinfection). Interpretation: We identified transmission-facilitating behaviours and atypical course(s) of infection as precipitating factors in shigellosis-affected MSM. The intercontinental spread of antimicrobial-resistant shigella through established transmission routes emphasises the need for new approaches to tackle the public health challenge of sexually transmitted infections in MSM. Funding: Wellcome Trust (grant number 098051).

Two lineages of Salmonella enterica serovar Typhimurium (S. Typhimurium) of multi-locus sequence type ST313 have been linked with the emergence of invasive Salmonella disease across sub-Saharan Africa. The expansion of these lineages has a temporal association with the HIV pandemic and antibiotic usage. We analysed the whole genome sequence of 129 ST313 isolates representative of the two lineages and found evidence of lineage-specific genome degradation, with some similarities to that observed in S. Typhi. Individual ST313 S. Typhimurium isolates exhibit a distinct metabolic signature and modified enteropathogenesis in both a murine and cattle model of colitis, compared to S. Typhimurium outside of the ST313 lineages. These data define phenotypes that distinguish ST313 isolates from other S. Typhimurium and may represent adaptation to a distinct pathogenesis and lifestyle linked to an-immuno-compromised human population., Author Summary Salmonella enterica is a diverse species, isolates of which can colonise or infect many different animals, including humans and can cause different disease syndromes. S. enterica can be sub-typed using serology into serovars. Isolates from some serovars, known as generalists, can infect multiple hosts (e.g. S. Typhimurium) and usually cause gastroenteritis. However, other serovars exhibit host adaptation or even restriction. Host-adapted serovars such as S. Dublin show preference for a particular host but can also infect other hosts, while host-restricted serovars are capable of infecting only a single host (e.g. S. Typhi in humans) and frequently cause febrile systemic disease (typhoid). In this study, we use genotypic and phenotypic methods to investigate clinical isolates representative of populations of two recently emerged S. Typhimurium lineages of type ST313 associated with invasive disease in sub-Saharan Africa. Our results identify potential characteristics in these isolates that may be associated with adaptation to invasive disease in humans with a compromised immunity.

Data produced by Biolog Phenotype MicroArrays are longitudinal measurements of cells' respiration on distinct substrates. We introduce a three-step pipeline to analyze phenotypic microarray data with novel procedures for grouping, normalization and effect identification. Grouping and normalization are standard problems in the analysis of phenotype microarrays defined as categorizing bacterial responses into active and non-active, and removing systematic errors from the experimental data, respectively. We expand existing solutions by introducing an important assumption that active and non-active bacteria manifest completely different metabolism and thus should be treated separately. Effect identification, in turn, provides new insights into detecting differing respiration patterns between experimental conditions, e.g. between different combinations of strains and temperatures, as not only the main effects but also their interactions can be evaluated. In the effect identification, the multilevel data are effectively processed by a hierarchical model in the Bayesian framework. The pipeline is tested on a data set of 12 phenotypic plates with bacterium Yersinia enterocolitica. Our pipeline is implemented in R language on the top of opm R package and is freely available for research purposes.

The emergence of new sequencing technologies has facilitated the use of bacterial whole genome alignments for evolutionary studies and outbreak analyses. These datasets, of increasing size, often include examples of multiple different mechanisms of horizontal sequence transfer resulting in substantial alterations to prokaryotic chromosomes. The impact of these processes demands rapid and flexible approaches able to account for recombination when reconstructing isolates’ recent diversification. Gubbins is an iterative algorithm that uses spatial scanning statistics to identify loci containing elevated densities of base substitutions suggestive of horizontal sequence transfer while concurrently constructing a maximum likelihood phylogeny based on the putative point mutations outside these regions of high sequence diversity. Simulations demonstrate the algorithm generates highly accurate reconstructions under realistically parameterized models of bacterial evolution, and achieves convergence in only a few hours on alignments of hundreds of bacterial genome sequences. Gubbins is appropriate for reconstructing the recent evolutionary history of a variety of haploid genotype alignments, as it makes no assumptions about the underlying mechanism of recombination. The software is freely available for download at github.com/sanger-pathogens/Gubbins, implemented in Python and C and supported on Linux and Mac OS X.

The bla NDM-1 gene is associated with extensive drug resistance in Gram-negative bacteria. This probably spread to Enterobacteriaceae from Acinetobacter spp., and we characterized plasmids associated with bla NDM-1 in Acinetobacter spp. to gain insight into their role in this dissemination. Four clinical NDM-1-producing Acinetobacter species strains from India and Pakistan were investigated. A plasmid harboring bla NDM-1 , pNDM-40-1, was characterized by whole-genome sequencing of Acinetobacter bereziniae CHI-40-1 and comparison with related plasmids. The presence of similar plasmids in strains from Pakistan was sought by PCR and sequencing of amplicons. Conjugation frequency was tested and stability of pNDM-40-1 investigated by real-time PCR of isolates passaged with and without antimicrobial selection pressure. A. bereziniae and Acinetobacter haemolyticus strains contained plasmids similar to the pNDM-BJ01-like plasmids identified in Acinetobacter spp. in China. The backbone of pNDM-40-1 was almost identical to that of pNDM-BJ01-like plasmids, but the transposon harboring bla NDM-1 , Tn 125 , contained two short deletions. Escherichia coli and Acinetobacter pittii transconjugants were readily obtained. Transconjugants retained pNDM-40-1 after a 14-day passage experiment, although stability was greater with meropenem selection. Fragments of pNDM-BJ01-like plasmid backbones are found near bla NDM-1 in some genetic contexts from Enterobacteriaceae , suggesting that cross-genus transfer has occurred. pNDM-BJ01-like plasmids have been described in isolates originating from a wide geographical region in southern Asia. In vitro data on plasmid transfer and stability suggest that these plasmids could have contributed to the spread of bla NDM-1 into Enterobacteriaceae .

Many bacterial pathogens are specialized, infecting one or few hosts, and this is often associated with more acute disease presentation. Specific genomes show markers of this specialization, which often reflect a balance between gene acquisition and functional gene loss. Within Salmonella enterica subspecies enterica, a single lineage exists that includes human and animal pathogens adapted to cause infection in different hosts, including S. enterica serovar Enteritidis (multiple hosts), S. Gallinarum (birds), and S. Dublin (cattle). This provides an excellent evolutionary context in which differences between these pathogen genomes can be related to host range. Genome sequences were obtained from ∼ 60 isolates selected to represent the known diversity of this lineage. Examination and comparison of the clades within the phylogeny of this lineage revealed signs of host restriction as well as evolutionary events that mark a path to host generalism. We have identified the nature and order of events for both evolutionary trajectories. The impact of functional gene loss was predicted based upon position within metabolic pathways and confirmed with phenotyping assays. The structure of S. Enteritidis is more complex than previously known, as a second clade of S. Enteritidis was revealed that is distinct from those commonly seen to cause disease in humans or animals, and that is more closely related to S. Gallinarum. Isolates from this second clade were tested in a chick model of infection and exhibited a reduced colonization phenotype, which we postulate represents an intermediate stage in pathogen-host adaptation.

Klebsiella pneumoniae is now recognized as an urgent threat to human health because of the emergence of multidrug-resistant strains associated with hospital outbreaks and hypervirulent strains associated with severe community-acquired infections. K. pneumoniae is ubiquitous in the environment and can colonize and infect both plants and animals. However, little is known about the population structure of K. pneumoniae, so it is difficult to recognize or understand the emergence of clinically important clones within this highly genetically diverse species. Here we present a detailed genomic framework for K. pneumoniae based on whole-genome sequencing of more than 300 human and animal isolates spanning four continents. Our data provide genome-wide support for the splitting of K. pneumoniae into three distinct species, KpI (K. pneumoniae), KpII (K. quasipneumoniae), and KpIII (K. variicola). Further, for K. pneumoniae (KpI), the entity most frequently associated with human infection, we show the existence of >150 deeply branching lineages including numerous multidrug-resistant or hypervirulent clones. We show K. pneumoniae has a large accessory genome approaching 30,000 protein-coding genes, including a number of virulence functions that are significantly associated with invasive community-acquired disease in humans. In our dataset, antimicrobial resistance genes were common among human carriage isolates and hospital-acquired infections, which generally lacked the genes associated with invasive disease. The convergence of virulence and resistance genes potentially could lead to the emergence of untreatable invasive K. pneumoniae infections; our data provide the whole-genome framework against which to track the emergence of such threats.

The genus Yersinia is a large and diverse bacterial genus consisting of human-pathogenic species, a fish-pathogenic species, and a large number of environmental species. Recently, the phylogenetic and population structure of the entire genus was elucidated through the genome sequence data of 241 strains encompassing every known species in the genus. Here we report the mining of this enormous data set to create a multilocus sequence typing-based scheme that can identify Yersinia strains to the species level to a level of resolution equal to that for whole-genome sequencing. Our assay is designed to be able to accurately subtype the important human-pathogenic species Yersinia enterocolitica to whole-genome resolution levels. We also report the validation of the scheme on 386 strains from reference laboratory collections across Europe. We propose that the scheme is an important molecular typing system to allow accurate and reproducible identification of Yersinia isolates to the species level, a process often inconsistent in nonspecialist laboratories. Additionally, our assay is the most phylogenetically informative typing scheme available for Y. enterocolitica .

Enterotoxigenic Escherichia coli (ETEC), a major cause of infectious diarrhea, produce heat-stable and/or heat-labile enterotoxins and at least 25 different colonization factors that target the intestinal mucosa. The genes encoding the enterotoxins and most of the colonization factors are located on plasmids found across diverse E. coli serogroups. Whole-genome sequencing of a representative collection of ETEC isolated between 1980 and 2011 identified globally distributed lineages characterized by distinct colonization factor and enterotoxin profiles. Contrary to current notions, these relatively recently emerged lineages might harbor chromosome and plasmid combinations that optimize fitness and transmissibility. These data have implications for understanding, tracking and possibly preventing ETEC disease.

In the multilocus Fisher-Wright model allele frequencies and the combinations of alleles (STs) fluctuate solely through genetic drift. Using multiple loci prevents recombination from obscuring the relationships between strains, because a change at one locus will not alter the information about relationships between strains present at the others. Computer simulation of bacterial populations is complementary to an analytic approach. Furthermore, simulation of the evolution of bacterial populations provides an alternative means of testing hypotheses: the effects of various evolutionary scenarios may be explored by introducing them into a simulation based on the neutral model. An example of this type of approach is discussed in this chapter, where the conditions under which an evolving bacterial population splits into separate clusters, mimicking speciation, are explored. In sum, while coalescent analyses may be the most appropriate for individual loci, and certainly have a much wider body of theoretical work to draw on than the methods discussed in the chapter, the authors are excited by the potential of methods such as multilocus sequence typing (MLST) that consider the genetics of populations of strains rather than of individual loci, to demonstrate forces other than neutrality influencing the structure of populations. The chapter also talks about the modeling of alleles and their association using recently developed techniques designed to simulate and analyze data sets obtained from MLST studies.

Salmonella enterica serovar Typhimurium definitive type 2 (DT2) is host restricted to Columba livia (rock or feral pigeon) but is also closely related to S. Typhimurium isolates that circulate in livestock and cause a zoonosis characterized by gastroenteritis in humans. DT2 isolates formed a distinct phylogenetic cluster within S. Typhimurium based on whole-genome-sequence polymorphisms. Comparative genome analysis of DT2 94-213 and S. Typhimurium SL1344, DT104, and D23580 identified few differences in gene content with the exception of variations within prophages. However, DT2 94-213 harbored 22 pseudogenes that were intact in other closely related S. Typhimurium strains. We report a novel in silico approach to identify single amino acid substitutions in proteins that have a high probability of a functional impact. One polymorphism identified using this method, a single-residue deletion in the Tar protein, abrogated chemotaxis to aspartate in vitro. DT2 94-213 also exhibited an altered transcriptional profile in response to culture at 42°C compared to that of SL1344. Such differentially regulated genes included a number involved in flagellum biosynthesis and motility., IMPORTANCE Whereas Salmonella enterica serovar Typhimurium can infect a wide range of animal species, some variants within this serovar exhibit a more limited host range and altered disease potential. Phylogenetic analysis based on whole-genome sequences can identify lineages associated with specific virulence traits, including host adaptation. This study represents one of the first to link pathogen-specific genetic signatures, including coding capacity, genome degradation, and transcriptional responses to host adaptation within a Salmonella serovar. We performed comparative genome analysis of reference and pigeon-adapted definitive type 2 (DT2) S. Typhimurium isolates alongside phenotypic and transcriptome analyses, to identify genetic signatures linked to host adaptation within the DT2 lineage.

Background:Staphylococcus aureus exhibits tropisms to many distinct animal hosts. While spillover events can occur wherever there is an interface between host species, changes in host tropism only occur with the establishment of sustained transmission in the new host species, leading to clonal expansion. Although the genomic variation underpinning adaptation in S. aureus genotypes infecting bovids and poultry has been well characterized the frequency of switches from one host to another remains obscure. We sought to identify sustained switches in host tropism in the S. aureus population, both anthroponotic and zoonotic, and their distribution over the species phylogeny. Methodologies/Results:We have used a sample of 3042 isolates, representing 696 distinct MLST genotypes, from a well-established database (www.mlst.net). Using an empirical parsimony approach (AdaptML) we have investigated the distribution of switches in host association between both human and non-human (henceforth referred to as animal) hosts. We reconstructed a credible description of past events in the form of a phylogenetic tree; the nodes and leaves of which are statistically associated with either human or animal habitats, estimated from extant host-association and the degree of sequence divergence between genotypes. We identified 15 likely historical switching events; 13 anthroponoses and two zoonoses. Importantly, we identified two human-associated clade candidates (CC25 and CC59) that have arisen from animal-associated ancestors; this demonstrates that a human-specific lineage can emerge from an animal host. We also highlight novel rabbit-associated genotypes arising from a human ancestor. Conclusions: S. aureus is an organism with the capacity to switch into and adapt to novel hosts, even after long periods of isolation in a single host species. Based on this evidence, animal-adapted S. aureus lineages exhibiting resistance to antibiotics must be considered a major threat to public health, as they can adapt to the human population.