- Volume 9, Issue 11, 2023

Carbapenem-resistant Enterobacterales (CRE) are an urgent public health threat. Genomic sequencing is an important tool for investigating CRE. Through the Division of Healthcare Quality Promotion Sentinel Surveillance system, we collected CRE and carbapenem-susceptible Enterobacterales (CSE) from nine clinical laboratories in the USA from 2013 to 2016 and analysed both phenotypic and genomic sequencing data for 680 isolates. We describe the molecular epidemiology and antimicrobial susceptibility testing (AST) data of this collection of isolates. We also performed a phenotype–genotype correlation for the carbapenems and evaluated the presence of virulence genes in Klebsiella pneumoniae complex isolates. These AST and genomic sequencing data can be used to compare and contrast CRE and CSE at these sites and serve as a resource for the antimicrobial resistance research community.

Long non-coding RNAs (lncRNAs) are regulatory molecules interacting in a wide array of biological processes. lncRNAs in fungal pathogens can be responsive to stress and play roles in regulating growth and nutrient acquisition. Recent evidence suggests that lncRNAs may also play roles in virulence, such as regulating pathogenicity-associated enzymes and on-host reproductive cycles. Despite the importance of lncRNAs, only a few model fungi have well-documented inventories of lncRNA. In this study, we apply a recent computational pipeline to predict high-confidence lncRNA candidates in Zymoseptoria tritici, an important global pathogen of wheat impacting global food production. We analyse genomic features of lncRNAs and the most likely associated processes through analyses of expression over a host infection cycle. We find that lncRNAs are frequently expressed during early infection, before the switch to necrotrophic growth. They are mostly located in facultative heterochromatic regions, which are known to contain many genes associated with pathogenicity. Furthermore, we find that lncRNAs are frequently co-expressed with genes that may be involved in responding to host defence signals, such as oxidative stress. Finally, we assess pangenome features of lncRNAs using four additional reference-quality genomes. We find evidence that the repertoire of expressed lncRNAs varies substantially between individuals, even though lncRNA loci tend to be shared at the genomic level. Overall, this study provides a repertoire and putative functions of lncRNAs in Z. tritici enabling future molecular genetics and functional analyses in an important pathogen.

Archamoebae comprises free-living or endobiotic amoebiform protists that inhabit anaerobic or microaerophilic environments and possess mitochondrion-related organelles (MROs) adapted to function anaerobically. We compared in silico reconstructed MRO proteomes of eight species (six genera) and found that the common ancestor of Archamoebae possessed very few typical components of the protein translocation machinery, electron transport chain and tricarboxylic acid cycle. On the other hand, it contained a sulphate activation pathway and bacterial iron–sulphur (Fe-S) assembly system of MIS-type. The metabolic capacity of the MROs, however, varies markedly within this clade. The glycine cleavage system is widely conserved among Archamoebae, except in Entamoeba, probably owing to its role in catabolic function or one-carbon metabolism. MRO-based pyruvate metabolism was dispensed within subgroups Entamoebidae and Rhizomastixidae, whereas sulphate activation could have been lost in isolated cases of Rhizomastix libera, Mastigamoeba abducta and Endolimax sp. The MIS (Fe-S) assembly system was duplicated in the common ancestor of Mastigamoebidae and Pelomyxidae, and one of the copies took over Fe-S assembly in their MRO. In Entamoebidae and Rhizomastixidae, we hypothesize that Fe-S cluster assembly in both compartments may be facilitated by dual localization of the single system. We could not find evidence for changes in metabolic functions of the MRO in response to changes in habitat; it appears that such environmental drivers do not strongly affect MRO reduction in this group of eukaryotes.

