- Volume 7, Issue 2, 2021

Neisseria gonorrhoeae , the bacterium responsible for the sexually transmitted disease gonorrhoea, has shown an extraordinary ability to develop antimicrobial resistance (AMR) to multiple classes of antimicrobials. With no available vaccine, managing N. gonorrhoeae infections demands effective preventive measures, antibiotic treatment and epidemiological surveillance. The latter two are progressively being supported by the generation of whole-genome sequencing (WGS) data on behalf of national and international surveillance programmes. In this context, this study aims to perform N. gonorrhoeae clustering into genogroups based on WGS data, for enhanced prospective laboratory surveillance. Particularly, it aims to identify the major circulating WGS-genogroups in Europe and to establish a relationship between these and AMR. Ultimately, it enriches public databases by contributing with WGS data from Portuguese isolates spanning 15 years of surveillance. A total of 3791 carefully inspected N. gonorrhoeae genomes from isolates collected across Europe were analysed using a gene-by-gene approach (i.e. using cgMLST). Analysis of cluster composition and stability allowed the classification of isolates into a two-step hierarchical genogroup level determined by two allelic distance thresholds revealing cluster stability. Genogroup clustering in general agreed with available N. gonorrhoeae typing methods [i.e. MLST (multilocus sequence typing), NG-MAST ( N. gonorrhoeae multi-antigen sequence typing) and PubMLST core-genome groups], highlighting the predominant genogroups circulating in Europe, and revealed that the vast majority of the genogroups present a dominant AMR profile. Additionally, a non-static gene-by-gene approach combined with a more discriminatory threshold for potential epidemiological linkage enabled us to match data with previous reports on outbreaks or transmission chains. In conclusion, this genogroup assignment allows a comprehensive analysis of N. gonorrhoeae genetic diversity and the identification of the WGS-based genogroups circulating in Europe, while facilitating the assessment (and continuous monitoring) of their frequency, geographical dispersion and potential association with specific AMR signatures. This strategy may benefit public-health actions through the prioritization of genogroups to be controlled, the identification of emerging resistance carriage, and the potential facilitation of data sharing and communication.

There is increased awareness of the worldwide spread of specific epidemic multidrug-resistant (MDR) lineages of the human commensal Staphylococcus epidermidis . Here, using bioinformatic analyses accounting for population structure, we determined genomic traits (genes, SNPs and k-mers) that distinguish S. epidermidis causing prosthetic-joint infections (PJIs) from commensal isolates from nares, by analysing whole-genome sequencing data from S. epidermidis from PJIs prospectively collected over 10 years in Sweden, and contemporary S. epidermidis from the nares of patients scheduled for arthroplasty surgery. Previously suggested virulence determinants and the presence of genes and mutations linked to antimicrobial resistance (AMR) were also investigated. Publicly available S. epidermidis sequences were used for international extrapolation and validation of findings. Our data show that S. epidermidis causing PJIs differed from nasal isolates not by virulence but by traits associated with resistance to compounds used in prevention of PJIs: β-lactams, aminoglycosides and chlorhexidine. Almost a quarter of the PJI isolates did not belong to any of the previously described major nosocomial lineages, but the AMR-related traits were also over-represented in these isolates, as well as in international S. epidermidis isolates originating from PJIs. Genes previously associated with virulence in S. epidermidis were over-represented in individual lineages, but failed to reach statistical significance when adjusted for population structure. Our findings suggest that the current strategies for prevention of PJIs select for nosocomial MDR S. epidermidis lineages that have arisen from horizontal gene transfer of AMR-related traits into multiple genetic backgrounds.

