- Volume 7, Issue 11, 2021

Staphylococcus aureus is a major bacterial pathogen in humans, and a dominant cause of severe bloodstream infections. Globally, antimicrobial resistance (AMR) in S. aureus remains challenging. While human risk factors for infection have been defined, contradictory evidence exists for the role of bacterial genomic variation in S. aureus disease. To investigate the contribution of bacterial lineage and genomic variation to the development of bloodstream infection, we undertook a genome-wide association study comparing bacteria from 1017 individuals with bacteraemia to 984 adults with asymptomatic S. aureus nasal carriage. Within 984 carriage isolates, we also compared healthcare-associated (HA) carriage with community-associated (CA) carriage. All major global lineages were represented in both bacteraemia and carriage, with no evidence for different infection rates. However, kmers tagging trimethoprim resistance-conferring mutation F99Y in dfrB were significantly associated with bacteraemia-vs-carriage (P=10-8.9-10-9.3). Pooling variation within genes, bacteraemia-vs-carriage was associated with the presence of mecA (HMP=10-5.3) as well as the presence of SCCmec (HMP=10-4.4). Among S. aureus carriers, no lineages were associated with HA-vs-CA carriage. However, we found a novel signal of HA-vs-CA carriage in the foldase protein prsA, where kmers representing conserved sequence allele were associated with CA carriage (P=10-7.1-10-19.4), while in gyrA, a ciprofloxacin resistance-conferring mutation, L84S, was associated with HA carriage (P=10-7.2). In an extensive study of S. aureus bacteraemia and nasal carriage in the UK, we found strong evidence that all S. aureus lineages are equally capable of causing bloodstream infection, and of being carried in the healthcare environment. Genomic variation in the foldase protein prsA is a novel genomic marker of healthcare origin in S. aureus but was not associated with bacteraemia. AMR determinants were associated with both bacteraemia and healthcare-associated carriage, suggesting that AMR increases the propensity not only to survive in healthcare environments, but also to cause invasive disease.

The increasing use of culture independent diagnostic testing for the diagnosis of   Neisseria gonorrhoeae   infection has led to gaps in surveillance of antimicrobial resistance (AMR) rates due to limited availability of cultures. Our study reports the findings of a second national survey of   N. gonorrhoeae   in New Zealand, utilizing whole-genome sequencing (WGS) to study the population structure, prevalence of AMR, epidemiology and transmission of gonorrhoea isolates. We analysed 314 isolates and found a strong correlation between carriage of acquired resistance genes or chromosomal point mutations and phenotypic susceptibility testing results. Overall, the New Zealand rates of azithromycin resistance and decreased susceptibility to ceftriaxone remain lower than in most countries, which are part of the World Health Organization (WHO) Global Gonococcal Antimicrobial Surveillance Programme (GASP). The phylogeny provides evidence of a diverse population significantly associated with sexual behaviour groups. Transmission clustering with a ten single nucleotide polymorphism (SNP) cut-off identified 49 clusters, of which ten were solely associated with men who have sex with men (MSM), whereas remaining clusters included heterosexual patients, as well as MSM, suggesting that bridging of sexual networks is occurring. Utilizing pairwise SNP differences between isolates of the same sequence types we determined genetic variation for the three typing schemes used in this study [Multi locus sequence typing (MLST), multi-antigen sequence typing (NG-MAST), and sequence typing for antimicrobial resistance (NG-STAR)]. A median of 0.0 to 52.5 pairwise SNP differences within a single NG-STAR sequence type underlines previous findings of the superiority of the NG-STAR typing scheme in terms of genomic inherency. With our analysis incorporating epidemiological and genomic data, we were able to show a comprehensive overview of the   N. gonorrhoeae   population circulating in New Zealand, focussing on AMR and transmission within sexual networks. Regular surveillance studies to understand the origin, evolution and spread of AMR for gonorrhoea remain necessary to make informed decisions about public health guidelines, as the internationally rising rates of ceftriaxone and azithromycin resistance have already led to adaptation of current treatment guidelines in the UK and the USA, highlighting the importance of regular surveillance in individual countries.