Horizontal gene transfer by plasmids can confer metabolic capabilities that expand a host cell’s niche. Yet, it is less understood whether the coalescence of specialized catabolic functions, antibiotic resistances and metal resistances on plasmids provides synergistic benefits. In this study, we report whole-genome assembly and phenotypic analysis of five Salmonella enterica strains isolated in the 1980s from milk powder in Munich, Germany. All strains exhibited the unusual phenotype of lactose-fermentation and encoded either of two variants of the lac operon. Surprisingly, all strains encoded the mobilized colistin resistance gene 9 (mcr-9), long before the first report of this gene in the literature. In two cases, the mcr-9 gene and the lac locus were linked within a large gene island that formed an IncHI2A-type plasmid in one strain but was chromosomally integrated in the other strain. In two other strains, the mcr-9 gene was found on a large IncHI1B/IncP-type plasmid, whereas the lac locus was encoded on a separate chromosomally integrated plasmidic island. The mcr-9 sequences were identical and genomic contexts could not explain the wide range of colistin resistances exhibited by the Salmonella strains. Nucleotide variants did explain phenotypic differences in motility and exopolysaccharide production. The observed linkage of mcr-9 to lactose metabolism, an array of heavy-metal detoxification systems, and other antibiotic resistance genes may reflect a coalescence of specialized phenotypes that improve the spread of colistin resistance in dairy facilities, much earlier than previously suspected.

The rapid warming of the Arctic is threatening the demise of its glaciers and their associated ecosystems. Therefore, there is an urgent need to explore and understand the diversity of genomes resident within glacial ecosystems endangered by human-induced climate change. In this study we use genome-resolved metagenomics to explore the taxonomic and functional diversity of different habitats within glacier-occupied catchments. Comparing different habitats within such catchments offers a natural experiment for understanding the effects of changing habitat extent or even loss upon Arctic microbiota. Through binning and annotation of metagenome-assembled genomes (MAGs) we describe the spatial differences in taxon distribution and their implications for glacier-associated biogeochemical cycling. Multiple taxa associated with carbon cycling included organisms with the potential for carbon monoxide oxidation. Meanwhile, nitrogen fixation was mediated by a single taxon, although diverse taxa contribute to other nitrogen conversions. Genes for sulphur oxidation were prevalent within MAGs implying the potential capacity for sulphur cycling. Finally, we focused on cyanobacterial MAGs, and those within cryoconite, a biodiverse microbe-mineral granular aggregate responsible for darkening glacier surfaces. Although the metagenome-assembled genome of Phormidesmis priestleyi, the cyanobacterium responsible for forming Arctic cryoconite was represented with high coverage, evidence for the biosynthesis of multiple vitamins and co-factors was absent from its MAG. Our results indicate the potential for cross-feeding to sustain P. priestleyi within granular cryoconite. Taken together, genome-resolved metagenomics reveals the vulnerability of glacier-associated microbiota to the deletion of glacial habitats through the rapid warming of the Arctic.

Each year, 15 million infants are born preterm (<37 weeks gestation), representing the leading cause of mortality for children under the age of five. Whilst there is no single cause, factors such as maternal genetics, environmental interactions, and the vaginal microbiome have been associated with an increased risk of preterm birth. Previous studies show that a vaginal microbiota dominated by Lactobacillus is, in contrast to communities containing a mixture of genera, associated with full-term birth. However, this binary principle does not fully consider more nuanced interactions between bacterial strains and the host. Here, through a combination of analyses involving genome-sequenced isolates and strain-resolved metagenomics, we identify that L. jensenii strains from preterm pregnancies are phylogenetically distinct from strains from full-term pregnancies. Detailed analysis reveals several genetic signatures that distinguish preterm birth strains, including genes predicted to be involved in cell wall synthesis, and lactate and acetate metabolism. Notably, we identify a distinct gene cluster involved in cell surface protein synthesis in our preterm strains, and profiling the prevalence of this gene cluster in publicly available genomes revealed it to be predominantly present in the preterm-associated clade. This study contributes to the ongoing search for molecular biomarkers linked to preterm birth and opens up new avenues for exploring strain-level variations and mechanisms that may contribute to preterm birth.