Pre-existing colonization with Staphylococcus aureus or Klebsiella pneumoniae has been found to increase the risk of infection in intensive care patients. We previously conducted a longitudinal study to characterize colonization of these two organisms in patients admitted to intensive care in a hospital in southern Vietnam. Here, using genomic and phylogenetic analyses, we aimed to assess the contribution these colonizing organisms made to infections. We found that in the majority of patients infected with S. aureus or K. pneumoniae , the sequence type of the disease-causing (infecting) isolate was identical to that of corresponding colonizing organisms in the respective patient. Further in-depth analysis revealed that in patients infected by S. aureus ST188 and by K. pneumoniae ST17, ST23, ST25 and ST86, the infecting isolate was closely related to and exhibited limited genetic variation relative to pre-infection colonizing isolates. Multidrug-resistant S. aureus ST188 was identified as the predominant agent of colonization and infection. Colonization and infection by K. pneumoniae were characterized by organisms with limited antimicrobial resistance profiles but extensive repertoires of virulence genes. Our findings augment the understanding of the link between bacterial colonization and infection in a low-resource setting, and could facilitate the development of novel evidence-based approaches to prevent and treat infections in high-risk patients in intensive care.

Salmonella enterica subspecies arizonae is frequently associated with animal reservoirs, particularly reptiles, and can cause illness in some mammals, including humans. Using whole-genome sequencing data, core genome phylogenetic analyses were performed using 112  S . enterica subsp. arizonae isolates, representing 46 of 102 described serovars. Nearly one-third of these are polyphyletic, including two serovars that appear in four and five distinct evolutionary lineages. Subspecies arizonae has a monophasic H antigen. Among the 46 serovars investigated, only 8 phase 1 H antigens were identified, demonstrating high conservation for this antigen. Prophages and plasmids were found throughout this subspecies, including five novel prophages. Polyphyly was also reflected in prophage content, although some clade-specific enrichment for some phages was observed. IncFII(S) was the most frequent plasmid replicon identified and was found in a quarter of S. enterica subsp. arizonae genomes. Salmonella pathogenicity islands (SPIs) 1 and 2 are present across all Salmonella , including this subspecies, although effectors sipA, sptP and arvA in SPI-1 and sseG and ssaI in SPI-2 appear to be lost in this lineage. SPI-20, encoding a type VI secretion system, is exclusive to this subspecies and is well maintained in all genomes sampled. A number of fimbral operons were identified, including the sas operon that appears to be a synapomorphy for this subspecies, while others exhibited more clade-specific patterns. This work reveals evolutionary patterns in S. enterica subsp. arizonae that make this subspecies a unique lineage within this very diverse species.

Streptococcus pneumoniae causes over one million deaths from lower respiratory infections per annum worldwide. Although mortality is very high in Southeast Asian countries, molecular epidemiological information remains unavailable for some countries. In this study, we report, for the first time, the whole-genome sequences and genetic profiles of pneumococcal strains isolated in Myanmar. We isolated 60 streptococcal strains from 300 children with acute respiratory infection at Yangon Children’s Hospital in Myanmar. We obtained whole-genome sequences and identified the species, serotypes, sequence types, antimicrobial resistance (AMR) profiles, virulence factor profiles and pangenome structure using sequencing-based analysis. Average nucleotide identity analysis indicated that 58 strains were S. pneumoniae and the other 2 strains were Streptococcus mitis . The major serotype was 19F (11 strains), followed by 6E (6B genetic variant; 7 strains) and 15 other serotypes; 5 untypable strains were also detected. Multilocus sequence typing analysis revealed 39 different sequence types, including 11 novel ones. In addition, genetic profiling indicated that AMR genes and mutations spread among pneumococcal strains in Myanmar. A minimum inhibitory concentration assay indicated that several pneumococcal strains had acquired azithromycin and tetracycline resistance, whereas no strains were found to be resistant against levofloxacin and high-dose penicillin G. Phylogenetic and pangenome analysis showed various pneumococcal lineages and that the pneumococcal strains contain a rich and mobile gene pool, providing them with the ability to adapt to selective pressures. This molecular epidemiological information can help in tracking global infection and supporting AMR control in addition to public health interventions in Myanmar.

