- Volume 9, Issue 4, 2023

Violacein is a water-insoluble violet pigment produced by various Gram-negative bacteria. The compound and the bacteria that produce it have been gaining attention due to the antimicrobial and proposed antitumour properties of violacein and the possibility that strains producing it may have broad industrial uses. Bacteria that produce violacein have been isolated from diverse environments including fresh and ocean waters, glaciers, tropical soils, trees, fish and the skin of amphibians. We report here the isolation and characterization of six violacein-producing bacterial strains and three non-violacein-producing close relatives, each isolated from either an aquatic environment or moist food materials in northern California, USA. For each isolate, we characterized traditional phenotypes, generated and analysed draft genome sequences, and carried out multiple types of taxonomic, phylogenetic and phylogenomic analyses. Based on these analyses we assign putative identifications to the nine isolates, which include representatives of the genera   Chromobacterium  ,   Aquitalea  ,   Iodobacter  ,   Duganella  ,   Massilia   and   Janthinobacterium  . In addition, we discuss the utility of various metrics for taxonomic assignment in these groups including average nucleotide identity, whole genome phylogenetic analysis and extent of recent homologous recombination using the software program PopCOGenT.

Genomic analyses are widely applied to epidemiological, population genetic and experimental studies of pathogenic fungi. A wide range of methods are employed to carry out these analyses, typically without including controls that gauge the accuracy of variant prediction. The importance of tracking outbreaks at a global scale has raised the urgency of establishing high-accuracy pipelines that generate consistent results between research groups. To evaluate currently employed methods for whole-genome variant detection and elaborate best practices for fungal pathogens, we compared how 14 independent variant calling pipelines performed across 35 Candida auris isolates from 4 distinct clades and evaluated the performance of variant calling, single-nucleotide polymorphism (SNP) counts and phylogenetic inference results. Although these pipelines used different variant callers and filtering criteria, we found high overall agreement of SNPs from each pipeline. This concordance correlated with site quality, as SNPs discovered by a few pipelines tended to show lower mapping quality scores and depth of coverage than those recovered by all pipelines. We observed that the major differences between pipelines were due to variation in read trimming strategies, SNP calling methods and parameters, and downstream filtration criteria. We calculated specificity and sensitivity for each pipeline by aligning three isolates with chromosomal level assemblies and found that the GATK-based pipelines were well balanced between these metrics. Selection of trimming methods had a greater impact on SAMtools-based pipelines than those using GATK. Phylogenetic trees inferred by each pipeline showed high consistency at the clade level, but there was more variability between isolates from a single outbreak, with pipelines that used more stringent cutoffs having lower resolution. This project generated two truth datasets useful for routine benchmarking of C. auris variant calling, a consensus VCF of genotypes discovered by 10 or more pipelines across these 35 diverse isolates and variants for 2 samples identified from whole-genome alignments. This study provides a foundation for evaluating SNP calling pipelines and developing best practices for future fungal genomic studies.

Sesquiterpene cyclases (STC) catalyse the cyclization of the C15 molecule farnesyl diphosphate into a vast variety of mono- or polycyclic hydrocarbons and, for a few enzymes, oxygenated structures, with diverse stereogenic centres. The huge diversity in sesquiterpene skeleton structures in nature is primarily the result of the type of cyclization driven by the STC. Despite the phenomenal impact of fungal sesquiterpenes on the ecology of fungi and their potentials for applications, the fungal sesquiterpenome is largely untapped. The identification of fungal STC is generally based on protein sequence similarity with characterized enzymes. This approach has improved our knowledge on STC in a few fungal species, but it has limited success for the discovery of distant sequences. Besides, the tools based on secondary metabolite biosynthesis gene clusters have shown poor performance for terpene cyclases. Here, we used four sets of sequences of fungal STC that catalyse four types of cyclization, and specific amino acid motives to identify phylogenetically related sequences in the genomes of basidiomycetes fungi from the order Polyporales. We validated that four STC genes newly identified from the genome sequence of Leiotrametes menziesii, each classified in a different phylogenetic clade, catalysed a predicted cyclization of farnesyl diphosphate. We built HMM models and searched STC genes in 656 fungal genomes genomes. We identified 5605 STC genes, which were classified in one of the four clades and had a predicted cyclization mechanism. We noticed that the HMM models were more accurate for the prediction of the type of cyclization catalysed by basidiomycete STC than for ascomycete STC.