The most common approach to sampling the bacterial populations within an infected or colonized host is to sequence genomes from a single colony obtained from a culture plate. However, it is recognized that this method does not capture the genetic diversity in the population. Sequencing a mixture of several colonies (pool-seq) is a better approach to detect population heterogeneity, but it is more complex to analyse due to different types of heterogeneity, such as within-clone polymorphisms, multi-strain mixtures, multi-species mixtures and contamination. Here, we compared 8 single-colony isolates (singles) and pool-seq on a set of 2286 Staphylococcus aureus culture samples to identify features that can distinguish pure samples, samples undergoing intraclonal variation and mixed strain samples. The samples were obtained by swabbing 3 body sites on 85 human participants quarterly for a year, who initially presented with a methicillin-resistant S. aureus skin and soft-tissue infection (SSTI). We compared parameters such as sequence quality, contamination, allele frequency, nucleotide diversity and pangenome diversity in each pool to those for the corresponding singles. Comparing singles from the same culture plate, we found that 18% of sample collections contained mixtures of multiple multilocus sequence types (MLSTs or STs). We showed that pool-seq data alone could predict the presence of multi-ST populations with 95% accuracy. We also showed that pool-seq could be used to estimate the number of intra-clonal polymorphic sites in the population. Additionally, we found that the pool may contain clinically relevant genes such as antimicrobial resistance markers that may be missed when only examining singles. These results highlight the potential advantage of analysing genome sequences of total populations obtained from clinical cultures rather than single colonies.

Due to the emergence of non-vaccine serotypes in vaccinated populations, Streptococcus pneumoniae remains a major global health challenge despite advances in vaccine development. Serotype 16F is among the predominant non-vaccine serotypes identified among vaccinated infants in South Africa (SA). To characterize lineages and antimicrobial resistance in 16F isolates obtained from South Africa and place the local findings in a global context, we analysed 10 923  S . pneumoniae carriage isolates obtained from infants recruited as part of a broader SA birth cohort. We inferred serotype, resistance profile for penicillin, chloramphenicol, cotrimoxazole, erythromycin and tetracycline, and global pneumococcal sequence clusters (GPSCs) from genomic data. To ensure global representation, we also included S. pneumoniae carriage and disease isolates from the Global Pneumococcal Sequencing (GPS) project database (n=19 607, collected from 49 countries across 5 continents, 1995–2018, accessed 17 March 2022). Nine per cent (934/10923) of isolates obtained from infants in the Drakenstein community in SA and 2 %(419/19607) of genomes in the GPS dataset were serotype 16F. Serotype 16F isolates were from 28 different lineages of S. pneumoniae, with GPSC33 and GPSC46 having the highest proportion of serotype 16F isolates at 26 % (346/1353) and 53 % (716/1353), respectively. Serotype 16F isolates were identified globally, but most isolates were collected from Africa. GPSC33 was associated with carriage [OR (95 % CI) 0.24 (0.09–0.66); P=0.003], while GPSC46 was associated with disease [OR (95 % CI) 19.9 (2.56–906.50); P=0.0004]. Ten per cent (37/346) and 15 % (53/346) of isolates within GPSC33 had genes associated with resistance to penicillin and co-trimoxazole, respectively, and 18 % (128/716) of isolates within GPSC46 had genes associated with resistance to co-trimoxazole. Resistant isolates formed genetic clusters, which may suggest emerging resistant lineages. Serotype 16F lineages were common in southern Africa. Some of these lineages were associated with disease and resistance to penicillin and cotrimoxazole. We recommend continuous genomic surveillance to determine the long-term impact of serotype 16F lineages on vaccine efficacy and antimicrobial therapy globally. Investing in vaccine strategies that offer protection over a wide range of serotypes/lineages remains essential. This paper contains data hosted by Microreact.