The medically relevant Trichophyton rubrum species complex has a variety of phenotypic presentations but shows relatively little genetic differences. Conventional barcodes, such as the internal transcribed spacer (ITS) region or the beta-tubulin gene, are not able to completely resolve the relationships between these closely related taxa. T. rubrum, T. soudanense and T. violaceum are currently accepted as separate species. However, the status of certain variants, including the T. rubrum morphotypes megninii and kuryangei and the T. violaceum morphotype yaoundei, remains to be deciphered. We conducted the first phylogenomic analysis of the T. rubrum species complex by studying 3105 core genes of 18 new strains from the BCCM/IHEM culture collection and nine publicly available genomes. Our analyses revealed a highly resolved phylogenomic tree with six separate clades. Trichophyton rubrum, T. violaceum and T. soudanense were confirmed in their status of species. The morphotypes T. megninii, T. kuryangei and T. yaoundei all grouped in their own respective clade with high support, suggesting that these morphotypes should be reinstituted to the species-level. Robinson-Foulds distance analyses showed that a combination of two markers (a ubiquitin-protein transferase and a MYB DNA-binding domain-containing protein) can mirror the phylogeny obtained using genomic data, and thus represent potential new markers to accurately distinguish the species belonging to the T. rubrum complex.

Mycobacterium abscessus comprises three subspecies: M. abscessus subsp. abscessus , M. abscessus subsp. bolletii , and M. abscessus subsp. massiliense . These closely related strains are typically multi-drug-resistant and can cause difficult-to-treat infections. Dominant clusters of isolates with increased pathogenic potential have been demonstrated in pulmonary infections in the global cystic fibrosis (CF) population. An investigation was performed on isolates cultured from an Asian, predominantly non-CF population to explore the phylogenomic relationships within our population and compare it to global M. abscessus isolates. Whole-genome-sequencing was performed on M. abscessus isolates between 2017 and 2019. Bioinformatic analysis was performed to determine multi-locus-sequence-type, to establish the phylogenetic relationships between isolates, and to identify virulence and resistance determinants in these isolates. A total of 210 isolates were included, of which 68.5 % (144/210) were respiratory samples. These isolates consisted of 140 (66.6 %) M . abscessus subsp. massiliense , 67 (31.9 %) M . abscessus subsp. abscessus, and three (1.4 %) M . abscessus subsp. bolletii . Dominant sequence-types in our population were similar to those of global CF isolates, but SNP differences in our population were comparatively wider despite the isolates being from the same geographical region. ESX (ESAT-6 secretory) cluster three appeared to occur most commonly in ST4 and ST6 M. abscessus subsp. massiliense , but other virulence factors did not demonstrate an association with isolate subspecies or sample source. We demonstrate that although similar predominant sequence-types are seen in our patient population, cross-transmission is absent. The risk of patient-to-patient transmission appears to be largely limited to the vulnerable CF population, indicating infection from environmental sources remains more common than human-to-human transmission. Resistance and virulence factors are largely consistent across the subspecies with the exception of clarithromycin susceptibility and ESX-3.

Microbial genomes are highly adaptable, with mobile genetic elements (MGEs) such as integrative conjugative elements (ICEs) mediating the dissemination of new genetic information throughout bacterial populations. This is countered by defence mechanisms such as CRISPR-Cas systems, which limit invading MGEs by sequence-specific targeting. Here we report the distribution of the pVir, pTet and PCC42 plasmids and a new 70–129 kb ICE (CampyICE1) in the foodborne bacterial pathogens   Campylobacter jejuni   and   Campylobacter coli  . CampyICE1 contains a degenerated Type II-C CRISPR system consisting of a sole Cas9 protein, which is distinct from the previously described Cas9 proteins from   C. jejuni   and   C. coli  . CampyICE1 is conserved in structure and gene order, containing blocks of genes predicted to be involved in recombination, regulation and conjugation. CampyICE1 was detected in 134/5829 (2.3 %)   C  .   jejuni   genomes and 92/1347 (6.8 %)   C  .   coli   genomes. Similar ICEs were detected in a number of non-jejuni/coli   Campylobacter   species, although these lacked a CRISPR-Cas system. CampyICE1 carries three separate short CRISPR spacer arrays containing a combination of 108 unique spacers and 16 spacer-variant families. A total of 69 spacers and 10 spacer-variant families (63.7 %) were predicted to target   Campylobacter   plasmids. The presence of a functional CampyICE1 Cas9 protein and matching anti-plasmid spacers was associated with the absence of the pVir, pTet and pCC42 plasmids (188/214 genomes, 87.9 %), suggesting that the CampyICE1-encoded CRISPR-Cas has contributed to the exclusion of competing plasmids. In conclusion, the characteristics of the CRISPR-Cas9 system on CampyICE1 suggests a history of plasmid warfare in   Campylobacter  .