Burkholderia pseudomallei , a Gram-negative pathogen, is the causative agent of melioidosis in humans. This bacterium can be isolated from the soil, stagnant and salt-water bodies, and human and animal clinical specimens. While extensive studies have contributed to our understanding of B. pseudomallei pathogenesis, little is known about how a harmless soil bacterium adapts when it shifts to a human host and exhibits its virulence. The bacterium’s large genome encodes an array of factors that support the pathogen’s ability to survive under stressful conditions, including the host’s internal milieu. In this study, we performed comparative transcriptome analysis of B. pseudomallei cultured in human plasma versus soil extract media to provide insights into B. pseudomallei gene expression that governs bacterial adaptation and infectivity in the host. A total of 455 genes were differentially regulated; genes upregulated in B. pseudomallei grown in human plasma are involved in energy metabolism and cellular processes, whilst the downregulated genes mostly include fatty acid and phospholipid metabolism, amino acid biosynthesis and regulatory function proteins. Further analysis identified a significant upregulation of biofilm-related genes in plasma, which was validated using the biofilm-forming assay and scanning electron microscopy. In addition, genes encoding known virulence factors such as capsular polysaccharide and flagella were also overexpressed, suggesting an overall enhancement of B. pseudomallei virulence potential when present in human plasma. This ex vivo gene expression profile provides comprehensive information on B. pseudomallei ’s adaptation when shifted from the environment to the host. The induction of biofilm formation under host conditions may explain the difficulty in treating septic melioidosis.

Complete Type VI Secretion Systems were identified in the genome sequence data of Neisseria subflava isolates sourced from throat swabs of human volunteers. The previous report was the first to describe two complete Type VI Secretion Systems in these isolates, both of which were distinct in terms of their gene organization and sequence homology. Since publication of the first report, Type VI Secretion System subtypes have been identified in Neisseria spp. The characteristics of each type in N. subflava are further investigated here and in the context of the other Neisseria spp., including identification of the lineages containing the different types and subtypes. Type VI Secretion Systems use VgrG for delivery of toxin effector proteins; several copies of vgrG and associated effector / immunity pairs are present in Neisseria spp. Based on sequence similarity between strains and species, these core Type VI Secretion System genes, vgrG, and effector / immunity genes may diversify via horizontal gene transfer, an instrument for gene acquisition and repair in Neisseria spp.

While the world is still recovering from the Covid-19 pandemic, monkeypox virus (MPXV) awaits to cause another global outbreak as a challenge to all of mankind. However, the Covid-19 pandemic has taught us a lesson to speed up the pace of viral genomic research for the implementation of preventive and treatment strategies. One of the important aspects of MPXV that needs immediate insight is its evolutionary lineage based on genomic studies. Utilizing high-quality isolates from the GISAID (Global Initiative on Sharing All Influenza Data) database, primarily sourced from Europe and North America, we employed a SNP-based whole-genome phylogeny method and identified four major clusters among 628 MPXV isolates. Our findings indicate a distinct evolutionary lineage for the first MPXV isolate, and a complex epidemiology and evolution of MPXV strains across various countries. Further analysis of the host–pathogen interaction network revealed key viral proteins, such as E3, SPI-2, K7 and CrmB, that play a significant role in regulating the network and inhibiting the host’s cellular innate immune system. Our structural analysis of proteins E3 and CrmB revealed potential disruption of stability due to certain mutations. While this study identified a large number of mutations within the new outbreak clade, it also reflected that we need to move fast with the genomic analysis of newly detected strains from around the world to develop better prevention and treatment methods.