Mobile genetic elements can innovate bacteria with new traits. In plant pathogenic Streptomyces, frequent and recent acquisition of integrative and conjugative or mobilizable genetic elements is predicted to lead to the emergence of new lineages that gained the capacity to synthesize Thaxtomin, a phytotoxin neccesary for induction of common scab disease on tuber and root crops. Here, we identified components of the Streptomyces -potato pathosystem implicated in virulence and investigated them as a nested and interacting system to reevaluate evolutionary models. We sequenced and analysed genomes of 166 strains isolated from over six decades of sampling primarily from field-grown potatoes. Virulence genes were associated to multiple subtypes of genetic elements differing in mechanisms of transmission and evolutionary histories. Evidence is consistent with few ancient acquisition events followed by recurrent loss or swaps of elements carrying Thaxtomin A-associated genes. Subtypes of another genetic element implicated in virulence are more distributed across Streptomyces . However, neither the subtype classification of genetic elements containing virulence genes nor taxonomic identity was predictive of pathogenicity on potato. Last, findings suggested that phytopathogenic strains are generally endemic to potato fields and some lineages were established by historical spread and further dispersed by few recent transmission events. Results from a hierarchical and system-wide characterization refine our understanding by revealing multiple mechanisms that gene and bacterial dispersion have had on shaping the evolution of a Gram-positive pathogen in agricultural settings.

The repeated emergence of multi-drug-resistant (MDR) Escherichia coli clones is a threat to public health globally. In recent work, drug-resistant E. coli were shown to be capable of displacing commensal E. coli in the human gut. Given the rapid colonization observed in travel studies, it is possible that the presence of a type VI secretion system (T6SS) may be responsible for the rapid competitive advantage of drug-resistant E. coli clones. We employed large-scale genomic approaches to investigate this hypothesis. First, we searched for T6SS genes across a curated dataset of over 20 000 genomes representing the full phylogenetic diversity of E. coli . This revealed large, non-phylogenetic variation in the presence of T6SS genes. No association was found between T6SS gene carriage and MDR lineages. However, multiple clades containing MDR clones have lost essential structural T6SS genes. We characterized the T6SS loci of ST410 and ST131 and identified specific recombination and insertion events responsible for the parallel loss of essential T6SS genes in two MDR clones.

Chlamydia abortus, an obligate intracellular bacterium, is a major causative agent of reproductive loss in ruminants, with zoonotic potential. Though this pathogen is primarily known to infect livestock, recent studies have detected and isolated genetically distinct avian strains of C. abortus from wild birds globally. Before this study, only five avian C. abortus genomes were publicly available. Therefore, we performed culture-independent probe-based whole-genome sequencing on clinical swabs positive for avian C. abortus obtained from Australian Torresian crows (Corvus orru) in 2019 and 2020. We successfully obtained draft genomes for three avian C. abortus strains (C1, C2 and C3), each comprising draft chromosomes with lengths of 1 115 667, 1 120 231 and 1 082 115 bp, and associated 7 553 bp plasmids, with a genome completeness exceeding 92 %. Molecular characterization revealed that these three strains comprise a novel sequence type (ST333), whilst phylogenetic analyses placed all three strains in a cluster with other avian C. abortus genomes. Interestingly, these three strains share a distant genomic relation (2693 single nucleotide variants) with the reference strain 15-58d/44 (ST152), isolated from a Eurasian magpie (Pica pica) in Poland, highlighting the need for more publicly available genomes. Broad comparative analyses with other avian C. abortus genomes revealed that the three draft genomes contain conserved Chlamydia genomic features, including genes coding for type III secretion system and polymorphic membrane proteins, and potential virulence factors such as the large chlamydial cytotoxin, warranting further studies. This research provides the first avian C. abortus draft genomes from Australian birds, highlighting Torresian crows as novel reservoir hosts for these potential pathogens, and demonstrates a practical methodology for sequencing novel Chlamydia genomes without relying on traditional cell culture.

