- Volume 9, Issue 10, 2023
Volume 9, Issue 10, 2023
- Reviews
-
- Pathogens and Epidemiology
-
-
Continuing genomic evolution of the Neisseria meningitidis cc11.2 urethritis clade, NmUC: a narrative review
More LessNeisseria meningitidis (Nm) is a bacterial pathogen responsible for invasive meningococcal disease. Though typically colonizing the nasopharynx, multiple outbreaks of meningococcal urethritis were first reported in 2015–2016; outbreaks originally presumed to be caused by Neisseria gonorrhoeae (Ng). Genomic analysis revealed that the Nm isolates causing these outbreaks were a distinct clade, and had integrated gonococcal DNA at multiple genomic sites, including the gonococcal denitrification apparatus aniA–norB, a partial gonococcal operon of five genes containing ispD, and the acetylglutamate kinase gene argB with the adjacent gonococcal locus NGO0843. The urethritis isolates had also deleted the group C capsule biosynthesis genes cssA/B/C and csc, resulting in loss of capsule. Collectively, these isolates form the N. meningitidis urethritis clade (NmUC). Genomic analysis of recent (2016–2022) NmUC isolates revealed that the genomic features have been maintained in the clade, implying that they are important for NmUC’s status as a urogenital pathogen. Furthermore, the analysis revealed the emergence of a sub-clade, designated NmUC-B, phylogenetically separated from the earlier NmUC-A. This sub-clade has integrated additional gonococcal alleles into the genome, including alleles associated with antimicrobial resistance. NmUC continues to adapt to a urethral niche and evolve as a urogenital pathogen.
-
- Research Articles
-
- Genomic Methodologies
-
-
The advantage of intergenic regions as genomic features for machine-learning-based host attribution of Salmonella Typhimurium from the USA
More LessSalmonella enterica is a taxonomically diverse pathogen with over 2600 serovars associated with a wide variety of animal hosts including humans, other mammals, birds and reptiles. Some serovars are host-specific or host-restricted and cause disease in distinct host species, while others, such as serovar S. Typhimurium (STm), are generalists and have the potential to colonize a wide variety of species. However, even within generalist serovars such as STm it is becoming clear that pathovariants exist that differ in tropism and virulence. Identifying the genetic factors underlying host specificity is complex, but the availability of thousands of genome sequences and advances in machine learning have made it possible to build specific host prediction models to aid outbreak control and predict the human pathogenic potential of isolates from animals and other reservoirs. We have advanced this area by building host-association prediction models trained on a wide range of genomic features and compared them with predictions based on nearest-neighbour phylogeny. SNPs, protein variants (PVs), antimicrobial resistance (AMR) profiles and intergenic regions (IGRs) were extracted from 3883 high-quality STm assemblies collected from humans, swine, bovine and poultry in the USA, and used to construct Random Forest (RF) machine learning models. An additional 244 recent STm assemblies from farm animals were used as a test set for further validation. The models based on PVs and IGRs had the best performance in terms of predicting the host of origin of isolates and outperformed nearest-neighbour phylogenetic host prediction as well as models based on SNPs or AMR data. However, the models did not yield reliable predictions when tested with isolates that were phylogenetically distinct from the training set. The IGR and PV models were often able to differentiate human isolates in clusters where the majority of isolates were from a single animal source. Notably, IGRs were the feature with the best performance across multiple models which may be due to IGRs acting as both a representation of their flanking genes, equivalent to PVs, while also capturing genomic regulatory variation, such as altered promoter regions. The IGR and PV models predict that ~45 % of the human infections with STm in the USA originate from bovine, ~40 % from poultry and ~14.5 % from swine, although sequences of isolates from other sources were not used for training. In summary, the research demonstrates a significant gain in accuracy for models with IGRs and PVs as features compared to SNP-based and core genome phylogeny predictions when applied within the existing population structure. This article contains data hosted by Microreact.
-
-
-
Genetic variants linked to the phenotypic outcome of invasive disease and carriage of Neisseria meningitidis
Neisseria meningitidis can be a human commensal in the upper respiratory tract but is also capable of causing invasive diseases such as meningococcal meningitis and septicaemia. No specific genetic markers have been detected to distinguish carriage from disease isolates. The aim here was to find genetic traits that could be linked to phenotypic outcomes associated with carriage versus invasive N. meningitidis disease through a bacterial genome-wide association study (GWAS). In this study, invasive N. meningitidis isolates collected in Sweden (n=103) and carriage isolates collected at Örebro University, Sweden (n=213) 2018–2019 were analysed. The GWAS analysis, treeWAS, was applied to single-nucleotide polymorphisms (SNPs), genes and k-mers. One gene and one non-synonymous SNP were associated with invasive disease and seven genes and one non-synonymous SNP were associated with carriage isolates. The gene associated with invasive disease encodes a phage transposase (NEIS1048), and the associated invasive SNP glmU S373C encodes the enzyme N-acetylglucosamine 1-phosphate (GlcNAC 1-P) uridyltransferase. Of the genes associated with carriage isolates, a gene variant of porB encoding PorB class 3, the genes pilE/pilS and tspB have known functions. The SNP associated with carriage was fkbp D33N, encoding a FK506-binding protein (FKBP). K-mers from PilS, tbpB and tspB were found to be associated with carriage, while k-mers from mtrD and tbpA were associated with invasiveness. In the genes fkbp, glmU, PilC and pilE, k-mers were found that were associated with both carriage and invasive isolates, indicating that specific variations within these genes could play a role in invasiveness. The data presented here highlight genetic traits that are significantly associated with invasive or carriage N. meningitidis across the species population. These traits could prove essential to our understanding of the pathogenicity of N. meningitidis and could help to identify future vaccine targets.
-
- Functional Genomics and Microbe–Niche Interactions
-
-
Comparative genomic analysis of Chryseobacterium species: deep insights into plant-growth-promoting and halotolerant capacities
More LessMembers of the genus Chryseobacterium have attracted great interest as beneficial bacteria that can promote plant growth and biocontrol. Given the recent risks of climate change, it is important to develop tolerance strategies for efficient applications of plant-beneficial bacteria in saline environments. However, the genetic determinants of plant-growth-promoting and halotolerance effects in Chryseobacterium have not yet been investigated at the genomic level. Here, a comparative genomic analysis was conducted with seven Chryseobacterium species. Phylogenetic and phylogenomic analyses revealed niche-specific evolutionary distances between soil and freshwater Chryseobacterium species, consistent with differences in genomic statistics, indicating that the freshwater bacteria have smaller genome sizes and fewer genes than the soil bacteria. Phosphorus- and zinc-cycling genes (required for nutrient acquisition in plants) were universally present in all species, whereas nitrification and sulphite reduction genes (required for nitrogen- and sulphur-cycling, respectively) were distributed only in soil bacteria. A pan-genome containing 6842 gene clusters was constructed, which reflected the general features of the core, accessory and unique genomes. Halotolerant species with an accessory genome shared a Kdp potassium transporter and biosynthetic pathways for branched-chain amino acids and the carotenoid lycopene, which are associated with countermeasures against salt stress. Protein–protein interaction network analysis was used to define the genetic determinants of Chryseobacterium salivictor NBC122 that reduce salt damage in bacteria and plants. Sixteen hub genes comprised the aromatic compound degradation and Por secretion systems, which are required to cope with complex stresses associated with saline environments. Horizontal gene transfer and CRISPR–Cas analyses indicated that C. salivictor NBC122 underwent more evolutionary events when interacting with different environments. These findings provide deep insights into genomic adaptation to dynamic interactions between plant-growth-promoting Chryseobacterium and salt stress.
-
-
-
Molecular basis of the phenotypic variants arising from a Pseudoalteromonas lipolytica mutator
More LessBacterial deficiencies in the DNA repair system can produce mutator strains that promote adaptive microevolution. However, the role of mutator strains in marine Pseudoalteromonas , capable of generating various gain-of-function genetic variants within biofilms, remains largely unknown. In this study, inactivation of mutS in Pseudoalteromonas lipolytica conferred an approximately 100-fold increased resistance to various antibiotics, including ciprofloxacin, rifampicin and aminoglycoside. Furthermore, the mutator of P. lipolytica generated variants that displayed enhanced biofilm formation but reduced swimming motility, indicating a high phenotypic diversity within the ΔmutS population. Additionally, we observed a significant production rate of approximately 50 % for the translucent variants, which play important roles in biofilm formation, when the ΔmutS strain was cultured on agar plates or under shaking conditions. Using whole-genome deep-sequencing combined with genetic manipulation, we demonstrated that point mutations in AT00_17115 within the capsular biosynthesis cluster were responsible for the generation of translucent variants in the ΔmutS subpopulation, while mutations in flagellar genes fliI and flgP led to a decrease in swimming motility. Collectively, this study reveals a specific mutator-driven evolution in P. lipolytica , characterized by substantial genetic and phenotypic diversification, thereby offering a reservoir of genetic attributes associated with microbial fitness.
-
- Pathogens and Epidemiology
-
-
Characterisation of neonatal Staphylococcus capitis NRCS-A isolates compared with non NRCS-A Staphylococcus capitis from neonates and adults
Staphylococcus capitis is a frequent cause of late-onset sepsis in neonates admitted to Neonatal Intensive Care Units (NICU). One clone of S. capitis , NRCS-A has been isolated from NICUs globally although the reasons for the global success of this clone are not well understood.
We analysed a collection of S. capitis colonising babies admitted to two NICUs, one in the UK and one in Germany as well as corresponding pathological clinical isolates. Genome analysis identified a population structure of three groups; non-NRCS-A isolates, NRCS-A isolates, and a group of ‘proto NRCS-A’ – isolates closely related to NRCS-A but not associated with neonatal infection. All bloodstream isolates belonged to the NRCS-A group and were indistinguishable from strains carried on the skin or in the gut. NRCS-A isolates showed increased tolerance to chlorhexidine and antibiotics relative to the other S. capitis as well as enhanced ability to grow at higher pH values. Analysis of the pangenome of 138 isolates identified characteristic nsr and tarJ genes in both the NRCS-A and proto groups. A CRISPR-cas system was only seen in NRCS-A isolates which also showed enrichment of genes for metal acquisition and transport.
We found evidence for transmission of S. capitis NRCS-A within NICU, with related isolates shared between babies and multiple acquisitions by some babies. Our data show NRCS-A strains commonly colonise uninfected babies in NICU representing a potential reservoir for potential infection. This work provides more evidence that adaptation to survive in the gut and on skin facilitates spread of NRCS-A, and that metal acquisition and tolerance may be important to the biology of NRCS-A. Understanding how NRCS-A survives in NICUs can help develop infection control procedures against this clone.
-
-
-
Comparative genomics of Mycoplasma feriruminatoris, a fast-growing pathogen of wild Caprinae
Mycoplasma feriruminatoris is a fast-growing Mycoplasma species isolated from wild Caprinae and first described in 2013. M. feriruminatoris isolates have been associated with arthritis, kerato conjunctivitis, pneumonia and septicemia, but were also recovered from apparently healthy animals. To better understand what defines this species, we performed a genomic survey on 14 strains collected from free-ranging or zoo-housed animals between 1987 and 2017, mostly in Europe. The average chromosome size of the M. feriruminatoris strains was 1,040±0,024 kbp, with 24 % G+C and 852±31 CDS. The core genome and pan-genome of the M. feriruminatoris species contained 628 and 1312 protein families, respectively. The M. feriruminatoris strains displayed a relatively closed pan-genome, with many features and putative virulence factors shared with species from the M. mycoides cluster, including the MIB-MIP Ig cleavage system, a repertoire of DUF285 surface proteins and a complete biosynthetic pathway for galactan. M. feriruminatoris genomes were found to be mostly syntenic, although repertoires of mobile genetic elements, including Mycoplasma Integrative and Conjugative Elements, insertion sequences, and a single plasmid varied. Phylogenetic- and gene content analyses confirmed that M. feriruminatoris was closer to the M. mycoides cluster than to the ruminant species M. yeatsii and M. putrefaciens . Ancestral genome reconstruction showed that the emergence of the M. feriruminatoris species was associated with the gain of 17 gene families, some of which encode defence enzymes and surface proteins, and the loss of 25 others, some of which are involved in sugar transport and metabolism. This comparative study suggests that the M. mycoides cluster could be extended to include M. feriruminatoris . We also find evidence that the specific organization and structure of the DnaA boxes around the oriC of M. feriruminatoris may contribute to drive the remarkable fast growth of this minimal bacterium.
-
-
-
Contrasting genes conferring short- and long-term biofilm adaptation in Listeria
Listeria monocytogenes is an opportunistic food-borne bacterium that is capable of infecting humans with high rates of hospitalization and mortality. Natural populations are genotypically and phenotypically variable, with some lineages being responsible for most human infections. The success of L. monocytogenes is linked to its capacity to persist on food and in the environment. Biofilms are an important feature that allow these bacteria to persist and infect humans, so understanding the genetic basis of biofilm formation is key to understanding transmission. We sought to investigate the biofilm-forming ability of L. monocytogenes by identifying genetic variation that underlies biofilm formation in natural populations using genome-wide association studies (GWAS). Changes in gene expression of specific strains during biofilm formation were then investigated using RNA sequencing (RNA-seq). Genetic variation associated with enhanced biofilm formation was identified in 273 genes by GWAS and differential expression in 220 genes by RNA-seq. Statistical analyses show that the number of overlapping genes flagged by either type of experiment is less than expected by random sampling. This novel finding is consistent with an evolutionary scenario where rapid adaptation is driven by variation in gene expression of pioneer genes, and this is followed by slower adaptation driven by nucleotide changes within the core genome.
-
-
-
Evidence of potential Campylobacter jejuni zooanthroponosis in captive macaque populations
Non-human primates share recent common ancestry with humans and exhibit comparable disease symptoms. Here, we explored the transmission potential of enteric bacterial pathogens in monkeys exhibiting symptoms of recurrent diarrhoea in a biomedical research facility in China. The common zoonotic bacterium Campylobacter jejuni was isolated from macaques (Macaca mulatta and Macaca fascicularis) and compared to isolates from humans and agricultural animals in Asia. Among the monkeys sampled, 5 % (44/973) tested positive for C. jejuni , 11 % (5/44) of which displayed diarrhoeal symptoms. Genomic analysis of monkey isolates, and 1254 genomes from various sources in Asia, were used to identify the most likely source of human infection. Monkey and human isolates shared high average nucleotide identity, common MLST clonal complexes and clustered together on a phylogeny. Furthermore, the profiles of putative antimicrobial resistance genes were similar between monkeys and humans. Taken together these findings suggest that housed macaques became infected with C. jejuni either directly from humans or via a common contamination source.
-
-
-
The clinical, genomic, and microbiological profile of invasive multi-drug resistant Escherichia coli in a major teaching hospital in the United Kingdom
Escherichia coli is a ubiquitous component of the human gut microbiome, but is also a common pathogen, causing around 40, 000 bloodstream infections (BSI) in the United Kingdom (UK) annually. The number of E. coli BSI has increased over the last decade in the UK, and emerging antimicrobial resistance (AMR) profiles threaten treatment options. Here, we combined clinical, epidemiological, and whole genome sequencing data with high content imaging to characterise over 300 E. coli isolates associated with BSI in a large teaching hospital in the East of England. Overall, only a limited number of sequence types (ST) were responsible for the majority of organisms causing invasive disease. The most abundant (20 % of all isolates) was ST131, of which around 90 % comprised the pandemic O25b:H4 group. ST131-O25b:H4 isolates were frequently multi-drug resistant (MDR), with a high prevalence of extended spectrum β-lactamases (ESBL) and fluoroquinolone resistance. There was no association between AMR phenotypes and the source of E. coli bacteraemia or whether the infection was healthcare-associated. Several clusters of ST131 were genetically similar, potentially suggesting a shared transmission network. However, there was no clear epidemiological associations between these cases, and they included organisms from both healthcare-associated and non-healthcare-associated origins. The majority of ST131 isolates exhibited strong binding with an anti-O25b antibody, raising the possibility of developing rapid diagnostics targeting this pathogen. In summary, our data suggest that a restricted set of MDR E. coli populations can be maintained and spread across both community and healthcare settings in this location, contributing disproportionately to invasive disease and AMR.
-
-
-
Phenotypic and genomic characterization of Pseudomonas aeruginosa isolates recovered from catheter-associated urinary tract infections in an Egyptian hospital
More LessCatheter-associated urinary tract infections (CAUTIs) represent one of the major healthcare-associated infections, and Pseudomonas aeruginosa is a common Gram-negative bacterium associated with catheter infections in Egyptian clinical settings. The present study describes the phenotypic and genotypic characteristics of 31 P . aeruginosa isolates recovered from CAUTIs in an Egyptian hospital over a 3 month period. Genomes of isolates were of good quality and were confirmed to be P. aeruginosa by comparison to the type strain (average nucleotide identity, phylogenetic analysis). Clonal diversity among the isolates was determined; eight different sequence types were found (STs 244, 357, 381, 621, 773, 1430, 1667 and 3765), of which ST357 and ST773 are considered to be high-risk clones. Antimicrobial resistance (AMR) testing according to European Committee on Antimicrobial Susceptibility Testing (EUCAST) guidelines showed that the isolates were highly resistant to quinolones [ciprofloxacin (12/31, 38.7 %) and levofloxacin (9/31, 29 %) followed by tobramycin (10/31, 32.5 %)] and cephalosporins (7/31, 22.5 %). Genotypic analysis of resistance determinants predicted all isolates to encode a range of AMR genes, including those conferring resistance to aminoglycosides, β-lactamases, fluoroquinolones, fosfomycin, sulfonamides, tetracyclines and chloramphenicol. One isolate was found to carry a 422 938 bp pBT2436-like megaplasmid encoding OXA-520, the first report from Egypt of this emerging family of clinically important mobile genetic elements. All isolates were able to form biofilms and were predicted to encode virulence genes associated with adherence, antimicrobial activity, anti-phagocytosis, phospholipase enzymes, iron uptake, proteases, secretion systems and toxins. The present study shows how phenotypic analysis alongside genomic analysis may help us understand the AMR and virulence profiles of P. aeruginosa contributing to CAUTIs in Egypt.
-
- Evolution and Responses to Interventions
-
-
Whole-genome sequencing and phylogenetic analysis of coxsackievirus-A16 strains causing hand, foot and mouth disease (HFMD) in India
Hand, foot and mouth disease (HFMD) is a common childhood infectious disease, caused by enteroviruses (EVs), which can present with typical or atypical lesions. The illness is self-limiting, but it can also have serious complications. Since 1997, HFMD infections have become endemic and have increased to epidemic proportions across the Asia Pacific region, including India. Coxsackievirus-A16 (CV-A16) outbreaks occurred in India from 2005 onwards, although the clinical symptoms were noticeably different during this period. Understanding the population dynamics of enteroviruses that cause HFMD is crucial in the post-polio era because one of the circulating strain may replace another as the dominant strain. The aim of this study is to describe the genetic features of the CV-A16 strains isolated from hand, foot and mouth disease (HFMD) patients in India. Reverse transcription PCR (RT-PCR) and cell-culture-based isolation of CV-A16 was done from the 55 clinical samples. The entire genome of the CV-A16 isolate was performed from the seven isolates. After the sequences were analysed, a phylogenetic tree was created using bioinformatics tools. The total genomic length obtained was 7411 base pairs (bp). Nucleotide similarity across various regions, including 5′UTR, P1, P2 and 3′UTR, ranged from 87.0–97.9 %, 77.0–95.4 %, 80.3–96.9 %, and 77.9–96.2 %, respectively. Correspondingly, similarities in the VP1 region’s nucleotide and amino acid sequences were 91.4–96.4 % and 99.3–99.7 %, respectively. Phylogenetic analysis highlighted that CV-A16 strains identified in Pune, Maharashtra, were grouped within the same cluster. The analysed CV-A16 isolates in this study aligned with subgenotype B1c. These findings have far-reaching implications for the surveillance, prevention and management of HFMD and CV-A16. Monitoring the dynamics of CV-A16 strains, informed by the genetic characteristics identified here, will significantly impact strategies aimed at tackling HFMD and its associated public health challenges.
-
- Methods
-
- Microbial Communities
-
-
Optimization of high-throughput 16S rRNA gene amplicon sequencing: an assessment of PCR pooling, mastermix use and contamination
16S rRNA gene sequencing is widely used to characterize human and environmental microbiomes. Sequencing at scale facilitates better powered studies but is limited by cost and time. We identified two areas in our 16S rRNA gene library preparation protocol where modifications could provide efficiency gains, including (1) pooling of multiple PCR amplifications per sample to reduce PCR drift and (2) manual preparation of mastermix to reduce liquid handling. Using nasal samples from healthy human participants and a serially diluted mock microbial community, we compared alpha and beta diversity, and compositional abundance where the PCR amplification was conducted in triplicate, duplicate or as a single reaction, and where manually prepared or premixed mastermix was used. One hundred and fifty-eight 16S rRNA gene sequencing libraries were prepared, including a replicate experiment. Comparing PCR pooling strategies, we found no significant difference in high-quality read counts and alpha diversity, and beta diversity by Bray–Curtis index clustered by replicate on principal coordinate analysis (PCoA) and non-metric dimensional scaling (NMDS) analysis. Choice of mastermix had no significant impact on high-quality read and alpha diversity, and beta diversity by Bray–Curtis index clustered by replicate in PCoA and NMDS analysis. Importantly, we observed contamination and variability of rare species (<0.01 %) across replicate experiments; the majority of contaminants were accounted for by removal of species present at <0.1 %, or were linked to reagents (including a primer stock). We demonstrate no requirement for pooling of PCR amplifications or manual preparation of PCR mastermix, resulting in a more efficient 16S rRNA gene PCR protocol.
-