Streptococcus suis colonizes the upper respiratory tract of healthy pigs at high abundance but can also cause opportunistic respiratory and systemic disease. Disease-associated S. suis reference strains are well studied, but less is known about commensal lineages. It is not known what mechanisms enable some S. suis lineages to cause disease while others persist as commensal colonizers, or to what extent gene expression in disease-associated and commensal lineages diverge. In this study we compared the transcriptomes of 21  S . suis strains grown in active porcine serum and Todd–Hewitt yeast broth. These strains included both commensal and pathogenic strains, including several strains of sequence type (ST) 1, which is responsible for most cases of human disease and is considered to be the most pathogenic S. suis lineage. We sampled the strains during their exponential growth phase and mapped RNA sequencing reads to the corresponding strain genomes. We found that the transcriptomes of pathogenic and commensal strains with large genomic divergence were unexpectedly conserved when grown in active porcine serum, but that regulation and expression of key pathways varied. Notably, we observed strong variation of expression across media of genes involved in capsule production in pathogens, and of the agmatine deiminase system in commensals. ST1 strains displayed large differences in gene expression between the two media compared to strains from other clades. Their capacity to regulate gene expression across different environmental conditions may be key to their success as zoonotic pathogens.

Human papillomavirus 52 (HPV52) infection is prevalent in the Chinese population, and variations in HPV52 show correlations with oncogenicity. However, no specific variation in HPV52 was reported to show relevancy to infection characteristics. In this study, we retrieved 222 isolates of E6 and L1 full-length genes from 197 Chinese women with HPV52 infection. After sequence alignment and phylogenetic tree construction, we found that 98.39 % of the collected variants belonged to the sublineage B2 and two variants displayed incongruence between the phylogenetic tree of E6 and L1. The analysis of the infection pattern showed that the presence of C6480A/T mutation in the L1 gene was associated with single infection (P=0.01) and persistent infection (P=0.047) of HPV52, while the A6516G nucleotide change was relevant to transient infection (P=0.018). Our data also indicated that variations T309C in the E6 gene and C6480T, C6600A in L1 were more commonly presented in patients with high-grade cytology (P<0.05). One HPV52 breakthrough infection after vaccination was identified, which hinted at the immune escape post-vaccination. Young coitarche age and non-condom usage were correlated to multiple infections. This study provided insight into the polymorphism of HPV52 and revealed the impact of variations in HPV52 on its infection characteristics.

Vibrio parahaemolyticus is the leading cause of seafood-borne gastroenteritis worldwide. A distinctive feature of the O3:K6 pandemic clone, and its derivatives, is the presence of a second, phylogenetically distinct, type III secretion system (T3SS2) encoded within the genomic island VPaI-7. The T3SS2 allows the delivery of effector proteins directly into the cytosol of infected eukaryotic cells to subvert key host-cell processes, critical for V. parahaemolyticus to colonize and cause disease. Furthermore, the T3SS2 also increases the environmental fitness of V. parahaemolyticus in its interaction with bacterivorous protists; hence, it has been proposed that it contributed to the global oceanic spread of the pandemic clone. Several reports have identified T3SS2-related genes in Vibrio and non- Vibrio species, suggesting that the T3SS2 gene cluster is not restricted to the Vibrionaceae and can mobilize through horizontal gene transfer events. In this work, we performed a large-scale genomic analysis to determine the phylogenetic distribution of the T3SS2 gene cluster and its repertoire of effector proteins. We identified putative T3SS2 gene clusters in 1130 bacterial genomes from 8 bacterial genera, 5 bacterial families and 47 bacterial species. A hierarchical clustering analysis allowed us to define six T3SS2 subgroups (I–VI) with different repertoires of effector proteins, redefining the concepts of T3SS2 core and accessory effector proteins. Finally, we identified a subset of the T3SS2 gene clusters (subgroup VI) that lacks most T3SS2 effector proteins described to date and provided a list of 10 novel effector candidates for this subgroup through bioinformatic analysis. Collectively, our findings indicate that the T3SS2 extends beyond the family Vibrionaceae and suggest that different effector protein repertories could have a differential impact on the pathogenic potential and environmental fitness of each bacterium that has acquired the Vibrio T3SS2 gene cluster.