The National Tuberculosis Programme (NTBP) monitors the occurrence and spread of tuberculosis (TB) and multidrug-resistant TB (MDR-TB) in Singapore. Since 2020, whole-genome sequencing (WGS) of Mycobacterium tuberculosis isolates has been performed at the National Public Health Laboratory (NPHL) for genomic surveillance, replacing spoligotyping and mycobacterial interspersed repetitive unit-variable number tandem repeats analysis (MIRU-VNTR). Four thousand three hundred and seven samples were sequenced from 2014 to January 2023, initially as research projects and later developed into a comprehensive public health surveillance programme. Currently, all newly diagnosed culture-positive cases of TB in Singapore are prospectively sent for WGS, which is used to perform lineage classification, predict drug resistance profiles and infer genetic relationships between TB isolates. This paper describes NPHL’s operational and technical experiences with implementing the WGS service in an urban TB-endemic setting, focusing on cluster detection and genomic drug susceptibility testing (DST). Cluster detection: WGS has been used to guide contact tracing by detecting clusters and discovering unknown transmission networks. Examples have been clusters in a daycare centre, housing apartment blocks and a horse-racing betting centre. Genomic DST: genomic DST prediction (gDST) identifies mutations in core genes known to be associated with TB drug resistance catalogued in the TBProfiler drug resistance mutation database. Mutations are reported with confidence scores according to a standardized approach referencing NPHL’s internal gDST confidence database, which is adapted from the World Health Organization (WHO) TB drug mutation catalogue. Phenotypic–genomic concordance was observed for the first-line drugs ranging from 2959/2998 (98.7 %) (ethambutol) to 2983/2996 (99.6 %) (rifampicin). Aspects of internal database management, reporting standards and caveats in results interpretation are discussed.

In the province of Alberta, Canada, invasive disease caused by Streptococcus pneumoniae serogroup 20 (serotypes 20A/20B) has been increasing in incidence. Here, we characterize provincial invasive serogroup 20 isolates collected from 1993 to 2019 alongside invasive and non-invasive serogroup 20 isolates from the Global Pneumococcal Sequencing (GPS) Project collected from 1998 to 2015. Trends in clinical metadata and geographic location were evaluated, and serogroup 20 isolate genomes were subjected to molecular sequence typing, virulence and antimicrobial resistance factor mining, phylogenetic analysis and pangenome calculation. Two hundred and seventy-four serogroup 20 isolates from Alberta were sequenced, and analysed along with 95 GPS Project genomes. The majority of invasive Alberta serogroup 20 isolates were identified after 2007 in primarily middle-aged adults and typed predominantly as ST235, a sequence type that was rare among GPS Project isolates. Most Alberta isolates carried a full-length whaF capsular gene, suggestive of serotype 20B. All Alberta and GPS Project genomes carried molecular resistance determinants implicated in fluoroquinolone and macrolide resistance, with a few Alberta isolates exhibiting phenotypic resistance to azithromycin, clindamycin, erythromycin, tetracycline and trimethoprim-sulfamethoxazole, as well as non-susceptibility to tigecycline. All isolates carried multiple virulence factors including those involved in adherence, immune modulation and nutrient uptake, as well as exotoxins and exoenzymes. Phylogenetically, Alberta serogroup 20 isolates clustered with predominantly invasive GPS Project isolates from the USA, Israel, Brazil and Nepal. Overall, this study highlights the increasing incidence of invasive S. pneumoniae serogroup 20 disease in Alberta, Canada, and provides insights into the genetic and clinical characteristics of these isolates within a global context.

Complete reference genomes, including correct feature annotations, are a fundamental aspect of genomic biology. In the case of protozoan species such as Giardia duodenalis, a major human and animal parasite worldwide, accurate genome annotation can deepen our understanding of the evolution of parasitism and pathogenicity by identifying genes underlying key traits and clinically relevant cellular mechanisms, and by extension, the development of improved prevention strategies and treatments. This study used bioinformatics analyses of Giardia mRNA libraries to characterize known introns and identify new intron candidates, working towards completion of the G. duodenalis assemblage A strain ‘WB’ genome and further elucidating Giardia’s gene expression. By using a set of experimentally validated positive control loci to calibrate our intron detection pipeline, we were able to detect evidence of previously missed candidate splice junctions directly from expressed transcript data. These intron candidates were further studied in silico using NMDS (non-metric multidimensional scaling) clustering to determine shared characteristics and their relative importance such as secondary structure, splicing efficiency and motif conservation, and thus to refine intron models. Results from this study identified 34 new intron candidates, with several potential introns showing evidence that secondary structure of the mRNA molecule might play a more significant role in splicing than previously reported eukaryotic splicing activity mediated by a reduced spliceosome present in G. duodenalis.

