- Volume 8, Issue 2, 2022
Volume 8, Issue 2, 2022
- Research Articles
-
- Genomic Methodologies
-
-
An inventory of early branch points in microbial phosphonate biosynthesis
More LessMicrobial phosphonate biosynthetic machinery has been identified in ~5 % of bacterial genomes and encodes natural products like fosfomycin as well as cell surface decorations. Almost all biological phosphonates originate from the rearrangement of phosphoenolpyruvate (PEP) to phosphonopyruvate (PnPy) catalysed by PEP mutase (Ppm), and PnPy is often converted to phosphonoacetaldehyde (PnAA) by PnPy decarboxylase (Ppd). Seven enzymes are known or likely to act on either PnPy or PnAA as early branch points en route to diverse biosynthetic outcomes, and these enzymes may be broadly classified into three reaction types: hydride transfer, aminotransfer, and carbon-carbon bond formation. However, the relative abundance of these branch points in microbial phosphonate biosynthesis is unknown. Also unknown is the proportion of ppm-containing gene neighbourhoods encoding new branch point enzymes and potentially novel phosphonates. In this study we computationally sorted 434 ppm-containing gene neighbourhoods based on these seven branch point enzymes. Unsurprisingly, the majority (56 %) of these pathways encode for production of the common naturally occurring compound 2-aminoethylphosphonate (AEP) or a hydroxylated derivative. The next most abundant genetically encoded intermediates were phosphonoalanine (PnAla, 9.2 %), 2-hydroxyethylphosphonate (HEP, 8.5 %), and phosphonoacetate (PnAc, 6 %). Significantly, about 13 % of the gene neighbourhoods could not be assigned to any of the seven branch points and may encode novel phosphonates. Sequence similarity network analysis revealed families of unusual gene neighbourhoods including possible production of phosphonoacrylate and phosphonofructose, the apparent biosynthetic use of the C-P lyase operon, and a virus-encoded phosphonate. Overall, these results highlight the utility of branch point inventories to identify novel gene neighbourhoods and guide future phosphonate discovery efforts.
-
- Functional Genomics and Microbe–Niche Interactions
-
-
Unveiling the genomic potential of Pseudomonas type strains for discovering new natural products
Microbes host a huge variety of biosynthetic gene clusters that produce an immeasurable array of secondary metabolites with many different biological activities such as antimicrobial, anticarcinogenic and antiviral. Despite the complex task of isolating and characterizing novel natural products, microbial genomic strategies can be useful for carrying out these types of studies. However, although genomic-based research on secondary metabolism is on the increase, there is still a lack of reports focusing specifically on the genus Pseudomonas . In this work, we aimed (i) to unveil the main biosynthetic systems related to secondary metabolism in Pseudomonas type strains, (ii) to study the evolutionary processes that drive the diversification of their coding regions and (iii) to select Pseudomonas strains showing promising results in the search for useful natural products. We performed a comparative genomic study on 194 Pseudomonas species, paying special attention to the evolution and distribution of different classes of biosynthetic gene clusters and the coding features of antimicrobial peptides. Using EvoMining, a bioinformatic approach for studying evolutionary processes related to secondary metabolism, we sought to decipher the protein expansion of enzymes related to the lipid metabolism, which may have evolved toward the biosynthesis of novel secondary metabolites in Pseudomonas . The types of metabolites encoded in Pseudomonas type strains were predominantly non-ribosomal peptide synthetases, bacteriocins, N-acetylglutaminylglutamine amides and ß-lactones. Also, the evolution of genes related to secondary metabolites was found to coincide with Pseudomonas species diversification. Interestingly, only a few Pseudomonas species encode polyketide synthases, which are related to the lipid metabolism broadly distributed among bacteria. Thus, our EvoMining-based search may help to discover new types of secondary metabolite gene clusters in which lipid-related enzymes are involved. This work provides information about uncharacterized metabolites produced by Pseudomonas type strains, whose gene clusters have evolved in a species-specific way. Our results provide novel insight into the secondary metabolism of Pseudomonas and will serve as a basis for the prioritization of the isolated strains. This article contains data hosted by Microreact.
-
-
-
Reciprocally rewiring and repositioning the Integration Host Factor (IHF) subunit genes in Salmonella enterica serovar Typhimurium: impacts on physiology and virulence
More LessThe Integration Host Factor (IHF) is a heterodimeric nucleoid-associated protein that plays roles in bacterial nucleoid architecture and genome-wide gene regulation. The ihfA and ihfB genes encode the subunits and are located 350 kbp apart, in the Right replichore of the Salmonella chromosome. IHF is composed of one IhfA and one IhfB subunit. Despite this 1 : 1 stoichiometry, MS revealed that IhfB is produced in 2-fold excess over IhfA. We re-engineered Salmonella to exchange reciprocally the protein-coding regions of ihfA and ihfB, such that each relocated protein-encoding region was driven by the expression signals of the other’s gene. MS showed that in this 'rewired' strain, IhfA is produced in excess over IhfB, correlating with enhanced stability of the hybrid ihfB–ihfA mRNA that was expressed from the ihfB promoter. Nevertheless, the rewired strain grew at a similar rate to the wild-type and was similar in competitive fitness. However, compared to the wild-type, it was less motile, had growth-phase-specific reductions in SPI-1 and SPI-2 gene expression, and was engulfed at a higher rate by RAW macrophage. Our data show that while exchanging the physical locations of its ihf genes and the rewiring of their regulatory circuitry are well tolerated in Salmonella , genes involved in the production of type 3 secretion systems exhibit dysregulation accompanied by altered phenotypes.
-
- Microbial Communities
-
-
A new perspective on ancient Mitis group streptococcal genetics
More LessMitis group Streptococcus are human obligate bacteria residing in the nasopharynx and oral cavity. They comprise both commensal and pathogenic species with the most well-known being Streptococcus pneumoniae – a leading cause of meningitis and pneumonia. A primary difference between the commensal and pathogenic species is the presence of the polysaccharide capsule – a major virulence factor in S. pneumoniae , also present in other commensal species. Our current understanding of the evolutionary divergence of the pathogenic and commensal species has been inferred from extant strains. Ancient genomes can further elucidate streptococcal evolutionary history. We extracted streptococcal genome reads from a 5700-year-old ancient metagenome and worked towards characterizing them. Due to excessive within- and between-species recombination common among streptococci we were unable to parse individual species. Further, the composite reads of the ancient metagenome do not fit within the diversity of any specific extant species. Using a capsular gene database and AT-content analysis we determined that this ancient metagenome is missing polysaccharide synthesis genes integral to streptococcal capsule formation. The presence of multiple zinc metalloproteases suggests that adaptation to host IgA1 had begun and the presence of other virulence factors further implies development of close host–microbe interactions, though the absence of a capsule suggests an inability to cause invasive disease. The presence of specific virulence factors such as pneumolysin implies stable maintenance of such genes through streptococcal evolution that may strengthen their value as anti-pneumococcal vaccine antigens, while maintaining awareness of their potential presence in commensal species. Following from Jensen et al.’s initial analysis we provide historical context for this long time human nasopharyngeal resident, the Mitis group Streptococcus .
-
- Pathogens and Epidemiology
-
-
Genomic differentiation within East Asian Helicobacter pylori
The East Asian region, including China, Japan and Korea, accounts for half of gastric cancer deaths. However, different areas have contrasting gastric cancer incidences and the population structure of Helicobacter pylori in this ethnically diverse region is yet unknown. We aimed to investigate genomic differences in H. pylori between these areas to identify sequence polymorphisms associated with increased cancer risk. We analysed 381 H . pylori genomes collected from different areas of the three countries using phylogenetic and population genetic tools to characterize population differentiation. The functional consequences of SNPs with a highest fixation index (Fst) between subpopulations were examined by mapping amino acid changes on 3D protein structure, solved or modelled. Overall, 329/381 genomes belonged to the previously identified hspEAsia population indicating that import of bacteria from other regions of the world has been uncommon. Seven subregional clusters were found within hspEAsia, related to subpopulations with various ethnicities, geographies and gastric cancer risks. Subpopulation-specific amino acid changes were found in multidrug exporters (hefC), transporters (frpB-4), outer membrane proteins (hopI) and several genes involved in host interaction, such as a catalase site, involved in H2O2 entrance, and a flagellin site mimicking host glycosylation. Several of the top hits, including frpB-4, hefC, alpB/hopB and hofC, have been found to be differentiated within the Americas in previous studies, indicating that a handful of genes may be key to local geographic adaptation. H. pylori within East Asia are not homogeneous but have become differentiated geographically at multiple loci that might have facilitated adaptation to local conditions and hosts. This has important implications for further evaluation of these changes in relation to the varying gastric cancer incidence between geographical areas in this region.
-
-
-
Widespread sharing of pneumococcal strains in a rural African setting: proximate villages are more likely to share similar strains that are carried at multiple timepoints
Madikay Senghore, Chrispin Chaguza, Ebrima Bojang, Peggy-Estelle Tientcheu, Rowan E. Bancroft, Stephanie W. Lo, Rebecca A. Gladstone, Lesley McGee, Archibald Worwui, Ebenezer Foster-Nyarko, Fatima Ceesay, Catherine Bi Okoi, Keith P. Klugman, Robert F. Breiman, Stephen D. Bentley, Richard Adegbola, Martin Antonio, William P. Hanage and Brenda A. Kwambana-AdamsThe transmission dynamics of Streptococcus pneumoniae in sub-Saharan Africa are poorly understood due to a lack of adequate epidemiological and genomic data. Here we leverage a longitudinal cohort from 21 neighbouring villages in rural Africa to study how closely related strains of S. pneumoniae are shared among infants. We analysed 1074 pneumococcal genomes isolated from 102 infants from 21 villages. Strains were designated for unique serotype and sequence-type combinations, and we arbitrarily defined strain sharing where the pairwise genetic distance between strains could be accounted for by the mean within host intra-strain diversity. We used non-parametric statistical tests to assess the role of spatial distance and prolonged carriage on strain sharing using a logistic regression model. We recorded 458 carriage episodes including 318 (69.4 %) where the carried strain was shared with at least one other infant. The odds of strain sharing varied significantly across villages (χ2=47.5, df=21, P-value <0.001). Infants in close proximity to each other were more likely to be involved in strain sharing, but we also show a considerable amount of strain sharing across longer distances. Close geographic proximity (<5 km) between shared strains was associated with a significantly lower pairwise SNP distance compared to strains shared over longer distances (P-value <0.005). Sustained carriage of a shared strain among the infants was significantly more likely to occur if they resided in villages within a 5 km radius of each other (P-value <0.005, OR 3.7). Conversely, where both infants were transiently colonized by the shared strain, they were more likely to reside in villages separated by over 15 km (P-value <0.05, OR 1.5). PCV7 serotypes were rare (13.5 %) and were significantly less likely to be shared (P-value <0.001, OR −1.07). Strain sharing was more likely to occur over short geographical distances, especially where accompanied by sustained colonization. Our results show that strain sharing is a useful proxy for studying transmission dynamics in an under-sampled population with limited genomic data. This article contains data hosted by Microreact.
-
-
-
Diversity of carbapenem-resistant Acinetobacter baumannii and bacteriophage-mediated spread of the Oxa23 carbapenemase
More LessCarbapenem-resistant Acinetobacter baumannii are prevalent in low- and middle-income countries such as Egypt, but little is known about the molecular epidemiology and mechanisms of resistance in these settings. Here, we characterize carbapenem-resistant A. baumannii from Alexandria, Egypt, and place it in a regional context. Fifty-four carbapenem-resistant isolates from Alexandria Main University Hospital (AMUH), Alexandria, Egypt, collected between 2010 and 2015 were genome sequenced using Illumina technology. Genomes were de novo assembled and annotated. Genomes for 36 isolates from the Middle East region were downloaded from GenBank. The core-gene compliment was determined using Roary, and analyses of recombination were performed in Gubbins. Multilocus sequence typing (MLST) sequence type (ST) and antibiotic-resistance genes were identified. The majority of Egyptian isolates belonged to one of three major clades, corresponding to Pasteur MLST clonal complex (CCPAS) 1, CCPAS2 and STPAS158. Strains belonging to STPAS158 have been reported almost exclusively from North Africa, the Middle East and Pakistan, and may represent a region-specific lineage. All isolates carried an oxa23 gene, six carried bla NDM-1 and one carried bla NDM-2. The oxa23 gene was located on a variety of different mobile elements, with Tn2006 predominant in CCPAS2 strains, and Tn2008 predominant in other lineages. Of particular concern, in 8 of the 13 CCPAS1 strains, the oxa23 gene was located in a temperate bacteriophage phiOXA, previously identified only once before in a CCPAS1 clone from the USA military. The carbapenem-resistant A. baumannii population in AMUH is very diverse, and indicates an endemic circulating population, including a region-specific lineage. A major mechanism for oxa23 dissemination in CCPAS1 isolates appears to be a bacteriophage, presenting new concerns about the ability of these carbapenemases to spread throughout the bacterial population.
-
-
-
Genomic diversity and antimicrobial resistance of Prevotella species isolated from chronic lung disease airways
Cystic fibrosis (CF) and chronic obstructive pulmonary disease (COPD) are characterized by increasingly frequent acute pulmonary exacerbations that reduce life quality and length. Human airways are home to a rich polymicrobial environment, which includes members of the obligately anaerobic genus Prevotella . Despite their commonness, surprisingly little is known about the prevalence, role, genomic diversity and antimicrobial resistance (AMR) potential of Prevotella species and strains in healthy and diseased airways. Here, we used comparative genomics to develop a real-time PCR assay to permit rapid Prevotella species identification and quantification from cultures and clinical specimens. Assay specificity was validated across a panel of Prevotella and non- Prevotella species, followed by PCR screening of CF and COPD respiratory-derived cultures. Next, 35 PCR-positive isolates were subjected to whole-genome sequencing. Of eight identified Prevotella species, P. histicola , P. melaninogenica , P. nanceiensis , P. salivae and P. denticola overlapped between participant cohorts. Phylogenomic analysis revealed considerable interhost but limited intrahost diversity, suggesting patient-specific lineages in the lower airways, probably from oral cavity aspirations. Correlation of phenotypic AMR profiles with AMR genes identified excellent correlation between tetQ presence and decreased doxycycline susceptibility, and ermF presence and decreased azithromycin susceptibility and clindamycin resistance. AMR rates were higher in the CF isolates, reflecting greater antibiotic use in this cohort. All tested Prevotella isolates were tobramycin-resistant, providing a potential selection method to improve Prevotella culture retrieval rates. Our addition of 35 airway-derived Prevotella genomes to public databases will enhance ongoing efforts to unravel the role of this diverse and enigmatic genus in both diseased and healthy lungs.
-
-
-
Population genomics of Bacillus anthracis from an anthrax hyperendemic area reveals transmission processes across spatial scales and unexpected within-host diversity
Genomic sequencing has revolutionized our understanding of bacterial disease epidemiology, but remains underutilized for zoonotic pathogens in remote endemic settings. Anthrax, caused by the spore-forming bacterium Bacillus anthracis , remains a threat to human and animal health and rural livelihoods in low- and middle-income countries. While the global genomic diversity of B. anthracis has been well-characterized, there is limited information on how its populations are genetically structured at the scale at which transmission occurs, critical for understanding the pathogen’s evolution and transmission dynamics. Using a uniquely rich dataset, we quantified genome-wide SNPs among 73 B. anthracis isolates derived from 33 livestock carcasses sampled over 1 year throughout the Ngorongoro Conservation Area, Tanzania, a region hyperendemic for anthrax. Genome-wide SNPs distinguished 22 unique B. anthracis genotypes (i.e. SNP profiles) within the study area. However, phylogeographical structure was lacking, as identical SNP profiles were found throughout the study area, likely the result of the long and variable periods of spore dormancy and long-distance livestock movements. Significantly, divergent genotypes were obtained from spatio-temporally linked cases and even individual carcasses. The high number of SNPs distinguishing isolates from the same host is unlikely to have arisen during infection, as supported by our simulation models. This points to an unexpectedly wide transmission bottleneck for B. anthracis , with an inoculum comprising multiple variants being the norm. Our work highlights that inferring transmission patterns of B. anthracis from genomic data will require analytical approaches that account for extended and variable environmental persistence, as well as co-infection.
-
-
-
Towards comprehensive understanding of bacterial genetic diversity: large-scale amplifications in Bordetella pertussis and Mycobacterium tuberculosis
Bacterial genetic diversity is often described solely using base-pair changes despite a wide variety of other mutation types likely being major contributors. Tandem duplication/amplifications are thought to be widespread among bacteria but due to their often-intractable size and instability, comprehensive studies of these mutations are rare. We define a methodology to investigate amplifications in bacterial genomes based on read depth of genome sequence data as a proxy for copy number. We demonstrate the approach with Bordetella pertussis , whose insertion sequence element-rich genome provides extensive scope for amplifications to occur. Analysis of data for 2430 B. pertussis isolates identified 272 putative amplifications, of which 94 % were located at 11 hotspot loci. We demonstrate limited phylogenetic connection for the occurrence of amplifications, suggesting unstable and sporadic characteristics. Genome instability was further described in vitro using long-read sequencing via the Nanopore platform, which revealed that clonally derived laboratory cultures produced heterogenous populations rapidly. We extended this research to analyse a population of 1000 isolates of another important pathogen, Mycobacterium tuberculosis . We found 590 amplifications in M. tuberculosis , and like B. pertussis , these occurred primarily at hotspots. Genes amplified in B. pertussis include those involved in motility and respiration, whilst in M. tuberuclosis, functions included intracellular growth and regulation of virulence. Using publicly available short-read data we predicted previously unrecognized, large amplifications in B. pertussis and M. tuberculosis . This reveals the unrecognized and dynamic genetic diversity of B. pertussis and M. tuberculosis , highlighting the need for a more holistic understanding of bacterial genetics.
-
-
-
Genomic and phenotypic analyses of diverse non-clinical Acinetobacter baumannii strains reveals strain-specific virulence and resistance capacity
Acinetobacter baumannii is a critically important pathogen known for its widespread antibiotic resistance and ability to persist in hospital-associated environments. Whilst the majority of A. baumannii infections are hospital-acquired, infections from outside the hospital have been reported with high mortality. Despite this, little is known about the natural environmental reservoir(s) of A. baumannii and the virulence potential underlying non-clinical strains. Here, we report the complete genome sequences of six diverse strains isolated from environments such as river, soil, and industrial sites around the world. Phylogenetic analyses showed that four of these strains were unrelated to representative nosocomial strains and do not share a monophyletic origin, whereas two had sequence types belonging to the global clone lineages GC1 and GC2. Further, the majority of these strains harboured genes linked to virulence and stress protection in nosocomial strains. These genotypic properties correlated well with in vitro virulence phenotypic assays testing resistance to abiotic stresses, serum survival, and capsule formation. Virulence potential was confirmed in vivo, with most environmental strains able to effectively kill Galleria mellonella greater wax moth larvae. Using phenomic arrays and antibiotic resistance profiling, environmental and nosocomial strains were shown to have similar substrate utilisation patterns although environmental strains were distinctly more sensitive to antibiotics. Taken together, these features of environmental A. baumannii strains suggest the existence of a strain-specific distinct gene pools for niche specific adaptation. Furthermore, environmental strains appear to be equally virulent as contemporary nosocomial strains but remain largely antibiotic sensitive.
-
-
-
Comparative genome analyses of five Vibrio penaeicida strains provide insights into their virulence-related factors
More LessVibrio penaeicida (family Vibrionaceae ) is an important bacterial pathogen that affects Japanese shrimp aquaculture. Only two whole-genome sequences of V. penaeicida are publicly available, which has hampered our understanding of the pathogenesis of shrimp vibriosis caused by this bacterium. To gain insight into the genetic features, evolution and pathogenicity of V. penaeicida, we sequenced five V. penaeicida strains (IFO 15640T, IFO 15641, IFO 15642, TUMSAT-OK1 and TUMSAT-OK2) and performed comparative genomic analyses. Virulence factors and mobile genetic elements were detected. Furthermore, average nucleotide identities (ANIs), clusters of orthologous groups and phylogenetic relationships were evaluated. The V. penaeicida genome consists of two circular chromosomes. Chromosome I sizes ranged from 4.1 to 4.3 Mb, the GC content ranged from 43.9 to 44.1 %, and the number of predicted protein-coding sequences (CDSs) ranged from 3620 to 3782. Chromosome II sizes ranged from 2.2 to 2.4 Mb, the GC content ranged from 43.5 to 43.8 %, and the number of predicted CDSs ranged from 1992 to 2273. All strains except IFO 15641 harboured one plasmid, having sizes that ranged from 150 to 285 kb. All five genomes had typical virulence factors, including adherence, anti-phagocytosis, flagella-related proteins and toxins (repeats-in-toxin and thermolabile haemolysin). The genomes also contained factors responsible for iron uptake and the type II, IV and VI secretion systems. The genome of strain TUMSAT-OK2 tended to encode more prophage regions than the other strains, whereas the genome of strain IFO 15640T had the highest number of regions encoding genomic islands. For comparative genome analysis, we used V. penaeicida (strain CAIM 285T) as a reference strain. ANIs between strain CAIM 285T and the five V. penaeicida strains were >95 %, which indicated that these strains belong to the same species. Orthology cluster analysis showed that strains TUMSAT-OK1 and TUMSAT-OK2 had the greatest number of shared gene clusters, followed by strains CAIM 285T and IFO 15640T. These strains were also the most closely related to each other in a phylogenetic analysis. This study presents the first comparative genome analysis of V. penaeicida and these results will be useful for understanding the pathogenesis of this bacterium.
-
-
-
Population genomics of the food-borne pathogen Vibrio fluvialis reveals lineage associated pathogenicity-related genetic elements
Vibrio fluvialis is a food-borne pathogen with epidemic potential that causes cholera-like acute gastroenteritis and sometimes extraintestinal infections in humans. However, research on its genetic diversity and pathogenicity-related genetic elements based on whole genome sequences is lacking. In this study, we collected and sequenced 130 strains of V. fluvialis from 14 provinces of China, and also determined the susceptibility of 35 of the strains to 30 different antibiotics. Combined with 52 publicly available V. fluvialis genomes, we inferred the population structure and investigated the characteristics of pathogenicity-related factors. The V. fluvialis strains exhibited high levels of homologous recombination and were assigned to two major populations, VflPop1 and VflPop2, according to the different compositions of their gene pools. VflPop2 was subdivided into groups 2.1 and 2.2. Except for VflPop2.2, which consisted only of Asian strains, the strains in VflPop1 and VflPop2.1 were distributed in the Americas, Asia and Europe. Analysis of the pathogenicity potential of V. fluvialis showed that most of the identified virulence-related genes or gene clusters showed high prevalence in V. fluvialis , except for three mobile genetic elements: pBD146, ICEVflInd1 and MGIVflInd1, which were scattered in only a few strains. A total of 21 antimicrobial resistance genes were identified in the genomes of the 182 strains analysed in this study, and 19 (90%) of them were exclusively present in VflPop2. Notably, the tetracycline resistance-related gene tet(35) was present in 150 (95%) of the strains in VflPop2, and in only one (4%) strain in VflPop1, indicating it was population-specific. In total, 91% of the 35 selected strains showed resistance to cefazolin, indicating V. fluvialis has a high resistance rate to cefazolin. Among the 15 genomes that carried the previously reported drug resistance-related plasmid pBD146, 11 (73%) showed resistance to trimethoprim-sulfamethoxazole, which we inferred was related to the presence of the dfr6 gene in the plasmid. On the basis of the population genomics analysis, the genetic diversity, population structure and distribution of pathogenicity-related factors of V. fluvialis were delineated in this study. The results will provide further clues regarding the evolution and pathogenic mechanisms of V. fluvialis , and improve our knowledge for the prevention and control of this pathogen.
-
-
-
Comparative genomics of 26 complete circular genomes of 18 different serotypes of Actinobacillus pleuropneumoniae
More LessActinobacillus pleuropneumoniae is a Gram-negative, rod-shaped bacterium of the family Pasteurellaceae causing pig pleuropneumonia associated with great economic losses worldwide. Nineteen serotypes with distinctive lipopolysaccharide (LPS) and capsular (CPS) compositions have been described so far, yet complete circular genomes are publicly available only for the reference strains of serotypes 1, 4 and 5b, and for field strains of serotypes 1, 3, 7 and 8. We aimed to complete this picture by sequencing the reference strains of 17 different serotypes with the MinION sequencer (Oxford Nanopore Technologies, ONT) and on an Illumina HiSeq (Illumina) platform. We also included two field isolates of serotypes 2 and 3 that were PacBio- and MinION-sequenced, respectively. Genome assemblies were performed following two different strategies, i.e. PacBio- or ONT-only de novo assemblies polished with Illumina reads or a hybrid assembly by directly combining ONT and Illumina reads. Both methods proved successful in obtaining accurate circular genomes with comparable qualities. blast-based genome comparisons and core-genome phylogeny based on core genes, SNP typing and multi-locus sequence typing (cgMLST) of the 26 circular genomes indicated well-conserved genomes across the 18 different serotypes, differing mainly in phage insertions, and CPS, LPS and RTX-toxin clusters, which, consistently, encode serotype-specific antigens. We also identified small antibiotic resistance plasmids, and complete subtype I-F and subtype II-C CRISPR-Cas systems. Of note, highly similar clusters encoding all those serotype-specific traits were also found in other pathogenic and commensal Actinobacillus species. Taken together with the presence of transposable elements surrounding these loci, we speculate a dynamic intra- and interspecies exchange of such virulence-related factors by horizontal gene transfer. In conclusion, our comprehensive genomics analysis provides useful information for diagnostic test and vaccine development, but also for whole-genome-based epidemiological studies, as well as for the surveillance of the evolution of antibiotic resistance and virulence genes in A. pleuropneumoniae .
-
-
-
Comparative genomic analysis of Staphylococcus aureus isolates associated with either bovine intramammary infections or human infections demonstrates the importance of restriction-modification systems in host adaptation
Staphylococcus aureus is a major etiological agent of clinical and subclinical bovine mastitis. The versatile and adaptative evolutionary strategies of this bacterium have challenged mastitis control and prevention globally, and the high incidence of S. aureus mastitis increases concerns about antimicrobial resistance (AMR) and zoonosis. This study aims to describe the evolutionary relationship between bovine intramammary infection (IMI)-associated S. aureus and human pathogenic S. aureus and further elucidate the specific genetic composition that leads to the emergence of successful bovine IMI-associated S. aureus lineages. We performed a phylogenomic analysis of 187 S . aureus isolates that originated from either dairy cattle or humans. Our results revealed that bovine IMI-associated S. aureus isolates showed distinct clades compared to human-originated S. aureus isolates. From a pan-genome analysis, 2070 core genes were identified. Host-specific genes and clonal complex (CC)-specific genes were also identified in bovine S. aureus isolates, mostly located in mobile genetic elements (MGEs). Additionally, the genome sequences of three apparent human-adapted isolates (two from CC97 and one from CC8), isolated from bovine mastitis samples, may provide an snapshot of the genomic characteristics in early host spillover events. Virulence and AMR genes were not conserved among bovine IMI-associated S. aureus isolates. Restriction-modification (R-M) genes in bovine IMI-associated S. aureus demonstrated that the Type I R-M system was lineage-specific and Type II R-M system was sequence type (ST)-specific. The distribution of exclusive, virulence, and AMR genes were closely correlated with the presence of R-M systems in S. aureus , suggesting that R-M systems may contribute to shaping clonal diversification by providing a genetic barrier to the horizontal gene transfer (HGT). Our findings indicate that the CC or ST lineage-specific R-M systems may limit genetic exchange between bovine-adapted S. aureus isolates from different lineages.
-
- Evolution and Responses to Interventions
-
-
Identification of novel, cryptic Clostridioides species isolates from environmental samples collected from diverse geographical locations
Clostridioides difficile is a pathogen often associated with hospital-acquired infection or antimicrobial-induced disease; however, increasing evidence indicates infections can result from community or environmental sources. Most genomic sequencing of C. difficile has focused on clinical strains, although evidence is growing that C. difficile spores are widespread in soil and water in the environment. In this study, we sequenced 38 genomes collected from soil and water isolates in Flagstaff (AZ, USA) and Slovenia in an effort targeted towards environmental surveillance of C. difficile . At the average nucleotide identity (ANI) level, the genomes were divergent to C. difficile at a threshold consistent with different species. A phylogenetic analysis of these divergent genomes together with Clostridioides genomes available in public repositories confirmed the presence of three previously described, cryptic Clostridioide s species and added two additional clades. One of the cryptic species (C-III) was almost entirely composed of Arizona and Slovenia genomes, and contained distinct sub-groups from each region (evidenced by SNP and gene-content differences). A comparative genomics analysis identified multiple unique coding sequences per clade, which can serve as markers for subsequent environmental surveys of these cryptic species. Homologues to the C. difficile toxin genes, tcdA and tcdB, were found in cryptic species genomes, although they were not part of the typical pathogenicity locus observed in C. difficile , and in silico PCR suggested that some would not amplify with widely used PCR diagnostic tests. We also identified gene homologues in the binary toxin cluster, including some present on phage and, for what is believed to be the first time, on a plasmid. All isolates were obtained from environmental samples, so the function and disease potential of these toxin homologues is currently unknown. Enzymatic profiles of a subset of cryptic isolates (n=5) demonstrated differences, suggesting that these isolates contain substantial metabolic diversity. Antimicrobial resistance (AMR) was observed across a subset of isolates (n=4), suggesting that AMR mechanisms are intrinsic to the genus, perhaps originating from a shared environmental origin. This study greatly expands our understanding of the genomic diversity of Clostridioides . These results have implications for C. difficile One Health research, for more sensitive C. difficile diagnostics, as well as for understanding the evolutionary history of C. difficile and the development of pathogenesis.
-
-
-
Comparative genomics and evolutionary analysis of Lactococcus garvieae isolated from human endocarditis
Lactococcus garvieae is a well-known pathogen of fish, but is rarely involved in infections in humans and other mammals. In humans, the main clinical manifestation of L. garvieae infections is endocarditis usually related to the ingestion of contaminated food, such as undercooked fish and shellfish. This study presents the first complete genomic sequence of a clinical L. garvieae strain isolated from a patient with endocarditis and its comparative analysis with other genomes. This human isolate contains a circular chromosome of 2 099 060 bp and one plasmid of 50 557 bp. In comparison with other fully sequenced L. garvieae strains, the chromosomal DNA of L. garvieae Lg-Granada carries a low proportion of insertion sequence elements and a higher number of putative prophages. Our results show that, in general, L. garvieae is a highly recombinogenic species with an open pangenome in which almost 30 % of its genome has undergone horizontal transfers. Within the genus Lactococcus, L. lactis is the main donor of genetic components to L. garvieae but, taking Lg-Granada as a representative, this bacterium tends to import more genes from Bacilli taxa than from other Lactococcus species.
-
-
-
Genomic analysis reveals high intra-species diversity of Shewanella algae
More LessShewanella algae is widely distributed in marine and freshwater habitats, and has been proved to be an emerging marine zoonotic and human pathogen. However, the genomic characteristics and pathogenicity of Shewanella algae are unclear. Here, the whole-genome features of 55 S . algae strains isolated from different sources were described. Pan-genome analysis yielded 2863 (19.4 %) genes shared among all strains. Functional annotation of the core genome showed that the main functions are focused on basic lifestyle such as metabolism and energy production. Meanwhile, the phylogenetic tree of the single nucleotide polymorphisms (SNPs) of core genome divided the 55 strains into three clades, with the majority of strains from China falling into the first two clades. As for the accessory genome, 167 genomic islands (GIs) and 65 phage-related elements were detected. The CRISPR-Cas system with a high degree of confidence was predicted in 23 strains. The GIs carried a suite of virulence genes and mobile genetic elements, while prophages contained several transposases and integrases. Horizontal genes transfer based on homology analysis indicated that these GIs and prophages were parts of major drivers for the evolution and the environmental adaptation of S. algae . In addition, a rich putative virulence-associated gene pool was found. Eight classes of antibiotic-associated resistance genes were detected, and the carriage rate of β-lactam resistance genes was 100 %. In conclusion, S. algae exhibits a high intra-species diversity in the aspects of population structure, virulence-associated genes and potential drug resistance, which is helpful for its evolution in pathogenesis and environmental adaptability.
-
- Short Communications
-
- Genomic Methodologies
-
-
Multi-laboratory evaluation of the Illumina iSeq platform for whole genome sequencing of Salmonella, Escherichia coli and Listeria
Patrick K. Mitchell, Leyi Wang, Bryce J. Stanhope, Brittany D. Cronk, Renee Anderson, Shipra Mohan, Lijuan Zhou, Susan Sanchez, Paula Bartlett, Carol Maddox, Vanessa DeShambo, Rinosh Mani, Lindsy M. Hengesbach, Sarah Gresch, Katie Wright, Sunil Mor, Shuping Zhang, Zhenyu Shen, Lifang Yan, Rebecca Mackey, Rebecca Franklin-Guild, Yan Zhang, Melanie Prarat, Katherine Shiplett, Akhilesh Ramachandran, Sai Narayanan, Justin Sanders, Andree A. Hunkapiller, Kevin Lahmers, Amanda A. Carbonello, Nicole Aulik, Ailam Lim, Jennifer Cooper, Angelica Jones, Jake Guag, Sarah M. Nemser, Gregory H. Tyson, Ruth Timme, Errol Strain, Renate Reimschuessel, Olgica Ceric and Laura B. GoodmanThere is a growing need for public health and veterinary laboratories to perform whole genome sequencing (WGS) for monitoring antimicrobial resistance (AMR) and protecting the safety of people and animals. With the availability of smaller and more affordable sequencing platforms coupled with well-defined bioinformatic protocols, the technological capability to incorporate this technique for real-time surveillance and genomic epidemiology has greatly expanded. There is a need, however, to ensure that data are of high quality. The goal of this study was to assess the utility of a small benchtop sequencing platform using a multi-laboratory verification approach. Thirteen laboratories were provided the same equipment, reagents, protocols and bacterial reference strains. The Illumina DNA Prep and Nextera XT library preparation kits were compared, and 2×150 bp iSeq i100 chemistry was used for sequencing. Analyses comparing the sequences produced from this study with closed genomes from the provided strains were performed using open-source programs. A detailed, step-by-step protocol is publicly available via protocols.io (https://www.protocols.io/view/iseq-bacterial-wgs-protocol-bij8kcrw). The throughput for this method is approximately 4–6 bacterial isolates per sequencing run (20–26 Mb total load). The Illumina DNA Prep library preparation kit produced high-quality assemblies and nearly complete AMR gene annotations. The Prep method produced more consistent coverage compared to XT, and when coverage benchmarks were met, nearly all AMR, virulence and subtyping gene targets were correctly identified. Because it reduces the technical and financial barriers to generating WGS data, the iSeq platform is a viable option for small laboratories interested in genomic surveillance of microbial pathogens.
-
- Pathogens and Epidemiology
-
-
Occurrence of genes encoding spore germination in Clostridium species that cause meat spoilage
More LessMembers of the genus Clostridium are frequently associated with meat spoilage. The ability for low numbers of spores of certain Clostridium species to germinate in cold-stored vacuum-packed meat can result in blown pack spoilage. However, little is known about the germination process of these clostridia, despite this characteristic being important for their ability to cause spoilage. This study sought to determine the genomic conditions for germination of 37 representative Clostridium strains from seven species ( C. estertheticum , C. tagluense , C. frigoris , C. gasigenes , C. putrefaciens , C. aligidicarnis and C. frigdicarnis) by comparison with previously characterized germination genes from C. perfringens , C. sporogenes and C. botulinum . All the genomes analysed contained at least one gerX operon. Seven different gerX operon configuration types were identified across genomes from C. estertheticum , C. tagluense and C. gasigenes . Differences arose between the C. gasigenes genomes and those belonging to C. tagluense / C. estertheticum in the number and type of genes coding for cortex lytic enzymes, suggesting the germination pathway of C. gasigenes is different. However, the core components of the germination pathway were conserved in all the Clostridium genomes analysed, suggesting that these species undergo the same major steps as Bacillus subtilis for germination to occur.
-
- Evolution and Responses to Interventions
-
-
Evidence for common ancestry and microevolution of passerine-adapted Salmonella enterica serovar Typhimurium in the UK and USA
More LessThe evolution of Salmonella enterica serovar Typhimurium (S. Typhimurium) within passerines has resulted in pathoadaptation of this serovar to the avian host in Europe. Recently, we identified an S. Typhimurium lineage from passerines in North America. The emergence of passerine-adapted S. Typhimurium in Europe and North America raises questions regarding its evolutionary origin. Here, we demonstrated that the UK and US passerine-adapted S. Typhimurium shared a common ancestor from ca. 1838, and larids played a key role in the clonal expansion by disseminating the common ancestor between North America and Europe. Further, we identified virulence gene signatures common in the passerine- and larid-adapted S. Typhimurium, including conserved pseudogenes in fimbrial gene lpfD and Type 3 Secretion System (T3SS) effector gene steC. However, the UK and US passerine-adapted S. Typhimurium also possessed unique virulence gene signatures (i.e. pseudogenes in fimbrial gene fimC and T3SS effector genes sspH2, gogB, sseJ and sseK2), and the majority of them (38/47) lost a virulence plasmid pSLT that was present in the larid-adapted S. Typhimurium. These results provide evidence that passerine-adapted S. Typhimurium share a common ancestor with those from larids, and the divergence of passerine- and larid-adapted S. Typhimurium might be due to pseudogenization or loss of specific virulence genes.
-
- Methods
-
- Genomic Methodologies
-
-
DNA isolation methods for Nanopore sequencing of the Streptococcus mitis genome
More LessStreptococcus mitis is a Gram-positive bacterium, member of the oral commensal microbiota, which can occasionally be the etiologic agent of diseases such as infective endocarditis, bacteraemia and septicaemia. The highly recombinogenic and repetitive nature of the S. mitis genome impairs the assembly of a complete genome relying only on short sequencing reads. Oxford Nanopore sequencing can overcome this limitation by generating long reads, enabling the resolution of genomic repeated regions and the assembly of a complete genome sequence. Since the output of a Nanopore sequencing run is strongly influenced by genomic DNA quality and molecular weight, the DNA isolation is the crucial step for an optimal sequencing run. In the present work, we have set up and compared three DNA isolation methods on two S. mitis strains, evaluating their capability of preserving genomic DNA integrity and purity. Sequencing of DNA isolated with a mechanical lysis-based method, despite being cheaper and quicker, did not generate ultra-long reads (maximum read length of 59516 bases) and did not allow the assembly of a circular complete genome. Two methods based on enzymatic lysis of the bacterial cell wall, followed by either (i) a modified CTAB DNA isolation procedure, or (ii) a DNA purification after osmotic lysis of the protoplasts allowed the sequencing of ultra-long reads up to 107294 and 181199 bases in length, respectively. The reconstruction of a circular complete genome was possible sequencing DNAs isolated using the enzymatic lysis-based methods.
-