Understanding the evolution of the SARS-CoV-2 virus in various regions of the world during the Covid-19 pandemic is essential to help mitigate the effects of this devastating disease. We describe the phylogenomic and population genetic patterns of the virus in Mexico during the pre-vaccination stage, including asymptomatic carriers. A real-time quantitative PCR screening and phylogenomic reconstructions directed at sequence/structure analysis of the spike glycoprotein revealed mutation of concern E484K in genomes from central Mexico, in addition to the nationwide prevalence of the imported variant 20C/S:452R (B.1.427/9). Overall, the detected variants in Mexico show spike protein mutations in the N-terminal domain (i.e. R190M), in the receptor-binding motif (i.e. T478K, E484K), within the S1–S2 subdomains (i.e. P681R/H, T732A), and at the basis of the protein, V1176F, raising concerns about the lack of phenotypic and clinical data available for the variants of interest we postulate: 20B/478K.V1 (B.1.1.222 or B.1.1.519) and 20B/P.4 (B.1.1.28.4). Moreover, the population patterns of single nucleotide variants from symptomatic and asymptomatic carriers obtained with a self-sampling scheme confirmed the presence of several fixed variants, and differences in allelic frequencies among localities. We identified the mutation N:S194L of the nucleocapsid protein associated with symptomatic patients. Phylogenetically, this mutation is frequent in Mexican sub-clades. Our results highlight the dual and complementary role of spike and nucleocapsid proteins in adaptive evolution of SARS-CoV-2 to their hosts and provide a baseline for specific follow-up of mutations of concern during the vaccination stage.

Here, we characterize an uncommon set of telomeres from   Streptomyces rimosus   ATCC 10970, the parental strain of a lineage of one of the earliest-discovered antibiotic producers. Following the closure of its genome sequence, we compared unusual telomeres from this organism with the other five classes of replicon ends found amongst streptomycetes. Closed replicons of streptomycete chromosomes were organized with respect to their phylogeny and physical orientation, which demonstrated that different telomeres were not associated with particular clades and are likely shared amongst different strains by plasmid-driven horizontal gene transfer. Furthermore, we identified a ~50 kb origin island with conserved synteny that is located at the core of all streptomycete chromosomes and forms an axis around which symmetrical chromosome inversions can take place. Despite this chromosomal bilateral symmetry, a bias in parS sites to the right of oriC is maintained across the family   Streptomycetaceae   and suggests that the formation of ParB/parS nucleoprotein complexes on the right replichore is a conserved feature in streptomycetes. Consequently, our studies reveal novel features of linear bacterial replicons that, through their manipulation, may lead to improvements in growth and productivity of this important industrial group of bacteria.

We report here the complete genome sequence of the widely studied Actinobacillus pleuropneumoniae serovar 8 reference strain 405, generated using the Pacific Biosciences (PacBio) RS II platform. Furthermore, we compared draft sequences generated by Illumina sequencing of six stocks of this strain, including the same original stock used to generate the PacBio sequence, held in different countries and found little genetic variation, with only three SNPs identified, all within the degS gene. However, sequences of two small plasmids, pARD3079 and p405tetH, detected by Illumina sequencing of the draft genomes were not identified in the PacBio sequence of the reference strain.

A special-care neonatal unit from a large public hospital in Malawi was noted as having more frequent, difficult-to-treat infections, and a suspected outbreak of multi-drug-resistant Klebsiella pneumoniae was investigated using genomic characterisation. All K. pneumoniae bloodstream infections (BSIs) from patients in the neonatal ward (n=62), and a subset of K. pneumoniae BSI isolates (n=38) from other paediatric wards in the hospital, collected over a 4 year period were studied. After whole genome sequencing, the strain sequence types (STs), plasmid types, virulence and resistance genes were identified. One ST340 clone, part of clonal complex 258 (CC258) and an ST that drives hospital outbreaks worldwide, harbouring numerous resistance genes and plasmids, was implicated as the likely cause of the outbreak. This study contributes molecular information necessary for tracking and characterizing this important hospital pathogen in sub-Saharan Africa.

Genomic epidemiology is a tool for tracing transmission of pathogens based on whole-genome sequencing. We introduce the mGEMS pipeline for genomic epidemiology with plate sweeps representing mixed samples of a target pathogen, opening the possibility to sequence all colonies on selective plates with a single DNA extraction and sequencing step. The pipeline includes the novel mGEMS read binner for probabilistic assignments of sequencing reads, and the scalable pseudoaligner Themisto. We demonstrate the effectiveness of our approach using closely related samples in a nosocomial setting, obtaining results that are comparable to those based on single-colony picks. Our results lend firm support to more widespread consideration of genomic epidemiology with mixed infection samples.