Distancing measures during the COVID-19 lockdown led to a temporary decrease of casual sex partners among clients of the Centre for Sexual Health (CSH) in Amsterdam, the Netherlands. We investigated the effect of this change on the genotypic and phenotypic distribution of Neisseria gonorrhoeae (Ng) isolates from CSH patients. From each Ng-positive patient we sequenced one isolate, resulting in 322 isolates which constituted two groups: 181 isolates cultured from 15 January to 29 February 2020 (before the first lockdown) and 141 cultured from 15 May to 30 June 2020 (during the first lockdown). Patient characteristics showed significantly more symptomatic patients and significantly fewer reported sex partners during the lockdown. Phenotypic data showed an increase in low-level azithromycin resistance and ceftriaxone susceptibility during the lockdown, and this remained after the study period. The diversity in sequence types (STs) decreased slightly during the lockdown. A shift occurred from ST 8156 being predominant before lockdown to ST 9362 during lockdown and a remarkably low median SNP distance of 17 SNPs was found between ST 9362 isolates obtained during lockdown. These findings reflect restricted travel and the change in sexual behaviour of CSH clients during the lockdown, with a potentially increased local transmission of the ST 9362 strain during this period, which led to genotypic and phenotypic changes in the Ng population. This shows that public health measures have far-reaching consequences and should be considered in the surveillance of other infectious diseases.

Acinetobacter baumannii is a common cause of multidrug-resistant (MDR) nosocomial infections around the world. However, little is known about the persistence and dynamics of A. baumannii in a healthy community. This study investigated the role of the community as a prospective reservoir for A. baumannii and explored possible links between hospital and community isolates. A total of 12 independent A. baumannii strains were isolated from human faecal samples from the community in Segamat, Malaysia, in 2018 and 2019. Another 15 were obtained in 2020 from patients at the co-located tertiary public hospital. The antimicrobial resistance profile and biofilm formation ability were analysed, and the relatedness of community and hospital isolates was determined using whole-genome sequencing (WGS). Antibiotic profile analysis revealed that 12 out of 15 hospital isolates were MDR, but none of the community isolates were MDR. However, phylogenetic analysis based on single-nucleotide polymorphisms (SNPs) and a pangenome analysis of core genes showed clustering between four community and two hospital strains. Such clustering of strains from two different settings based on their genomes suggests that these strains could persist in both. WGS revealed 41 potential resistance genes on average in the hospital strains, but fewer (n=32) were detected in the community strains. In contrast, 68 virulence genes were commonly seen in strains from both sources. This study highlights the possible transmission threat to public health posed by virulent A. baumannii present in the gut of asymptomatic individuals in the community.

In the Netherlands, more than half of domestic shigellosis cases are among men who have sex with men (MSM), particularly in the Amsterdam region. However, there is limited insight into which Shigella strains circulate in the Netherlands. Our objective was to assess the added value of whole-genome sequencing (WGS)-based surveillance for Shigella . To this end, we determined the relatedness among Shigella spp. isolates from patients in the Amsterdam region, as well as in an international context, including antimicrobial resistance markers, using WGS. The following criteria were used: it should provide insight into (1) clustering of shigellosis cases and the affected population, (2) the extent of admixture of MSM-associated isolates with those from the broader population and (3) the presence of antimicrobial resistance. It should then lead to more opportunities for targeted control measures. For this study, Shigella isolates from three laboratories in the Amsterdam region obtained between February 2019 and October 2021 were subjected to Illumina WGS at the National Institute for Public Health and the Environment (RIVM). Raw data were quality-checked and assembled, the Shigella serotype was determined with ShigaTyper, and antimicrobial resistance markers were detected using ResFinder and PointFinder. For Shigella sonnei , subclades were determined using Mykrobe. Relatedness of isolates, including 21 international reference genomes, was assessed with core genome multilocus sequence typing. In total, 109 isolates were included, of which 27 were from females (25 %) and 66 were from males (61 %), with which the majority (n=48, 73 %) being from MSM. No information on sex was available for the remaining 16 cases. The WGS data for all isolates, comprising 55  S . sonnei , 52  Shigella flexneri , 1  Shigella boydii and 1  Shigella dysenteriae , met the quality criteria. In total, 14 clusters containing 51 isolates (49 %) were identified, with a median cluster size of 2.5 cases (range: 2–15). Nine out of 14 clusters were MSM-associated, and 8 clusters (57 %) were travel-related. Six of the MSM clusters were related to international reference genomes. The prevalence of antimicrobial resistance markers was higher among isolates from MSM than non-MSM patients, particularly for ciprofloxacin (89 vs 33 %) and azithromycin (58 vs 17 %). In conclusion, about half of Shigella spp. patients were part of a cluster, of which a substantial part were related to international reference genomes, particularly among MSM, and a high prevalence of antimicrobial resistance markers was found. These findings indicate widespread international circulation of Shigella spp., particularly among MSM, with multidrug resistance that hampers treatment of patients. Moreover, the results of this study led to the implementation of a national WGS-based laboratory surveillance programme for Shigella spp. that started in April 2022.