Chronic rhinosinusitis (CRS) is a common chronic sinonasal mucosal inflammation associated with Staphylococcus aureus biofilm and relapsing infections. This study aimed to determine rates of S. aureus persistence and pathoadaptation in CRS patients by investigating the genomic relatedness and antibiotic resistance/tolerance in longitudinally collected S. aureus clinical isolates. A total of 68 S . aureus paired isolates (34 pairs) were sourced from 34 CRS patients at least 6 months apart. Isolates were grown into 48 h biofilms and tested for tolerance to antibiotics. A hybrid sequencing strategy was used to obtain high-quality reference-grade assemblies of all isolates. Single nucleotide variants (SNV) divergence in the core genome and sequence type clustering were used to analyse the relatedness of the isolate pairs. Single nucleotide and structural genome variations, plasmid similarity, and plasmid copy numbers between pairs were examined. Our analysis revealed that 41 % (14/34 pairs) of S. aureus isolates were persistent, while 59 % (20/34 pairs) were non-persistent. Persistent isolates showed episode-specific mutational changes over time with a bias towards events in genes involved in adhesion to the host and mobile genetic elements such as plasmids, prophages, and insertion sequences. Furthermore, a significant increase in the copy number of conserved plasmids of persistent strains was observed. This was accompanied by a significant increase in biofilm tolerance against all tested antibiotics, which was linked to a significant increase in biofilm biomass over time, indicating a potential biofilm pathoadaptive process in persistent isolates. In conclusion, our study provides important insights into the mutational changes during S. aureus persistence in CRS patients highlighting potential pathoadaptive mechanisms in S. aureus persistent isolates culminating in increased biofilm biomass.

Plasmids play important roles in bacterial genome diversification. In the Serratia marcescens complex (SMC), a notable contribution of plasmids to genome diversification was also suggested by our recent analysis of >600 draft genomes. As accurate analyses of plasmids in draft genomes are difficult, in this study we analysed 142 closed genomes covering the entire complex, 67 of which were obtained in this study, and identified 132 plasmids (1.9–244.4 kb in length) in 77 strains. While the average numbers of plasmids in clinical and non-clinical strains showed no significant difference, strains belonging to clade 2 (one of the two hospital-adapted lineages) contained more plasmids than the others. Pangenome analysis revealed that of the 28 954 genes identified, 12.8 % were plasmid-specific, and 1.4 % were present in plasmids or chromosomes depending on the strain. In the latter group, while transposon-related genes were most prevalent (31.4 % of the function-predicted genes), genes related to antimicrobial resistance and heavy metal resistance accounted for a notable proportion (22.7 %). Mash distance-based clustering separated the 132 plasmids into 23 clusters and 50 singletons. Most clusters/singletons showed notably different GC contents compared to those of host chromosomes, suggesting their recent or relatively recent appearance in the SMC. Among the 23 clusters, 17 were found in only clinical or only non-clinical strains, suggesting the possible preference of their distribution on the environmental niches of host strains. Regarding the host strain phylogeny, 16 clusters were distributed in two or more clades, suggesting their interclade transmission. Moreover, for many plasmids, highly homologous plasmids were found in other species, indicating the broadness of their potential host ranges, beyond the genus, family, order, class or even phylum level. Importantly, highly homologous plasmids were most frequently found in Klebsiella pneumoniae and other species in the family Enterobacteriaceae , suggesting that this family, particularly K. pneumoniae , is the main source for plasmid exchanges with the SMC. These results highlight the power of closed genome-based analysis in the investigation of plasmids and provide important insights into the nature of plasmids distributed in the SMC.