Mycobacterium microti is an animal-adapted member of the Mycobacterium tuberculosis complex (MTBC), which was originally isolated from voles, but has more recently also been isolated from other selected mammalian hosts, including occasionally from humans. Here, we have generated and analysed the complete genome sequences of five representative vole and clinical M. microti isolates using PacBio- and Illumina-based technologies, and have tested their virulence and vaccine potential in SCID (severe combined immune deficient) mouse and/or guinea pig infection models. We show that the clinical isolates studied here cluster separately in the phylogenetic tree from vole isolates and other clades from publicly available M. microti genome sequences. These data also confirm that the vole and clinical M. microti isolates were all lacking the specific RD1mic region, which in other tubercle bacilli encodes the ESX-1 type VII secretion system. Biochemical analysis further revealed marked phenotypic differences between isolates in type VII-mediated secretion of selected PE and PPE proteins, which in part were attributed to specific genetic polymorphisms. Infection experiments in the highly susceptible SCID mouse model showed that the clinical isolates were significantly more virulent than the tested vole isolates, but still much less virulent than the M. tuberculosis H37Rv control strain. The strong attenuation of the ATCC 35872 vole isolate in immunocompromised mice, even compared to the attenuated BCG (bacillus Calmette–Guérin) vaccine, and its historic use in human vaccine trials encouraged us to test this strain’s vaccine potential in a guinea pig model, where it demonstrated similar protective efficacy as a BCG control, making it a strong candidate for vaccination of immunocompromised individuals in whom BCG vaccination is contra-indicated. Overall, we provide new insights into the genomic and phenotypic variabilities and particularities of members of an understudied clade of the MTBC, which all share a recent common ancestor that is characterized by the deletion of the RD1mic region.

Poly(A) polymerases (PAPs) and tRNA nucleotidyltransferases belong to a superfamily of nucleotidyltransferases and modify RNA 3′-ends. The product of the pcnB gene, PAP I, has been characterized in a few β-, γ- and δ- Proteobacteria . Using the PAP I signature sequence, putative PAPs were identified in bacterial species from the α- and ε- Proteobacteria and from four other bacterial phyla ( Firmicutes , Actinobacteria , Bacteroidetes and Aquificae ). Phylogenetic analysis, alien index and G+C content calculations strongly suggest that the PAPs in the species identified in this study arose by horizontal gene transfer from the β- and γ- Proteobacteria .

Burkholderia pseudomallei , a soil-dwelling Gram-negative bacterium, is the causative agent of the endemic tropical disease melioidosis. Clinical manifestations of B. pseudomallei infection range from acute or chronic localized infection in a single organ to fulminant septicaemia in multiple organs. The diverse clinical manifestations are attributed to various factors, including the genome plasticity across B. pseudomallei strains. We previously characterized B. pseudomallei strains isolated in Malaysia and noted different levels of virulence in model hosts. We hypothesized that the difference in virulence might be a result of variance at the genome level. In this study, we sequenced and assembled four Malaysian clinical B. pseudomallei isolates, UKMR15, UKMPMC2000, UKMD286 and UKMH10. Phylogenomic analysis showed that Malaysian subclades emerged from the Asian subclade, suggesting that the Malaysian strains originated from the Asian region. Interestingly, the low-virulence strain, UKMH10, was the most distantly related compared to the other Malaysian isolates. Genomic island (GI) prediction analysis identified a new island of 23 kb, GI9c, which is present in B. pseudomallei and Burkholderia mallei , but not Burkholderia thailandensis . Genes encoding known B. pseudomallei virulence factors were present across all four genomes, but comparative analysis of the total gene content across the Malaysian strains identified 104 genes that are absent in UKMH10. We propose that these genes may encode novel virulence factors, which may explain the reduced virulence of this strain. Further investigation on the identity and role of these 104 proteins may aid in understanding B. pseudomallei pathogenicity to guide the design of new therapeutics for treating melioidosis.