The severity and progression of lung disease are highly variable across individuals with cystic fibrosis (CF) and are imperfectly predicted by mutations in the human gene CFTR, lung microbiome variation or other clinical factors. The opportunistic pathogen Pseudomonas aeruginosa (Pa) dominates airway infections in most CF adults. Here we hypothesized that within–host genetic variation of Pa populations would be associated with lung disease severity. To quantify Pa genetic variation within CF sputum samples, we used deep amplicon sequencing (AmpliSeq) of 209 Pa genes previously associated with pathogenesis or adaptation to the CF lung. We trained machine learning models using Pa single nucleotide variants (SNVs), microbiome diversity data and clinical factors to classify lung disease severity at the time of sputum sampling, and to predict lung function decline after 5 years in a cohort of 54 adult CF patients with chronic Pa infection. Models using Pa SNVs alone classified lung disease severity with good sensitivity and specificity (area under the receiver operating characteristic curve: AUROC=0.87). Models were less predictive of lung function decline after 5 years (AUROC=0.74) but still significantly better than random. The addition of clinical data, but not sputum microbiome diversity data, yielded only modest improvements in classifying baseline lung function (AUROC=0.92) and predicting lung function decline (AUROC=0.79), suggesting that Pa AmpliSeq data account for most of the predictive value. Our work provides a proof of principle that Pa genetic variation in sputum tracks lung disease severity, moderately predicts lung function decline and could serve as a disease biomarker among CF patients with chronic Pa infections.

Limited information is available regarding the genomic characteristics of P. aeruginosa causing ear infections. Our aim is to characterize the genotypic features of an emerging ST316 sublineage causing aural infections in Shanghai. A total of 199 ear swab isolates were subjected to whole genome sequencing (WGS). Complete genomes for two isolates were resolved. We showed this recently emerged sublineage exhibited high-level resistance to fluoroquinolones (FQs) primarily by accumulation of known mutations in quinolone resistance determining regions (QRDRs). Loss-of-function mutations in mexR and mexCD were frequently detected. Mutations in fusA1 (P166S) and parE (S492F) were resident in this sublinage about 2 years after its emergence. Recombination events might be a key driver of genomic diversity in this sublineage. Convergent evolution events on Multidrug-resistant (MDR) determinants were also observed. We generated predictive machine models and identified biomarkers of resistance to gentamicin, fosfomycin, and cefoperazone-sulbactam in this sublineage. This sublineage tended to be less virulent by loss of a series virulence genes represented by ppkA, rhlI, and iron uptake- and antimicrobial activity-related genes. Specific mutations were detected in pilU and lpxB genes that related to surface structures. Moreover, this sublineage differed from non-ST316 isolates in several ways, including virulence genes related to cell surface structure. Our analysis suggested acquisition of a roughly 390 kbp MDR plasmid carrying qnrVC1 might play an important role in the success of this sublinage. Clonal expansion of this sublineage exhibiting enhanced adaptation to cause ear infections is concerning, which requires urgent control measures to be implemented.

The ARTIC protocol uses a multiplexed PCR approach with two primer pools tiling the entire SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) genome. Primer pool updates are necessary for accurate amplicon sequencing of evolving SARS-CoV-2 variants with novel mutations. The suitability of the ARTIC V4 and updated V4.1 primer scheme was assessed using whole genome sequencing of Omicron from clinical samples using Oxford Nanopore Technology. Analysis of Omicron BA.1 genomes revealed that 93.22 % of clinical samples generated improved genome coverage at 50× read depth with V4.1 primers when compared to V4 primers. Additionally, the V4.1 primers improved coverage of BA.1 across amplicons 76 and 88, which resulted in the detection of the variant-defining mutations G22898A, A26530G and C26577G. The Omicron BA.2 sub-variant (VUI-22JAN-01) replaced BA.1 as the dominant variant by March 2022, and analysis of 168 clinical samples showed reduced coverage across amplicons 15 and 75. Upon further interrogation of primer binding sites, a mutation at C4321T [present in 163/168 (97 %) of samples] was identified as a possible cause of complete dropout of amplicon 15. Furthermore, two mutations were identified within the primer binding regions for amplicon 75: A22786C (present in 90 % of samples) and C22792T (present in 12.5 % of samples). Together, these mutations may result in reduced coverage of amplicon 75, and further primer updates would allow the identification of the two BA.2-defining mutations present in amplicon 75: A22688G and T22679C. This work highlights the need for ongoing surveillance of primer matches as circulating variants evolve and change.

Aggregation of children in schools has been established to be a key driver of transmission of infectious diseases. Mathematical models of transmission used to predict the impact of control measures, such as vaccination and testing, commonly depend on self-reported contact data. However, the link between self-reported social contacts and pathogen transmission has not been well described. To address this, we used Staphylococcus aureus as a model organism to track transmission within two secondary schools in England and test for associations between self-reported social contacts, test positivity and the bacterial strain collected from the same students. Students filled out a social contact survey and their S. aureus colonization status was ascertained through self-administered swabs from which isolates were sequenced. Isolates from the local community were also sequenced to assess the representativeness of school isolates. A low frequency of genome-linked transmission precluded a formal analysis of links between genomic and social networks, suggesting that S. aureus transmission within schools is too rare to make it a viable tool for this purpose. Whilst we found no evidence that schools are an important route of transmission, increased colonization rates found within schools imply that school-age children may be an important source of community transmission.

Streptococcus pyogenes genotype emm1 is a successful, globally distributed epidemic clone that is regarded as inherently virulent. An emm1 sublineage, M1UK, that produces increased levels of SpeA toxin was associated with increased scarlet fever and invasive infections in England in 2015/2016. Defined by 27 SNPs in the core genome, M1UK is now dominant in England. To more fully characterize M1UK, we undertook comparative transcriptomic and proteomic analyses of M1UK and contemporary non-M1UK emm1 strains (M1global). Just seven genes were differentially expressed by M1UK compared with contemporary M1global strains. In addition to speA, five genes in the operon that includes glycerol dehydrogenase were upregulated in M1UK (gldA, mipB/talC, pflD, and phosphotransferase system IIC and IIB components), while aquaporin (glpF2) was downregulated. M1UK strains have a stop codon in gldA. Deletion of gldA in M1global abrogated glycerol dehydrogenase activity, and recapitulated upregulation of gene expression within the operon that includes gldA, consistent with a feedback effect. Phylogenetic analysis identified two intermediate emm1 sublineages in England comprising 13/27 (M113SNPs) and 23/27 SNPs (M123SNPs), respectively, that had failed to expand in the population. Proteomic analysis of invasive strains from the four phylogenetic emm1 groups highlighted sublineage-specific changes in carbohydrate metabolism, protein synthesis and protein processing; upregulation of SpeA was not observed in chemically defined medium. In rich broth, however, expression of SpeA was upregulated ~10-fold in both M123SNPs and M1UK sublineages, compared with M113SNPs and M1global. We conclude that stepwise accumulation of SNPs led to the emergence of M1UK. While increased expression of SpeA is a key indicator of M1UK and undoubtedly important, M1UK strains have outcompeted M123SNPs and other emm types that produce similar or more superantigen toxin. We speculate that an accumulation of adaptive SNPs has contributed to a wider fitness advantage in M1UK on an inherently successful emm1 streptococcal background.